title: Effective low-energy description of the two impurity Anderson model: RKKY interaction and quantum criticality

abstract: We show that the RKKY interaction in the two-impurity Anderson model comprise two contributions: a ferromagnetic part stemming from the symmetrized hybridization functions and an anti-ferromagnetic part. We demonstrate that this anti-ferromagnetic contribution can also be generated by an effective local tunneling term between the two impurities. This tunneling can be analytically calculated for particle-hole symmetric impurities. Replacing the full hybridization functions by the symmetric part and this tunneling term leads to the identical low-temperature fixed point spectrum in the numerical renormalization group. Compensating this tunneling term is used to restore the Varma-Jones quantum critical point between a strong coupling phase and a local singlet phase even in the absence of particle-hole symmetry in the hybridization functions. We analytically investigate the spatial frequencies of the effective tunneling term based on the combination of the band dispersion and the shape of the Fermi surface. Numerical renormalization group calculations provide a comparison of the distance dependent tunneling term and the local spin-spin correlation function. Derivations between the spatial dependency of the full spin-spin correlation function and the textbook RKKY interaction are reported.


title: High-Efficiency Ultra-Violet Dielectric Meta-Holograms with Antiferromagnetic Resonances

abstract: Metasurfaces with spatially varying subwavelength structures enable full control of electromagnetic waves over a wide spectrum. High-efficiency metasurfaces, especially in a transmission mode, are of practical significance in optical elements and systems, hitherto their operating frequencies have been expanded down to visible-wavelength ranges. Challenges of developing shorter-wavelength metasurfaces originate from electromagnetic loss caused by strong absorption for most high-refractive-index materials. Here we introduce a large-bandgap semiconductor material-niobium pentoxide (Nb2O5)-to engineer a ultraviolet meta-hologram with a total efficiency of 81% at 355nm wavelength. This meta-hologram modulates the geometric phase of transmitted circular-polarization light via orientation-varying high-aspect-ratio nano-bricks that are elaborately designed to excite antiferromagnetic resonances. We reveal that the induced antiferromagnetic modes maintain one component (e.g., Ex-component) of incident light through even-numbered antiparallel magnetic dipoles (AMDs) but reverse the other component (e.g., Ey-component) via odd-numbered AMDs, thus realizing high-efficiency polarization conversion. Our approach might open the door towards high-performance ultraviolet-light nanophotonics and meta-optics.


title: Abundance of zinc in the red giants of Galactic globular cluster 47 Tucanae

abstract: We investigate possible relations between the abundances of zinc and the light elements sodium, magnesium, and potassium in the atmospheres ofred giant branch (RGB) stars of the Galactic globular cluster 47~Tuc and study connections between the chemical composition and dynamical properties of the cluster RGB stars. The abundance of zinc was determined in 27 RGB stars of 47~Tuc using 1D local thermal equilibrium (LTE) synthetic line profile fitting to the high-resolution 2dF spectra obtained with the Anglo-Australian Telescope (AAT). Synthetic spectra used in the fitting procedure were computed with the SYNTHE code and 1D ATLAS9 stellar model atmospheres. The average 1D~LTE zinc-to-iron abundance ratio and its RMS variations due to star-to-star abundance spread determined in the sample of 27 RGB stars is $\langle{\rm[Zn/Fe]}\rangle^{\rm 1D~LTE}=0.11\pm0.09$. We did not detect any statistically significant relations between the abundances of zinc and those of light elements. Neither did we find any significant correlation or anticorrelation between the zinc abundance in individual stars and their projected distance from the cluster center. Finally, no statistically significant relation between the absolute radial velocities of individual stars and the abundance of zinc in their atmospheres was detected. The determined average [Zn/Fe]$^{\rm~1DLTE}$ ratio agrees well with those determined in this cluster in earlier studies and nearly coincides with that of Galactic field stars at this metallicity. All these results suggest that nucleosynthesis of zinc and light elements proceeded in separate, unrelated pathways in 47~Tuc.


title: Periodic P\'olya urns and an application to Young tableaux

