title: The point-charge self-energy in a nonminimal Lorentz violating Maxwell Electrodynamics

abstract: In this letter we study the self-energy of a point-like charge for the electromagnetic field in a non minimal Lorentz symmetry breaking scenario in a $n+1$ dimensional space time. We consider two variations of a model where the Lorentz violation is caused by a background vector $d^{\nu}$ that appears in a higher derivative interaction. We restrict our attention to the case where $d^{\mu}$ is a time-like background vector, namely $d^{2}=d^{\mu}d_{\mu}>0$, and we verify that the classical self-energy is finite for any odd spatial dimension $n$ and diverges for even $n$. We also make some comments regarding obstacles in the quantization of the proposed model.


title: Lexicase selection in Learning Classifier Systems

abstract: The lexicase parent selection method selects parents by considering performance on individual data points in random order instead of using a fitness function based on an aggregated data accuracy. While the method has demonstrated promise in genetic programming and more recently in genetic algorithms, its applications in other forms of evolutionary machine learning have not been explored. In this paper, we investigate the use of lexicase parent selection in Learning Classifier Systems (LCS) and study its effect on classification problems in a supervised setting. We further introduce a new variant of lexicase selection, called batch-lexicase selection, which allows for the tuning of selection pressure. We compare the two lexicase selection methods with tournament and fitness proportionate selection methods on binary classification problems. We show that batch-lexicase selection results in the creation of more generic rules which is favorable for generalization on future data. We further show that batch-lexicase selection results in better generalization in situations of partial or missing data.


title: Integrable 2d sigma models: quantum corrections to geometry from RG flow

abstract: Classically integrable $\sigma$-models are known to be solutions of the 1-loop RG equations, or "Ricci flow", with only a few couplings running. In some of the simplest examples of integrable deformations we find that in order to preserve this property at 2 (and higher) loops the classical $\sigma$-model should be corrected by quantum counterterms. The pattern is similar to that of effective $\sigma$-models associated to gauged WZW theories. We consider in detail the examples of the $\eta$-deformation of $S^2$ ("sausage model") and $H^2$, as well as the closely related $\lambda$-deformation of the $SO(1,2)/SO(2)$ coset. We also point out that similar counterterms are required in order for non-abelian duality to commute with RG flow beyond the 1-loop order.


title: Metric characterisation of unitaries in JB$^*$-algebras

abstract: Let $M$ be a unital JB$^*$-algebra whose closed unit ball is denoted by $\mathcal{B}_M$. Let $\partial_e(\mathcal{B}_M)$ denote the set of all extreme points of $\mathcal{B}_M$. We prove that an element $u\in \partial_e(\mathcal{B}_M)$ is a unitary if and only if the set $$\mathcal{M}_{u} = \{e\in \partial_e(\mathcal{B}_M) : \|u\pm e\|\leq \sqrt{2} \}$$ contains an isolated point. This is a new geometric characterisation of unitaries in $M$ in terms of the set of extreme points of $\mathcal{B}_M$.


title: The Impact of Private and Work-Related Smartphone Usage on Interruptibility

abstract: In the last decade, the effects of interruptions through mobile notifications have been extensively researched in the field of Human-Computer Interaction. Breakpoints in tasks and activities, cognitive load, and personality traits have all been shown to correlate with individuals' interruptibility. However, concepts that explain interruptibility in a broader sense are needed to provide a holistic understanding of its characteristics. In this paper, we build upon the theory of social roles to conceptualize and investigate the correlation between individuals' private and work-related smartphone usage and their interruptibility. Through our preliminary study with four participants over 11 weeks, we found that application sequences on smartphones correlate with individuals' private and work roles. We observed that participants engaged in these roles tend to follow specific interruptibility strategies - integrating, combining, or segmenting private and work-related engagements. Understanding these strategies breaks new ground for attention and interruption management systems in ubiquitous computing.


title: Optimal mechanisms with budget for user generated contents

abstract: In this paper, we design gross product maximization mechanisms which incentivize users to upload high-quality contents on user-generated-content (UGC) websites. We show that, the proportional division mechanism, which is widely used in practice, can perform arbitrarily bad in the worst case. The problem can be formulated using a linear program with bounded and increasing variables. We then present an $O(n\log n)$ algorithm to find the optimal mechanism, where n is the number of players.


title: Matroid Bases with Cardinality Constraints on the Intersection

abstract: Given two matroids $\mathcal{M}_{1} = (E, \mathcal{B}_{1})$ and $\mathcal{M}_{2} = (E, \mathcal{B}_{2})$ on a common ground set $E$ with base sets $\mathcal{B}_{1}$ and $\mathcal{B}_{2}$, some integer $k \in \mathbb{N}$, and two cost functions $c_{1}, c_{2} \colon E \rightarrow \mathbb{R}$, we consider the optimization problem to find a basis $X \in \mathcal{B}_{1}$ and a basis $Y \in \mathcal{B}_{2}$ minimizing cost $\sum_{e\in X} c_1(e)+\sum_{e\in Y} c_2(e)$ subject to either a lower bound constraint $|X \cap Y| \le k$, an upper bound constraint $|X \cap Y| \ge k$, or an equality constraint $|X \cap Y| = k$ on the size of the intersection of the two bases $X$ and $Y$. The problem with lower bound constraint turns out to be a generalization of the Recoverable Robust Matroid problem under interval uncertainty representation for which the question for a strongly polynomial-time algorithm was left as an open question by Hradovich et al. We show that the two problems with lower and upper bound constraints on the size of the intersection can be reduced to weighted matroid intersection, and thus be solved with a strongly polynomial-time primal-dual algorithm. The question whether the problem with equality constraint can also be solved efficiently turned out to be a lot harder. As our main result, we present a strongly-polynomial, primal-dual algorithm for the problem with equality constraint on the size of the intersection. Additionally, we discuss generalizations of the problems from matroids to polymatroids, and from two to three or more matroids.


title: Degeneration of spectral sequences and complex Lagrangian submanifolds

abstract: There is a local-to-global $\mathrm{Ext}$ spectral sequence $\mathrm{E}_2^{p,q} = \mathrm{H}^p(\mathrm{L}, \Omega^q_\mathrm{L}) \Rightarrow \mathrm{Ext}^{p+q}(i_*\mathscr{O}_\mathrm{L}, i_*\mathscr{O}_{\mathrm{L}})$ for a smooth Lagrangian subvariety in a hyperk\"ahler variety. We prove its degeneration on $\mathrm{E}_{2}$, and various generalisations thereof.


title: Interpretable Deep Learning Model for the Detection and Reconstruction of Dysarthric Speech

abstract: This paper proposed a novel approach for the detection and reconstruction of dysarthric speech. The encoder-decoder model factorizes speech into a low-dimensional latent space and encoding of the input text. We showed that the latent space conveys interpretable characteristics of dysarthria, such as intelligibility and fluency of speech. MUSHRA perceptual test demonstrated that the adaptation of the latent space let the model generate speech of improved fluency. The multi-task supervised approach for predicting both the probability of dysarthric speech and the mel-spectrogram helps improve the detection of dysarthria with higher accuracy. This is thanks to a low-dimensional latent space of the auto-encoder as opposed to directly predicting dysarthria from a highly dimensional mel-spectrogram.


title: Modeling Semantic Compositionality with Sememe Knowledge

abstract: Semantic compositionality (SC) refers to the phenomenon that the meaning of a complex linguistic unit can be composed of the meanings of its constituents. Most related works focus on using complicated compositionality functions to model SC while few works consider external knowledge in models. In this paper, we verify the effectiveness of sememes, the minimum semantic units of human languages, in modeling SC by a confirmatory experiment. Furthermore, we make the first attempt to incorporate sememe knowledge into SC models, and employ the sememeincorporated models in learning representations of multiword expressions, a typical task of SC. In experiments, we implement our models by incorporating knowledge from a famous sememe knowledge base HowNet and perform both intrinsic and extrinsic evaluations. Experimental results show that our models achieve significant performance boost as compared to the baseline methods without considering sememe knowledge. We further conduct quantitative analysis and case studies to demonstrate the effectiveness of applying sememe knowledge in modeling SC. All the code and data of this paper can be obtained on https://github.com/thunlp/Sememe-SC.


title: Constant-Time Dynamic $(\Delta+1)$-Coloring and Weight Approximation for Minimum Spanning Forest: Dynamic Algorithms Meet Property Testing

abstract: With few exceptions (namely, algorithms for maximal matching, $2$-approximate vertex cover, and certain constant-stretch spanners), all known fully dynamic algorithms in general graphs require (amortized) $\Omega(\log n)$ update/query time. Showing for the first time that techniques from property testing can lead to constant-time fully dynamic graph algorithms we prove the following results: (1) We give a fully dynamic (Las-Vegas style) algorithm with constant expected amortized time per update that maintains a proper $(\Delta+1)$-vertex coloring of a graph with maximum degree at most $\Delta$. This improves upon the previous $O(\log \Delta)$-time algorithm by Bhattacharya et al. (SODA 2018). We show that our result does not only have optimal running time, but is also optimal in the sense that already deciding whether a $\Delta$-coloring exists in a dynamically changing graph with maximum degree at most $\Delta$ takes $\Omega(\log n)$ time per operation. (2) We give two fully dynamic algorithms that maintain a $(1+\varepsilon)$-approximation of the weight $M$ of the minimum spanning forest of a graph $G$ with edges weights in $[1,W]$. Our deterministic algorithm takes $O({W^2 \log W}/{\varepsilon^3})$ worst-case time, which is constant if both $W$ and $\varepsilon$ are constant. This is somewhat surprising as a lower bound by Patrascu and Demaine (SIAM J. Comput. 2006) shows that it takes $\Omega(\log n)$ time per operation to maintain the exact weight of the MSF that holds even for $W=1$. Our randomized (Monte-Carlo style) algorithm works with high probability and runs in worst-case $O(\frac{1}{\varepsilon^4}\log^2(\frac{1}{\varepsilon}))$ time if $W= O({(m^*)^{1/3}}/{\log^3 n})$, where $m^*$ is the minimum number of edges in the graph throughout all the updates. It works even against an adaptive adversary.


