title: On the Estimation of Network Complexity: Dimension of Graphons

abstract: Network complexity has been studied for over half a century and has found a wide range of applications. Many methods have been developed to characterize and estimate the complexity of networks. However, there has been little research with statistical guarantees. In this paper, we develop a statistical theory of graph complexity in a general model of random graphs, the so-called graphon model. Given a graphon, we endow the latent space of the nodes with the neighborhood distance that measures the propensity of two nodes to be connected with similar nodes. Our complexity index is then based on the covering number and the Minkowski dimension of (a purified version of) this metric space. Although the latent space is not identifiable, these indices turn out to be identifiable. This notion of complexity has simple interpretations on popular examples of random graphs: it matches the number of communities in stochastic block models; the dimension of the Euclidean space in random geometric graphs; the regularity of the link function in H\"older graphon models. From a single observation of the graph, we construct an estimator of the neighborhood-distance and show universal non-asymptotic bounds for its risk, matching minimax lower bounds. Based on this estimated distance, we compute the corresponding covering number and Minkowski dimension and we provide optimal non-asymptotic error bounds for these two plug-in estimators.


title: On the vanishing of discrete singular cubical homology for graphs

abstract: We prove that if G is a graph without 3-cycles and 4-cycles, then the discrete cubical homology of G is trivial in dimension d, for all d\ge 2. We also construct a sequence { G_d } of graphs such that this homology is non-trivial in dimension d for d\ge 1. Finally, we show that the discrete cubical homology induced by certain coverings of G equals the ordinary singular homology of a 2-dimensional cell complex built from G, although in general it differs from the discrete cubical homology of the graph as a whole.


title: Dynamic Spatial-Temporal Representation Learning for Traffic Flow Prediction