abstract: P{\'o}lya urns are urns where at each unit of time a ball is drawn and is replaced with some other balls according to its colour. We introduce a more general model: The replacement rule depends on the colour of the drawn ball and the value of the time (mod p). We discuss some intriguing properties of the differential operators associated to the generating functions encoding the evolution of these urns. The initial partial differential equation indeed leads to ordinary linear differential equations and we prove that the moment generating functions are D-finite. For a subclass, we exhibit a closed form for the corresponding generating functions (giving the exact state of the urns at time n). When the time goes to infinity, we show that these periodic P{\'o}lya urns follow a rich variety of behaviours: their asymptotic fluctuations are described by a family of distributions, the generalized Gamma distributions, which can also be seen as powers of Gamma distributions. En passant, we establish some enumerative links with other combinatorial objects, and we give an application for a new result on the asymptotics of Young tableaux: This approach allows us to prove that the law of the lower right corner in a triangular Young tableau follows asymptotically a product of generalized Gamma distributions.


title: On Minimal Sets to Destroy the $k$-Core in Random Networks

abstract: We study the problem of finding the smallest set of nodes in a network whose removal results in an empty $k$-core; where the $k$-core is the sub-network obtained after the iterative removal of all nodes of degree smaller than $k$. This problem is also known in the literature as finding the minimal contagious set. The main contribution of our work is an analysis of the performance of the recently introduced corehd algorithm [Scientific Reports, 6, 37954 (2016)] on random networks taken from the configuration model via a set of deterministic differential equations. Our analyses provides upper bounds on the size of the minimal contagious set that improve over previously known bounds. Our second contribution is a new heuristic called the weak-neighbor algorithm that outperforms all currently known local methods in the regimes considered.

categories: physics.soc-ph cond-mat.dis-nn cs.SI math.PR


title: Semi-parametric estimation of the variogram of a Gaussian process with stationary increments

abstract: We consider the semi-parametric estimation of a scale parameter of a one-dimensional Gaussian process with known smoothness. We suggest an estimator based on quadratic variations and on the moment method. We provide asymptotic approximations of the mean and variance of this estimator, together with asymptotic normality results, for a large class of Gaussian processes. We allow for general mean functions and study the aggregation of several estimators based on various variation sequences. In extensive simulation studies, we show that the asymptotic results accurately depict thefinite-sample situations already for small to moderate sample sizes. We also compare various variation sequences and highlight the efficiency of the aggregation procedure.


title: An upper-limit on the linear polarization fraction of the GW170817 radio continuum

abstract: We present late-time radio observations of GW170817, the first binary neutron star merger discovered through gravitational waves by the advanced LIGO and Virgo detectors. Our observations, carried out with the Karl G. Jansky Very Large Array, were optimized to detect polarized radio emission, and thus to constrain the linear polarization fraction of GW170817. At an epoch of ~244 days after the merger, we rule out linearly polarized emission above a fraction of ~12% at a frequency of 2.8 GHz (99% confidence). Within the structured jet scenario (a.k.a. successful jet plus cocoon system) for GW170817, the derived upper-limit on the radio continuum linear polarization fraction strongly constrains the magnetic field configuration in the shocked ejecta. We show that our results for GW170817 are compatible with the low level of linear polarization found in afterglows of cosmological long gamma-ray bursts. Finally, we discuss our findings in the context of future expectations for the study of radio counterparts of binary neutron star mergers identified by ground-based gravitational-wave detectors.


title: Annihilation of tor\_p(G\_S^ab) for real abelian extensions

abstract: Preprint of a paper to appear in "Communications in Advanced Mathematical Sciences". Let K be a real abelian extension of Q. Let p be a prime number, S the set of p-places of K and G\_K,S the Galois group of the maximal S-ramified pro-p-extension of K (i.e., unramified outside p and infinity). We revisit the problem of annihilation of the p-torsion group T\_K:=tor\_Z\_p(G\_K,S^ab) initiated by us and Oriat then systematized in our paper on the construction of p-adic L-functions in which we obtained a canonical ideal annihilator of T\_K in full generality (1978--1981). Afterwards (1992--2014) some annihilators, using cyclotomic units, were proposed by Solomon, Belliard--Nguyen Quang Do, Nguyen Quang Do--Nicolas, All, Belliard--Martin.In this text, we improve our original papers and show that, in general, the Solomon elements are not optimal and/or partly degenerated. We obtain, whatever K and p, an universal non-degenerated annihilator in terms of p-adic logarithms of cyclotomic numbers related to L\_p-functions at s=1 of its primitive characters of K (Theorem 9.4). Some computations are given with PARI programs; the case p=2 is analyzed and illustrated in degrees 2, 3, 4 to test a conjecture.


title: The T-dual symmetries of a bosonic string

abstract: We investigate whether the symmetry transformations of a bosonic string are connected by T-duality. We start with a standard closed string theory. We continue with a modified open string theory, modified to preserve the symmetry transformations possessed by the closed string theory. Because the string theory is conformally invariant world sheet field theory, in order to find the transformations which preserve the physics, one has to demand the isomorphism between the conformal field theories corresponding to the initial and the transformed field configurations. We find the symmetry transformations corresponding to the similarity transformation of the energy-momentum tensor, and find that their generators are T-dual. Particularly, we find that the general coordinate and local gauge transformations are T-dual, so we conclude that T-duality in addition to the well known exchanges, transforms symmetries of the initial and its T-dual theory into each other.


title: Single-photon characteristics of superposed weak coherent states

abstract: We study a superposed weak coherent state that can fundamentally mimic an ideal single photon not only with respect to the number of photons but also in terms of an indeterminate phase. It is close to the single-photon state with high fidelity and exhibits fundamental features of single photons such as antibunching and Hong-Ou-Mandel interference. The emergence and vanishing of single-photon characteristics can be directly observed by changing two parameters, i.e., the mean photon number and number of phases. Our result shows that the uncertainty between the photon number and phase indeed constitutes the characteristics of single photons. Finally, we apply the superposed weak coherent state to quantum key distribution and demonstrate that it outperforms the typical approach using phase-randomized weak coherent states.


title: The invasion speed of cell migration models with realistic cell cycle time distributions

abstract: Cell proliferation is typically incorporated into stochastic mathematical models of cell migration by assuming that cell divisions occur after an exponentially distributed waiting time. Experimental observations, however, show that this assumption is often far from the real cell cycle time distribution (CCTD). Recent studies have suggested an alternative approach to modelling cell proliferation based on a multi-stage representation of the CCTD. In order to validate and parametrise these models, it is important to connect them to experimentally measurable quantities. In this paper we investigate the connection between the CCTD and the speed of the collective invasion. We first state a result for a general CCTD, which allows the computation of the invasion speed using the Laplace transform of the CCTD. We use this to deduce the range of speeds for the general case. We then focus on the more realistic case of multi-stage models, using both a stochastic agent-based model and a set of reaction-diffusion equations for the cells' average density. By studying the corresponding travelling wave solutions, we obtain an analytical expression for the speed of invasion for a general N-stage model with identical transition rates, in which case the resulting cell cycle times are Erlang distributed. We show that, for a general N-stage model, the Erlang distribution and the exponential distribution lead to the minimum and maximum invasion speed, respectively. This result allows us to determine the range of possible invasion speeds in terms of the average proliferation time for any multi-stage model.


title: Parameter free quantitative analysis of atom probe data by correlation functions: Application to the precipitation in Al-Zn-Mg-Cu

abstract: Atom probe tomography enables precise quantification of the composition of second phase particles from their early stages, leading to improved understanding of the thermodynamic and kinetic mechanisms of phase formation and quantify structure-property relationships. Here we demonstrate how approaches developed for small-angle scattering can be adapted to atom probe tomography. By exploiting nearest-neighbor distributions and radial distribution function, we introduce a parameter free methodology to efficiently extract information such as particle size, composition, volume fraction, number density and inter-particle distance. We demonstrate the strength of this approach in the analysis of a precipitation-hardened model Al-Zn-Mg-Cu high-strength lightweight alloy.


title: CLEAR: Coverage-based Limiting-cell Experiment Analysis for RNA-seq

abstract: Direct cDNA preamplification protocols developed for single-cell RNA-seq have enabled transcriptome profiling of precious clinical samples and rare cells without sample pooling or RNA extraction. Currently, there is no algorithm optimized to reveal and remove noisy transcripts in limiting-cell RNA-seq (lcRNA-seq) data for downstream analyses. Herein, we present CLEAR, a workflow that identifies reliably quantifiable transcripts in lcRNA-seq data for differentially expressed gene (DEG) analysis. Libraries at three input amounts of FACS-derived CD5+ and CD5- cells from a chronic lymphocytic leukemia patient were used to develop CLEAR. When using CLEAR transcripts vs. using all transcripts, downstream analyses revealed more shared transcripts across different input RNA amounts, improved Principal Component Analysis (PCA) separation, and yielded more DEGs between cell types. As proof-of-principle, CLEAR was applied to an in-house lcRNA-seq dataset and two public datasets. When imputation is used, CLEAR is also adaptable to large clinical studies and for single cell analyses.


title: Black Box FDR

abstract: Analyzing large-scale, multi-experiment studies requires scientists to test each experimental outcome for statistical significance and then assess the results as a whole. We present Black Box FDR (BB-FDR), an empirical-Bayes method for analyzing multi-experiment studies when many covariates are gathered per experiment. BB-FDR learns a series of black box predictive models to boost power and control the false discovery rate (FDR) at two stages of study analysis. In Stage 1, it uses a deep neural network prior to report which experiments yielded significant outcomes. In Stage 2, a separate black box model of each covariate is used to select features that have significant predictive power across all experiments. In benchmarks, BB-FDR outperforms competing state-of-the-art methods in both stages of analysis. We apply BB-FDR to two real studies on cancer drug efficacy. For both studies, BB-FDR increases the proportion of significant outcomes discovered and selects variables that reveal key genomic drivers of drug sensitivity and resistance in cancer.


title: Automatic Identification of Research Fields in Scientific Papers

abstract: The TERRE-ISTEX project aims to identify scientific research dealing with specific geographical territories areas based on heterogeneous digital content available in scientific papers. The project is divided into three main work packages: (1) identification of the periods and places of empirical studies, and which reflect the publications resulting from the analyzed text samples, (2) identification of the themes which appear in these documents, and (3) development of a web-based geographical information retrieval tool (GIR). The first two actions combine Natural Language Processing patterns with text mining methods. The integration of the spatial, thematic and temporal dimensions in a GIR contributes to a better understanding of what kind of research has been carried out, of its topics and its geographical and historical coverage. Another originality of the TERRE-ISTEX project is the heterogeneous character of the corpus, including PhD theses and scientific articles from the ISTEX digital libraries and the CIRAD research center.


title: Fidelity-based Probabilistic Q-learning for Control of Quantum Systems

abstract: The balance between exploration and exploitation is a key problem for reinforcement learning methods, especially for Q-learning. In this paper, a fidelity-based probabilistic Q-learning (FPQL) approach is presented to naturally solve this problem and applied for learning control of quantum systems. In this approach, fidelity is adopted to help direct the learning process and the probability of each action to be selected at a certain state is updated iteratively along with the learning process, which leads to a natural exploration strategy instead of a pointed one with configured parameters. A probabilistic Q-learning (PQL) algorithm is first presented to demonstrate the basic idea of probabilistic action selection. Then the FPQL algorithm is presented for learning control of quantum systems. Two examples (a spin- 1/2 system and a lamda-type atomic system) are demonstrated to test the performance of the FPQL algorithm. The results show that FPQL algorithms attain a better balance between exploration and exploitation, and can also avoid local optimal policies and accelerate the learning process.


title: Neural Message Passing with Edge Updates for Predicting Properties of Molecules and Materials

abstract: Neural message passing on molecular graphs is one of the most promising methods for predicting formation energy and other properties of molecules and materials. In this work we extend the neural message passing model with an edge update network which allows the information exchanged between atoms to depend on the hidden state of the receiving atom. We benchmark the proposed model on three publicly available datasets (QM9, The Materials Project and OQMD) and show that the proposed model yields superior prediction of formation energies and other properties on all three datasets in comparison with the best published results. Furthermore we investigate different methods for constructing the graph used to represent crystalline structures and we find that using a graph based on K-nearest neighbors achieves better prediction accuracy than using maximum distance cutoff or the Voronoi tessellation graph.


title: A direct linear inversion for discontinuous elastic parameters recovery from internal displacement information only

abstract: The aim of this paper is to present and analyze a new direct method for solving the linear elasticity inverse problem. Given measurements of some displacement fields inside a medium, we show that a stable reconstruction of elastic parameters is possible, even for discontinuous parameters and without boundary information. We provide a general approach based on the weak definition of the stiffness-to-force operator which conduces to see the problem as a linear system. We prove that in the case of shear modulus reconstruction, we have an $L^2$-stability with only one measurement under minimal smoothness assumptions. This stability result is obtained though the proof that the linear operator to invert has closed range. We then describe a direct discretization which provides stable reconstructions of both isotropic and anisotropic stiffness tensors.


title: Real-Space Magnetic Imaging of the Multiferroic Spinels MnV2O4 and Mn3O4

abstract: Controlling multiferroic behavior in materials will enable the development of a wide variety of technological applications. However, the exact mechanisms driving multiferroic behavior are not well understood in most materials. Two such materials are the spinels MnV2O4 and Mn3O4, where mechanical strain is thought to play a role in determining magnetic behavior. Bulk studies of MnV2O4 have yielded conflicting and inconclusive results, due in part to the presence of mesoscale magnetic inhomogeneity, which complicates the interpretation of bulk measurements. To study the sub-micron-scale magnetic properties of Mn-based spinel materials, we performed magnetic force microscopy (MFM) on MnV2O4 samples subject to different levels of mechanical strain. We also used a crystal grain mapping technique to perform spatially registered MFM on Mn3O4. These local investigations revealed 100-nm-scale "stripe" modulations in the magnetic structure of both materials. In MnV2O4, the magnetization of these stripes is estimated to be Mz $\approx$ 105 A/m, which is on the order of the saturation magnetization reported previously. Cooling in a strong magnetic field eliminated the stripe patterning only in the low-strain sample of MnV2O4. The discovery of nanoscale magnetostructural inhomogeneity that is highly susceptible to magnetic field control in these materials necessitates both a revision of theoretical proposals and a reinterpretation of experimental data regarding the low-temperature phases and magnetic-field-tunable properties of these Mn-based spinels.


title: Several recent developments in estimation and robust control of quantum systems

abstract: This paper summarizes several recent developments in the area of estimation and robust control of quantum systems and outlines several directions for future research. Quantum state tomography via linear regression estimation and adaptive quantum state estimation are introduced and a Hamiltonian identification algorithm is outlined. Two quantum robust control approaches including sliding mode control and sampling-based learning control are illustrated.


title: A User Evaluation of Automated Process Discovery Algorithms

abstract: Process mining methods allow analysts to use logs of historical executions of business processes in order to gain knowledge about the actual behavior of these processes. One of the most widely studied process mining operations is automated process discovery. An event log is taken as input by an automated process discovery method and produces a business process model as output that captures the control-flow relations between tasks that are described by the event log. In this setting, this paper provides a systematic comparative evaluation of existing implementations of automated process discovery methods with domain experts by using a real-life event log extracted from an international software engineering company and four quality metrics. The evaluation results highlight gaps and unexplored trade-offs in the field and allow researchers to improve the lacks in the automated process discovery methods in terms of usability of process discovery techniques in industry.


title: A superconducting absolute spin valve

abstract: A superconductor with a spin-split excitation spectrum behaves as an ideal ferromagnetic spin-injector in a tunneling junction. It was theoretical predicted that the combination of two such spin-split superconductors with independently tunable magnetizations, may be used as an ideal $absolute$ spin-valve. Here we report on the first switchable superconducting spin-valve based on two EuS/Al bilayers coupled through an aluminum oxide tunnel barrier. The spin-valve shows a relative resistance change between the parallel and antiparallel configuration of the EuS layers up to 900% that demonstrates a highly spin-polarized currents through the junction. Our device may be pivotal for realization of thermoelectric radiation detectors, logical element for a memory cell in cryogenics superconductor-based computers and superconducting spintronics in general.


title: Wrinkles and creases in the bending, unbending and eversion of soft sectors

abstract: We study what is clearly one of the most common modes of deformation found in nature, science and engineering, namely the large elastic bending of curved structures, as well as its inverse, unbending, which can be brought beyond complete straightening to turn into eversion. We find that the suggested mathematical solution to these problems always exists and is unique when the solid is modelled as a homogeneous, isotropic, incompressible hyperelastic material with a strain-energy satisfying the strong ellipticity condition. We also provide explicit asympto-tic solutions for thin sectors. When the deformations are severe enough, the compressed side of the elastic material may buckle and wrinkles could then develop. We analyse in detail the onset of this instability for the Mooney-Rivlin strain energy, which covers the cases of the neo-Hookean model in exact non-linear elasticity and of third-order elastic materials in weakly non-linear elasticity. In particular the associated theoretical and numerical treatment allows us to predict the number and wavelength of the wrinkles. Guided by experimental observations we finally look at the development of creases, which we simulate through advanced finite element computations. In some cases the linearised analysis allows us to predict correctly the number and the wavelength of the creases, which turn out to occur only a few percent of strain earlier than the wrinkles.


title: On the Relation Between Linearity-Generating Processes and Linear-Rational Models

abstract: We review the notion of a linearity-generating (LG) process introduced by Gabaix (2007) and relate LG processes to linear-rational (LR) models studied by Filipovic, Larsson, and Trolle (2017). We show that every LR model can be represented as an LG process and vice versa. We find that LR models have two basic properties which make them an important representation of LG processes. First, LR models can be easily specified and made consistent with nonnegative interest rates. Second, LR models go naturally with the long-term risk factorization due to Alvarez and Jermann (2005), Hansen and Scheinkman (2009), and Qin and Linetsky (2017). Every LG process under the long forward measure can be represented as a lower dimensional LR model.


title: Asymptotics to all orders of the Euler--Darboux equation in a triangle

abstract: In Einstein's theory of relativity, the interaction of two collinearly polarized plane gravitational waves can be described by a Goursat problem for the Euler--Darboux equation in a triangular domain. In this paper, using a representation of the solution in terms of Abel integrals, we give a full asymptotic expansion of the solution near the diagonal of the triangle. The expansion is related to the formation of a curvature singularity of the spacetime. In particular, our framework allows for boundary data with derivatives which are singular at the corners. This level of generality is crucial for the application to gravitational waves.


title: Evaluating CBR Similarity Functions for BAM Switching in Networks with Dynamic Traffic Profile

abstract: In an increasingly complex scenario for network management, a solution that allows configuration in more autonomous way with less intervention of the network manager is expected. This paper presents an evaluation of similarity functions that are necessary in the context of using a learning strategy for finding solutions. The learning approach considered is based on Case-Based Reasoning (CBR) and is applied to a network scenario where different Bandwidth Allocation Models (BAMs) behaviors are used and must be eventually switched looking for the best possible network operation. In this context, it is required to identify and configure an adequate similarity function that will be used in the learning process to recover similar solutions previously considered. This paper introduces the similarity functions, explains the relevant aspects of the learning process in which the similarity function plays a role and, finally, presents a proof of concept for a specific similarity function adopted. Results show that the similarity function was capable to get similar results from the existing use case database. As such, the use of similarity functions with CBR technique has proved to be potentially satisfactory for supporting BAM switching decisions mostly driven by the dynamics of input traffic profile.


title: A SDN-based Flexible System for On-the-Fly Monitoring and Treatment of Security Events

abstract: The Software Defined Networking (SDN) paradigm decouples control and data planes, offering high programmability and a global view of the network. However, it is a challenge not only provide security in these next generation networks as well as allow that network attacks could be subjected to an incident and forensic treatment procedure. This paper proposes the implementation of flexible mechanisms of monitoring and treatment of security events categorized per type of attack and associated with whitelist and blacklist resources by means of the SDN controller programmability. The resources to perform intrusion and attack analysis are validated by means of a real SDN/OpenFlow testbed.


title: A Publish/Subscribe QoS-aware Framework for Massive IoT Traffic Orchestration

abstract: Internet of Things (IoT) application deployment requires the allocation of resources such as virtual machines, storage, and network elements that must be deployed over distinct infrastructures such as cloud computing, Cloud of Things (CoT), datacenters and backbone networks. For massive IoT data acquisition, a gateway-based data aggregation approach is commonly used featuring sensor/ actuator seamless access and providing cache/ buffering and preprocessing functionality. In this perspective , gateways acting as producers need to allocate network resources to send IoT data to consumers. In this paper, it is proposed a Publish/-Subscribe (PubSub) quality of service (QoS) aware framework (PSIoT-Orch) that orchestrates IoT traffic and allocates network resources between aggregates and consumers for massive IoT traffic. PSIoT-Orch schedules IoT data flows based on its configured QoS requirements. Additionally , the framework allocates network resources (LSP/ bandwidth) over a controlled backbone network with limited and constrained resources between IoT data users and consumers. Network resources are allocated using a Bandwidth Allocation Model (BAM) to achieve efficient network resource allocation for scheduled IoT data streams. PSIoT-Orch adopts an ICN (Information-Centric Network) PubSub architecture approach to handle IoT data transfers requests among framework components. The proposed framework aims at gathering the inherent advantages of an ICN-centric approach using a PubSub message scheme while allocating resources efficiently keeping QoS awareness and handling restricted network resources (bandwidth) for massive IoT traffic.


title: A Topological Invariant for Modular Fusion Categories

abstract: The modular data of a modular category $\mathcal{C}$, consisting of the $S$-matrix and the $T$-matrix, is known to be an incomplete invariant of $\mathcal{C}$. More generally, the invariants of framed links and knots defined by a modular category as part of a topological quantum field theory can be viewed as numerical invariants of the category. Among these invariants, we study the invariant defined by the Borromean link colored by three objects. Thus we obtain a tensor that we call $B$. We derive a formula for the Borromean tensor for the twisted Drinfeld doubles of finite groups. Along with $T$, it distinguishes the $p$ non-equivalent modular categories of the form $\mathcal{Z}({\rm Vec}_G^\omega)$ for $G$ the non-abelian group $\mathbb{Z}/q\mathbb{Z} \rtimes \mathbb{Z}/p\mathbb{Z}$, which are not distinguished by the modular data.

categories: math.QA math-ph math.CT math.MP math.RA math.RT


title: Multiphase Partitions of Lattice Random Walks

abstract: Considering the dynamics of non-interacting particles randomly moving on a lattice, the occurrence of a discontinuous transition in the values of the lattice parameters (lattice spacing and hopping times) determines the uprisal of two lattice phases. In this Letter we show that the hyperbolic hydrodynamic model obtained by enforcing the boundedness of lattice velocities derived by Giona (2018) correctly describes the dynamics of the system and permits to derive easily the boundary condition at the interface, which, contrarily to the common belief, involves the lattice velocities in the two phases and not the phase diffusivities. The dispersion properties of independent particles moving on an infinite lattice composed by the periodic repetition of a multiphase unit cell are investigated. It is shown that the hyperbolic transport theory correctly predicts the effective diffusion coefficient over all the range of parameter values, while the corresponding continuous parabolic models deriving from Langevin equations for particle motion fail. The failure of parabolic transport models is shown via a simple numerical experiment.


title: Reducing Metadata Leakage from Encrypted Files and Communication with PURBs

abstract: Most encrypted data formats leak metadata via their plaintext headers, such as format version, encryption schemes used, number of recipients who can decrypt the data, and even the recipients' identities. This leakage can pose security and privacy risks to users, e.g., by revealing the full membership of a group of collaborators from a single encrypted e-mail, or by enabling an eavesdropper to fingerprint the precise encryption software version and configuration the sender used. We propose that future encrypted data formats improve security and privacy hygiene by producing $\textit{Padded Uniform Random Blobs}$ or PURBs: ciphertexts indistinguishable from random bit strings to anyone without a decryption key. A PURB's content leaks $\textit{nothing at all}$, even the application that created it, and is padded such that even its length leaks as little as possible. Encoding and decoding ciphertexts with $\textit{no}$ cleartext markers presents efficiency challenges, however. We present cryptographically agile encodings enabling legitimate recipients to decrypt a PURB efficiently, even when encrypted for any number of recipients' public keys and/or passwords, and when these public keys are from different cryptographic suites. PURBs employ Padm\'e, a~novel padding scheme that limits information leakage via ciphertexts of maximum length $M$ to a practical optimum of $O(\log \log M)$ bits, comparable to padding to a power of two, but with lower overhead of at most $12\%$ and decreasing with larger payloads.


title: Electromagnetic duality and central charge

abstract: We provide a full realization of the electromagnetic duality at the boundary by extending the phase space of Maxwell's theory through the introduction of edge modes and their conjugate momenta. We show how such extension, which follows from a boundary action, is necessary in order to have well defined canonical generators of the boundary magnetic symmetries. In this way, both electric and magnetic soft modes are encoded in a boundary gauge field and its conjugate dual. This implementation of the electromagnetic duality has striking consequences. In particular, we show first how the electric charge quantization follows straightforwardly from the topological properties of the $U(1)$-bundle of the boundary dual potential. Moreover, having a well defined canonical action of the electric and magnetic symmetry generators on the phase space, we can compute their algebra and reveal the presence of a central charge between them. We conclude with possible implications of these results in the quantum theory.


title: Variable binaries and variables in binaries in the Binary Star Database BDB

abstract: The BDB, Binary star DataBase http://bdb.inasan.ru combines data of the catalogues of binary and multiple stars of all observational types. There is a number of ways for variable stars to form or to be a part of binary or multiple systems. We describe how such stars are represented in the database.


title: A Visual Proof that $\pi^e < e^\pi$

abstract: In this article, we give another visual proof of $\pi^e < e^\pi$.


title: The characterization of perfect Roman domination stable trees

abstract: A \emph{perfect Roman dominating function} (PRDF) on a graph $G = (V, E)$ is a function $f : V \rightarrow \{0, 1, 2\}$ satisfying the condition that every vertex $u$ for which $f(u) = 0$ is adjacent to exactly one vertex $v$ for which $f(v) = 2$. The weight of a PRDF is the value $w(f) = \sum_{u \in V}f(u)$. The minimum weight of a PRDF on a graph $G$ is called the \emph{perfect Roman domination number $\gamma_R^p(G)$} of $G$. A graph $G$ is perfect Roman domination domination stable if the perfect Roman domination number of $G$ remains unchanged under the removal of any vertex. In this paper, we characterize all trees that are perfect Roman domination stable.


title: A Generalized Matrix Splitting Algorithm

abstract: Composite function minimization captures a wide spectrum of applications in both computer vision and machine learning. It includes bound constrained optimization, $\ell_1$ norm regularized optimization, and $\ell_0$ norm regularized optimization as special cases. This paper proposes and analyzes a new Generalized Matrix Splitting Algorithm (GMSA) for minimizing composite functions. It can be viewed as a generalization of the classical Gauss-Seidel method and the Successive Over-Relaxation method for solving linear systems in the literature. Our algorithm is derived from a novel triangle operator mapping, which can be computed exactly using a new generalized Gaussian elimination procedure. We establish the global convergence, convergence rate, and iteration complexity of GMSA for convex problems. In addition, we also discuss several important extensions of GMSA. Finally, we validate the performance of our proposed method on three particular applications: nonnegative matrix factorization, $\ell_0$ norm regularized sparse coding, and $\ell_1$ norm regularized Dantzig selector problem. Extensive experiments show that our method achieves state-of-the-art performance in term of both efficiency and efficacy.


title: Observer design and practical stability of nonlinear systems under unknown time-delay

abstract: In the present paper, we study observer design and we establish some sufficient conditions for practical exponential stability for a class of time-delay nonlinear systems written in triangular form. In case of delay, the exponential convergence of the observer was confirmed. Based on the Lyapunov-Krasovskii functionals, the practical stability of the proposed observer is achieved. Finally, a physical model and simulation findings show the feasibility of the suggested strategy.


title: A spatial likelihood analysis for MAGIC telescope data

abstract: Context. The increase in sensitivity of Imaging Atmospheric Cherenkov Telescopes (IACTs) has lead to numerous detections of extended $\gamma$-ray sources at TeV energies, sometimes of sizes comparable to the instrument's field of view (FoV). This creates a demand for advanced and flexible data analysis methods, able to extract source information by utilising the photon counts in the entire FoV. Aims. We present a new software package, "SkyPrism", aimed at performing 2D (3D if energy is considered) fits of IACT data, possibly containing multiple and extended sources, based on sky images binned in energy. Though the development of this package was focused on the analysis of data collected with the MAGIC telescopes, it can further be adapted to other instruments, such as the future Cherenkov Telescope Array (CTA). Methods. We have developed a set of tools that, apart from sky images (count maps), compute the instrument response functions (IRFs) of MAGIC (effective exposure throughout the FoV, point spread function (PSF), energy resolution and background shape), based on the input data, Monte-Carlo simulations and the pointing track of the telescopes. With this information, the presented package can perform a simultaneous maximum likelihood fit of source models of arbitrary morphology to the sky images providing energy spectra, detection significances, and upper limits. Results. We demonstrate that the SkyPrism tool accurately reconstructs the MAGIC PSF, on and off-axis performance as well as the underlying background. We further show that for a point source analysis with MAGIC's default observational settings, SkyPrism gives results compatible with those of the standard tools while being more flexible and widely applicable.


title: Data-driven Analytics for Business Architectures: Proposed Use of Graph Theory

abstract: Business Architecture (BA) plays a significant role in helping organizations understand enterprise structures and processes, and align them with strategic objectives. However, traditional BAs are represented in fixed structure with static model elements and fail to dynamically capture business insights based on internal and external data. To solve this problem, this paper introduces the graph theory into BAs with aim of building extensible data-driven analytics and automatically generating business insights. We use IBM's Component Business Model (CBM) as an example to illustrate various ways in which graph theory can be leveraged for data-driven analytics, including what and how business insights can be obtained. Future directions for applying graph theory to business architecture analytics are discussed.


title: On-chip optical isolator and nonreciprocal parity-time symmetry induced by stimulated Brillouin scattering

abstract: Realization of chip-scale nonreciprocal optics such as isolators and circulators is highly demanding for all-optical signal routing and protection with standard photonics foundry process. Owing to the significant challenge for incorporating magneto-optical materials on chip, the exploration of magnetic-free alternatives has become exceedingly imperative in integrated photonics. Here, we demonstrate a chip-based, tunable all-optical isolator at the telecommunication band based upon bulk stimulated Brillouin scattering (SBS) in a high-Q silica microtoroid resonator. This device exhibits remarkable characteristics over most state-of-the-art implements, including high isolation ratio, no insertion loss, and large working power range. Thanks to the guided acoustic wave and accompanying momentum-conservation condition, SBS also enables us to realize the first nonreciprocal parity-time symmetry in two directly-coupled microresonators. The breach of time-reversal symmetry further makes the design a versatile arena for developing many formidable ultra-compact devices such as unidirectional single-mode Brillouin lasers and supersensitive photonic sensors.


title: On Critical Threshold Value for Simple Games

abstract: In this note, we show that for every simple game with n players the critical threshold value is at most n/4. This verifies the conjecture of Freixas and Kurz.


title: A compact ultrahigh vacuum scanning tunneling microscope with dilution refrigeration

abstract: We have designed and built a scanning tunneling microscope (STM) setup for operation at millikelvin temperatures in ultra high vacuum. A compact cryostat with an integrated dilution refrigerator has been built, that allows measurements at a base temperature of 25 mK in magnetic field up to 7.5 T with low mechanical and electronic noise. The cryostat is not larger than conventional helium bath cryostats (23 and 13 liters of nitrogen and helium respectively), so that the setup does not require a large experimental hall and fits easily into a standard lab space. Mechanical vibrations with running dilution circulation were kept below 300 fm/$\sqrt{\mathrm{Hz}}$ by mechanically decoupling the STM from the cryostat and the pumping system. All electronic input lines were low-pass filtered, reducing the electronic temperature to below 100 mK, as deduced from the quasiparticle peaks of superconducting aluminium. The microscope is optically accessible in the parked position, making sample and tip exchange fast and user-friendly. For measurement the STM is lowered 60 mm down so that the sample ends in the middle of a wet superconducting magnetic coil.


title: Brownian motion of magnetic domain walls and skyrmions, and their diffusion constants

abstract: Extended numerical simulations enable to ascertain the diffusive behavior at finite temperatures of chiral walls and skyrmions in ultra-thin model Co layers exhibiting symmetric - Heisenberg - as well as antisymmetric - Dzyaloshinskii-Moriya - exchange interactions. The Brownian motion of walls and skyrmions is shown to obey markedly different diffusion laws as a function of the damping parameter. Topology related skyrmion diffusion suppression with vanishing damping parameter, albeit already documented, is shown to be restricted to ultra-small skyrmion sizes or, equivalently, to ultra-low damping coefficients, possibly hampering observation.


title: The Computational 2D Materials Database: High-Throughput Modeling and Discovery of Atomically Thin Crystals

abstract: We introduce the Computational 2D Materials Database (C2DB), which organises a variety of structural, thermodynamic, elastic, electronic, magnetic, and optical properties of around 1500 two-dimensional materials distributed over more than 30 different crystal structures. Material properties are systematically calculated by state-of-the art density functional theory and many-body perturbation theory (G$_0\!$W$\!_0$ and the Bethe-Salpeter Equation for $\sim$200 materials) following a semi-automated workflow for maximal consistency and transparency. The C2DB is fully open and can be browsed online or downloaded in its entirety. In this paper, we describe the workflow behind the database, present an overview of the properties and materials currently available, and explore trends and correlations in the data. Moreover, we identify a large number of new potentially synthesisable 2D materials with interesting properties targeting applications within spintronics, (opto-)electronics, and plasmonics. The C2DB offers a comprehensive and easily accessible overview of the rapidly expanding family of 2D materials and forms an ideal platform for computational modeling and design of new 2D materials and van der Waals heterostructures.


title: A Markov Variation Approach to Smooth Graph Signal Interpolation

abstract: In this paper we present the Markov variation, a smoothness measure which offers a probabilistic interpretation of graph signal smoothness. This measure is then used to develop an optimization framework for graph signal interpolation. Our approach is based on diffusion embedding vectors and the connection between diffusion maps and signal processing on graphs. As diffusion embedding vectors may be expensive to compute for large graphs, we present a computationally efficient method, based on the Nystr\"{o}m extension, for interpolation of signals over a graph. We demonstrate our approach on the MNIST dataset and a dataset of daily average temperatures around the US. We show that our method outperforms state of the art graph signal interpolation techniques on both datasets, and that our computationally efficient reconstruction achieves slightly reduced accuracy with a large computational speedup.


title: A Multi-task Network to Detect Junctions in Retinal Vasculature

abstract: Junctions in the retinal vasculature are key points to be able to extract its topology, but they vary in appearance, depending on vessel density, width and branching/crossing angles. The complexity of junction patterns is usually accompanied by a scarcity of labels, which discourages the usage of very deep networks for their detection. We propose a multi-task network, generating labels for vessel interior, centerline, edges and junction patterns, to provide additional information to facilitate junction detection. After the initial detection of potential junctions in junction-selective probability maps, candidate locations are re-examined in centerline probability maps to verify if they connect at least 3 branches. The experiments on the DRIVE and IOSTAR showed that our method outperformed a recent study in which a popular deep network was trained as a classifier to find junctions. Moreover, the proposed approach is applicable to unseen datasets with the same degree of success, after training it only once.


title: Ensemble-LES Analysis of Perturbation Response of Turbulent Partially-Premixed Flames

abstract: The response to small perturbations of the Sandia CH4/air flame series is studied using a dynamical systems formula- tion. Here, the Sandia D and E turbulent partially-premixed flames are computed using large eddy simulation (LES) with a flamelet/progress variable approach. Using 300 simultaneous LES computations of each flame, the partial Lyapunov spectrum of the flow is obtained, which provides a set of Lyapunov exponents (LEs) and Lyapunov vectors (LVs). Special numerical procedures to handle such large data sets and perform matrix manipulations on-the-fly are discussed. From the LEs, the dimension of the attractor for the flames was estimated, for the first time, to be at least 5000 (flame D) and 10500 (flame E). This is a significant result indicating that the chaoticity of the turbulent flame is confined to a volume of phase-space that is quantified. The growth of perturbations occurs at time-scales that are much larger than the smallest strain-rate based times scales, but smaller than the integral scales. Further, the Lya- punov vectors that determine the perturbation response show that the strongest perturbations are aligned with the jet breakdown region, while other LVs are aligned with the shear layer. It is found that both the extinction/re-ignition process and the turbulence-induced jet motion strongly contribute to the flame chaoticity but are not equally important for all the LVs. These results provide new insights into the perturbation response of turbulent flames. More broadly, this Lyapunov approach allows a full characterization of the chaotic dynamics associated with turbulent flames, which have applications for both modeling and control of such complex flows.


title: Singularity methods for meromorphic solutions of differential equations

abstract: This chapter is mainly a tutorial introduction to methods recently developed in order to find all (as opposed to some) meromorphic particular solutions of given nonintegrable, autonomous, algebraic ordinary differential equations of any order. The examples are taken from physics and include Kuramoto-Sivashinsky and complex Ginzburg-Landau equations.


title: Chiral properties of graphene h-BN hybrid systems

abstract: The application of the chiral decomposition procedure to hybrid graphene h-BN systems revealed rules for the partition of the system into effective subsystems being bilayers plus monolayer in case the number of layers is odd. Three types of subsystems have been detected namely purely graphene bilayers and monolayers, mixed bilayers and pure h-BN monolayers depending on the hybrid composition. The effective parameters characterizing these chiral subsystems consist of the interlayer couplings and on-site potentials which shows the mechanism of compensation of the asymmetry introduced into the system by h-BN layers. For illustration, we provide a pedagogical overview about chiral tunneling in graphene subsystems (MLG, BLG) present in hybrid with one h-BN layer. We have established the parameter ranges for which the characteristic features in the spectrum are observed, such as Fabry-P\'erot resonances in the case of MLG and "magic angles" in the case of effective BLG. We also consider different hybrid stacking in order to indicate effective systems with the desired properties required in the electronic and spintronic applications.


title: Quasi-monotonicity formulas for classical obstacle problems with Sobolev coefficients and applications

abstract: We establish Weiss' and Monneau's type quasi-monotonicity formulas for quadratic energies having matrix of coefficients in a Sobolev space $W^{1,p}$, $p>n$, and provide an application to the corresponding free boundary analysis for the related classical obstacle problems.


title: The arithmetic puncturing problem and integral points

abstract: In this paper, we provide evidence to support a positive answer to a question of Hassett and Tschinkel. In particular, if an algebraic variety V has a dense set of rational points, they ask whether or not the set of D-integral points is potentially dense, where D is a set of codimension at least two. We give a positive answer to this question in many cases, including varieties whose generic linear section is a smooth rational curve, and certain K3 surfaces. We also discuss some stronger notions of integrality of points, and give some positive answers to some cases of the analogous question in the stronger context.


title: Equivalent partial differential equations of a lattice Boltzmann scheme

abstract: We show that when we formulate the lattice Boltzmann equation with a small time step $\Delta$t and an associated space scale $\Delta$x, a Taylor expansion joined with the so-called equivalent equation methodology leads to establish macroscopic fluid equations as a formal limit. We recover the Euler equations of gas dynamics at the first order and the compressible Navier-Stokes equations at the second order. 1) Discrete geometry $\bullet$ We denote by d the dimension of space and by L a regular d-dimensional lattice. Such a lattice is composed by a set L 0 of nodes or vertices and a set L 1 of links or edges between two vertices. From a practical point of view, given a vertex x, there exists a set V (x) of neighbouring nodes, including the node x itself. We consider here that the lattice L is parametrized by a space step $\Delta$x > 0. For the fundamental example called D2Q9 (see e.g. Lallemand and Luo, 2000), the set V (x) is given with the help of the family of vectors (e j) 0$\le$j$\le$J defined by J = 8, (1.1) e j = 0 0 , 1 0 , 0 1 , --1 0 , 0 --1 , 1 1 , --1 1 , --1 --1 , 1 --1 and the vicinity (1.2) V (x) = { x + $\Delta$x e j , 0 $\le$ j $\le$ J } .