title: On-axis scattering of scalar fields by charged rotating black holes

abstract: We investigate the scattering of a massless scalar field by a Kerr-Newman black hole, considering the case of on-axis incidence. We use the partial wave method to find numerical results for the scattering cross section, which we compare with classical and semiclassical analytical results, obtaining excellent agreement. We present a selection of plots for different values of the black hole parameters.


title: Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method

abstract: By listening to gravity in the low frequency band, between 0.1 mHz and 1 Hz, the future space-based gravitational-wave observatory LISA will be able to detect tens of thousands of astrophysical sources from cosmic dawn to the present. The detection and characterization of all resolvable sources is a challenge in itself, but LISA data analysis will be further complicated by interruptions occurring in the interferometric measurements. These interruptions will be due to various causes occurring at various rates, such as laser frequency switches, high-gain antenna re-pointing, orbit corrections, or even unplanned random events. Extracting long-lasting gravitational-wave signals from gapped data raises problems such as noise leakage and increased computational complexity. We address these issues by using Bayesian data augmentation, a method that reintroduces the missing data as auxiliary variables in the sampling of the posterior distribution of astrophysical parameters. This provides a statistically consistent way to handle gaps while improving the sampling efficiency and mitigating leakage effects. We apply the method to the estimation of galactic binaries parameters with different gap patterns, and we compare the results to the case of complete data.


title: Spin-spin correlations in central rows of Ising models with holes

abstract: In our previous works on infinite horizontal Ising strips of width $m$ alternating with layers of strings of Ising chains of length $n$, we found the surprising result that the specific heats are not much different for different values of $N$, the separation of the strings. For this reason, we study here for $N=1$ the spin-spin correlation in the central row of each strip, and also the central row of a strings layer. We show that these can be written as a Toeplitz determinants. Their generating functions are ratios of two polynomials, which in the limit of infinite vertical size become square roots of polynomials whose degrees are $m+1$ where $m$ is the size of the strips. We find the asymptotic behaviors near the critical temperature to be two-dimensional Ising-like. But in regions not very close to criticality the behavior may be different for different $m$ and $n$. Finally, in the appendix we shall present results for generating functions in more general models.


title: Dense Peelable Random Uniform Hypergraphs

abstract: We describe a new family of $k$-uniform hypergraphs with independent random edges. The hypergraphs have a high probability of being peelable, i.e. to admit no sub-hypergraph of minimum degree $2$, even when the edge density (number of edges over vertices) is close to $1$. In our construction, the vertex set is partitioned into linearly arranged segments and each edge is incident to random vertices of $k$ consecutive segments. Quite surprisingly, the linear geometry allows our graphs to be peeled "from the outside in". The density thresholds $f_k$ for peelability of our hypergraphs ($f_3 \approx 0.918$, $f_4 \approx 0.977$, $f_5 \approx 0.992$, ...) are well beyond the corresponding thresholds ($c_3 \approx 0.818$, $c_4 \approx 0.772$, $c_5 \approx 0.702$, ...) of standard $k$-uniform random hypergraphs. To get a grip on $f_k$, we analyse an idealised peeling process on the random weak limit of our hypergraph family. The process can be described in terms of an operator on functions and $f_k$ can be linked to thresholds relating to the operator. These thresholds are then tractable with numerical methods. Random hypergraphs underlie the construction of various data structures based on hashing. These data structures frequently rely on peelability of the hypergraph or peelability allows for simple linear time algorithms. To demonstrate the usefulness of our construction, we used our $3$-uniform hypergraphs as a drop-in replacement for the standard $3$-uniform hypergraphs in a retrieval data structure by Botelho et al. This reduces memory usage from $1.23m$ bits to $1.12m$ bits ($m$ being the input size) with almost no change in running time.


title: Efficient Gauss Elimination for Near-Quadratic Matrices with One Short Random Block per Row, with Applications

abstract: In this paper we identify a new class of sparse near-quadratic random Boolean matrices that have full row rank over $\mathbb{F}_2=\{0,1\}$ with high probability and can be transformed into echelon form in almost linear time by a simple version of Gauss elimination. The random matrix with dimensions $n(1-\varepsilon) \times n$ is generated as follows: In each row, identify a block of length $L=O((\log n)/\varepsilon)$ at a random position. The entries outside the block are 0, the entries inside the block are given by fair coin tosses. Sorting the rows according to the positions of the blocks transforms the matrix into a kind of band matrix, on which, as it turns out, Gauss elimination works very efficiently with high probability. For the proof, the effects of Gauss elimination are interpreted as a ("coin-flipping") variant of Robin Hood hashing, whose behaviour can be captured in terms of a simple Markov model from queuing theory. Bounds for expected construction time and high success probability follow from results in this area. By employing hashing, this matrix family leads to a new implementation of a retrieval data structure, which represents an arbitrary function $f\colon S \to \{0,1\}$ for some set $S$ of $m=(1-\varepsilon)n$ keys. It requires $m/(1-\varepsilon)$ bits of space, construction takes $O(m/\varepsilon^2$) expected time on a word RAM, while queries take $O(1/\varepsilon)$ time and access only one contiguous segment of $O((\log m)/\varepsilon)$ bits in the representation. The method is competitive with state-of-the-art methods. By well-established methods the retrieval data structure leads to efficient constructions of (static) perfect hash functions and (static) Bloom filters with almost optimal space and very local storage access patterns for queries.


title: Julian Schwinger and the Semiclassical Atom

abstract: In the early 1980s, Schwinger made seminal contributions to the semiclassical theory of atoms. There had, of course, been earlier attempts at improving upon the Thomas--Fermi model of the 1920s. Yet, a consistent derivation of the leading and next-to-leading corrections to the formula for the total binding energy of neutral atoms, $$-\frac{E}{e^2/a_0} = 0.768745Z^{7/3} - \frac{1}{2}Z^2+0.269900Z^{5/3} + \cdots\,,$$ had not been accomplished before Schwinger got interested in the matter; here, $Z$ is the atomic number and $e^2/a_0$ is the Rydberg unit of energy. The corresponding improvements upon the Thomas--Fermi density were next on his agenda with, perhaps, less satisfactory results. Schwinger's work not only triggered extensive investigations by mathematicians, who eventually convinced themselves that Schwinger got it right, but also laid the ground, in passing, for later refinements --- some of them very recent.


title: Sparse Regular Expression Matching

abstract: A regular expression specifies a set of strings formed by single characters combined with concatenation, union, and Kleene star operators. Given a regular expression $R$ and a string $Q$, the regular expression matching problem is to decide if $Q$ matches any of the strings specified by $R$. Regular expressions are a fundamental concept in formal languages and regular expression matching is a basic primitive for searching and processing data. A standard textbook solution [Thompson, CACM 1968] constructs and simulates a nondeterministic finite automaton, leading to an $O(nm)$ time algorithm, where $n$ is the length of $Q$ and $m$ is the length of $R$. Despite considerable research efforts only polylogarithmic improvements of this bound are known. Recently, conditional lower bounds provided evidence for this lack of progress when Backurs and Indyk [FOCS 2016] proved that, assuming the strong exponential time hypothesis (SETH), regular expression matching cannot be solved in $O((nm)^{1-\epsilon})$, for any constant $\epsilon > 0$. Hence, the complexity of regular expression matching is essentially settled in terms of $n$ and $m$. In this paper, we take a new approach and introduce a \emph{density} parameter, $\Delta$, that captures the amount of nondeterminism in the NFA simulation on $Q$. The density is at most $nm+1$ but can be significantly smaller. Our main result is a new algorithm that solves regular expression matching in $$O\left(\Delta \log \log \frac{nm}{\Delta} +n + m\right)$$ time. This essentially replaces $nm$ with $\Delta$ in the complexity of regular expression matching. We complement our upper bound by a matching conditional lower bound that proves that we cannot solve regular expression matching in time $O(\Delta^{1-\epsilon})$ for any constant $\epsilon > 0$ assuming SETH.


title: Endpoint estimates for the maximal function over prime numbers

abstract: Given an ergodic dynamical system $(X, \mathcal{B}, \mu, T)$, we prove that for each function $f$ belonging to the Orlicz space $L(\log L)^2(\log \log L)(X, \mu)$, the ergodic averages \[ \frac{1}{\pi(N)} \sum_{p \in \mathbb{P}_N} f\big(T^p x\big), \] converge for $\mu$-almost all $x \in X$, where $\mathbb{P}_N$ is the set of prime numbers not larger that $N$ and $\pi(N) = \# \mathbb{P}_N$.


title: The dual volume of quasi-Fuchsian manifolds and the Weil-Petersson distance

abstract: Making use of the dual Bonahon-Schl\"afli formula, we prove that the dual volume of the convex core of a quasi-Fuchsian manifold $M$ is bounded by an explicit constant, depending only on the topology of $M$, times the Weil-Petersson distance between the hyperbolic structures on the upper and lower boundary components of the convex core of $M$.


title: Ergodicity and partial hyperbolicity on Seifert manifolds

abstract: We show that conservative partially hyperbolic diffeomorphism isotopic to the identity on Seifert 3-manifolds are ergodic.


title: The Atacama Large Aperture Submillimeter Telescope (AtLAST)

abstract: The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.


title: Thin Film Liquid Crystals with Oblique Anchoring and Boojums

abstract: We study a two-dimensional variational problem which arises as a thin-film limit of the Landau-de Gennes energy of nematic liquid crystals. We impose an oblique angle condition for the nematic director on the boundary, via boundary penalization (weak anchoring.) We show that for strong anchoring strength (relative to the usual Ginzburg-Landau length scale parameter,) defects will occur in the interior, as in the case of strong (Dirichlet) anchoring, but for weaker anchoring strength all defects will occur on the boundary. These defects will each carry a fractional winding number; such boundary defects are known as "boojums". The boojums will occur in ordered pairs along the boundary; for nematic director with angle oblique to the normal vector, they serve to reduce the winding of the phase in two steps, in order to avoid the formation of interior defects. We determine the number and location of the defects via a Renormalized Energy and numerical simulations.