abstract: As a crucial component in intelligent transportation systems, traffic flow prediction has recently attracted widespread research interest in the field of artificial intelligence (AI) with the increasing availability of massive traffic mobility data. Its key challenge lies in how to integrate diverse factors (such as temporal rules and spatial dependencies) to infer the evolution trend of traffic flow. To address this problem, we propose a unified neural network called Attentive Traffic Flow Machine (ATFM), which can effectively learn the spatial-temporal feature representations of traffic flow with an attention mechanism. In particular, our ATFM is composed of two progressive Convolutional Long Short-Term Memory (ConvLSTM \cite{xingjian2015convolutional}) units connected with a convolutional layer. Specifically, the first ConvLSTM unit takes normal traffic flow features as input and generates a hidden state at each time-step, which is further fed into the connected convolutional layer for spatial attention map inference. The second ConvLSTM unit aims at learning the dynamic spatial-temporal representations from the attentionally weighted traffic flow features. Further, we develop two deep learning frameworks based on ATFM to predict citywide short-term/long-term traffic flow by adaptively incorporating the sequential and periodic data as well as other external influences. Extensive experiments on two standard benchmarks well demonstrate the superiority of the proposed method for traffic flow prediction. Moreover, to verify the generalization of our method, we also apply the customized framework to forecast the passenger pickup/dropoff demands in traffic prediction and show its superior performance. Our code and data are available at {\color{blue}\url{https://github.com/liulingbo918/ATFM}}.


title: Betti realization of varieties defined by formal Laurent series

abstract: We give two constructions of functorial topological realizations for schemes of finite type over the field $\mathbb{C}(\!(t)\!)$ of formal Laurent series with complex coefficients, with values in the homotopy category of spaces over the circle. The problem of constructing such a realization was stated by D. Treumann, motivated by certain questions in mirror symmetry. The first construction uses spreading out and the usual Betti realization over $\mathbb{C}$. The second uses generalized semistable models and log Betti realization defined by Kato and Nakayama, and applies to smooth rigid analytic spaces as well. We provide comparison theorems between the two constructions and relate them to the \'etale homotopy type and de Rham cohomology. As an illustration of the second construction, we treat two examples, the Tate curve and the non-archimedean Hopf surface.


title: Coherence-variance uncertainty relation and coherence cost for quantum measurement under conservation laws

abstract: Uncertainty relations are one of the fundamental principles in physics. It began as a fundamental limitation in quantum mechanics, and today the word {\it uncertainty relation} is a generic term for various trade-off relations in nature. In this letter, we improve the Kennard-Robertson uncertainty relation, and clarify how much coherence we need to implement quantum measurement under conservation laws. Our approach systematically improves and reproduces the previous various refinements of the Kennard-Robertson inequality. As a direct consequence of our inequalities, we improve a well-known limitation of quantum measurements, the Wigner-Araki-Yanase-Ozawa theorem. This improvement gives an asymptotic equality for the necessary and sufficient amount of coherence to implement a quantum measurement with the desired accuracy under conservation laws.


title: Hidden walls: STEM course barriers identified by students with disabilities

abstract: Historically, non-disabled individuals have viewed disability as a personal deficit requiring change to the disabled individual. However, models have emerged from disability activists and disabled intellectuals that emphasize the role of disabling social structures in preventing or hindering equal access across the ability continuum. We used the social relational proposition, which situates disability within the interaction of impairments and particular social structures, to identify disabling structures in introductory STEM courses. We conducted interviews with nine students who identified with a range of impairments about their experiences in introductory STEM courses. We assembled a diverse research team and analyzed the interviews through phenomenological analysis. Participants reported course barriers that prevented effective engagement with course content. These barriers resulted in challenges with time management as well as feelings of stress and anxiety. We discuss recommendations for supporting students to more effectively engage with introductory STEM courses.


title: Analysis of quasielastic $(e,e')$ electromagnetic responses and scaling for nuclear matter in the relativistic mean field model

abstract: A detailed study of the electromagnetic responses for quasielastic $(e,e')$ reactions from nuclear matter is presented within the framework of a relativistic model including momentum dependent scalar and vector mean field potentials in both the initial and final nucleon wave function states. The effects ascribed to the use of energy-dependent potentials are carefully analyzed by comparing their predictions with the responses obtained in the case of constant, i.e., energy-independent, potentials, as well as with the plane wave limit in the final state. The study is extended to the scaling phenomenon. Results are provided for the scaling functions corresponding to different nuclear systems evaluated at several kinematics ranging from intermediate to high values of the transfer momentum. Emphasis is placed on scaling of the first and second kinds and the role played by the relativistic scalar and vector mean field potentials.


title: Analysis of the dissipative range of the energy spectrum in grid turbulence and in direct numerical simulations

abstract: We present a statistical analysis of the behavior of the kinetic energy spectrum in the dissipative range of scales of fully developed three-dimensional turbulence. We analyze spectra recorded in experiments of grid turbulence generated in the Modane wind tunnel, and spectra obtained from high-resolution direct numerical simulations of the forced Navier-Stokes equation. We focus on wave-number scales in the near-dissipative range. We find that in this domain, the spectra behave as a stretched exponential with an exponent $\alpha = 0.68 \pm 0.19$. This result corroborates previous DNS studies which found that the spectrum in the near-dissipative range is best modeled by a stretched exponential with $\alpha<1$. Moreover, we here give an estimate for the actual value of $\alpha$, which turns out to be in precise agreement with the theoretical prediction $\alpha=2/3$ from the non-perturbative renormalisation group.


title: Two Numerical Methods for the 3D Anisotropic Propagation of Galactic Cosmic Rays

abstract: Conventional cosmic-ray propagation models usually assume an isotropic diffusion coefficient to account for the random deflection of cosmic rays by the turbulent interstellar magnetic field. Such a picture is very successful in explaining many observational phenomena related to the propagation of Galactic cosmic rays, such as broken power-law energy spectra, secondary-to-primary ratios, etc. However, the isotropic diffusion presupposition is facing severe challenges from recent observations. In particular, such observations on the large-scale anisotropy of TeV cosmic rays show that the dipole direction differs from the prediction of the conventional model. One possible reason is that the large-scale regular magnetic field, which leads to an anisotropic diffusion of cosmic rays, has not been included in the model provided by the public numerical packages. In this work, we propose two numerical schemes to solve the $3$-dimensional anisotropic transport equation: the pseudo source method and Hundsdorfer-Verwer scheme. Both methods are verified by reproducing the measured B/C and proton spectrum and the radial variation of spectral index expected by former 2D simulation. As a demonstration of the prediction capability, dipole anisotropy is also calculated by a toy simulation with a rough magnetic field.


title: Game Theoretical Approach to Sequential Hypothesis Test with Byzantine Sensors

abstract: In this paper, we consider the problem of sequential binary hypothesis test in adversary environment based on observations from s sensors, with the caveat that a subset of c sensors is compromised by an adversary, whose observations can be manipulated arbitrarily. We choose the asymptotic Average Sample Number (ASN) required to reach a certain level of error probability as the performance metric of the system. The problem is cast as a game between the detector and the adversary, where the detector aims to optimize the system performance while the adversary tries to deteriorate it. We propose a pair of flip attack strategy and voting hypothesis testing rule and prove that they form an equilibrium strategy pair for the game. We further investigate the performance of our proposed detection scheme with unknown number of compromised sensors and corroborate our result with simulation.


title: Accelerating equilibrium isotope effect calculations: II. Stochastic implementation of direct estimators

abstract: Path integral calculations of equilibrium isotope effects and isotopic fractionation are expensive due to the presence of path integral discretization errors, statistical errors, and thermodynamic integration errors. Whereas the discretization errors can be reduced by high-order factorization of the path integral and statistical errors by using centroid virial estimators, two recent papers proposed alternative ways to completely remove the thermodynamic integration errors: Cheng and Ceriotti [J. Chem. Phys. 141, 244112 (2015)] employed a variant of free-energy perturbation called "direct estimators," while Karandashev and Van\'{\i}\v{c}ek [J. Chem. Phys. 143, 194104 (2017)] combined the thermodynamic integration with a stochastic change of mass and piecewise-linear umbrella biasing potential. Here we combine the former approach with the stochastic change of mass in order to decrease its statistical errors when applied to larger isotope effects, and perform a thorough comparison of different methods by computing isotope effects first on a harmonic model, and then on methane and methanium, where we evaluate all isotope effects of the form $\mathrm{CH}_{\mathrm{4-x}}\mathrm{D}_{\mathrm{x}}/\mathrm{CH}_{4}$ and $\mathrm{CH}_{\mathrm{5-x}}\mathrm{D}^{+}_{\mathrm{x}}/\mathrm{CH}^{+}_{5}$, respectively. We discuss thoroughly the reasons for a surprising behavior of the original method of direct estimators, which performed well for a much larger range of isotope effects than what had been expected previously.


title: A graphon counter example

abstract: We give an example of a graphon such that there is no equivalent graphon with a degree function that is (weakly) increasing.


title: Parameter identification for the Landau-Lifshitz-Gilbert equation in Magnetic Particle Imaging

abstract: Magnetic particle imaging (MPI) is a tracer-based technique for medical imaging where the tracer consists of ironoxide nanoparticles. The key idea is to measure the particle response to a temporally changing external magnetic field to compute the spatial concentration of the tracer inside the object. A decent mathematical model demands for a data-driven computation of the system function which does not only describe the measurement geometry but also encodes the interaction of the particles with the external magnetic field. The physical model of this interaction is given by the Landau-Lifshitz-Gilbert (LLG) equation. The determination of the system function can be seen as an inverse problem of its own which can be interpreted as a calibration problem for MPI. In this contribution the calibration problem is formulated as an inverse parameter identification problem for the LLG equation. We give a detailed analysis of the direct as well as the inverse problem in an all-at-once as well as in a reduced setting. The analytical results yield a deeper understanding of inverse problems connected to the LLG equation and provide a starting point for the development of robust numerical solution methods in MPI.


title: A simulation study of methods for handling disease progression in dose-finding clinical trials

abstract: In traditional dose-finding studies, dose-limiting toxicity (DLT) is determined within a fixed time observation window where DLT is often defined as a binary outcome. In the setting of oncology dose-finding trials, often patients in advanced stage of diseases are enrolled. Therefore, disease progression may occur within the DLT observation window leading to treatment discontinuation and rendering the patient unevaluable for DLT assessment. As a result, additional patients have to be enrolled, increasing the sample size. We propose and compare several practical methods for handling disease progression which occurs within the DLT observation window, in the context of the time-to-event continual reassessment method (TITE-CRM) which allows using partial observations. The methods differ on the way they define an evaluable patient and in the way incomplete observations are included. The methods are illustrated and contrasted in the context of a single simulated trial, and compared via simulations under various scenarios of dose-progression relationship, in the setting of advanced soft-tissue sarcoma.


title: Blockchain Technologies for Smart Energy Systems: Fundamentals, Challenges and Solutions

abstract: In this paper, we discuss the integration of blockchain in smart energy systems. We present various blockchain technology solutions, review important blockchain platforms, and several blockchain based smart energy projects in different smart energy domains. The majority of blockchain platforms with embedded combination of blockchain technology solutions are computing- and resource- intensive, and hence not entirely suitable for smart energy applications. We consider the requirements of smart energy systems and accordingly identify appropriate blockchain technology solutions for smart energy applications. Our analysis can help in the development of flexible blockchain platforms for smart energy systems.


title: Strong suppression of the resistivity near the transition to superconductivity in narrow micro-bridges in external magnetic fields

abstract: We have investigated a series of superconducting bridges based on homogeneous amorphous WSi and MoSi films, with bridge widths w ranging from 2 um to 1000 um and film thicknesses d ~ 4-6 nm and 100 nm. Upon decreasing the bridge widths below the respective Pearl lengths, we observe in all cases distinct changes in the characteristics of the resistive transitions to superconductivity. For each of the films, the resistivity curves R(B,T) separate at a well-defined and field-dependent temperature T*(B) with decreasing the temperature, resulting in a dramatic suppression of the resistivity and a sharpening of the transitions with decreasing bridge width w. The associated excess conductivity in all the bridges scales as 1/w, which may suggest the presence of a highly conducting region that is dominating the electric transport in narrow bridges. We argue that this effect can only be observed in materials with sufficiently weak vortex pinning.


title: Optimal stopping times for a class of Ito diffusion bridges

abstract: The scope of this paper is to study the optimal stopping problems associated to a stochastic process, which may represent the gain of an investment, for which information on the final value is available a priori. This information may proceed, for example, from insider trading or from pinning at expiration of stock options. We solve and provide explicit solutions to these optimization problems. As special case, we discuss different processes whose optimal barrier has the same shape as the optimal barrier of the Brownian bridge. So doing we provide a catalogue of alternatives to the Brownian bridge which in practice could be better adapted to the data. Moreover, we investigate if, for any given (decreasing) curve, there exists a process with this curve as optimal barrier. This provides a model for the optimal liquidation time, i.e. the optimal time at which the investor should liquidate a position in order to maximize the gain.


title: Une construction d'extensions faiblement non ramifi\'ees d'un anneau de valuation

abstract: \'Etant donn\'e un anneau de valuation $V$, de corps r\'esiduel $F$ et de groupe des valeurs $\Gamma$, on donne une condition suffisante pour qu'un anneau local dominant $V$ soit un anneau de valuation de groupe $\Gamma$. Lorsque $V$ contient un corps $k$, ce r\'esultat est appliqu\'e \`a la construction d'un anneau de valuation contenant $V$ et une extension donn\'ee $k'$ de $k$, de groupe $\Gamma$ et de corps r\'esiduel engendr\'e par $k'$ et $F$. Cela s'av\`ere possible, notamment, lorsque $k'$ ou $F$ est s\'eparable sur $k$. Given a valuation ring $V$, with residue field $F$ and value group $\Gamma$, we give a sufficient condition for a local ring dominating $V$ to be a valuation ring with the same value group. When $V$ contains a field $k$, we apply this result to the problem of constructing a valuation ring $W$ containing $V$ and a prescribed extension $k'$ of $k$, with value group $\Gamma$ and residue field generated by $k'$ and $F$; this is possible in particular when either $k'$ or $F$ is separable over $k$.


title: Blackbox Attacks on Reinforcement Learning Agents Using Approximated Temporal Information

abstract: Recent research on reinforcement learning (RL) has suggested that trained agents are vulnerable to maliciously crafted adversarial samples. In this work, we show how such samples can be generalised from White-box and Grey-box attacks to a strong Black-box case, where the attacker has no knowledge of the agents, their training parameters and their training methods. We use sequence-to-sequence models to predict a single action or a sequence of future actions that a trained agent will make. First, we show our approximation model, based on time-series information from the agent, consistently predicts RL agents' future actions with high accuracy in a Black-box setup on a wide range of games and RL algorithms. Second, we find that although adversarial samples are transferable from the target model to our RL agents, they often outperform random Gaussian noise only marginally. This highlights a serious methodological deficiency in previous work on such agents; random jamming should have been taken as the baseline for evaluation. Third, we propose a novel use for adversarial samplesin Black-box attacks of RL agents: they can be used to trigger a trained agent to misbehave after a specific time delay. This appears to be a genuinely new type of attack. It potentially enables an attacker to use devices controlled by RL agents as time bombs.


title: Active binary mixtures of fast and slow hard spheres

abstract: We computationally studied the phase behavior and dynamics of binary mixtures of active particles, where each 'species' had distinct activities leading to distinct velocities, fast and slow. We obtained phase diagrams demonstrating motility-induced phase separation (MIPS) upon varying the activity and concentration of each species, and extended current kinetic theory of active/passive mixtures to active/active mixtures. We discovered two regimes of behavior quantified through the participation of each species in the dense phase compared to their monodisperse counterparts. In regime I (active/passive and active/weakly-active), we found that the dense phase was segregated by particle type into domains of fast and slow particles. Moreover, fast particles were suppressed from entering the dense phase while slow particles were enhanced entering the dense phase, compared to monodisperse systems of all-fast or all-slow particles. These effects decayed asymptotically as the activity of the slow species increased, approaching the activity of the fast species until they were negligible (regime II). In regime II, the dense phase was homogeneously mixed and each species participated in the dense phase as if it were it a monodisperse system; each species behaved as if it weren't mixed at all. Finally, we collapsed our data defining two quantities that measure the total activity of the system. We thus found expressions that predict, a priori, the percentage of each particle type that participates in the dense phase through MIPS.


title: IB-homogeneous graphs

abstract: The Lachlan-Woodrow Theorem identifies ultrahomogeneous graphs up to isomorphism. Recently, the present author and D. Hartman classified MB-homogeneous graphs up to bimorphism-equivalence. We extend those results in this paper, showing that every IB-homogeneous graph is either ultrahomogeneous or MB-homogeneous, and thus all the IB-homogeneous graphs are known up to bimorphism-equivalence.


title: Odd-frequency superconductivity and Meissner effect in the doped topological insulator Bi$_2$Se$_3$

abstract: Doped Bi$_2$Se$_3$ is proposed to be a nematic superconductor originating from unusual inter-orbital pairing. We classify all induced superconducting pair correlations in Bi$_2$Se$_3$ and discover that intra-orbital odd-frequency pairing dominates over a significant range of frequencies. Moreover, we explore the contributions of even- and odd-frequency pairing to the Meissner effect, including separating intra- and inter-band processes in the response function. Contrary to expectations, and due to inter-band contributions, we find a diamagnetic odd-frequency Meissner effect.


title: Dynamics for discrete subgroups of $\operatorname{SL}_2(\mathbb C)$

abstract: Margulis wrote in the preface of his book Discrete subgroups of semisimple Lie groups that "A number of important topics have been omitted. The most significant of these is the theory of Kleinian groups and Thurston's theory of 3-dimensional manifolds: these two theories can be united under the common title of Theory of discrete subgroups of $\operatorname{SL}_2(\mathbb C)$." In this article, we will discuss a few recent advances regarding this missing topic from his book, which were influenced by his earlier works.


title: Data Driven Vulnerability Exploration for Design Phase System Analysis

abstract: Applying security as a lifecycle practice is becoming increasingly important to combat targeted attacks in safety-critical systems. Among others there are two significant challenges in this area: (1) the need for models that can characterize a realistic system in the absence of an implementation and (2) an automated way to associate attack vector information; that is, historical data, to such system models. We propose the cybersecurity body of knowledge (CYBOK), which takes in sufficiently characteristic models of systems and acts as a search engine for potential attack vectors. CYBOK is fundamentally an algorithmic approach to vulnerability exploration, which is a significant extension to the body of knowledge it builds upon. By using CYBOK, security analysts and system designers can work together to assess the overall security posture of systems early in their lifecycle, during major design decisions and before final product designs. Consequently, assisting in applying security earlier and throughout the systems lifecycle.


title: An Edge Computing Robot Experience for Automatic Elderly Mental Health Care Based on Voice

abstract: We need open platforms driven by specialists, in which queries can be created and collected for long periods and the diagnosis made, based on a rigorous clinical follow-up. In this work, we developed a multi-language robot interface helping to evaluate the mental health of seniors by interacting through questions. The specialist can propose questions, as well as to receive users' answers, in text form. The robot can automatically interact with the user using the appropriate language. It can process the answers and under the guidance of a specialist, questions and answers can be oriented towards the desired therapy direction. The prototype, was implemented on an embedded device meant for edge computing, thus it is able to filter environmental noise and can be placed anywhere at home. The experience is now available for specialists to create queries and answers through a Web-based interface.


title: Conditional propagation of chaos for mean field systems of interacting neurons

abstract: We study the stochastic system of interacting neurons introduced in De Masi et al. (2015) and in Fournier and L\"ocherbach (2016) in a diffusive scaling. The system consists of $N$ neurons, each spiking randomly with rate depending on its membrane potential. At its spiking time, the potential of the spiking neuron is reset to $0$ and all other neurons receive an additional amount of potential which is a centred random variable of order $ 1 / \sqrt{N}.$ In between successive spikes, each neuron's potential follows a deterministic flow. We prove the convergence of the system, as $N \to \infty$, to a limit nonlinear jumping stochastic differential equation driven by Poisson random measure and an additional Brownian motion $W$ which is created by the central limit theorem. This Brownian motion is underlying each particle's motion and induces a common noise factor for all neurons in the limit system. Conditionally on $W,$ the different neurons are independent in the limit system. This is the {\it conditional propagation of chaos} property. We prove the well-posedness of the limit equation by adapting the ideas of Graham (1992) to our frame. To prove the convergence in distribution of the finite system to the limit system, we introduce a new martingale problem that is well suited for our framework. The uniqueness of the limit is deduced from the exchangeability of the underlying system.


title: Macrospin analysis of RF excitations within fully perpendicular magnetic tunnel junctions with second order easy-axis magnetic anisotropy contribution

abstract: The conditions of field and voltage for inducing steady state excitations in fully perpendicular magnetic tunnel junctions (pMTJs), adapted for memory applications, were numerically investigated by the resolution of the Landau-Lifshitz-Gilbert equation in the macrospin approach. Both damping-like and the field-like spin transfer torque terms were taken into account in the simulations, as well as the contribution of the second order uniaxial anisotropy term (K2), which has been recently revealed in MgO-based pMTJs. An in-plane applied magnetic field balances the out of plane symmetry of the pMTJ and allows the signal detection. Using this model, we assessed the states of the free layer magnetization as a function of strength of K2 and polar theta_H angle of the applied field (varied from 90 to 60 deg.). There are two stable states, with the magnetization in-plane or out of plane of the layer, and two dynamic states with self-sustained oscillations, called in-plane precession state (IPP) or out of plane precession state (OPP). The IPP mode, with oscillation frequencies up to 7 GHz, appears only for positive voltages if theta_H = 90 deg. However, it shows a more complex distribution when the field is slightly tilted out of plane. The OPP mode is excited only if K2 is considered and reaches a maximum oscillation frequency of 15 GHz. Large areas of dynamic states with high frequencies are obtained for strong values of the field-like torque and K2, when applying a slightly tilted external field toward the out of plane direction. The non-zero temperature does not modify the phase diagrams, but reduces drastically the power spectral density peak amplitudes.


title: Discrete Element Method Model of Elastic Fiber Uniaxial Compression

abstract: A flexible fiber model based on the discrete element method (DEM) is presented and validated for the simulation of uniaxial compression of flexible fibers in a cylindrical container. It is found that the contact force models in the DEM simulations have a significant impact on compressive forces exerted on the fiber bed. Only when the geometry-dependent normal contact force model and the static friction model are employed, the simulation results are in good agreement with experimental results. Systematic simulation studies show that the compressive force initially increases and eventually saturates with an increase in the fiber-fiber friction coefficient, and the fiber-fiber contact forces follow a similar trend. The compressive force and lateral shear-to-normal stress ratio increase linearly with increasing fiber-wall friction coefficient. In uniaxial compression of frictional fibers, more static friction contacts occur than dynamic friction contacts with static friction becoming more predominant as the fiber-fiber friction coefficient increases.


title: Oceanographic Considerations for Exoplanet Life Detection

abstract: Liquid water oceans are at the center of our search for life on exoplanets because water is a strict requirement for life as we know it. However, oceans are dynamic habitats---and some oceans may be better hosts for life than others. In Earth's ocean, circulation transports essential nutrients such as phosphate and is a first-order control on the distribution and productivity of life. Of particular importance is upward flow from the dark depths of the ocean in response to wind-driven divergence in surface layers. This `upwelling' returns essential nutrients that tend to accumulate at depth via sinking of organic particulates back to the sunlit regions where photosynthetic life thrives. Ocean dynamics are likely to impose constraints on the activity and atmospheric expression of photosynthetic life in exo-oceans as well, but we lack an understanding of how ocean dynamics may differ on other planets. We address this issue by exploring the sensitivity of ocean dynamics to a suite of planetary parameters using ROCKE-3D, a fully coupled ocean-atmosphere GCM. Our results suggest that planets that rotate slower and have higher surface pressure than Earth may be the most attractive targets for remote life detection because upwelling is enhanced under these conditions, resulting in greater nutrient supply to the surface biosphere. Seasonal deepening of the mixed layer on high obliquity planets may also enhance nutrient replenishment from depth into the surface mixed layer. Efficient nutrient recycling favors greater biological activity, more biosignature production, and thus more detectable life. More generally, our results demonstrate the importance of considering oceanographic phenomena for exoplanet life detection and motivate future interdisciplinary contributions to the emerging field of exo-oceanography.


title: BNB autoregressions for modeling integer-valued time series with extreme observations

abstract: This article introduces a general class of heavy-tailed autoregressions for modeling integer-valued time series with outliers. The proposed specification is based on a heavy-tailed mixture of negative binomial distributions that features an observation-driven dynamic equation for the conditional expectation. The existence of a unique stationary and ergodic solution for the class of autoregressive processes is shown under a general contraction condition. The estimation of the model can be easily performed by Maximum Likelihood given the closed form of the likelihood function. The strong consistency and the asymptotic normality of the estimator are formally derived. Two examples of specifications illustrate the flexibility of the approach and the relevance of the theoretical results. In particular, a linear dynamic equation and a score-driven equation for the conditional expectation are considered. The score-driven specification is shown to be particularly appealing as it delivers a robust filtering method that attenuates the impact of outliers. An empirical application to the time series of narcotics trafficking reports in Sydney illustrates the effectiveness of the method in handling extreme observations.


title: Structured Query Construction via Knowledge Graph Embedding

abstract: In order to facilitate the accesses of general users to knowledge graphs, an increasing effort is being exerted to construct graph-structured queries of given natural language questions. At the core of the construction is to deduce the structure of the target query and determine the vertices/edges which constitute the query. Existing query construction methods rely on question understanding and conventional graph-based algorithms which lead to inefficient and degraded performances facing complex natural language questions over knowledge graphs with large scales. In this paper, we focus on this problem and propose a novel framework standing on recent knowledge graph embedding techniques. Our framework first encodes the underlying knowledge graph into a low-dimensional embedding space by leveraging generalized local knowledge graphs. Given a natural language question, the learned embedding representations of the knowledge graph are utilized to compute the query structure and assemble vertices/edges into the target query. Extensive experiments were conducted on the benchmark dataset, and the results demonstrate that our framework outperforms state-of-the-art baseline models regarding effectiveness and efficiency.


title: Integrated Information in the Spiking-Bursting Stochastic Model

abstract: This study presents a comprehensive analytic description in terms of the empirical "whole minus sum" version of Integrated Information in comparison to the "decoder based" version for the "spiking-bursting" discrete-time, discrete-state stochastic model, which was recently introduced to describe a specific type of dynamics in a neuron-astrocyte network. The "whole minus sum" information may change sign, and an interpretation of this transition in terms of "net synergy" is available in the literature. This motivates our particular interest to the sign of the "whole minus sum" information in our analytical consideration. The behavior of the "whole minus sum" and "decoder based" information measures are found to bear a lot of similarity, showing their mutual asymptotic convergence as time-uncorrelated activity is increased, with the sign transition of the "whole minus sum" information associated to a rapid growth in the "decoder based" information. The study aims at creating a theoretical base for using the spiking-bursting model as a well understood reference point for applying Integrated Information concepts to systems exhibiting similar bursting behavior (in particular, to neuron-astrocyte networks). The model can also be of interest as a new discrete-state test bench for different formulations of Integrated Information.


title: A Note on Lipschitz Continuiy of the Solutions of a Class of Elliptic Free Boundary Problems

abstract: We provide a new and simple proof based on Harnack's inequality to the Lipschitz continuity of the solutions of a class of free boundary problems.


title: AR-based interaction for safe human-robot collaborative manufacturing

abstract: Industrial standards define safety requirements for Human-Robot Collaboration (HRC) in industrial manufacturing. The standards particularly require real-time monitoring and securing of the minimum protective distance between a robot and an operator. In this work, we propose a depth-sensor based model for workspace monitoring and an interactive Augmented Reality (AR) User Interface (UI) for safe HRC. The AR UI is implemented on two different hardware: a projector-mirror setup anda wearable AR gear (HoloLens). We experiment the workspace model and UIs for a realistic diesel motor assembly task. The AR-based interactive UIs provide 21-24% and 57-64% reduction in the task completion and robot idle time, respectively, as compared to a baseline without interaction and workspace sharing. However, subjective evaluations reveal that HoloLens based AR is not yet suitable for industrial manufacturing while the projector-mirror setup shows clear improvements in safety and work ergonomics.


title: Elliptic quasi-variational inequalities under a smallness assumption: Uniqueness, differential stability and optimal control

abstract: We consider a quasi-variational inequality governed by a moving set. We employ the assumption that the movement of the set has a small Lipschitz constant. Under this requirement, we show that the quasi-variational inequality has a unique solution which depends Lipschitz-continuously on the source term. If the data of the problem is (directionally) differentiable, the solution map is directionally differentiable as well. We also study the optimal control of the quasi-variational inequality and provide necessary optimality conditions of strongly stationary type.


title: Franck-Condon factors via compressive sensing

abstract: Probabilities of vibronic transitions in molecules are referred to as Franck-Condon factors (FCFs). Although several approaches for calculating FCFs have been developed, such calculations are still challenging. Recently it was shown that there exists a correspondence between the problem of calculating FCFs and boson sampling. However, if the output photon number distribution of boson sampling is sparse then it can be classically simulated. Exploiting these results, we develop a method to approximately reconstruct the distribution of FCFs of certain molecules. We demonstrate this method by applying it to formic acid and thymine at $0$ K. In our method, we first obtain the marginal photon number distributions for pairs of modes of a Gaussian state associated with the molecular transition. We then apply a compressive sensing method called polynomial-time matching pursuit to recover FCFs.


title: Schwarz type lemmas for pseudo-Hermitian manifolds

abstract: In this paper, we consider some generalized holomorphic maps between pseudo-Hermitian manifolds. These maps include the \emph{CR} maps and the transversally holomorphic maps. In terms of some sub-Laplacian or Hessian type Bochner formulas, and comparison theorems in the pseudo-Hermitian version, we are able to establish several Schwarz type results for both the \emph{CR} maps and the transversally holomorphic maps between pseudo-Hermitian manifolds. Finally, we also discuss the \emph{CR} hyperbolicity problem for pseudo-Hermitian manifolds.


title: Introducing Open boundary conditions in modeling nonperiodic materials and interfaces: the impact of the periodic assumption

abstract: Simulations are essential to accelerate the discovery of new materials and to gain full understanding of known ones. Although hard to realize experimentally, periodic boundary conditions are omnipresent in material simulations. In this work, we intro-duce ROBIN (recursive open boundary and interfaces), the first method allowing open boundary conditions in material and interface modeling. The computational costs are limited to solving quantum properties in a focus area which allows explicitly discretizing millions of atoms in real space and to consider virtually any type of environment (be it periodic, regular, or ran-dom). The impact of the periodicity assumption is assessed in detail with silicon dopants in graphene. Graphene was con-firmed to produce a band gap with periodic substitution of 3% carbon with silicon in agreement with published periodic boundary condition calculations. Instead, 3% randomly distributed silicon in graphene only shifts the energy spectrum. The predicted shift agrees quantitatively with published experimental data. Key insight of this assessment is, assuming periodici-ty elevates a small perturbation of a periodic cell into a strong impact on the material property prediction. Periodic boundary conditions can be applied on truly periodic systems only. More general systems should apply an open boundary method for reliable predictions.


title: Dynamics and escape of active particles in a harmonic trap

abstract: The dynamics of active particles is of interest at many levels and is the focus of theoretical and experimental research. There have been many attempts to describe the dynamics of particles affected by random active forces in terms of an effective temperature. This kind of description is tempting due to the similarities (or lack thereof) with systems in or near thermal equilibrium. However, the generality and validity of the effective temperature is not yet fully understood. Here, we studied the dynamics of trapped particles subjected to both thermal and active forces. The particles were not overdamped. Expressions for the effective temperature due to the potential and kinetic energies were derived, and they differ from each other. A third possible effective temperature can be derived from the escape time of the particle from the trap, using a Kramers-like expression for the mean escape time. We found that over a large fraction of the parameter space, the potential energy effective temperature is in agreement with the escape temperature, while the kinetic effective temperature only agrees with the former two in the overdamped limit. Moreover, we show that the specific implementation of the random active force, and not only its first two moments and the two point auto-correlation function, affects the escape time distribution.


title: Optimizing Generalized Rate Metrics through Game Equilibrium

abstract: We present a general framework for solving a large class of learning problems with non-linear functions of classification rates. This includes problems where one wishes to optimize a non-decomposable performance metric such as the F-measure or G-mean, and constrained training problems where the classifier needs to satisfy non-linear rate constraints such as predictive parity fairness, distribution divergences or churn ratios. We extend previous two-player game approaches for constrained optimization to a game between three players to decouple the classifier rates from the non-linear objective, and seek to find an equilibrium of the game. Our approach generalizes many existing algorithms, and makes possible new algorithms with more flexibility and tighter handling of non-linear rate constraints. We provide convergence guarantees for convex functions of rates, and show how our methodology can be extended to handle sums of ratios of rates. Experiments on different fairness tasks confirm the efficacy of our approach.


title: Reinforcement Learning for Joint Optimization of Multiple Rewards

abstract: Finding optimal policies which maximize long term rewards of Markov Decision Processes requires the use of dynamic programming and backward induction to solve the Bellman optimality equation. However, many real-world problems require optimization of an objective that is non-linear in cumulative rewards for which dynamic programming cannot be applied directly. For example, in a resource allocation problem, one of the objectives is to maximize long-term fairness among the users. We notice that when an agent aim to optimize some function of the sum of rewards is considered, the problem loses its Markov nature. This paper addresses and formalizes the problem of optimizing a non-linear function of the long term average of rewards. We propose model-based and model-free algorithms to learn the policy, where the model-based policy is shown to achieve a regret of $\Tilde{O}\left(LKDS\sqrt{\frac{A}{T}}\right)$ for $K$ objectives combined with a concave $L$-Lipschitz function. Further, using the fairness in cellular base-station scheduling, and queueing system scheduling as examples, the proposed algorithm is shown to significantly outperform the conventional RL approaches.

categories: cs.LG cs.AI cs.GT cs.IT cs.MA math.IT stat.ML


title: Quantum marginal problem and incompatibility

abstract: One of the basic distinctions between classical and quantum mechanics is the existence of fundamentally incompatible quantities. Such quantities are present on all levels of quantum objects: states, measurements, quantum channels, and even higher order dynamics. In this manuscript, we show that two seemingly different aspects of quantum incompatibility: the quantum marginal problem of states and the incompatibility on the level of quantum channels are in many-to-one correspondence. Importantly, as incompatibility of measurements is a special case of the latter, it also forms an instance of the quantum marginal problem. The generality of the connection is harnessed by solving the marginal problem for Gaussian and Bell diagonal states, as well as for pure states under depolarizing noise. Furthermore, we derive entropic criteria for channel compatibility, and develop a converging hierarchy of semi-definite programs for quantifying the strength of quantum memories.


title: A refinement of Christol's theorem for algebraic power series

abstract: A famous result of Christol gives that a power series $F(t)=\sum_{n\ge 0} f(n)t^n$ with coefficients in a finite field $\mathbb{F}_q$ of characteristic $p$ is algebraic over the field of rational functions in $t$ if and only if there is a finite-state automaton accepting the base-$p$ digits of $n$ as input and giving $f(n)$ as output for every $n\ge 0$. An extension of Christol's theorem, giving a complete description of the algebraic closure of $\mathbb{F}_q(t)$, was later given by Kedlaya. When one looks at the support of an algebraic power series, that is the set of $n$ for which $f(n)\neq 0$, a well-known dichotomy for sets generated by finite-state automata shows that the support set is either sparse---with the number of $n\le x$ for which $f(n)\neq 0$ bounded by a polynomial in $\log(x)$---or it is reasonably large in the sense that the number of $n\le x$ with $f(n)\neq 0$ grows faster than $x^{\alpha}$ for some positive $\alpha$. The collection of algebraic power series with sparse supports forms a ring and we give a purely algebraic characterization of this ring in terms of Artin-Schreier extensions and we extend this to the context of Kedlaya's work on generalized power series.


title: Alternating multiple zeta values, and explicit formulas of some Euler-Apery-type series

abstract: In this paper, we study some Euler-Ap\'ery-type series which involve central binomial coefficients and (generalized) harmonic numbers. In particular, we establish elegant explicit formulas of some series by iterated integrals and alternating multiple zeta values. Based on these formulas, we further show that some other series are reducible to ln(2), zeta values, and alternating multiple zeta values by considering the contour integrals related to gamma functions, polygamma functions and trigonometric functions. The evaluations of a large number of special Euler-Ap\'ery-type series are presented as examples.


title: Cooperative regulation of cellular identity in systems with intercellular communication defects

abstract: The cooperative dynamics of cellular populations emerging from the underlying interactions determines cellular functions, and thereby their identity in tissues. Global deviations from this dynamics on the other hand reflects pathological conditions. However, how these conditions are stabilized from dis-regulation on the level of the single entities is still unclear. Here we tackle this question using generic Hodgkin-Huxley type of models that describe physiological bursting dynamics of pancreatic beta-cells, and introduce channel dis-function to mimic pathological silent dynamics. The probability for pathological behavior in beta-cell populations is $\sim100\%$ when all cell have these defects, despite the negligible size of the silent state basin of attraction for single cells. In stark contrast, in a more realistic scenario for a heterogeneous population, stabilization of the pathological state depends on the size of the sub-population which acquired the defects. However, the probability to exhibit stable pathological dynamics in this case is less than $\sim10\%$. These results therefore suggest that the physiological bursting dynamics of a population of beta-cells is cooperatively regulated, even under intercellular communication defects induced by dis-functional channels of single-cells.


title: Encoders and Decoders for Quantum Expander Codes Using Machine Learning

abstract: Quantum key distribution (QKD) allows two distant parties to share encryption keys with security based on laws of quantum mechanics. In order to share the keys, the quantum bits have to be transmitted from the sender to the receiver over a noisy quantum channel. In order to transmit this information, efficient encoders and decoders need to be designed. However, large-scale design of quantum encoders and decoders have to depend on the channel characteristics and require look-up tables which require memory that is exponential in the number of qubits. In order to alleviate that, this paper aims to design the quantum encoders and decoders for expander codes by adapting techniques from machine learning including reinforcement learning and neural networks to the quantum domain. The proposed quantum decoder trains a neural network which is trained using the maximum aposteriori error for the syndromes, eliminating the use of large lookup tables. The quantum encoder uses deep Q-learning based techniques to optimize the generator matrices in the quantum Calderbank-Shor-Steane (CSS) codes. The evaluation results demonstrate improved performance of the proposed quantum encoder and decoder designs as compared to the quantum expander codes.


title: Astroalign: A Python module for astronomical image registration

abstract: We present an algorithm implemented in the astroalign Python module for image registration in astronomy. Our module does not rely on WCS information and instead matches 3-point asterisms (triangles) on the images to find the most accurate linear transformation between the two. It is especially useful in the context of aligning images prior to stacking or performing difference image analysis. Astroalign can match images of different point-spread functions, seeing, and atmospheric conditions.


title: Stable fixed points of combinatorial threshold-linear networks

abstract: Combinatorial threshold-linear networks (CTLNs) are a special class of recurrent neural networks whose dynamics are tightly controlled by an underlying directed graph. Recurrent networks have long been used as models for associative memory and pattern completion, with stable fixed points playing the role of stored memory patterns in the network. In prior work, we showed that target-free cliques of the graph correspond to stable fixed points of the dynamics, and we conjectured that these are the only stable fixed points possible. In this paper, we prove that the conjecture holds in a variety of special cases, including for networks with very strong inhibition and graphs of size $n \leq 4$. We also provide further evidence for the conjecture by showing that sparse graphs and graphs that are nearly cliques can never support stable fixed points. Finally, we translate some results from extremal combinatorics to obtain an upper bound on the number of stable fixed points of CTLNs in cases where the conjecture holds.


title: Elliptic solutions of Boussinesq type lattice equations and the elliptic Nth root of unity

abstract: We establish an infinite family of solutions in terms of elliptic functions of the lattice Boussinesq systems by setting up a direct linearisation scheme, which provides the solution structure for those equations in the elliptic case. The latter, which contains as main structural element a Cauchy kernel on the torus, is obtained from a dimensional reduction of the elliptic direct linearisation scheme of the lattice Kadomtsev-Petviashvili equation, which requires the introduction of a novel technical concept, namely the "elliptic cube root of unity". Thus, in order to implement the reduction we define, more generally, the notion of {\em elliptic $N^{\rm th}$ root of unity}, and discuss some of its properties in connection with a special class of elliptic addition formulae. As a particular concrete application we present the class of elliptic $N$-soliton solutions of the lattice Boussinesq systems.


title: Species subsets and embedded networks of S-systems

abstract: Magombedze and Mulder (2013) studied the gene regulatory system of \textit{Mycobacterium Tuberculosis} (\textit{Mtb}) by partitioning this into three subsystems based on putative gene function and role in dormancy/latency development. Each subsystem, in the form of $S$-system, is represented by an embedded chemical reaction network (CRN), defined by a species subset and a reaction subset induced by the set of digraph vertices of the subsystem. Based on the network decomposition theory initiated by Feinberg in 1987, we have introduced the concept of incidence-independent and developed the theory of $\mathscr{C}$- and $\mathscr{C}^*$-decompositions including their structure theorems in terms of linkage classes. With the $S$-system CRN $\mathscr{N}$ of Magombedze and Mulder's \textit{Mtb} model, its reaction set partition induced decomposition of subnetworks that are not CRNs of $S$-system but constitute independent decomposition of $\mathscr{N}$. We have also constructed a new $S$-system CRN $\mathscr{N}^*$ for which the embedded networks are $\mathscr{C}^*$-decomposition. We have shown that subnetworks of $\mathscr{N}$ and the embedded networks (subnetworks of $\mathscr{N}^*$) are digraph homomorphisms. Lastly, we attempted to explore modularity in the context of CRN.


title: Supervised Multimodal Bitransformers for Classifying Images and Text

abstract: Self-supervised bidirectional transformer models such as BERT have led to dramatic improvements in a wide variety of textual classification tasks. The modern digital world is increasingly multimodal, however, and textual information is often accompanied by other modalities such as images. We introduce a supervised multimodal bitransformer model that fuses information from text and image encoders, and obtain state-of-the-art performance on various multimodal classification benchmark tasks, outperforming strong baselines, including on hard test sets specifically designed to measure multimodal performance.


title: Collinear Factorization in Wide-Angle Hadron Pair Production in $e^+e^-$ Annihilation

abstract: We compute the inclusive unpolarized dihadron production cross section in the far from back-to-back region of $e^+ e^-$ annihilation in leading order pQCD using existing fragmentation function fits and standard collinear factorization, focusing on the large transverse momentum region where transverse momentum is comparable to the hard scale (the center-of-mass energy). We compare with standard transverse-momentum-dependent (TMD) fragmentation function-based predictions intended for the small transverse momentum region with the aim of testing the expectation that the two types of calculation roughly coincide at intermediate transverse momentum. We find significant tension, within the intermediate transverse momentum region, between calculations done with existing non-perturbative TMD fragmentation functions and collinear factorization calculations if the center-of-mass energy is not extremely large. We argue that $e^+ e^-$ measurements are ideal for resolving this tension and exploring the large-to-small transverse momentum transition, given the typically larger hard scales ($\gtrsim 10$ GeV) of the process as compared with similar scenarios that arise in semi-inclusive deep inelastic scattering and fixed-target Drell-Yan measurements.


title: Transcendental Galois Theory and Noether's Problem

abstract: In 1918, Noether published a paper where she studied such a problem, now called Noether's problem on rationality: Let $L=K\left( t_{1},t_{2},\cdots ,t_{n}\right) $ be a purely transcendental extension over a field $K$ and $G$ a finite subgroup acting transitively on $t_{1},t_{2},\cdots ,t_{n}$ in an evident manner. Is it true that the invariant subfield $L^{G}$ of $L$ under $% G$ is still purely transcendental over $K$? The problem has been open in general except for minor particular cases. In this paper we will attempt to understand a general theory for Noether's problem on rationality by transcendental Galois theory. Then new particular cases will be obtained. We will also give a generalization for the remarkable counter-example given by Swan in 1969.


title: Unsupervised Clustering of Quantitative Imaging Phenotypes using Autoencoder and Gaussian Mixture Model

abstract: Quantitative medical image computing (radiomics) has been widely applied to build prediction models from medical images. However, overfitting is a significant issue in conventional radiomics, where a large number of radiomic features are directly used to train and test models that predict genotypes or clinical outcomes. In order to tackle this problem, we propose an unsupervised learning pipeline composed of an autoencoder for representation learning of radiomic features and a Gaussian mixture model based on minimum message length criterion for clustering. By incorporating probabilistic modeling, disease heterogeneity has been taken into account. The performance of the proposed pipeline was evaluated on an institutional MRI cohort of 108 patients with colorectal cancer liver metastases. Our approach is capable of automatically selecting the optimal number of clusters and assigns patients into clusters (imaging subtypes) with significantly different survival rates. Our method outperforms other unsupervised clustering methods that have been used for radiomics analysis and has comparable performance to a state-of-the-art imaging biomarker.


title: Mendeley Reader Counts for US Computer Science Conference Papers and Journal articles

abstract: Although bibliometrics are normally applied to journal articles when used to support research evaluations, conference papers are at least as important in fast-moving computing-related fields. It is therefore important to assess the relative advantages of citations and altmetrics for computing conference papers to make an informed decision about which, if any, to use. This paper compares Scopus citations with Mendeley reader counts for conference papers and journal articles that were published between 1996 and 2018 in 11 computing fields and had at least one US author. The data showed high correlations between Scopus citation counts and Mendeley reader counts in all fields and most years, but with few Mendeley readers for older conference papers and few Scopus citations for new conference papers and journal articles. The results therefore suggest that Mendeley reader counts have a substantial advantage over citation counts for recently-published conference papers due to their greater speed, but are unsuitable for older conference papers.


title: Extracting and Learning a Dependency-Enhanced Type Lexicon for Dutch

abstract: This thesis is concerned with type-logical grammars and their practical applicability as tools of reasoning about sentence syntax and semantics. The focal point is narrowed to Dutch, a language exhibiting a large degree of word order variability. In order to overcome difficulties arising as a result of that variability, the thesis explores and expands upon a type grammar based on Multiplicative Intuitionistic Linear Logic, agnostic to word order but enriched with decorations that aim to reduce its proof-theoretic complexity. An algorithm for the conversion of dependency-annotated sentences into type sequences is then implemented, populating the type logic with concrete, data-driven lexical types. Two experiments are ran on the resulting grammar instantiation. The first pertains to the learnability of the type-assignment process by a neural architecture. A novel application of a self-attentive sequence transduction model is proposed; contrary to established practices, it constructs types inductively by internalizing the type-formation syntax, thus exhibiting generalizability beyond a pre-specified type vocabulary. The second revolves around a deductive parsing system that can resolve structural ambiguities by consulting both word and type information; preliminary results suggest both excellent computational efficiency and performance.


title: Archetypical "push the band critical point" mechanism for peaking of the density of states in three-dimensional crystals: Theory and case study of cubic H$_{3}$S

abstract: The point of zero gradient of the electronic band structure--critical point--generally induces the singularity in the density of states (DOS), but no isolated critical point yields strict divergence of the DOS in three dimensions, differently from the lower dimensional cases. In view of the band structure as a smooth hypersurface on the reciprocal space, we discuss the minimal deformation of the band structure that yields non-divergent but large sharp DOS peaks in three dimensions. By "pushing down" the energy level at the second order saddle point (maximum), a continuous closed loop of saddle points (sphere of maxima) encircling the original position of the saddle point (maximum) emerges, with which the DOS peak is formed. Being high dimensional features, the saddle loop and extremum shell thus formed are difficult to locate with standard band structure analysis on linear ${\bf k}$-point paths. The Lifshitz transition occurring over a linear or planar manifold is discussed as an indicator of such features. We also find that the celebrated DOS peak in the recently discovered superconducting hydride H$_{3}$S originates from the saddle loop. On this basis, we successfully extract the minimal model that explains how the DOS peak is formed. Our theory characterizes a large class of DOS peaks sometimes found in the three dimensional electronic structures, building a basis for profound understanding of their origins.


title: Measurement of the electron reconstruction efficiency at LHCb

abstract: The single electron track-reconstruction efficiency is calibrated using a sample corresponding to $1.3~\mathrm{fb}^{-1}$of $pp$ collision data recorded with the LHCb detector in 2017. This measurement exploits $B^+\to J/\psi (e^+e^-)K^+$ decays, where one of the electrons is fully reconstructed and paired with the kaon, while the other electron is reconstructed using only the information of the vertex detector. Despite this partial reconstruction, kinematic and geometric constraints allow the $B^+$-meson mass to be reconstructed and the signal to be well separated from backgrounds. This in turn allows the electron reconstruction efficiency to be measured by matching the partial track segment found in the vertex detector to tracks found by LHCb's regular reconstruction algorithms. The agreement between data and simulation is evaluated, and corrections are derived for simulated electrons in bins of kinematics. These correction factors allow LHCb to measure branching fractions involving single electrons with a systematic uncertainty below $1\%$.


title: LOFAR early-time search for coherent radio emission from GRB 180706A

abstract: The nature of the central engines of gamma-ray bursts (GRBs) and the composition of their relativistic jets are still under debate. If the jets are Poynting flux dominated rather than baryon dominated, a coherent radio flare from magnetic re-connection events might be expected with the prompt gamma-ray emission. There are two competing models for the central engines of GRBs; a black hole or a newly formed milli-second magnetar. If the central engine is a magnetar it is predicted to produce coherent radio emission as persistent or flaring activity. In this paper, we present the deepest limits to date for this emission following LOFAR rapid response observations of GRB 180706A. No emission is detected to a 3$\sigma$ limit of 1.7 mJy beam$^{-1}$ at 144 MHz in a two-hour LOFAR observation starting 4.5 minutes after the gamma-ray trigger. A forced source extraction at the position of GRB 180706A provides a marginally positive (1 sigma) peak flux density of $1.1 \pm 0.9$ mJy. The data were time-sliced into different sets of snapshot durations to search for FRB like emission. No short duration emission was detected at the location of the GRB. We compare these results to theoretical models and discuss the implications of a non-detection.


title: Discriminative and Robust Online Learning for Siamese Visual Tracking

abstract: The problem of visual object tracking has traditionally been handled by variant tracking paradigms, either learning a model of the object's appearance exclusively online or matching the object with the target in an offline-trained embedding space. Despite the recent success, each method agonizes over its intrinsic constraint. The online-only approaches suffer from a lack of generalization of the model they learn thus are inferior in target regression, while the offline-only approaches (e.g., convolutional siamese trackers) lack the target-specific context information thus are not discriminative enough to handle distractors, and robust enough to deformation. Therefore, we propose an online module with an attention mechanism for offline siamese networks to extract target-specific features under L2 error. We further propose a filter update strategy adaptive to treacherous background noises for discriminative learning, and a template update strategy to handle large target deformations for robust learning. Effectiveness can be validated in the consistent improvement over three siamese baselines: SiamFC, SiamRPN++, and SiamMask. Beyond that, our model based on SiamRPN++ obtains the best results over six popular tracking benchmarks and can operate beyond real-time.


title: Real-time stock analysis for blending recipes in industrial plants

abstract: Many companies use Excel spreadsheets to keep stock records and to calculate process-specific data. These spreadsheets are often hard to understand and track. And if the user does not protect them, there is a risk that the user randomly changes or erase formulas. The paper focuses on the stocks of products used in a blending process with a known recipe. Developing an application that can bring this data in a centralized form and that can assist the operator in decide is a necessity. When a programmer implements an application that uses data from plants he needs to consider one fundamental aspect as reading real-time data from the process. The real-time stock analysis application takes into account all the above elements. The application is easy to use by an operator in the command room of installation because of the planning algorithms integrated into it. The algorithms proposed and implemented in this paper have well-defined goals: identifying the ingredients needed to achieve the blending process for required quantities, determine the quantities of the finished product that can be made with the existing ingredients and determine the optimum quantities of the finished product. The application implemented in C# intensively uses these algorithms and gives the user the ability to build the result step by step.


title: Full Convergence of the Iterative Bayesian Update and Applications to Mechanisms for Privacy Protection

abstract: The iterative Bayesian update (IBU) and the matrix inversion (INV) are the main methods to retrieve the original distribution from noisy data resulting from the application of privacy protection mechanisms. We show that the theoretical foundations of the IBU established in the literature are flawed, as they rely on an assumption that in general is not satisfied in typical real datasets. We then fix the theory of the IBU, by providing a general convergence result for the underlying Expectation-Maximization method. Our framework does not rely on the above assumption, and also covers a more general local privacy model. Finally we evaluate the precision of the IBU on data sanitized with the Geometric, $k$-RR, and RAPPOR mechanisms, and we show that it outperforms INV in the first case, while it is comparable to INV in the other two cases.


title: Probing High-Energy Hadronic Interactions with Extensive Air Showers

abstract: In this article, I will review how observations performed with extensive air showers are being used to probe hadronic interactions at high energy. I will briefly overview the new studies exploring the connection between the dynamics of air showers and multiparticle production, and how this knowledge can be translated into constraints on high energy hadronic models. I will also overview direct measurements, complementary to, and beyond the reach of, accelerator experiments.


title: Regression Under Human Assistance

abstract: Decisions are increasingly taken by both humans and machine learning models. However, machine learning models are currently trained for full automation -- they are not aware that some of the decisions may still be taken by humans. In this paper, we take a first step towards the development of machine learning models that are optimized to operate under different automation levels. More specifically, we first introduce the problem of ridge regression under human assistance and show that it is NP-hard. Then, we derive an alternative representation of the corresponding objective function as a difference of nondecreasing submodular functions. Building on this representation, we further show that the objective is nondecreasing and satisfies $\alpha$-submodularity, a recently introduced notion of approximate submodularity. These properties allow a simple and efficient greedy algorithm to enjoy approximation guarantees at solving the problem. Experiments on synthetic and real-world data from two important applications -- medical diagnosis and content moderation-demonstrate that our algorithm outsources to humans those samples in which the prediction error of the ridge regression model would have been the highest if it had to make a prediction, it outperforms several competitive baselines, and its performance is robust with respect to several design choices and hyperparameters used in the experiments.


title: Multi-Objective Multi-Agent Decision Making: A Utility-based Analysis and Survey

abstract: The majority of multi-agent system (MAS) implementations aim to optimise agents' policies with respect to a single objective, despite the fact that many real-world problem domains are inherently multi-objective in nature. Multi-objective multi-agent systems (MOMAS) explicitly consider the possible trade-offs between conflicting objective functions. We argue that, in MOMAS, such compromises should be analysed on the basis of the utility that these compromises have for the users of a system. As is standard in multi-objective optimisation, we model the user utility using utility functions that map value or return vectors to scalar values. This approach naturally leads to two different optimisation criteria: expected scalarised returns (ESR) and scalarised expected returns (SER). We develop a new taxonomy which classifies multi-objective multi-agent decision making settings, on the basis of the reward structures, and which and how utility functions are applied. This allows us to offer a structured view of the field, to clearly delineate the current state-of-the-art in multi-objective multi-agent decision making approaches and to identify promising directions for future research. Starting from the execution phase, in which the selected policies are applied and the utility for the users is attained, we analyse which solution concepts apply to the different settings in our taxonomy. Furthermore, we define and discuss these solution concepts under both ESR and SER optimisation criteria. We conclude with a summary of our main findings and a discussion of many promising future research directions in multi-objective multi-agent systems.


title: User Evaluation of a Multi-dimensional Statistical Dialogue System

abstract: We present the first complete spoken dialogue system driven by a multi-dimensional statistical dialogue manager. This framework has been shown to substantially reduce data needs by leveraging domain-independent dimensions, such as social obligations or feedback, which (as we show) can be transferred between domains. In this paper, we conduct a user study and show that the performance of a multi-dimensional system, which can be adapted from a source domain, is equivalent to that of a one-dimensional baseline, which can only be trained from scratch.


title: Robust Barrier Functions for a Fully Autonomous, Remotely Accessible Swarm-Robotics Testbed

abstract: The Robotarium, a remotely accessible swarm-robotics testbed, has provided free, open access to robotics and controls research for hundreds of users in thousands of experiments. This high level of usage requires autonomy in the system, which mainly corresponds to constraint satisfaction in the context of users' submissions. In other words, in case that the users' inputs to the robots may lead to collisions, these inputs must be altered to avoid these collisions automatically. However, these alterations must be minimal so as to preserve the users' objective in the experiment. Toward this end, the system has utilized barrier functions, which admit a minimally invasive controller-synthesis procedure. However, barrier functions are yet to be robustified with respect to unmodeled disturbances (e.g., wheel slip or packet loss) in a manner conducive to real-time synthesis. As such, this paper formulates robust barrier functions for a general class of disturbed control-affine systems that, in turn, is key for the Robotarium to operate fully autonomously (i.e., without human supervision). Experimental results showcase the effectiveness of this robust formulation in a long-term experiment in the Robotarium.


title: Explicit Facial Expression Transfer via Fine-Grained Representations

abstract: Facial expression transfer between two unpaired images is a challenging problem, as fine-grained expression is typically tangled with other facial attributes. Most existing methods treat expression transfer as an application of expression manipulation, and use predicted global expression, landmarks or action units (AUs) as a guidance. However, the prediction may be inaccurate, which limits the performance of transferring fine-grained expression. Instead of using an intermediate estimated guidance, we propose to explicitly transfer facial expression by directly mapping two unpaired input images to two synthesized images with swapped expressions. Specifically, considering AUs semantically describe fine-grained expression details, we propose a novel multi-class adversarial training method to disentangle input images into two types of fine-grained representations: AU-related feature and AU-free feature. Then, we can synthesize new images with preserved identities and swapped expressions by combining AU-free features with swapped AU-related features. Moreover, to obtain reliable expression transfer results of the unpaired input, we introduce a swap consistency loss to make the synthesized images and self-reconstructed images indistinguishable. Extensive experiments show that our approach outperforms the state-of-the-art expression manipulation methods for transferring fine-grained expressions while preserving other attributes including identity and pose.


title: Limit theorems for Bajraktarevi\'c and Cauchy quotient means of independent identically distributed random variables

abstract: We derive strong laws of large numbers and central limit theorems for Bajraktarevi\'c, Gini and exponential- (also called Beta-type) and logarithmic Cauchy quotient means of independent identically distributed (i.i.d.) random variables. The exponential- and logarithmic Cauchy quotient means of a sequence of i.i.d. random variables behave asymptotically normal with the usual square root scaling just like the geometric means of the given random variables. Somewhat surprisingly, the multiplicative Cauchy quotient means of i.i.d. random variables behave asymptotically in a rather different way: in order to get a non-trivial normal limit distribution a time dependent centering is needed.


title: Activity time series of old stars from late F to early K. II. Radial velocity jitter and exoplanet detectability

abstract: The effect of stellar activity on RV appears to be a limiting factor in detecting Earth-mass planets in the habitable zone of a star similar to the Sun in spectral type and activity level. It is crucial to estimate if this conclusion remain true for other stars with current correction methods. We built realistic time series in RV and chromospheric emission for old main-sequence F6-K4 stars. The stellar parameters are spectral type, activity level, rotation period, cycle period and amplitude, latitude coverage, and spot constrast, which we chose to be in ranges that are compatible with our current knowledge. This very large set of synthetic time series allowed us to study the effect of the parameters on the RV jitter and how the different contributions to the RV are affected in this first analysis of the data set. The RV jitter was used to provide a first-order detection limit for each time series and different temporal samplings. We find that the coverage in latitude of the activity pattern and the cycle amplitudes have a strong effect on the RV jitter, as has stellar inclination. RV jitter trends with B-V and LogR'HK are similar to observations, but activity cannot be responsible for RV jitter larger than 2-3 m/s for very quiet stars: this observed jitter is therefore likely to be due to other causes. We show that based on the RV jitter that is associated with each time series and using a simple criterion, a planet with one Earth mass and a period of one to two years probably cannot be detected with current analysis techniques, except for the lower mass stars in our sample, but many observations would be required. The effect of inclination is critical. The results are very important in the context of future RV follow-ups of transit detections of such planets. A significant improvement of analysis techniques and/or observing strategies must be made to reach such low detection limits.


title: Scalings Pertaining to Current Sheet Disruption Mediated by the Plasmoid Instability

abstract: Analytic scaling relations are derived for a phenomenological model of the plasmoid instability in an evolving current sheet, including the effects of reconnection outflow. Two scenarios are considered, where the plasmoid instability can be triggered either by an injected initial perturbation or by the natural noise of the system (here referred to as the system noise). The two scenarios lead to different scaling relations because the initial noise decays when the linear growth of the plasmoid instability is not sufficiently fast to overcome the advection loss caused by the reconnection outflow, whereas the system noise represents the lowest level of fluctuations in the system. The leading order approximation for the current sheet width at disruption takes the form of a power-law multiplied by a logarithmic factor, and from that, the scaling relations for the wavenumber and the linear growth rate of the dominant mode are obtained. When the effects of the outflow are neglected, the scaling relations agree, up to the leading order approximation, with previously derived scaling relations based on a principle of least time. The analytic scaling relations are validated with numerical solutions of the model.


title: Automated Polysomnography Analysis for Detection of Non-Apneic and Non-Hypopneic Arousals using Feature Engineering and a Bidirectional LSTM Network

abstract: Objective: The aim of this study is to develop an automated classification algorithm for polysomnography (PSG) recordings to detect non-apneic and non-hypopneic arousals. Our particular focus is on detecting the respiratory effort-related arousals (RERAs) which are very subtle respiratory events that do not meet the criteria for apnea or hypopnea, and are more challenging to detect. Methods: The proposed algorithm is based on a bidirectional long short-term memory (BiLSTM) classifier and 465 multi-domain features, extracted from multimodal clinical time series. The features consist of a set of physiology-inspired features (n = 75), obtained by multiple steps of feature selection and expert analysis, and a set of physiology-agnostic features (n = 390), derived from scattering transform. Results: The proposed algorithm is validated on the 2018 PhysioNet challenge dataset. The overall performance in terms of the area under the precision-recall curve (AUPRC) is 0.50 on the hidden test dataset. This result is tied for the second-best score during the follow-up and official phases of the 2018 PhysioNet challenge. Conclusions: The results demonstrate that it is possible to automatically detect subtle non-apneic/non-hypopneic arousal events from PSG recordings. Significance: Automatic detection of subtle respiratory events such as RERAs together with other non-apneic/non-hypopneic arousals will allow detailed annotations of large PSG databases. This contributes to a better retrospective analysis of sleep data, which may also improve the quality of treatment.


title: Portfolio optimisation under rough Heston models

abstract: This thesis investigates Merton's portfolio problem under two different rough Heston models, which have a non-Markovian structure. The motivation behind this choice of problem is due to the recent discovery and success of rough volatility processes. The optimisation problem is solved from two different approaches: firstly by considering an auxiliary random process, which solves the optimisation problem with the martingale optimality principle, and secondly, by a finite dimensional approximation of the volatility process which casts the problem into its classical stochastic control framework. In addition, we show how classical results from Merton's portfolio optimisation problem can be used to help motivate the construction of the solution in both cases. The optimal strategy under both approaches is then derived in a semi-closed form, and comparisons between the results made. The approaches discussed in this thesis, combined with the historical works on the distortion transformation, provide a strong foundation to build models capable of handling increasing complexity demanded by the ever growing financial market.


title: $I = 3$ three-pion scattering amplitude from lattice QCD

abstract: We analyze the spectrum of two- and three-pion states of maximal isospin obtained recently for isosymmetric QCD with pion mass $M\approx 200\;$MeV in Ref. [1]. Using the relativistic three-particle quantization condition, we find $\sim 2 \sigma$ evidence for a nonzero value for the contact part of the three-$\pi^+$ ($I=3$) scattering amplitude. We also compare our results to leading-order chiral perturbation theory. We find good agreement at threshold, and some tension in the energy dependent part of the three-$\pi^+$scattering amplitude. We also find that the two-$\pi^+$ ($I=2$) spectrum is fit well by an $s$-wave phase shift that incorporates the expected Adler zero.


title: Sharp asymptotics of the first eigenvalue on some degenerating surfaces

abstract: We study sharp asymptotics of the first eigenvalue on Riemannian surfaces obtained from a fixed Riemannian surface by attaching a collapsing flat handle or cross cap to it. Through a careful choice of parameters this construction can be used to strictly increase the first eigenvalue normalized by area if the initial surface has some symmetries. If these symmetries are not present we show that the first eigenvalue normalized by area strictly decreases for the same range of parameters. These results are the main motivation for the construction in \cite{MS3}, where we show a monotonicity result for the normalized first eigenvalue without any symmetry assumptions.


title: Testing multi-photon interference on a silicon chip

abstract: Multi-photon interference in large multi-port interferometers is key to linear optical quantum computing and in particular to boson sampling. Silicon photonics enables complex interferometric circuits with many components in a small footprint and has the potential to extend these experiments to larger numbers of interfering modes. However, loss has generally limited the implementation of multi-photon experiments in this platform. Here, we make use of high-efficiency grating couplers to combine bright and pure quantum light sources based on ppKTP waveguides with silicon circuits. We interfere up to 5 photons in up to 15 modes, verifying genuine multi-photon interference by comparing the results against various models including partial distinguishability between photons.


title: SystemDS: A Declarative Machine Learning System for the End-to-End Data Science Lifecycle

abstract: Machine learning (ML) applications become increasingly common in many domains. ML systems to execute these workloads include numerical computing frameworks and libraries, ML algorithm libraries, and specialized systems for deep neural networks and distributed ML. These systems focus primarily on efficient model training and scoring. However, the data science process is exploratory, and deals with underspecified objectives and a wide variety of heterogeneous data sources. Therefore, additional tools are employed for data engineering and debugging, which requires boundary crossing, unnecessary manual effort, and lacks optimization across the lifecycle. In this paper, we introduce SystemDS, an open source ML system for the end-to-end data science lifecycle from data integration, cleaning, and preparation, over local, distributed, and federated ML model training, to debugging and serving. To this end, we aim to provide a stack of declarative language abstractions for the different lifecycle tasks, and users with different expertise. We describe the overall system architecture, explain major design decisions (motivated by lessons learned from Apache SystemML), and discuss key features and research directions. Finally, we provide preliminary results that show the potential of end-to-end lifecycle optimization.


title: Parallel Computation of Graph Embeddings

abstract: Graph embedding aims at learning a vector-based representation of vertices that incorporates the structure of the graph. This representation then enables inference of graph properties. Existing graph embedding techniques, however, do not scale well to large graphs. We therefore propose a framework for parallel computation of a graph embedding using a cluster of compute nodes with resource constraints. We show how to distribute any existing embedding technique by first splitting a graph for any given set of constrained compute nodes and then reconciling the embedding spaces derived for these subgraphs. We also propose a new way to evaluate the quality of graph embeddings that is independent of a specific inference task. Based thereon, we give a formal bound on the difference between the embeddings derived by centralised and parallel computation. Experimental results illustrate that our approach for parallel computation scales well, while largely maintaining the embedding quality.


title: Gravitational wave and collider signals in complex two-Higgs doublet model with dynamical CP-violation at finite temperature

abstract: Extra CP-violating source for electroweak baryogenesis can dynamically appear at finite temperature in the complex two-Higgs doublet model, which might help to alleviate the strong constraints from the electric dipole moment experiments. In this scenario, we study the detailed phase transition dynamics and the corresponding gravitational wave signals in synergy with the collider signals at future lepton colliders. For some parameter spaces, various phase transition patterns can occur, such as the multi-step phase transition and supercooling. Gravitational waves complementary to collider signals can help to pin down the underlying phase transition dynamics or different phase transition patterns.


title: On the Gamma convergence of functionals defined over pairs of measures and energy-measures

abstract: A novel general framework for the study of $\Gamma$-convergence of functionals defined over pairs of measures and energy-measures is introduced. This theory allows us to identify the $\Gamma$-limit of these kind of functionals by knowing the $\Gamma$-limit of the underlining energies. In particular, the interaction between the functionals and the underlining energies results, in the case these latter converge to a non continuous energy, in an additional effect in the relaxation process. This study was motivated by a question in the context of epitaxial growth evolution with adatoms. Interesting cases of application of the general theory are also presented.


title: Multipartite entanglement structure in the Eigenstate Thermalization Hypothesis

abstract: We study the quantum Fisher information (QFI) and, thus, the multipartite entanglement structure of thermal pure states in the context of the Eigenstate Thermalization Hypothesis (ETH). In both the canonical ensemble and the ETH, the quantum Fisher information may be explicitly calculated from the response functions. In the case of ETH, we find that the expression of the QFI bounds the corresponding canonical expression from above. This implies that although average values and fluctuations of local observables are indistinguishable from their canonical counterpart, the entanglement structure of the state is starkly different; with the difference amplified, e.g., in the proximity of a thermal phase transition. We also provide a state-of-the-art numerical example of a situation where the quantum Fisher information in a quantum many-body system is extensive while the corresponding quantity in the canonical ensemble vanishes. Our findings have direct relevance for the entanglement structure in the asymptotic states of quenched many-body dynamics.


title: Stationary states of weak coupling limit type Markov generators and quantum transport models

abstract: We prove that every stationary state in the annihilator of all Kraus operators of a weak coupling limit type Markov generator consists of two pieces, one of them supported on the interaction-free subspace and the second one on its orthogonal complement. In particular, we apply the previous result to describe in detail the structure of a slightly modified quantum transport model due to Arefeva, Kozyrev and Volovich (modified AKV's model) studied first in Ref.\cite{gggq}, in terms of generalized annihilation and creation operators.


title: DRLViz: Understanding Decisions and Memory in Deep Reinforcement Learning

abstract: We present DRLViz, a visual analytics interface to interpret the internal memory of an agent (e.g. a robot) trained using deep reinforcement learning. This memory is composed of large temporal vectors updated when the agent moves in an environment and is not trivial to understand due to the number of dimensions, dependencies to past vectors, spatial/temporal correlations, and co-correlation between dimensions. It is often referred to as a black box as only inputs (images) and outputs (actions) are intelligible for humans. Using DRLViz, experts are assisted to interpret decisions using memory reduction interactions, and to investigate the role of parts of the memory when errors have been made (e.g. wrong direction). We report on DRLViz applied in the context of video games simulators (ViZDoom) for a navigation scenario with item gathering tasks. We also report on experts evaluation using DRLViz, and applicability of DRLViz to other scenarios and navigation problems beyond simulation games, as well as its contribution to black box models interpretability and explainability in the field of visual analytics.


title: Common and not so common high-energy theory methods for condensed matter physics

abstract: This chapter is a collection of techniques, warnings, facts and ideas that are sometimes regarded as theoretical curiosities in high-energy physics but have important consequences in condensed matter physics. In particular, we describe theories that have the property of having finite but undetermined radiative corrections that also happen to describe topological semi-metallic phases in condensed matter. In the process, we describe typical methods in high-energy physics that illustrate the working principles to describe a given phase of matter and its response to external fields.


title: Dissecting the hot bubbles in LMC-N57 with XMM-Newton

abstract: We present a study of the diffuse X-ray emission from the star forming region LMC-N 57 in the Large Magellanic Cloud (LMC). We use archival XMM-Newton observations to unveil in detail the distribution of hot bubbles in this complex. X-ray emission is detected from the central superbubble (SB) DEM L 229, the supernova remnant (SNR) 0532$-$675 and the Wolf-Rayet (WR) bubble DEM L 231 around the WR star Br 48. Comparison with infrared images unveils the powerful effect of massive stars in destroying their nurseries. The distribution of the hot gas in the SNR and the SB display their maxima in regions in contact with the filamentary cold material detected by IR images. Our observations do not reveal extended X-ray emission filling DEM L 231, although several point-like sources are detected in the field of view of this WR nebula. The X-ray properties of Br 48 are consistent with a binary WN4$+$O as proposed by other authors. We modeled the X-ray emission from the SB and found that its X-ray emission can be simply explained by pressure-driven wind model, that is, there is no need to invoke the presence of a SN explosion as previously suggested. The pressure calculations of the hot gas confirms that the dynamical evolution of the SB DEM L 229 is dominated by the stellar winds from the star cluster LH 76.


title: Scattering diagrams, stability conditions, and coherent sheaves on $\mathbb{P}^2$

abstract: We show that a purely algebraic structure, a two-dimensional scattering diagram, describes a large part of the wall-crossing behavior of moduli spaces of Bridgeland semistable objects in the derived category of coherent sheaves on $\mathbb{P}^2$. This gives a new algorithm computing the Hodge numbers of the intersection cohomology of the classical moduli spaces of Gieseker semistable sheaves on $\mathbb{P}^2$, or equivalently the refined Donaldson-Thomas invariants for compactly supported sheaves on local $\mathbb{P}^2$. As applications, we prove that the intersection cohomology of moduli spaces of Gieseker semistable sheaves on $\mathbb{P}^2$ is Hodge-Tate, and we give the first non-trivial numerical checks of the general $\chi$-independence conjecture for refined Donaldson-Thomas invariants of one-dimensional sheaves on local $\mathbb{P}^2$.


title: A Proposed Framework for Interactive Virtual Reality In Situ Visualization of Parallel Numerical Simulations

abstract: As computer simulations progress to increasingly complex, non-linear, and three-dimensional systems and phenomena, intuitive and immediate visualization of their results is becoming crucial. While Virtual Reality (VR) and Natural User Interfaces (NUIs) have been shown to improve understanding of complex 3D data, their application to live in situ visualization and computational steering is hampered by performance requirements. Here, we present the design of a software framework for interactive VR in situ visualization of parallel numerical simulations, as well as a working prototype implementation. Our design is targeted towards meeting the performance requirements for VR, and our work is packaged in a framework that allows for easy instrumentation of simulations. Our preliminary results inform about the technical feasibility of the architecture, as well as the challenges that remain.


title: Local-to-global frames and applications to dynamical sampling problem

abstract: In this paper we consider systems of vectors in a Hilbert space $\mathcal{H}$ of the form $\{g_{jk}: j \in J, \, k\in K\}\subset \mathcal{H}$ where $J$ and $K$ are countable sets of indices. We find conditions under which the local reconstruction properties of such a system extend to global stable recovery properties on the whole space. As a particular case, we obtain new local-to-global results for systems of type $\{A^ng\}_{g\in\mathcal{G},0\leq n\leq L }$ arising in the dynamical sampling problem.