title: Deep learning based inverse method for layout design

abstract: Layout design with complex constraints is a challenging problem to solve due to the non-uniqueness of the solution and the difficulties in incorporating the constraints into the conventional optimization-based methods. In this paper, we propose a design method based on the recently developed machine learning technique, Variational Autoencoder (VAE). We utilize the learning capability of the VAE to learn the constraints and the generative capability of the VAE to generate design candidates that automatically satisfy all the constraints. As such, no constraints need to be imposed during the design stage. In addition, we show that the VAE network is also capable of learning the underlying physics of the design problem, leading to an efficient design tool that does not need any physical simulation once the network is constructed. We demonstrated the performance of the method on two cases: inverse design of surface diffusion induced morphology change and mask design for optical microlithography.


title: Base Station Switch-off with Mutual Repulsion in 5G Massive MIMO Networks

abstract: When small cells are densely deployed in the fifth generation (5G) cellular networks, switching off a part of base stations (BSs) is a practical approach for saving energy consumption considering the variation of traffic load. The small cell network with the massive multi-input multi-output (massive MIMO) system is analyzed in this paper due to the dense deployment and low power consumption. Based on the BS switch-off strategy with distance constraints, the energy and coverage efficiency are investigated to illustrate the performance of the BS switch-off strategy. Simulation results indicate that the energy efficiency and coverage efficiency of the proposed strategy are better than the random strategy. The energy efficiency increases with the BS intensity and the minimal distance, and a maximum coverage efficiency can be achieved with the increase of the BS intensity and the minimum distance. In this case, the optimal BS switch-off strategy can be designed under this work in the actual scene.