title: SynthCity: A large scale synthetic point cloud

abstract: With deep learning becoming a more prominent approach for automatic classification of three-dimensional point cloud data, a key bottleneck is the amount of high quality training data, especially when compared to that available for two-dimensional images. One potential solution is the use of synthetic data for pre-training networks, however the ability for models to generalise from synthetic data to real world data has been poorly studied for point clouds. Despite this, a huge wealth of 3D virtual environments exist which, if proved effective can be exploited. We therefore argue that research in this domain would be of significant use. In this paper we present SynthCity an open dataset to help aid research. SynthCity is a 367.9M point synthetic full colour Mobile Laser Scanning point cloud. Every point is assigned a label from one of nine categories. We generate our point cloud in a typical Urban/Suburban environment using the Blensor plugin for Blender.


title: Toward a Procedural Fruit Tree Rendering Framework for Image Analysis

abstract: We propose a procedural fruit tree rendering framework, based on Blender and Python scripts allowing to generate quickly labeled dataset (i.e. including ground truth semantic segmentation). It is designed to train image analysis deep learning methods (e.g. in a robotic fruit harvesting context), where real labeled training datasets are usually scarce and existing synthetic ones are too specialized. Moreover, the framework includes the possibility to introduce parametrized variations in the model (e.g. lightning conditions, background), producing a dataset with embedded Domain Randomization aspect.


title: Bound state solutions of the Klein Gordon equation with energy-dependent potentials

abstract: In this manuscript, we investigate the exact bound state solution of the Klein-Gordon equation for an energy-dependent Coulomb-like vector plus scalar potential energies. To the best of our knowledge, this problem is examined in literature with a constant and position dependent mass functions. As a novelty, we assume a mass-function that depends on energy and position and revisit the problem with the following cases: First, we examine the case where the mixed vector and scalar potential energy possess equal magnitude and equal sign as well as an opposite sign. Then, we study pure scalar and pure vector cases. In each case, we derive an analytic expression of the energy spectrum by employing the asymptotic iteration method. We obtain a non-trivial relation among the tuning parameters which lead the examined problem to a constant mass one. Finally, we calculate the energy spectrum by the Secant method and show that the corresponding unnormalized wave functions satisfy the boundary conditions. We conclude the manuscript with a comparison of the calculated energy spectra versus tuning parameters.


title: Towards Affordance Prediction with Vision via Task Oriented Grasp Quality Metrics

abstract: While many quality metrics exist to evaluate the quality of a grasp by itself, no clear quantification of the quality of a grasp relatively to the task the grasp is used for has been defined yet. In this paper we propose a framework to extend the concept of grasp quality metric to task-oriented grasping by defining affordance functions via basic grasp metrics for an open set of task affordances. We evaluate both the effectivity of the proposed task oriented metrics and their practical applicability by learning to infer them from vision. Indeed, we assess the validity of our novel framework both in the context of perfect information, i.e., known object model, and in the partial information context, i.e., inferring task oriented metrics from vision, underlining advantages and limitations of both situations. In the former, physical metrics of grasp hypotheses on an object are defined and computed in known object model simulation, in the latter deep models are trained to infer such properties from partial information in the form of synthesized range images.


title: Electronic Transport Evidence for Topological Nodal-Line Semimetals of ZrGeSe single crystals

abstract: Although the band topology of ZrGeSe has been studied via magnetic torque technique, the electronic transport behaviors related to the relativistic Fermions in ZrGeSe are still unknown. Here, we first report systematic electronic transport properties of high-quality ZrGeSe single crystals under magnetic fields up to 14 T. Resistivity plateaus of temperature dependent resistivity curves both in the presence and absence of magnetic fields as well as large, non-saturating magnetoresistance in low-temperature region were observed. By analyzing the temperature- and angular-dependent Shubnikov-de Haas oscillations and fitting it via the Lifshitz-Kosevich (LK) formula with the Berry phase being taken into account, we proved that Dirac fermions dominate the electronic transport behaviors of ZrGeSe and the presence of non-trivial Berry phase. First principles calculations demonstrate that ZrGeSe possesses Dirac bands and normal bands near Fermi surface, resulting in the observed magnetotransport phenomena. These results demonstrate that ZrGeSe is a topological nodal-line semimetal, which provides a fundamentally important platform to study the quantum physics of topological semimetals.


title: Approximate Bayesian inference for analysis of spatio-temporal flood frequency data

abstract: Extreme floods cause casualties, and widespread damage to property and vital civil infrastructure. We here propose a Bayesian approach for predicting extreme floods using the generalized extreme-value (GEV) distribution within gauged and ungauged catchments. A major methodological challenge is to find a suitable parametrization for the GEV distribution when covariates or latent spatial effects are involved. Other challenges involve balancing model complexity and parsimony using an appropriate model selection procedure, and making inference using a reliable and computationally efficient approach. Our approach relies on a latent Gaussian modeling framework with a novel multivariate link function designed to separate the interpretation of the parameters at the latent level and to avoid unreasonable estimates of the shape and time trend parameters. Structured additive regression models are proposed for the four parameters at the latent level. For computational efficiency with large datasets and richly parametrized models, we exploit an accurate and fast approximate Bayesian inference approach. We applied our proposed methodology to annual peak river flow data from 554 catchments across the United Kingdom (UK). Our model performed well in terms of flood predictions for both gauged and ungauged catchments. The results show that the spatial model components for the transformed location and scale parameters, and the time trend, are all important. Posterior estimates of the time trend parameters correspond to an average increase of about $1.5\%$ per decade and reveal a spatial structure across the UK. To estimate return levels for spatial aggregates, we further develop a novel copula-based post-processing approach of posterior predictive samples, in order to mitigate the effect of the conditional independence assumption at the data level, and we show that our approach provides accurate results.


title: Spontaneous Radiofrequency Emission from Electron Spins within Drosophila: a preliminary report on a novel biological signal

abstract: Using simple radiofrequency (RF) instrumentation, we detect spontaneous RF emission from Drosophila immersed in a magnetic field within a waveguide or a shielded RF resonator. Remarkably, the RF emissions are abolished by chloroform anesthesia. Conversely, activation of the nervous system by temperature-sensitive cation channels causes bursts of radiofrequency emission. Both observations suggest that the radiofrequency emissions come from the nervous system. RF frequency and magnetic field dependence are consistent with RF emission being related to electron Zeeman energy. Since the RF emission occurs without external energy input, its energy must come from cell metabolism.We propose that RF emissions are due to in vivo spin-polarised cellular electron currents. We suggest that spin-polarised currents relax radiatively to equilibrium at the Larmor frequency and are detected under our experimental conditions (resonator cavity and lock-in detection). We invite scrutiny of our methods and results, aim to elicit theorists interest in this novel phenomenon, and encourage replication in this and other preparations.


title: Evaluating bird collision risk of a high-speed railway crossing the habitat of the crested ibis (Nipponia nippon) in Qinling Mountains, China

abstract: Bird collisions with high-speed transport modes is a vital topic on vehicle safety and wildlife protection, especially when high-speed trains, with an average speed of 250km/h, have to run across the habitat of an endangered bird species. This paper evaluates the bird-train collision risk associated with a recent high-speed railway project in Qinling Mountains, China, for the crested ibis (Nipponia nippon) and other local bird species. Using line transect surveys and walking monitoring techniques, we surveyed the population abundance, spatial-temporal distributions, and bridge-crossing behaviors of the birds in the study area. The results show that: (1) The crested ibis and the egret were the two most abundant waterfowl species in the study area. The RAI of these two species were about 43.69% and 42.91%, respectively; (2) Crested ibises overall habitat closer to the railway bridge. 91.63% of them were firstly detected within the range of 0m to 25m of the vicinity of the bridge; (3) the ratio between crossing over and under the railway bridge was about 7:3. Crested ibises were found to prefer to fly over the railway bridge (89.29% of the total crossing activities observed for this species). Egrets were more likely to cross the railway below the bridge, and they accounted for 60.27% of the total observations of crossing under the bridge. We recommend that, while the collision risk of crested ibises could be low, barrier-like structures, such as fences, should still be considered to promote the conservation of multiple bird species in the area. This paper provides a practical case for railway ecology studies in China. To our best knowledge, this is the first high-speed railway project that takes protecting crested ibises as one of the top priorities, and exemplifies the recent nationwide initiative towards the construction of "eco-civilization" in the country.


title: Broken adiabaticity induced by Lifshitz transition in MoS$_2$ and WS$_2$ single layers

abstract: The breakdown of the adiabatic Born-Oppenheimer approximation is striking dynamical phenomenon, however, it occurs only in a handful of layered materials. Here, I show that adiabaticity breaks down in doped single-layer transition metal dichalcogenides in a quite intriguing manner. Namely, significant nonadiabatic coupling, which acts on frequencies of the Raman-active modes, is prompted by a Lifshitz transition due to depopulation and population of multiple valence and conduction valleys, respectively. The outset of the latter event is shown to be dictated by the interplay of highly non-local electron-electron interaction and spin-orbit coupling. In addition, intense electron-hole pair scatterings due to electron-phonon coupling are inducing phonon linewidth modifications as a function of doping. Comprehending these intricate dynamical effects turns out to be a key for mastering characterization of electron doping in two-dimensional nano-devices by means of Raman spectroscopy.


title: The role of topological defects in the two-stage melting and elastic behavior of active Brownian particles