title: Surgical Data Science: A Consensus Perspective

abstract: Surgical data science is a scientific discipline with the objective of improving the quality of interventional healthcare and its value through capturing, organization, analysis, and modeling of data. The goal of the 1st workshop on Surgical Data Science was to bring together researchers working on diverse topics in surgical data science in order to discuss existing challenges, potential standards and new research directions in the field. Inspired by current open space and think tank formats, it was organized in June 2016 in Heidelberg. While the first day of the workshop, which was dominated by interactive sessions, was open to the public, the second day was reserved for a board meeting on which the information gathered on the public day was processed by (1) discussing remaining open issues, (2) deriving a joint definition for surgical data science and (3) proposing potential strategies for advancing the field. This document summarizes the key findings.


title: Wave-U-Net: A Multi-Scale Neural Network for End-to-End Audio Source Separation

abstract: Models for audio source separation usually operate on the magnitude spectrum, which ignores phase information and makes separation performance dependant on hyper-parameters for the spectral front-end. Therefore, we investigate end-to-end source separation in the time-domain, which allows modelling phase information and avoids fixed spectral transformations. Due to high sampling rates for audio, employing a long temporal input context on the sample level is difficult, but required for high quality separation results because of long-range temporal correlations. In this context, we propose the Wave-U-Net, an adaptation of the U-Net to the one-dimensional time domain, which repeatedly resamples feature maps to compute and combine features at different time scales. We introduce further architectural improvements, including an output layer that enforces source additivity, an upsampling technique and a context-aware prediction framework to reduce output artifacts. Experiments for singing voice separation indicate that our architecture yields a performance comparable to a state-of-the-art spectrogram-based U-Net architecture, given the same data. Finally, we reveal a problem with outliers in the currently used SDR evaluation metrics and suggest reporting rank-based statistics to alleviate this problem.


title: Local law and Tracy-Widom limit for sparse sample covariance matrices

abstract: We consider spectral properties of sparse sample covariance matrices, which includes biadjacency matrices of the bipartite Erd\H{o}s-R\'enyi graph model. We prove a local law for the eigenvalue density up to the upper spectral edge. Under a suitable condition on the sparsity, we also prove that the limiting distribution of the rescaled, shifted extremal eigenvalues is given by the GOE Tracy-Widom law with an explicit formula on the deterministic shift of the spectral edge. For the biadjacency matrix of an Erd\H{o}s-R\'enyi graph with two vertex sets of comparable sizes $M$ and $N$, this establishes Tracy-Widom fluctuations of the second largest eigenvalue when the connection probability $p$ is much larger than $N^{-2/3}$ with a deterministic shift of order $(Np)^{-1}$.


title: Strategies to cure numerical shock instability in HLLEM Riemann solver

abstract: The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver for cheap and accurate computation of high speed gasdynamical flows. Unfortunately this scheme is known to be plagued by various forms of numerical shock instability. In this paper we present various strategies to save the HLLEM scheme from developing such spurious solutions. A linear scale analysis of its mass and interface-normal momentum flux discretizations reveal that its antidiffusive terms, which are primarily responsible for resolution of linear wavefields, are inadvertently activated along a normal shock front due to numerical perturbations. These erroneously activated terms counteract the favourable damping mechanism provided by its inherent HLL-type diffusive terms and trigger the shock instability. To avoid this, two different strategies are proposed for discretization of these critical flux components in the vicinity of a shock: one that deals with increasing the magnitude of inherent HLL-type dissipation through careful manipulation of specific non-linear wave speed estimates while the other deals with reducing the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and estimate von-Neumann type stability bounds on the CFL number arising from their use. Results from classic numerical test cases show that both types of modified HLLEM schemes are able to provide excellent shock stable solutions while retaining commendable accuracy on shear dominated viscous flows.