abstract: We find that crystalline states of repulsive active Brownian particles at high activity melt into a hexatic state but this transition is not driven by an unbinding of bound dislocation pairs as suggested by the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) theory. Upon reducing the density, the crystalline state melts into a high-density hexatic state devoid of any defects. Decreasing the density further, the dislocations proliferate and introduce plasticity in the system, nevertheless maintaining the hexatic state, but eventually melting into a fluid state. Remarkably, the elastic constants of active solids are equal to those of their passive counterparts, as the swim contribution to the stress tensor is negligible in the solid state. The sole effect of activity is that the stable solid regime shifts to higher densities. Furthermore, discontinuities in the elastic constants as a function of density correspond to changes in the defect concentrations rather than to the solid-hexatic transition.


title: Kippenhahn's Theorem for joint numerical ranges and quantum states

abstract: Kippenhahn's Theorem asserts that the numerical range of a matrix is the convex hull of a certain algebraic curve. Here, we show that the joint numerical range of finitely many Hermitian matrices is similarly the convex hull of a semi-algebraic set. We discuss an analogous statement regarding the dual convex cone to a hyperbolicity cone and prove that the class of bases of these dual cones is closed under linear operations. The result offers a new geometric method to analyze quantum states.


title: Improving Variational Quantum Optimization using CVaR

abstract: Hybrid quantum/classical variational algorithms can be implemented on noisy intermediate-scale quantum computers and can be used to find solutions for combinatorial optimization problems. Approaches discussed in the literature minimize the expectation of the problem Hamiltonian for a parameterized trial quantum state. The expectation is estimated as the sample mean of a set of measurement outcomes, while the parameters of the trial state are optimized classically. This procedure is fully justified for quantum mechanical observables such as molecular energies. In the case of classical optimization problems, which yield diagonal Hamiltonians, we argue that aggregating the samples in a different way than the expected value is more natural. In this paper we propose the Conditional Value-at-Risk as an aggregation function. We empirically show -- using classical simulation as well as quantum hardware -- that this leads to faster convergence to better solutions for all combinatorial optimization problems tested in our study. We also provide analytical results to explain the observed difference in performance between different variational algorithms.


title: Hard topological versus soft geometrical magnetic particle transport

abstract: The transport on top of a periodic two-dimensional hexagonal magnetic pattern of (i) a single macroscopic steel sphere, (ii) a doublet of wax/magnetite composite spheres, and (iii) an immiscible mixture of ferrofluid droplets with a perfluorinated liquid is studied experimentally and analyzed theoretically. The transport of all these magnetic objects is achieved by moving an external permanent magnet on a closed modulation loop around the two-dimensional magnetic pattern. The transport of one and also that of two objects per unit cell is topologically protected and characterized by discrete displacements of the particles as we continuously scan through a family of modulation loops. The direction and the type of transport is characterized by the winding numbers of the modulation loops around special objects in control space, which is the space of possible directions of the external magnetic field. The number of winding numbers necessary for characterizing the topological transport increases with the number of particles per unit cell. The topological character of the transport is destroyed when transporting a large collection of particles per unit cell like it is the case for a macroscopic assembly of magnetic nano particles in a ferrofluid droplet for which the transport is geometrical and no longer topological. To characterize the change in the transport from topological to geometrical, we perform computer simulations of the transport of an increasing number of particles per unit cell.


title: How cold gas continuously entrains mass and momentum from a hot wind

abstract: The existence of fast moving, cold gas ubiquitously observed in galactic winds is theoretically puzzling, since the destruction time of cold gas is much smaller than its acceleration time. In previous work, we showed that cold gas can accelerate to wind speeds and grow in mass if the radiative cooling time of mixed gas is shorter than the cloud destruction time. Here, we study this process in much more detail, and find remarkably robust cloud acceleration and growth in a wide variety of scenarios. Radiative cooling, rather than the Kelvin-Helmholtz instability, enables self-sustaining entrainment of hot gas onto the cloud via cooling-induced pressure gradients. Indeed, growth peaks when the cloud is almost co-moving. The entrainment velocity is of order the cold gas sound speed, and growth is accompanied by cloud pulsations. Growth is also robust to the background wind and initial cloud geometry. In an adiabatic Chevalier-Clegg type wind, for instance, the mass growth rate is constant. Although growth rates are similar with magnetic fields, cloud morphology changes dramatically, with low density, magnetically supported filaments which have a small mass fraction but dominate by volume. This could bias absorption line observations. Cloud growth from entraining and cooling hot gas can potentially account for the cold gas content of the CGM. It can also fuel star formation in the disk as cold gas recycled in a galactic fountain accretes and cools halo gas. We speculate that galaxy-scale simulations should converge in cold gas mass once cloud column densities of ${\rm N} \sim 10^{18} \ {\rm cm^{-2}}$ are resolved.


title: Coulomb Energy Density Functionals for Nuclear Systems: Recent Studies of Coulomb Exchange and Correlation Functionals

abstract: The Coulomb exchange and correlation energy density functionals for electron systems are applied to nuclear systems. It is found that the exchange functionals in the generalized gradient approximation provide agreements with the exact-Fock energy with one adjustable parameter within a few dozen $ \mathrm{keV} $ accuracy, whereas the correlation functionals are not directly applicable to nuclear systems due to the existence of the nuclear force.


title: Eppur si riscalda -- and yet, it (just) heats up: Further Comments on "Quantifying hot carrier and thermal contributions in plasmonic photocatalysis"

abstract: Our Comment [Sivan et al., Science 2019] (as well as its longer version [Dubi, Un, & Sivan, ArXiv 2019], and the supporting theoretical studies [Dubi and Sivan, ArXiv 2018]) on recent attempts to distinguish thermal and non-thermal ("hot carrier") contributions to plasmon-assisted photocatalysis [Zhou et al., Science 2018] initiated a re-evaluation process of previous literature on the topic within the nano-plasmonics and chemistry communities. The Response of Zhou et al. attempts to defend the claims of the original paper. In this manuscript, we show that the Response presents additional data that further validates our central criticism: inaccurately measured temperatures (that are lower than the actual temperature of the catalyst) led Zhou et al. to incorrectly claim conclusive evidence of non-thermal effects. We identify flaws in the experimental setup (e.g. the use of the default settings for the thermal camera and incorrect positioning of the thermometer) that may have led Zhou et al. to make such claims. We further show that the Response contains several factual errors and does not address the technical problems we identified with the data acquisition in [Zhou et al., Science 2018]. We demonstrate that both the Response and the original paper contain additional faults, for example, in the power determination and in the normalization of the rate to the catalyst volume, and exhibit misconceptions regarding the thermo-optic response of metal nanostructures. The burden of proof required by the proposal of a novel physical mechanism has simply not been met, especially when the existing data can be modeled exquisitely by conventional theory.


title: Metamorphic Detection of Adversarial Examples in Deep Learning Models With Affine Transformations

abstract: Adversarial attacks are small, carefully crafted perturbations, imperceptible to the naked eye; that when added to an image cause deep learning models to misclassify the image with potentially detrimental outcomes. With the rise of artificial intelligence models in consumer safety and security intensive industries such as self-driving cars, camera surveillance and face recognition, there is a growing need for guarding against adversarial attacks. In this paper, we present an approach that uses metamorphic testing principles to automatically detect such adversarial attacks. The approach can detect image manipulations that are so small, that they are impossible to detect by a human through visual inspection. By applying metamorphic relations based on distance ratio preserving affine image transformations which compare the behavior of the original and transformed image; we show that our proposed approach can determine whether or not the input image is adversarial with a high degree of accuracy.


title: Robust Transmission Network Expansion Planning Problem Considering Storage Units

abstract: This paper addresses the transmission network expansion planning problem considering storage units under uncertain demand and generation capacity. A two-stage adaptive robust optimization framework is adopted whereby short- and long-term uncertainties are accounted for. This work differs from previously reported solutions in an important aspect, namely, we include binary recourse variables to avoid the simultaneous charging and discharging of storage units once uncertainty is revealed. Two-stage robust optimization with discrete recourse problems is a challenging task, so we propose using a nested column-and-constraint generation algorithm to solve the resulting problem. This algorithm guarantees convergence to the global optimum in a finite number of iterations. The performance of the proposed algorithm is illustrated using the Garver's test system.


title: The EL Theorem

abstract: The combined universal probability $\mathbf{m}(D)$ of strings $x$ in sets $D$ is close to max $\mathbf{m}(x)$ over $x$ in $D$: their logs differ by at most $D$'s information $\mathbf{I}(D:\mathcal{H})$ about the halting sequence $\mathcal{H}$.


title: The cohomology rings of homogeneous spaces

abstract: Let $G$ be a compact connected Lie group and $K$ a closed connected subgroup. Assume that the order of any torsion element in the integral cohomology of $G$ and $K$ is invertible in a given principal ideal domain $k$. It is known that in this case the cohomology of the homogeneous space $G/K$ with coefficients in $k$ and the torsion product of $H^{*}(BK)$ and $k$ over $H^{*}(BG)$ are isomorphic as $k$-modules. We show that this isomorphism is multiplicative and natural in the pair $(G,K)$ provided that 2 is invertible in $k$. The proof uses homotopy Gerstenhaber algebras in an essential way. In particular, we show that the normalized singular cochains on the classifying space of a torus are formal as a homotopy Gerstenhaber algebra.


title: Homotopy Gerstenhaber algebras are strongly homotopy commutative

abstract: We show that any homotopy Gerstenhaber algebra is naturally a strongly homotopy commutative (shc) algebra in the sense of Stasheff-Halperin with a homotopy associative structure map. In the presence of certain additional operations corresponding to a cup-1 product on the bar construction, the structure map becomes homotopy commutative, so that one obtains an shc algebra in the sense of Munkholm.


title: Ultracontractive Properties for Directed Graph Semigroups with Applications to Coupled Oscillators

abstract: It is now well known that ultracontractive properties of semigroups with infinitesimal generator given by an undirected graph Laplacian operator can be obtained through an understanding of the geometry of the underlying infinite weighted graph. The aim of this work is to extend these results to semigroups with infinitesimal generator given by a directed graph Laplacian operator, through an analogous inspection of the geometry of the underlying directed graph. In particular, we introduce appropriate nomenclature to discuss the geometry of an infinite directed graph, as well as provide sufficient conditions to extend ultracontractive properties of undirected graph Laplacians to those of the directed variety. Such directed graph Laplacians can often be observed in the study of coupled oscillators, where recent work made explicit the link between synchronous patterns to systems of identically coupled oscillators and ultracontractive properties of undirected graph semigroups. Therefore, in this work we demonstrate the applicability of our results on directed graph semigroups by extending the aforementioned investigation beyond the idealized case of identically coupled oscillators.


title: ReQA: An Evaluation for End-to-End Answer Retrieval Models

abstract: Popular QA benchmarks like SQuAD have driven progress on the task of identifying answer spans within a specific passage, with models now surpassing human performance. However, retrieving relevant answers from a huge corpus of documents is still a challenging problem, and places different requirements on the model architecture. There is growing interest in developing scalable answer retrieval models trained end-to-end, bypassing the typical document retrieval step. In this paper, we introduce Retrieval Question-Answering (ReQA), a benchmark for evaluating large-scale sentence-level answer retrieval models. We establish baselines using both neural encoding models as well as classical information retrieval techniques. We release our evaluation code to encourage further work on this challenging task.


title: Efficient quasiparticle traps with low dissipation through gap engineering

abstract: Quasiparticles represent an intrinsic source of perturbation for superconducting qubits, leading to both dissipation of the qubit energy and dephasing. Recently, it has been shown that normal-metal traps may efficiently reduce the quasiparticle population and improve the qubit lifetime, provided the trap surpasses a certain characteristic size. Moreover, while the trap itself introduces new relaxation mechanisms, they are not expected to harm state-of-the-art transmon qubits under the condition that the traps are not placed too close to extremal positions where electric fields are high. Here, we study a different type of trap, realized through gap engineering. We find that gap-engineered traps relax the remaining constraints imposed on normal metal traps. Firstly, the characteristic trap size, above which the trap is efficient, is reduced with respect to normal metal traps, such that here, strong traps are possible in smaller devices. Secondly, the losses caused by the trap are now greatly reduced, providing more flexibility in trap placement. The latter point is of particular importance, since for efficient protection from quasiparticles, the traps ideally should be placed close to the active parts of the qubit device, where electric fields are typically high.


title: Homotopy Gerstenhaber formality of Davis-Januszkiewicz spaces

abstract: A homotopy Gerstenhaber structure on a differential graded algebra is essentially a family of operations defining a multiplication on its bar construction. We prove that the normalized singular cochain algebra of a Davis-Januszkiewicz space is formal as a homotopy Gerstenhaber algebra, for any coefficient ring. This generalizes a recent result by the author about classifying spaces of tori and also strengthens the well-known dga formality result for Davis-Januszkiewicz spaces due to the author and Notbohm-Ray. As an application, we determine the cohomology rings of free and based loop spaces of Davis-Januszkiewicz spaces.


title: Partition of Optical Properties Into Orbital Contributions

abstract: Nonlinear optical properties (NLOPs) play a major role in photonics, electro-optics and optoelectronics, and other fields of modern optics. The design of new NLO molecules and materials has benefited from the development of computational tools to analyze the relationship between the electronic structure of molecules and its optical response. In this paper, we present a new means to analyze the response property through the partition of NLOPs in terms of orbital contributions (PNOC). This tool can be used to obtain a local representation of the NLOPs, providing a powerful visualization aid to connect the magnitude of the optical property with some parts of the molecule. Unlike other methods to analyze NLOPs, the PNOC decomposes the optical property into orbitals of the unperturbed system, furnishing this method with the ability to assess the performance of single- and multi-determinant electronic structure methods. PNOC can be also used to design small basis sets for an accurate description of large systems, saving a substantial amount of computer time for the calculation of optical properties.


title: AWG-based Nonblocking Shuffle-Exchange Networks

abstract: Optical shuffle-exchange networks (SENs) have wide application in different kinds of interconnection networks. This paper proposes an approach to construct modular optical SENs, using a set of arrayed waveguide gratings (AWGs) and tunable wavelength converters (TWCs). According to the wavelength routing property of AWGs, we demonstrate for the first time that an AWG is functionally equivalent to a classical shuffle network by nature. Based on this result, we devise a systematic method to design a large-scale wavelength-division-multiplexing (WDM) shuffle network using a set of small-size AWGs associated with the same wavelength set. Combining the AWG-based WDM shuffle networks and the TWCs with small conversion range, we finally obtain an AWG-based WDM SEN, which not only is scalable in several ways, but also can achieve 100% utilization when the input wavelength channels are all busy. We also study the routing and wavelength assignment (RWA) problem of the AWG-based WDM SEN, and prove that the self-routing property and the nonblocking routing conditions of classical SENs are preserved in such AWG-based WDM SEN.


title: Symmetry reduction of the three-body problem based on Euler angles

abstract: We consider the classical three-body problem with an arbitrary pair potential which depends on the inter-body distance. A general three-body configuration is set by three "radial" and three angular variables, which determine the shape and orientation, respectively, of a triangle with the three bodies located at the vertices. The radial variables are given by the distances between a reference body and the other two, and by the angle at the reference body between the other two. Such radial variables set the potential energy of the system, and they are reminiscent of the inter-body distance in the two-body problem. On the other hand, the angular variables are the Euler angles relative to a rigid rotation of the triangle, and they are analogous to the polar and azimuthal angle of the vector between the two bodies in the two-body problem. We show that the rotational symmetry allows us to obtain a closed set of eight Hamilton equations of motion, whose generalized coordinates are the thee radial variables and one additional angle, for which we provide the following geometrical interpretation. Given a reference body, we consider the plane through it which is orthogonal to the line between the reference and a second body. We show that the angular variable above is the angle between the plane projection of the angular-momentum vector, and the projection of the radius between the reference and the third body.


title: Fast Haar Transforms for Graph Neural Networks

abstract: Graph Neural Networks (GNNs) have become a topic of intense research recently due to their powerful capability in high-dimensional classification and regression tasks for graph-structured data. However, as GNNs typically define the graph convolution by the orthonormal basis for the graph Laplacian, they suffer from high computational cost when the graph size is large. This paper introduces Haar basis which is a sparse and localized orthonormal system for a coarse-grained chain on graph. The graph convolution under Haar basis, called Haar convolution, can be defined accordingly for GNNs. The sparsity and locality of the Haar basis allow Fast Haar Transforms (FHTs) on graph, by which a fast evaluation of Haar convolution between graph data and filters can be achieved. We conduct experiments on GNNs equipped with Haar convolution, which demonstrates state-of-the-art results on graph-based regression and node classification tasks.


title: Accurate Frequency Domain Identification of ODEs with Arbitrary Signals

abstract: The difficulty in frequency domain identification is that frequency components of arbitrary inputs and outputs are not related by the system's transfer function if signals are windowed. When rectangular windows are used, it is well known that this difference is related to transient effects that can be estimated alongside the systems' parameters windows. In this work, we generalize the approach for arbitrary windows, showing that signal windowing introduces additional terms in the system's equations. The formalism is useful for frequency-domain input-output analysis of a system, and also for system identification. For the latter application, the approach considerably reduces aliasing effects and allows the computation of the associated correction terms, reducing the number of parameters that need to be estimated. The system identification approach has features of the modulating-function technique, filtering out the effects of initial conditions while retaining the spectral interpretation of frequency-domain methods.


title: A Mobile Cloud Collaboration Fall Detection System Based on Ensemble Learning

abstract: Falls are one of the important causes of accidental or unintentional injury death worldwide. Therefore, this paper presents a reliable fall detection algorithm and a mobile cloud collaboration system for fall detection. The algorithm is an ensemble learning method based on decision tree, named Falldetection Ensemble Decision Tree (FEDT). The mobile cloud collaboration system can be divided into three stages: 1) mobile stage: use a light-weighted threshold method to filter out the activities of daily livings (ADLs), 2) collaboration stage: transmit data to cloud and meanwhile extract features in the cloud, 3) cloud stage: deploy the model trained by FEDT to give the final detection result with the extracted features. Experiments show that the performance of the proposed FEDT outperforms the others' over 1-3% both on sensitivity and specificity, and more importantly, the system can provide reliable fall detection in practical scenario.


title: A Survey of Intelligent Reflecting Surfaces (IRSs): Towards 6G Wireless Communication Networks

abstract: Intelligent reflecting surfaces (IRSs) tune wireless environments to increase spectrum and energy efficiencies. In view of much recent attention to the IRS concept as a promising technology for 6G wireless communications, we present a survey of IRSs in this paper. Specifically, we categorize recent research studies of IRSs as follows. For IRS-aided communications, the summary includes capacity/data rate analyses, power/spectral optimizations, channel estimation, deep learning-based design, and reliability analysis. Then we review IRSs implementations as well as the use of IRSs in secure communications, terminal-positioning, and other novel applications. We further identify future research directions for IRSs, with an envision of the IRS technology playing a critical role in 6G communication networks similar to that of massive MIMO in 5G networks. As a timely summary of IRSs, our work will be of interest to both researchers and practitioners working on IRSs for 6G networks.


title: Dynamical mean-field study of Vanadium diselenide monolayer ferromagnetism

abstract: To understand the magnetism of VSe$_2$, whose monolayer form has recently been reported to be a room temperature ferromagnet, Here, the combined method of conventional density functional theory with dynamical mean-field theory has been adopted. This higher-level computation method enables us to resolve many of existing controversies and contradictions in between theory and experiment. First of all, this new approach is shown to give the correct magnetic properties of both bulk and two-dimensional limit of VSe$_2$ which demonstrates its superiority to the conventional methods. The results demonstrate that monolayer VSe$_2$ without charge density waves is a ferromagnet with ordering temperature of 250K. From the direct simulation of temperature-dependent magnetic susceptibility and ordered moment, it is shown that its ferromagnetism is clearly two-dimensional in nature. Further, it is shown that this ferromagnetic order is vulnerable to extra charge dopings which provides the important insight to elucidate recent experimental controversies.