title: Mixed integer nonlinear programming for Joint Coordination of Plug-in Electrical Vehicles Charging and Smart Grid Operations

abstract: The problem of joint coordination of plug-in electric vehicles (PEVs) charging and grid power control is to minimize both PEVs charging cost and energy generation cost while meeting both residential and PEVs' power demands and suppressing the potential impact of PEVs integration. A bang-bang PEV charging strategy is adopted to exploit its simple online implementation, which requires computation of a mixed integer nonlinear programming problem (MINP) in binary variables of the PEV charging strategy and continuous variables of the grid voltages. A new solver for this MINP is proposed. Its efficiency is shown by numerical simulations.


title: Local Lie derivations on von Neumann algebras and algebras of locally measurable operators

abstract: Let $\mathcal{A}$ be a unital associative algebra and $\mathcal{M}$ be an $\mathcal{A}$-bimodule. A linear mapping $\varphi$ from $\mathcal{A}$ into an $\mathcal{A}$-bimodule $\mathcal{M}$ is called a Lie derivation if $\varphi[A,B]=[\varphi(A),B]+[A,\varphi(B)]$ for each $A,B$ in $\mathcal{A}$, and $\varphi$ is called a \emph{local Lie derivation} if for every $A$ in $\mathcal{A}$, there exists a Lie derivation $\varphi_{A}$ (depending on $A$) from $\mathcal{A}$ into $\mathcal{M}$ such that $\varphi(A)=\varphi_{A}(A)$. In this paper, we prove that every local Lie derivation on von Neumann algebras is a Lie derivation; and we show that if $\mathcal M$ is a type I von Neumann algebra with atomic lattice of projections, then every local Lie derivation on $LS(\mathcal M)$ is a Lie derivation.


title: The Well Tempered Lasso

abstract: We study the complexity of the entire regularization path for least squares regression with 1-norm penalty, known as the Lasso. Every regression parameter in the Lasso changes linearly as a function of the regularization value. The number of changes is regarded as the Lasso's complexity. Experimental results using exact path following exhibit polynomial complexity of the Lasso in the problem size. Alas, the path complexity of the Lasso on artificially designed regression problems is exponential. We use smoothed analysis as a mechanism for bridging the gap between worst case settings and the de facto low complexity. Our analysis assumes that the observed data has a tiny amount of intrinsic noise. We then prove that the Lasso's complexity is polynomial in the problem size. While building upon the seminal work of Spielman and Teng on smoothed complexity, our analysis is morally different as it is divorced from specific path following algorithms. We verify the validity of our analysis in experiments with both worst case settings and real datasets. The empirical results we obtain closely match our analysis.


title: Hearst Patterns Revisited: Automatic Hypernym Detection from Large Text Corpora

abstract: Methods for unsupervised hypernym detection may broadly be categorized according to two paradigms: pattern-based and distributional methods. In this paper, we study the performance of both approaches on several hypernymy tasks and find that simple pattern-based methods consistently outperform distributional methods on common benchmark datasets. Our results show that pattern-based models provide important contextual constraints which are not yet captured in distributional methods.


title: Assessing the impact of machine intelligence on human behaviour: an interdisciplinary endeavour

abstract: This document contains the outcome of the first Human behaviour and machine intelligence (HUMAINT) workshop that took place 5-6 March 2018 in Barcelona, Spain. The workshop was organized in the context of a new research programme at the Centre for Advanced Studies, Joint Research Centre of the European Commission, which focuses on studying the potential impact of artificial intelligence on human behaviour. The workshop gathered an interdisciplinary group of experts to establish the state of the art research in the field and a list of future research challenges to be addressed on the topic of human and machine intelligence, algorithm's potential impact on human cognitive capabilities and decision making, and evaluation and regulation needs. The document is made of short position statements and identification of challenges provided by each expert, and incorporates the result of the discussions carried out during the workshop. In the conclusion section, we provide a list of emerging research topics and strategies to be addressed in the near future.


title: Dynamical typicality approach to eigenstate thermalization

abstract: We consider the set of all initial states within a microcanonical energy shell of an isolated many-body quantum system, which exhibit the same, arbitrary but fixed non-equilibrium expectation value for some given observable $A$. On condition that this set is not too small, it is shown by means of a dynamical typicality approach that most such initial states exhibit thermalization if and only if $A$ satisfies the so-called weak eigenstate thermalization hypothesis (wETH). Here, thermalization means that the expectation value of $A$ spends most of its time close to the microcanonical value after initial transients have died out. The wETH means that, within the energy shell, most eigenstates of the pertinent system Hamiltonian exhibit very similar expectation values of $A$.


title: Evolutionary Game Dynamics for Two Interacting Populations under Environmental Feedback

abstract: We study the evolutionary dynamics of games under environmental feedback using replicator equations for two interacting populations. One key feature is to consider jointly the co-evolution of the dynamic payoff matrices and the state of the environment: the payoff matrix varies with the changing environment and at the same time, the state of the environment is affected indirectly by the changing payoff matrix through the evolving population profiles. For such co-evolutionary dynamics, we investigate whether convergence will take place, and if so, how. In particular, we identify the scenarios where oscillation offers the best predictions of long-run behavior by using reversible system theory. The obtained results are useful to describe the evolution of multi-community societies in which individuals' payoffs and societal feedback interact.


title: Obtaining fairness using optimal transport theory

abstract: Statistical algorithms are usually helping in making decisions in many aspects of our lives. But, how do we know if these algorithms are biased and commit unfair discrimination of a particular group of people, typically a minority? \textit{Fairness} is generally studied in a probabilistic framework where it is assumed that there exists a protected variable, whose use as an input of the algorithm may imply discrimination. There are different definitions of Fairness in the literature. In this paper we focus on two of them which are called Disparate Impact (DI) and Balanced Error Rate (BER). Both are based on the outcome of the algorithm across the different groups determined by the protected variable. The relationship between these two notions is also studied. The goals of this paper are to detect when a binary classification rule lacks fairness and to try to fight against the potential discrimination attributable to it. This can be done by modifying either the classifiers or the data itself. Our work falls into the second category and modifies the input data using optimal transport theory.


title: Approximation of Hermitian Matrices by Positive Semidefinite Matrices using Modified Cholesky Decompositions

abstract: A new algorithm to approximate Hermitian matrices by positive semidefinite Hermitian matrices based on modified Cholesky decompositions is presented. In contrast to existing algorithms, this algorithm allows to specify bounds on the diagonal values of the approximation. It has no significant runtime and memory overhead compared to the computation of a classical Cholesky decomposition. Hence it is suitable for large matrices as well as sparse matrices since it preserves the sparsity pattern of the original matrix. The algorithm tries to minimize the approximation error in the Frobenius norm as well as the condition number of the approximation. Since these two objectives often contradict each other, it is possible to weight these two objectives by parameters of the algorithm. In numerical experiments, the algorithm outperforms existing algorithms regarding these two objectives. A Cholesky decomposition of the approximation is calculated as a byproduct. This is useful, for example, if a corresponding linear equation should be solved. A fully documented and extensively tested implementation is available. Numerical optimization and statistics are two fields of application in which the algorithm can be of particular interest.


title: Relation modules for finite W-algebras and tensor products of highest weight evaluation modules for Yangians of type A

abstract: We construct explicitly a large family of Gelfand-Tsetlin modules for an arbitrary finite W-algebra of type A and establish their irreducibility. A basis of these modules is formed by the Gelfand-Tsetlin tableaux whose entries satisfy certain admissible sets of relations. Characterization and an effective method of constructing such admissible relations are given. In the case of the Yangian of gl_n we prove the sufficient condition for the irreducibility of the tensor product of two highest weight relation modules and establish irreducibility of any number of highest weight relation modules with generic highest weights. This extends the results of Molev to infinite dimensional highest modules.


title: Spreading vectors for similarity search

abstract: Discretizing multi-dimensional data distributions is a fundamental step of modern indexing methods. State-of-the-art techniques learn parameters of quantizers on training data for optimal performance, thus adapting quantizers to the data. In this work, we propose to reverse this paradigm and adapt the data to the quantizer: we train a neural net which last layer forms a fixed parameter-free quantizer, such as pre-defined points of a hyper-sphere. As a proxy objective, we design and train a neural network that favors uniformity in the spherical latent space, while preserving the neighborhood structure after the mapping. We propose a new regularizer derived from the Kozachenko--Leonenko differential entropy estimator to enforce uniformity and combine it with a locality-aware triplet loss. Experiments show that our end-to-end approach outperforms most learned quantization methods, and is competitive with the state of the art on widely adopted benchmarks. Furthermore, we show that training without the quantization step results in almost no difference in accuracy, but yields a generic catalyzer that can be applied with any subsequent quantizer.