title: The cohomology rings of smooth toric varieties and quotients of moment-angle complexes

abstract: Partial quotients of moment-angle complexes are topological analogues of smooth, not necessarily compact toric varieties. In 1998, Buchstaber and Panov proposed a formula for the cohomology ring of such a partial quotient in terms of a torsion product involving the corresponding Stanley-Reisner ring. We show that their formula gives the correct cup product if 2 is invertible in the chosen coefficient ring, but not in general. We rectify this by defining an explicit deformation of the canonical multiplication on the torsion product.


title: Deformation classes of real Cayley M-octads

abstract: We study 8-point configurations in the real projective space forming an intersection locus of three quadrics and containing no coplanar quadruples. We found that there exists precisely 8 mirror-pairs of deformation classes of such configurations. We describe also the mutual position of these 8 pairs and find the real monodromy groups acting on the 8-point configurations, for each deformation class.


title: Uniform stability of some large-scale parallel server networks

abstract: In this paper we study the uniform stability properties of two classes of parallel server networks with multiple classes of jobs and multiple server pools of a tree topology. These include a class of networks with a single non-leaf server pool, such as the 'N' and 'M' models, and networks of any tree topology with class-dependent service rates. We show that with $\sqrt{n}$ safety staffing, and no abandonment, in the Halfin--Whitt regime, the diffusion-scaled controlled queueing processes are exponentially ergodic and their invariant probability distributions are tight, uniformly over all stationary Markov controls. We use a unified approach in which the same Lyapunov function is used in the study of the prelimit and diffusion limit. A parameter called the spare capacity (safety staffing) of the network plays a central role in characterizing the stability results: the parameter being positive is necessary and sufficient that the limiting diffusion is uniformly exponentially ergodic over all stationary Markov controls. We introduce the concept of "system-wide work conserving policies", which are defined as policies that minimize the number of idle servers at all times. This is stronger than the so-called joint work conservation. We show that, provided the spare capacity parameter is positive, the diffusion-scaled processes are geometrically ergodic and the invariant distributions are tight, uniformly over all "system-wide work conserving policies". In addition, when the spare capacity is negative we show that the diffusion-scaled processes are transient under any stationary Markov control, and when it is zero, they cannot be positive recurrent.


title: Millimeter and far-IR study of the IRCD SDC 341.232-0.268

abstract: We analyze the molecular gas and dust associated with the infrared dark cloud SDC 341.232-0.268 in order to investigate the characteristics and parameters of the gas, determine the evolutionary status of four embedded EGO candidates and establish possible infall or outflow gas motions. We based our study on $^{12}$CO(2-1), $^{13}$CO(2-1), and C$^{18}$O(2-1) data obtained with the APEX telescope, molecular data of high density tracers from the MALT90 survey and IR images from Spitzer, Herschel and ATLASGAL. The study reveals two clumps at -44 km/s towards the IRDC, with densities of > 10$^4$ cm$^{-3}$ , typical of IRDCs, while high density tracers show H$_2$ densities > 10$^5$ . FIR images reveals the presence of cold dust linked to the molecular clumps and EGOs. A comparison of spectra of the optically thin and optically thick molecular lines towards the EGOs suggests the existence of infall and outflow motions.


title: Characterizations of nested GVZ-groups by central series

abstract: Many properties of groups can be defined by the existence of a particular normal series. The classic examples being solvability, supersolvability and nilpotence. Among the nilpotent groups are the so-called nested GVZ-groups --- groups where the centers of the irreducible characters form a chain, and where every irreducible character vanishes off of its center. In this paper, we show that nested GVZ-groups can be characterized by the existence of a certain ascending central series, or by the existence of a certain descending central series.


title: Bayesian Optimization in Variational Latent Spaces with Dynamic Compression

abstract: Data-efficiency is crucial for autonomous robots to adapt to new tasks and environments. In this work we focus on robotics problems with a budget of only 10-20 trials. This is a very challenging setting even for data-efficient approaches like Bayesian optimization (BO), especially when optimizing higher-dimensional controllers. Simulated trajectories can be used to construct informed kernels for BO. However, previous work employed supervised ways of extracting low-dimensional features for these. We propose a model and architecture for a sequential variational autoencoder that embeds the space of simulated trajectories into a lower-dimensional space of latent paths in an unsupervised way. We further compress the search space for BO by reducing exploration in parts of the state space that are undesirable, without requiring explicit constraints on controller parameters. We validate our approach with hardware experiments on a Daisy hexapod robot and an ABB Yumi manipulator. We also present simulation experiments with further comparisons to several baselines on Daisy and two manipulators. Our experiments indicate the proposed trajectory-based kernel with dynamic compression can offer ultra data-efficient optimization.


title: Do Design Metrics Capture Developers Perception of Quality? An Empirical Study on Self-Affirmed Refactoring Activities

abstract: Background. Refactoring is a critical task in software maintenance and is generally performed to enforce the best design and implementation practices or to cope with design defects. Several studies attempted to detect refactoring activities through mining software repositories allowing to collect, analyze and get actionable data-driven insights about refactoring practices within software projects. Aim. We aim at identifying, among the various quality models presented in the literature, the ones that are more in-line with the developer's vision of quality optimization, when they explicitly mention that they are refactoring to improve them. Method. We extract a large corpus of design-related refactoring activities that are applied and documented by developers during their daily changes from 3,795 curated open source Java projects. In particular, we extract a large-scale corpus of structural metrics and anti-pattern enhancement changes, from which we identify 1,245 quality improvement commits with their corresponding refactoring operations, as perceived by software engineers. Thereafter, we empirically analyze the impact of these refactoring operations on a set of common state-of-the-art design quality metrics. Results. The statistical analysis of the obtained results shows that (i) a few state-of-the-art metrics are more popular than others; and (ii) some metrics are being more emphasized than others. Conclusions. We verify that there are a variety of structural metrics that can represent the internal quality attributes with different degrees of improvement and degradation of software quality. Most of the metrics that are mapped to the main quality attributes do capture developer intentions of quality improvement reported in the commit messages.


title: The least signless Laplacian eigenvalue of the complements of bicyclic graphs

abstract: Suppose that $G$ is a connected simple graph with the vertex set $V(G)=\{v_1, v_2,\cdots,v_n\}$. Then the adjacency matrix of $G$ is $A(G)=(a_{ij})_{n\times n}$, where $a_{ij}=1$ if $v_i$ is adjacent to $v_j$, and otherwise $a_{ij}=0$. The degree matrix $D(G)=diag(d_{G}(v_1), d_{G}(v_2), \dots, d_{G}(v_n)),$ where $d_{G}(v_i)$ denotes the degree of $v_i$ in the graph $G$ ($1\leq i\leq n$). The matrix $Q(G)=D(G)+A(G)$ is called the signless Laplacian matrix of $G$. The least eigenvalue of $Q(G)$ is also called the least signless Laplacian eigenvalue of $G$. In this paper we give two graft transformations and then use them to characterize the unique connected graph whose least signless Laplacian eigenvalue is minimum among the complements of all bicyclic graphs.


title: RL-RRT: Kinodynamic Motion Planning via Learning Reachability Estimators from RL Policies

abstract: This paper addresses two challenges facing sampling-based kinodynamic motion planning: a way to identify good candidate states for local transitions and the subsequent computationally intractable steering between these candidate states. Through the combination of sampling-based planning, a Rapidly Exploring Randomized Tree (RRT) and an efficient kinodynamic motion planner through machine learning, we propose an efficient solution to long-range planning for kinodynamic motion planning. First, we use deep reinforcement learning to learn an obstacle-avoiding policy that maps a robot's sensor observations to actions, which is used as a local planner during planning and as a controller during execution. Second, we train a reachability estimator in a supervised manner, which predicts the RL policy's time to reach a state in the presence of obstacles. Lastly, we introduce RL-RRT that uses the RL policy as a local planner, and the reachability estimator as the distance function to bias tree-growth towards promising regions. We evaluate our method on three kinodynamic systems, including physical robot experiments. Results across all three robots tested indicate that RL-RRT outperforms state of the art kinodynamic planners in efficiency, and also provides a shorter path finish time than a steering function free method. The learned local planner policy and accompanying reachability estimator demonstrate transferability to the previously unseen experimental environments, making RL-RRT fast because the expensive computations are replaced with simple neural network inference. Video: https://youtu.be/dDMVMTOI8KY


title: A statistically significant lack of debris discs in medium separation binary systems

abstract: We compile a sample of 341 binary and multiple star systems with the aim of searching for and characterising Kuiper belt-like debris discs. The sample is assembled by combining several smaller samples studied in previously published work with targets from two unpublished Herschel surveys. We find that 38 systems show excess emission at 70 or 100 $\mu$m suggestive of a debris disc. While nine of the discs appear to be unstable to perturbations from their host binary based on a simple analysis of their inferred radii, we argue that the evidence for genuine instability is not strong, primarily because of uncertainty in the true disc radii, uncertainty in the boundaries of the unstable regions, and orbital projection effects. The binary separation distributions of the disc-bearing and disc-free systems are different at a confidence level of $99.4\%$, indicating that binary separation strongly influences the presence of detectable levels of debris. No discs are detected for separations between $\sim$25 and 135 au; this is likely a result of binaries whose separations are comparable with typical disc radii clearing out their primordial circumstellar or circumbinary material via dynamical perturbations. The disc detection rate is $19^{+5}_{-3}\%$ for binaries wider than 135 au, similar to published results for single stars. Only $8^{+2}_{-1}\%$ of systems with separations below 25 au host a detectable disc, which may suggest that planetesimal formation is inhibited in binaries closer than a few tens of au, similar to the conclusions of studies of known planet-hosting binaries.


title: An equilibrium problem on the sphere with two equal charges