title: Galaxy Cluster Mass Reconstruction Project - IV. Understanding the effects of imperfect membership on cluster mass estimation

abstract: The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.


title: Quantum advantage of unitary Clifford circuits with magic state inputs

abstract: We study the computational power of unitary Clifford circuits with solely magic state inputs (CM circuits), supplemented by classical efficient computation. We show that CM circuits are hard to classically simulate up to multiplicative error (assuming PH non-collapse), and also up to additive error under plausible average-case hardness conjectures. Unlike other such known classes, a broad variety of possible conjectures apply. Along the way we give an extension of the Gottesman-Knill theorem that applies to universal computation, showing that for Clifford circuits with joint stabiliser and non-stabiliser inputs, the stabiliser part can be eliminated in favour of classical simulation, leaving a Clifford circuit on only the non-stabiliser part. Finally we discuss implementational advantages of CM circuits.


title: How long does the surplus stay close to its historical high?

abstract: In this paper we find the Laplace transforms of the weighted occupation times for a spectrally negative L\'evy surplus process to spend below its running maximum up to the first exit times. The results are expressed in terms of generalized scale functions. For step weight functions, the Laplace transforms can be further expressed in terms of scale functions.


title: No correlation between Solar flares and the decay rate of several $\beta$-decaying isotopes

abstract: We report on finding no correlation between the two strongest observed Solar flares in September 2017 and the decay rates of $^{60}$Co, $^{44}$Ti and $^{137}$Cs sources, which are continuously measured by two independent NaI(Tl) detector setups. We test for variations in the number of observed counts with respect to the number of expected counts over multiple periods with timescales varying from 1 to 109 hours around the Solar flare. No excess or deficit exceeds the 2$\sigma$ global significance. We set a conservative lower limit on the decay rate deviation over an 84-hour period around the two correlated Solar flares in September 2017 to $0.044\%$ with 2$\sigma$ confidence. A fractional change of $~0.1\%$ in the decay rate of $^{54}$Mn over a period of 84 hours was claimed with 7$\sigma$ significance during multiple Solar flares in December 2006. We exclude such an effect at 4.7$\sigma$ significance.


title: The $10^{-3}$ eV Frontier in Neutrinoless Double Beta Decay

abstract: The observation of neutrinoless double beta decay would allow to establish lepton number violation and the Majorana nature of neutrinos. The rate of this process in the case of 3-neutrino mixing is controlled by the neutrinoless double beta decay effective Majorana mass $|\langle m \rangle|$. For a neutrino mass spectrum with normal ordering, which is favoured over the spectrum with inverted ordering by recent global fits, $|\langle m \rangle|$ can be significantly suppressed. Taking into account updated data on the neutrino oscillation parameters, we investigate the conditions under which $|\langle m \rangle|$ in the case of spectrum with normal ordering exceeds $10^{-3}~(5\times 10^{-3})$ eV: $|\langle m \rangle|_\text{NO} > 10^{-3}~(5\times 10^{-3})$ eV. We analyse first the generic case with unconstrained leptonic CP violation Majorana phases. We show, in particular, that if the sum of neutrino masses is found to satisfy $\Sigma > 0.10$ eV, then $|\langle m \rangle|_\text{NO} > 5\times 10^{-3}$ eV for any values of the Majorana phases. We consider also cases where the values for these phases are either CP conserving or are in line with predictive schemes combining flavour and generalised CP symmetries.


title: Low-Complexity Multiuser QAM Detection for Uplink 1-bit Massive MIMO

abstract: This work studies multiuser detection for one-bit massive multiple-input multiple-output (MIMO) systems in order to diminish the power consumption at the base station (BS). A low-complexity near-maximum-likelihood (nML) multiuser detection algorithm is designed, assuming that each BS antenna port is connected with a pair of single-bit resolution analog-to-digital converters (ADCs) and each user equipment (UE) transmits symbols from a quadrature amplitude modulation (QAM) constellation. First, a novel convex program is formulated as a convex surrogate of the ML detector and subsequently solved through an accelerated first-order method. Then, the solution of the convex optimization problem is harnessed to solve a refined combinatorial problem with reduced search space, requiring non-exponential complexity on the number of the UEs. Judicious simulation study corroborates the efficacy of the resulting two-phase detection algorithm. The proposed two-phase algorithm can achieve symbol error rate (SER) performance close to the ML detector, with significantly reduced computation cost compared to the nML detection schemes in prior art.


title: Formal Small-step Verification of a Call-by-value Lambda Calculus Machine

abstract: We formally verify an abstract machine for a call-by-value lambda-calculus with de Bruijn terms, simple substitution, and small-step semantics. We follow a stepwise refinement approach starting with a naive stack machine with substitution. We then refine to a machine with closures, and finally to a machine with a heap providing structure sharing for closures. We prove the correctness of the three refinement steps with compositional small-step bottom-up simulations. There is an accompanying Coq development verifying all results.


title: Cantor-Bendixson type ranks on Polish spaces

abstract: For any Polish space $X$ it is well-known that the Cantor-Bendixson rank provides a co-analytic rank on $F_{\aleph_0}(X)$ if and only if $X$ is a $\sigma$-compact. In the case of $\omega^\omega$ one may recover a co-analytic rank on $F_{\aleph_0}(\omega^\omega)$ by considering the Cantor-Bendixson rank of the induced trees instead. In this paper we will generalize this idea to arbitrary Polish spaces and thereby construct a family of co-analytic ranks on $F_{\aleph_0}(X)$ for any Polish space $X$. We study the behaviour of this family and compare the ranks to the original Cantor-Bendixson rank. The main results are characterizations of the compact and $\sigma$-compact Polish spaces in terms of this behaviour.


title: Learning in Integer Latent Variable Models with Nested Automatic Differentiation

abstract: We develop nested automatic differentiation (AD) algorithms for exact inference and learning in integer latent variable models. Recently, Winner, Sujono, and Sheldon showed how to reduce marginalization in a class of integer latent variable models to evaluating a probability generating function which contains many levels of nested high-order derivatives. We contribute faster and more stable AD algorithms for this challenging problem and a novel algorithm to compute exact gradients for learning. These contributions lead to significantly faster and more accurate learning algorithms, and are the first AD algorithms whose running time is polynomial in the number of levels of nesting.


title: Prediction of the FIFA World Cup 2018 - A random forest approach with an emphasis on estimated team ability parameters

abstract: In this work, we compare three different modeling approaches for the scores of soccer matches with regard to their predictive performances based on all matches from the four previous FIFA World Cups 2002 - 2014: Poisson regression models, random forests and ranking methods. While the former two are based on the teams' covariate information, the latter method estimates adequate ability parameters that reflect the current strength of the teams best. Within this comparison the best-performing prediction methods on the training data turn out to be the ranking methods and the random forests. However, we show that by combining the random forest with the team ability parameters from the ranking methods as an additional covariate we can improve the predictive power substantially. Finally, this combination of methods is chosen as the final model and based on its estimates, the FIFA World Cup 2018 is simulated repeatedly and winning probabilities are obtained for all teams. The model slightly favors Spain before the defending champion Germany. Additionally, we provide survival probabilities for all teams and at all tournament stages as well as the most probable tournament outcome.


title: Analysis of Length Normalization in End-to-End Speaker Verification System

abstract: The classical i-vectors and the latest end-to-end deep speaker embeddings are the two representative categories of utterance-level representations in automatic speaker verification systems. Traditionally, once i-vectors or deep speaker embeddings are extracted, we rely on an extra length normalization step to normalize the representations into unit-length hyperspace before back-end modeling. In this paper, we explore how the neural network learns length-normalized deep speaker embeddings in an end-to-end manner. To this end, we add a length normalization layer followed by a scale layer before the output layer of the common classification network. We conducted experiments on the verification task of the Voxceleb1 dataset. The results show that integrating this simple step in the end-to-end training pipeline significantly boosts the performance of speaker verification. In the testing stage of our L2-normalized end-to-end system, a simple inner-product can achieve the state-of-the-art.


title: Does The Cloud Need Stabilizing?

abstract: The last decade has witnessed rapid proliferation of cloud computing. While even the smallest distributed programs (with 3-5 actions) produce many unanticipated error cases due to concurrency involved, it seems short of a miracle these web-services are able to operate at those vast scales. In this paper, we explore the factors that contribute most to the high-availability of cloud computing services and examine where self-stabilization could fit in that picture.


title: The landscape of NeuroImage-ing research

abstract: As the field of neuroimaging grows, it can be difficult for scientists within the field to gain and maintain a detailed understanding of its ever-changing landscape. While collaboration and citation networks highlight important contributions within the field, the roles of and relations among specific areas of study can remain quite opaque. Here, we apply techniques from network science to map the landscape of neuroimaging research documented in the journal NeuroImage over the past decade. We create a network in which nodes represent research topics, and edges give the degree to which these topics tend to be covered in tandem. The network displays small-world architecture, with communities characterized by common imaging modalities and medical applications, and with bridges that integrate these distinct subfields. Using node-level analysis, we quantify the structural roles of individual topics within the neuroimaging landscape, and find high levels of clustering within the structural MRI subfield as well as increasing participation among topics related to psychiatry. The overall prevalence of a topic is unrelated to the prevalence of its neighbors, but the degree to which a topic becomes more or less popular over time is strongly related to changes in the prevalence of its neighbors. Broadly, this work presents a cohesive model for understanding the landscape of neuroimaging research across the field, in broad subfields, and within specific topic areas.


title: Nonadiabatic coupling effects in MgB$_2$ reexamined

abstract: The unusual Raman spectrum of MgB$_2$ and its formidable temperature dependence are successfully reproduced by means of a parameter-free \emph{ab initio} nonadiabatic theory that accounts for the electron-hole pair scattering mechanisms with the system phonons. This example turns out to be a prototypical case where a strong nonadiabatic renormalization of the phonon frequency is partially washed out by the aforementioned scattering events, bringing along a characteristic temperature dependence. Both electron-hole pair lifetime and energy renormalization effects due to dynamical electron-phonon coupling turn out to play a crucial role. This theory could aid in comprehending other Raman spectra characterized with unconventionally strong electron-phonon interaction.


title: Impact of End-User Behavior on User/Network Association in HetNets

abstract: We study the impact of end-user behavior on user/network association in a HetNet with multiple service providers (SPs). Specifically, we consider the uncertainty in the service guarantees offered by SPs in a HetNet, and use Prospect Theory (PT) to model end-user decision making. We formulate user association with SPs as a multiple leader Stackelberg game where each SP offers a data rate to each user with a certain service guarantee and at a certain price, while the user chooses the best offer among multiple such bids. Using the specific example of a HetNet with one cellular base station and one WiFi access point, we show that when the end users underestimate the advertised service guarantees, then some of the Nash Equilibrium strategies under the Expected Utility Theory (EUT) model become infeasible under PT, and for those Nash Equilibria that are feasible under both EUT and PT, the resulting user utilities are less under PT. We propose resource allocation and bidding mechanisms for the SPs to mitigate these effects.


title: Rupture of amorphous graphene via void formation

abstract: Apart from its unique and exciting electronic properties, many sensor based applications of graphene are purely based on its mechanical and structural properties. Here we report a numerical and analytical study of a void in amorphous (small domain polycrystalline) graphene, and show that the energetics of a void is a balance between the line tension cost versus the increased area gain. Using the concepts of classical nucleation theory, we show that the critical radius of a void formed in amorphous graphene at constant pressure is simply the ratio of line tension at the void and the applied pressure. The values of the critical radius of the void for flat and buckled graphene are 3.48\AA ~and 3.31\AA, respectively at 2 eV/\AA${^2}$ pressure. We also show that the dominant finite size correction to the line tension is inversely proportional to the radius of the void in both flat and buckled cases. Contrary to conventional wisdom, with the help of a simple analytical model we find that the shear modulus sets the lower limit of the line tension in the samples. This makes our study relevant for other two-dimensional amorphous materials such as h-BN, phosphorene, borophene, and transition metal dichalcogenides. Our results are useful for the better understanding of polycrystalline graphene under tension and therefore have direct implications on the very fascinating field of strain engineering known as "straintronics" to manipulate or improve graphene's properties.


title: (In)Secure Configuration Practices of WPA2 Enterprise Supplicants

abstract: WPA2 Enterprise is a fundamental technology for secure communication in enterprise wireless networks. A key requirement of this technology is that WiFi-enabled devices (i.e., supplicants) be correctly configured before connecting to the enterprise wireless network. Supplicants that are not configured correctly may fall prey of attacks aimed at stealing the network credentials very easily. Such credentials have an enormous value because they usually unlock access to all enterprise services. In this work we investigate whether users and technicians are aware of these important and widespread risks. We conducted two extensive analyses: a survey among approximately 1000 users about how they configured their WiFi devices for enterprise network access; and, a review of approximately 310 network configuration guides made available by enterprise network administrators. The results provide strong indications that the key requirement of WPA2 Enterprise is violated systematically and thus can no longer be considered realistic.


title: Standard conjectures for abelian fourfolds

abstract: Let A be an abelian fourfold. We prove the Standard Conjecture of Hodge type for A. By combining this result with a theorem of Clozel we deduce that numerical equivalence on A coincides with l-adic homological equivalence on A for infinitely many l. The approach consists in reformulating this question into a p-adic problem and then using p-adic Hodge Theory to solve it.


title: Peak positions of strongly unimodal sequences

abstract: We study combinatorial and asymptotic properties of the rank of strongly unimodal sequences. We find a generating function for the rank enumeration function, and give a new combinatorial interpretation of the ospt-function introduced by Andrews, Chan, and Kim. We conjecture that the enumeration function for the number of unimodal sequences of a fixed size and varying rank is log-concave, and prove an asymptotic result in support of this conjecture. Finally, we determine the asymptotic behavior of the rank for strongly unimodal sequences, and prove that its values (when appropriately renormalized) are normally distributed with mean zero in the asymptotic limit.


title: Data-driven model for the identification of the rock type at a drilling bit

abstract: Directional oil well drilling requires high precision of the wellbore positioning inside the productive area. However, due to specifics of engineering design, sensors that explicitly determine the type of the drilled rock are located farther than 15m from the drilling bit. As a result, the target area runaways can be detected only after this distance, which in turn, leads to a loss in well productivity and the risk of the need for an expensive re-boring operation. We present a novel approach for identifying rock type at the drilling bit based on machine learning classification methods and data mining on sensors readings. We compare various machine-learning algorithms, examine extra features coming from mathematical modeling of drilling mechanics, and show that the real-time rock type classification error can be reduced from 13.5 % to 9 %. The approach is applicable for precise directional drilling in relatively thin target intervals of complex shapes and generalizes appropriately to new wells that are different from the ones used for training the machine learning model.


title: Comment on "Measurement of the $Q^2$ Dependence of the Deuteron Spin Structure Function $g_1$ and its Moments at Low $Q^2$ with CLAS"

abstract: We argue that the recently published CLAS results on the deuteron spin polarizability $\gamma_0$ [Adhikari {\it et al.}, Phys.\ Rev.\ Lett.\ {\bf 120}, 062501 (2018)], as well as their comparisons with chiral perturbation theory ($\chi$PT), are misleading. In reality, the deuteron polarizability is larger by 4 orders of magnitude, as demonstrated here by a novel calculation in pionless EFT. The CLAS paper, on the other hand, presents only a tiny correction to it, based on a partial evaluation of a sum rule for $\gamma_0$. The sum rule is assumed to have the same form as for the nucleon; we argue it does not. Moreover, their "test of $\chi$PT" tacitly involves assumptions which, as we demonstrate, may not be valid at the claimed accuracy.


title: A Scenario Decomposition Algorithm for Strategic Time Window Assignment Vehicle Routing Problems

abstract: We study the strategic decision-making problem of assigning time windows to customers in the context of vehicle routing applications that are affected by operational uncertainty. This problem, known as the Time Window Assignment Vehicle Routing Problem, can be viewed as a two-stage stochastic optimization problem, where time window assignments constitute first-stage decisions, vehicle routes adhering to the assigned time windows constitute second-stage decisions, and the objective is to minimize the expected routing costs. We prove that a sampled deterministic equivalent of this stochastic model can be reduced to a variant of the Consistent Vehicle Routing Problem, and we leverage this result to develop a new scenario decomposition algorithm to solve it. From a modeling viewpoint, our approach can accommodate both continuous and discrete sets of feasible time window assignments as well as general scenario-based models of uncertainty for several routing-specific parameters, including customer demands and travel times, among others. From an algorithmic viewpoint, our approach can be easily parallelized, can utilize any available vehicle routing solver as a black box, and can be readily modified as a heuristic for large-scale instances. We perform a comprehensive computational study to demonstrate that our algorithm strongly outperforms all existing solution methods, as well as to quantify the trade-off between computational tractability and expected cost savings when considering a larger number of future scenarios during strategic time window assignment.


title: A Proof the Functional Equation Conjecture

abstract: In the early 2000's the first and second named authors worked for a period of six years in an attempt of proving the Compositional Shuffle Conjecture [1]. Their approach was based on the discovery that all the Combinatorial properties predicted by the Compositional Shuffle Conjecture remain valid for each family of Parking Functions with prescribed diagonal cars. The validity of this property was reduced to the proof of a functional equation satisfied by a Catalan family of univariate polynomials. The main result in this paper is a proof of this functional equation. The Compositional Shuffle Conjecture was proved in 2015 by Eric Carlsson and Anton Mellit [3]. Our proof of the Functional Equation removes one of the main obstacles in the completion of the Garsia-Hicks approach to the proof of the Compositional Shuffle Conjecture. At the end of this writing we formulate a few further conjectures including what remains to be proved to complete this approach.


title: Tunneling spectroscopy of gate-induced superconductivity in MoS$_2$

abstract: The ability to gate-induce superconductivity by electrostatic charge accumulation is a recent breakthrough in physics and nano-electronics. With the exception of LaAlO$_3$/SrTiO$_3$ interfaces, experiments on gate-induced superconductors have been largely confined to resistance measurements, which provide very limited information about the superconducting state. Here, we explore gate-induced superconductivity in MoS$_2$ by performing tunneling spectroscopy to determine the energy-dependent density of states (DOS) for different levels of electron density $\textit{n}$. In the superconducting state, the DOS is strongly suppressed at energy smaller than the gap, \Delta, which is maximum (\Delta ~ 2 meV) for $\textit{n}$ of ~ 10$^{14}$ cm$^{-2}$ and decreases monotonously for larger $\textit{n}$. A perpendicular magnetic field $\textit{B}$ generates states at $E<\Delta$ that fill the gap, but a 20% DOS suppression of superconducting origin unexpectedly persists much above the transport critical field. Conversely, an in-plane field up to 10 T leaves the DOS entirely unchanged. Our measurements exclude that the superconducting state in MoS$_2$ is fully gapped and reveal the presence of a DOS that vanishes linearly with energy, the explanation of which requires going beyond a conventional, purely phonon-driven Bardeen-Cooper-Schrieffer mechanism.