abstract: We study the equilibrium measure on the two dimensional sphere in the presence of an external field generated by two equal point charges. The support of the equilibrium measure is known as the droplet. Brauchart et al. showed that the complement of the droplet consists of two spherical caps when the charges are small. When the charges are bigger the droplet becomes simply connected and we prove that the boundary of the droplet is mapped by stereographic projection to an ellipse in the plane. Moreover, we compute a mother body for the droplet that we derive from an equilibrium problem with a weakly admissible external field on the real line.


title: Towards Ultimate Parton Distributions from p+p and e+p Collisions

abstract: We present results from a detailed assessment of the ultimate constraining power of LHC data on the PDFs that can be expected from the complete dataset, in particular after the High-Luminosity (HL) phase. To achieve this, HL-LHC pseudo-data for different projections of the experimental uncertainties are generated, and the resulting constraints on the PDF4LHC15 set are quantified by means of the Hessian profiling method. We find that HL-LHC measurements can reduce PDF uncertainties by up to a factor of 2 to 4 in comparison to state-of-the-art fits, leading to few-percent uncertainties for important observables such as the Higgs boson transverse momentum distribution via gluon-fusion. Our results illustrate the significant improvement in the precision of PDF fits achievable from hadron collider data alone. In addition, we apply the same methodology to the final anticipated data sample from the proposed LHeC, and compare these with the HL-LHC projections, demonstrating an encouraging complementarity between the projected HL-LHC and LHeC constraints.


title: New fifth and seventh order mock theta function identities

abstract: We give simple proofs of Hecke-Rogers indefinite binary theta series identities for the two Ramanujan fifth order mock theta functions $\chi_0(q)$ and $\chi_1(q)$ and all three of Ramanujan's seventh order mock theta functions. We find that the coefficients of the three mock theta functions of order 7 are surprisingly related.


title: The random strategy in Maker-Breaker graph minor games

abstract: In a $(1:b)$ biased Maker-Breaker game, how good a strategy is for a player can be measured by the bias range for which its rival can win, choosing an appropriate counterstrategy. Bednarska and {\L}uczak proved that, in the $H$-subgraph game, the uniformly random strategy for Maker is essentially optimal with high probability. Here we prove an analogous result for the $H$-graph minor game, and we study for which choices of $H$ the random strategy is within a factor of $1+o(1)$ of being optimal.


title: Quantifying Error in the Presence of Confounders for Causal Inference

abstract: Estimating average causal effect (ACE) is useful whenever we want to know the effect of an intervention on a given outcome. In the absence of a randomized experiment, many methods such as stratification and inverse propensity weighting have been proposed to estimate ACE. However, it is hard to know which method is optimal for a given dataset or which hyperparameters to use for a chosen method. To this end, we provide a framework to characterize the loss of a causal inference method against the true ACE, by framing causal inference as a representation learning problem. We show that many popular methods, including back-door methods can be considered as weighting or representation learning algorithms, and provide general error bounds for their causal estimates. In addition, we consider the case when unobserved variables can confound the causal estimate and extend proposed bounds using principles of robust statistics, considering confounding as contamination under the Huber contamination model. These bounds are also estimable; as an example, we provide empirical bounds for the Inverse Propensity Weighting (IPW) estimator and show how the bounds can be used to optimize the threshold of clipping extreme propensity scores. Our work provides a new way to reason about competing estimators, and opens up the potential of deriving new methods by minimizing the proposed error bounds.


title: Environmental Stability of Bismuthene: Oxidation Mechanism and Structural Stability of 2D Pnictogens

abstract: Recently a new group of two dimensional (2D) materials, originating from the group V elements (pnictogens), has gained global attention owing to their outstanding properties.


title: Hacking VMAF with Video Color and Contrast Distortion

abstract: Video quality measurement takes an important role in many applications. Full-reference quality metrics which are usually used in video codecs comparisons are expected to reflect any changes in videos. In this article, we consider different color corrections of compressed videos which increase the values of full-reference metric VMAF and almost don't decrease other widely-used metric SSIM. The proposed video contrast enhancement approach shows the metric inapplicability in some cases for video codecs comparisons, as it may be used for cheating in the comparisons via tuning to improve this metric values.


title: A measure for characterizing heavy-tailed networks

abstract: Heavy-tailed networks are often characterized in the literature by their degree distribution's similarity to a power law. However, many heavy-tailed networks in real life do not have power-law degree distributions, and in many applications the scale-free nature of the network is irrelevant so long as the network possesses hubs. Here we present the Cooke-Nieboer index (CNI), a non-asymptotic measure of the heavy-tailedness of a network's degree distribution which does not presume a power-law form. The CNI is easy to calculate, and clearly distinguishes between networks with power-law, exponential, and symmetric degree distributions.


title: Variational Autoencoders and Nonlinear ICA: A Unifying Framework

abstract: The framework of variational autoencoders allows us to efficiently learn deep latent-variable models, such that the model's marginal distribution over observed variables fits the data. Often, we're interested in going a step further, and want to approximate the true joint distribution over observed and latent variables, including the true prior and posterior distributions over latent variables. This is known to be generally impossible due to unidentifiability of the model. We address this issue by showing that for a broad family of deep latent-variable models, identification of the true joint distribution over observed and latent variables is actually possible up to very simple transformations, thus achieving a principled and powerful form of disentanglement. Our result requires a factorized prior distribution over the latent variables that is conditioned on an additionally observed variable, such as a class label or almost any other observation. We build on recent developments in nonlinear ICA, which we extend to the case with noisy, undercomplete or discrete observations, integrated in a maximum likelihood framework. The result also trivially contains identifiable flow-based generative models as a special case.


title: Cake cutting: Explicit examples for impossibility results

abstract: In this article we suggest a model of computation for the cake cutting problem. In this model the mediator can ask the same queries as in the Robertson-Webb model but he or she can only perform algebraic operations as in the Blum-Shub-Smale model. All existing algorithms described in the Robertson-Webb model can be described in this new model.We show that in this model there exist explicit couples of measures for which no algorithm outputs an equitable fair division with connected parts.We also show that there exist explicit set of measures for which no algorithm in this model outputs a fair division which maximizes the utilitarian social welfare function.The main tool of our approach is Galois theory.


title: ABSense: Sensing Electromagnetic Waves on Metasurfaces via Ambient Compilation of Full Absorption

abstract: Metasurfaces constitute effective media for manipulating and transforming impinging EM waves. Related studies have explored a series of impactful MS capabilities and applications in sectors such as wireless communications, medical imaging and energy harvesting. A key-gap in the existing body of work is that the attributes of the EM waves to-be-controlled (e.g., direction, polarity, phase) are known in advance. The present work proposes a practical solution to the EM wave sensing problem using the intelligent and networked MS counterparts-the HyperSurfaces (HSFs), without requiring dedicated field sensors. An nano-network embedded within the HSF iterates over the possible MS configurations, finding the one that fully absorbs the impinging EM wave, hence maximizing the energy distribution within the HSF. Using a distributed consensus approach, the nano-network then matches the found configuration to the most probable EM wave traits, via a static lookup table that can be created during the HSF manufacturing. Realistic simulations demonstrate the potential of the proposed scheme. Moreover, we show that the proposed workflow is the first-of-its-kind embedded EM compiler, i.e., an autonomic HSF that can translate high-level EM behavior objectives to the corresponding, low-level EM actuation commands.


title: Mobile phone data's potential for informing infrastructure planning in developing countries

abstract: High quality census data are not always available in developing countries. Instead, mobile phone data are becoming a trending proxy to evaluate population density, activity and social characteristics. They offer additional advantages for infrastructure planning such as being updated in real-time, including mobility information and recording visitors' activity. We combine various data sets from Senegal to evaluate mobile phone data's potential to replace insufficient census data for infrastructure planning in developing countries. As an applied case, we test their ability at predicting domestic electricity consumption. We show that, contrary to common belief, average mobile phone activity does not correlate well with population density. However, it can provide better electricity consumption estimates than basic census data. More importantly, we use curve and network clustering techniques to enhance the accuracy of the predictions, to recover good population mapping potential and to reduce the collection of required data to substantially smaller samples.


title: Helical Antiferromagnetic Ordering in EuNi{1.95}As2

abstract: The Eu{+2} spins-7/2 in EuNi2As2 with the body-centered tetragonal ThCr2Si2 structure order antiferromagnetically below the Neel temperature TN =15 K into a helical antiferromagnetic (AFM) structure with the helix axis aligned along the tetragonal c axis and the Eu ordered moments aligned ferromagnetically within the ab plane as previously reported from neutron diffraction measurements [T. Jin, et al., Phys. Rev. B 99, 014425 (2019)]. Here we study the crystallographic, magnetic, thermal, and electronic transport properties of Bi-flux-grown single crystals using single-crystal x-ray diffraction, anisotropic magnetic susceptibility chi, isothermal magnetization M, heat capacity Cp, and electrical resistivity rho measurements versus applied magnetic field H and temperature T. Vacancies are found on the Ni sites corresponding to the composition EuNi{1.95}As2. The chi{ab}(T) data below TN are fitted well by molecular field theory (MFT), and the helix turn angle kd and the Eu-Eu Heisenberg exchange constants are extracted from the fit parameters. The kd value is in good agreement with the neutron-diffraction result. The magnetic contribution to the zero-field heat capacity below TN is also fitted by MFT. The isothermal in-plane magnetization Mab exhibits two metamagnetic transitions versus H, whereas Mc(T = 2 K) is nearly linear up to H =14 T, both behaviors being consistent with MFT. The Mc(H,T), rho(Hc,T), and Cp(Hc,T) data yielded a Hc-T phase diagram separating the AFM and paramagnetic phases in good agreement with MFT. Anisotropic chi(T) literature data for the ThCr2Si2-type helical antiferromagnet EuRh2As2 are also fitted well by MFT. A comparison is made between the crystallographic and magnetic properties of ThCr2Si2- type EuM2Pn2 compounds with M = Fe, Co, Ni, Cu, or Rh, and Pn = P or As, where only ferromagnetic and c-axis helical AFM structures are found.


title: Discrepancy of minimal Riesz energy points

abstract: We find upper bounds for the spherical cap discrepancy of the set of minimizers of the Riesz $s$-energy on the sphere $\mathbb S^d.$ Our results are based in bounds for a Sobolev discrepancy introduced by Thomas Wolff in an unpublished manuscript where estimates for the spherical cap discrepancy of the logarithmic energy minimizers in $\mathbb S^2$ were obtained. Our result improves previously known bounds for $0\le s<2$ and $s\neq 1$ in $\mathbb S^2,$ where $s=0$ is Wolff's result, and for $d-t_0<s<d$ with $t_0\approx 2.5$ when $d\ge 3$ and $s\neq d-1.$


title: The evolution of electron dispersion in the series of rare-earth tritelluride compounds obtained from their charge-density-wave properties and susceptibility calculations

abstract: We calculated electron susceptibility of rare-earth tritelluride compounds RTe$_3$ as a function of temperature, wave vector and electron-dispersion parameters. Comparison of results obtained with the available experimental data on the transition temperature and on the wave vector of a charge-density wave in these compounds allowed us to predict values and the evolution of electron-dispersion parameters with the variation of atomic number of rare-earth element R.


title: A Data-Driven Approach for Accurate Rainfall Prediction

abstract: In recent years, there has been growing interest in using Precipitable Water Vapor (PWV) derived from Global Positioning System (GPS) signal delays to predict rainfall. However, the occurrence of rainfall is dependent on a myriad of atmospheric parameters. This paper proposes a systematic approach to analyze various parameters that affect precipitation in the atmosphere. Different ground-based weather features like Temperature, Relative Humidity, Dew Point, Solar Radiation, PWV along with Seasonal and Diurnal variables are identified, and a detailed feature correlation study is presented. While all features play a significant role in rainfall classification, only a few of them, such as PWV, Solar Radiation, Seasonal and Diurnal features, stand out for rainfall prediction. Based on these findings, an optimum set of features are used in a data-driven machine learning algorithm for rainfall prediction. The experimental evaluation using a four-year (2012-2015) database shows a true detection rate of 80.4%, a false alarm rate of 20.3%, and an overall accuracy of 79.6%. Compared to the existing literature, our method significantly reduces the false alarm rates.


title: Area Law Behaviour of Mutual Information at Finite Temperature

abstract: Entanglement entropy in free scalar field theory at its ground state is dominated by an area law term. However, when mixed states are considered this property ceases to exist. We show that in such cases the mutual information obeys an "area law". The proportionality constant connecting the area to the mutual information has an interesting dependence on the temperature. At infinite temperature it tends to a finite value which coincides with the classical calculation.


title: Time-domain characterization and correction of on-chip distortion of control pulses in a quantum processor

abstract: We introduce Cryoscope, a method for sampling on-chip baseband pulses used to dynamically control qubit frequency in a quantum processor. We specifically use Cryoscope to measure the step response of the dedicated flux control lines of two-junction transmon qubits in circuit QED processors with the temporal resolution of the room-temperature arbitrary waveform generator producing the control pulses. As a first application, we iteratively improve this step response using optimized real-time digital filters to counter the linear-dynamical distortion in the control line, as needed for high-fidelity, repeatable one- and two-qubit gates based on dynamical control of qubit frequency.


title: Hydrodynamically interrupted droplet growth in scalar active matter

abstract: Suspensions of spherical active particles often show microphase separation. At a continuum level, coupling their scalar density to fluid flow, there are two distinct explanations. Each involves an effective interfacial tension: the first mechanical (causing flow) and the second diffusive (causing Ostwald ripening). Here we show how the negative mechanical tension of contractile swimmers creates, via a self-shearing instability, a steady-state life cycle of droplet growth interrupted by division whose scaling behavior we predict. When the diffusive tension is also negative, this is replaced by an arrested regime (mechanistically distinct, but with similar scaling) where division of small droplets is prevented by reverse Ostwald ripening.


title: Optimizing defence, counter-defence and counter-counter defence in parasitic and trophic interactions -- A modelling study

abstract: In host-pathogen interactions, often the host (attacked organism) defends itself by some toxic compound and the parasite, in turn, responds by producing an enzyme that inactivates that compound. In some cases, the host can respond by producing an inhibitor of that enzyme, which can be considered as a counter-counter defence. An example is provided by cephalosporins, beta-lactamases and clavulanic acid (an inhibitor of beta-lactamases). Here, we tackle the question under which conditions it pays, during evolution, to establish a counter-counter defence rather than to intensify or widen the defence mechanisms. We establish a mathematical model describing this phenomenon, based on enzyme kinetics for competitive inhibition. We use an objective function based on Haber's rule, which says that the toxic effect is proportional to the time integral of toxin concentration. The optimal allocation of defence and counter-counter defence can be calculated in an analytical way despite the nonlinearity in the underlying differential equation. The calculation provides a threshold value for the dissociation constant of the inhibitor. Only if the inhibition constant is below that threshold, that is, in the case of strong binding of the inhibitor, it pays to have a counter-counter defence. This theoretical prediction accounts for the observation that not for all defence mechanisms, a counter-counter defence exists. Our results should be of interest for computing optimal mixtures of beta-lactam antibiotics and beta-lactamase inhibitors such as sulbactam, as well as for plant-herbivore and other molecular-ecological interactions and to fight antibiotic resistance in general.


title: On Characters of a Class of P-polynomial table algebras and applications

abstract: In this paper, we study the characters of homogeneous monotonic P-polynomial table algebras with finite dimension d>=5. We then apply them to association schemes. To this end, we develop some methods using tridiagonal matrices and Z-transform. Moreover, we calculate the eigenvalues of a special tridiagonal matrix which is found through the first intersection matrix of P-polynomial table algebras.


title: Improving Prognostic Performance in Resectable Pancreatic Ductal Adenocarcinoma using Radiomics and Deep Learning Features Fusion in CT Images

abstract: As an analytic pipeline for quantitative imaging feature extraction and analysis, radiomics has grown rapidly in the past a few years. Recent studies in radiomics aim to investigate the relationship between tumors imaging features and clinical outcomes. Open source radiomics feature banks enable the extraction and analysis of thousands of predefined features. On the other hand, recent advances in deep learning have shown significant potential in the quantitative medical imaging field, raising the research question of whether predefined radiomics features have predictive information in addition to deep learning features. In this study, we propose a feature fusion method and investigate whether a combined feature bank of deep learning and predefined radiomics features can improve the prognostics performance. CT images from resectable Pancreatic Adenocarcinoma (PDAC) patients were used to compare the prognosis performance of common feature reduction and fusion methods and the proposed risk-score based feature fusion method for overall survival. It was shown that the proposed feature fusion method significantly improves the prognosis performance for overall survival in resectable PDAC cohorts, elevating the area under ROC curve by 51% compared to predefined radiomics features alone, by 16% compared to deep learning features alone, and by 32% compared to existing feature fusion and reduction methods for a combination of deep learning and predefined radiomics features.


title: Efficient simulation of filament elastohydrodynamics in three dimensions

abstract: Fluid-structure simulations of slender inextensible filaments in a viscous fluid are often plagued by numerical stiffness. Recent coarse-graining studies have reduced the computational requirements of simulating such systems, though have thus far been limited to the motion of planar filaments. In this work we extend such frameworks to filament motion in three dimensions, identifying and circumventing coordinate-system singularities introduced by filament parameterisation via repeated changes of basis. The resulting methodology enables efficient and rapid study of the motion of flexible filaments in three dimensions, and is readily extensible to a wide range of problems, including filament motion in confined geometries, large-scale active matter simulations, and the motility of mammalian spermatozoa.


title: Scheduling With Inexact Job Sizes: The Merits of Shortest Processing Time First

abstract: It is well known that size-based scheduling policies, which take into account job size (i.e., the time it takes to run them), can perform very desirably in terms of both response time and fairness. Unfortunately, the requirement of knowing a priori the exact job size is a major obstacle which is frequently insurmountable in practice. Often, it is possible to get a coarse estimation of job size, but unfortunately analytical results with inexact job sizes are challenging to obtain, and simulation-based studies show that several size-based algorithm are severely impacted by job estimation errors. For example, Shortest Remaining Processing Time (SRPT), which yields optimal mean sojourn time when job sizes are known exactly, can drastically underperform when it is fed inexact job sizes. Some algorithms have been proposed to better handle size estimation errors, but they are somewhat complex and this makes their analysis challenging. We consider Shortest Processing Time (SPT), a simplification of SRPT that skips the update of "remaining" job size and results in a preemptive algorithm that simply schedules the job with the shortest estimated processing time. When job size is inexact, SPT performs comparably to the best known algorithms in the presence of errors, while being definitely simpler. In this work, SPT is evaluated through simulation, showing near-optimal performance in many cases, with the hope that its simplicity can open the way to analytical evaluation even when inexact inputs are considered.


title: Superdiffusive limits for deterministic fast-slow dynamical systems

abstract: We consider deterministic fast-slow dynamical systems on $\mathbb{R}^m\times Y$ of the form \[ \begin{cases} x_{k+1}^{(n)} = x_k^{(n)} + n^{-1} a(x_k^{(n)}) + n^{-1/\alpha} b(x_k^{(n)}) v(y_k)\;,\quad y_{k+1} = f(y_k)\;, \end{cases} \] where $\alpha\in(1,2)$. Under certain assumptions we prove convergence of the $m$-dimensional process $X_n(t)= x_{\lfloor nt \rfloor}^{(n)}$ to the solution of the stochastic differential equation \[ \mathop{}\!\mathrm{d} X = a(X)\mathop{}\!\mathrm{d} t + b(X) \diamond \mathop{}\!\mathrm{d} L_\alpha \; , \] where $L_\alpha$ is an $\alpha$-stable L\'evy process and $\diamond$ indicates that the stochastic integral is in the Marcus sense. In addition, we show that our assumptions are satisfied for intermittent maps $f$ of Pomeau-Manneville type.