title: Structure-Feature based Graph Self-adaptive Pooling

abstract: Various methods to deal with graph data have been proposed in recent years. However, most of these methods focus on graph feature aggregation rather than graph pooling. Besides, the existing top-k selection graph pooling methods have a few problems. First, to construct the pooled graph topology, current top-k selection methods evaluate the importance of the node from a single perspective only, which is simplistic and unobjective. Second, the feature information of unselected nodes is directly lost during the pooling process, which inevitably leads to a massive loss of graph feature information. To solve these problems mentioned above, we propose a novel graph self-adaptive pooling method with the following objectives: (1) to construct a reasonable pooled graph topology, structure and feature information of the graph are considered simultaneously, which provide additional veracity and objectivity in node selection; and (2) to make the pooled nodes contain sufficiently effective graph information, node feature information is aggregated before discarding the unimportant nodes; thus, the selected nodes contain information from neighbor nodes, which can enhance the use of features of the unselected nodes. Experimental results on four different datasets demonstrate that our method is effective in graph classification and outperforms state-of-the-art graph pooling methods.


title: Using Sampled Network Data With The Autologistic Actor Attribute Model

abstract: Social science research increasingly benefits from statistical methods for understanding the structured nature of social life, including for social network data. However, the application of statistical network models within large-scale community research is hindered by too little understanding of the validity of their inferences under realistic data collection conditions, including sampled or missing network data. The autologistic actor attribute model (ALAAM) is a statistical model based on the well-established exponential random graph model (ERGM) for social networks. ALAAMs can be regarded as a social influence model, predicting an individual-level outcome based on the actor's network ties, concurrent outcomes of his/her network partners, and attributes of the actor and his/her network partners. In particular, an ALAAM can be used to measure contagion effects, that is, the propensity of two actors connected by a social network tie to both have the same value of an attribute. We investigate the effect of using simple random samples and snowball samples of network data on ALAAM parameter inference, and find that parameter inference can still work well even with a nontrivial fraction of missing nodes. However it is safer to take a snowball sample of the network and estimate conditional on the snowball sampling structure.


title: A Kernel of Truth: Determining Rumor Veracity on Twitter by Diffusion Pattern Alone

abstract: Recent work in the domain of misinformation detection has leveraged rich signals in the text and user identities associated with content on social media. But text can be strategically manipulated and accounts reopened under different aliases, suggesting that these approaches are inherently brittle. In this work, we investigate an alternative modality that is naturally robust: the pattern in which information propagates. Can the veracity of an unverified rumor spreading online be discerned solely on the basis of its pattern of diffusion through the social network? Using graph kernels to extract complex topological information from Twitter cascade structures, we train accurate predictive models that are blind to language, user identities, and time, demonstrating for the first time that such "sanitized" diffusion patterns are highly informative of veracity. Our results indicate that, with proper aggregation, the collective sharing pattern of the crowd may reveal powerful signals of rumor truth or falsehood, even in the early stages of propagation.


title: Atomic-scale insights into semiconductor heterostructures: from experimental three-dimensional analysis of the interface to a generalized theory of interface roughness scattering

abstract: We develop a generalized theory for the scattering process produced by interface roughness on charge carriers and which is suitable for any semiconductor heterostructure. By exploiting our experimental insights into the three-dimensional atomic landscape obtained on Ge/GeSi heterointerfaces obtained by atom probe tomography, we have been able to define the full set of interface parameters relevant to the scattering potential, including both the in-plane and axial correlation inside real diffuse interfaces. Our experimental findings indicate a partial coherence of the interface roughness along the growth direction within the interfaces. We show that it is necessary to include this feature, previously neglected by theoretical models, when heterointerfaces characterized by finite interface widths are taken into consideration.


title: The exponentially weighted average forecaster in geodesic spaces of non-positive curvature

abstract: This paper addresses the problem of prediction with expert advice for outcomes in a geodesic space with non-positive curvature in the sense of Alexandrov. Via geometric considerations, and in particular the notion of barycenters, we extend to this setting the definition and analysis of the classical exponentially weighted average forecaster. We also adapt the principle of online to batch conversion to this setting. We shortly discuss the application of these results in the context of aggregation and for the problem of barycenter estimation.

categories: math.ST cs.GT cs.LG math.MG stat.ML stat.TH


title: The topological dimension of radial Julia sets

abstract: We prove that the meandering set for $f_a(z)=e^z+a$ is homeomorphic to the space of irrational numbers whenever $a$ belongs to the Fatou set of $f_a$. This extends recent results by Vasiliki Evdoridou and Lasse Rempe. It implies that the radial Julia set of $f_a$ has topological dimension zero for all attracting and parabolic parameters, including all $a\in (-\infty,-1]$. Similar results are obtained for Fatou's function $f(z)=z+1+e^{-z}$.


title: Measuring relative opinion from location-based social media: A case study of the 2016 U.S. presidential election

abstract: Social media has become an emerging alternative to opinion polls for public opinion collection, while it is still posing many challenges as a passive data source, such as structurelessness, quantifiability, and representativeness. Social media data with geotags provide new opportunities to unveil the geographic locations of users expressing their opinions. This paper aims to answer two questions: 1) whether quantifiable measurement of public opinion can be obtained from social media and 2) whether it can produce better or complementary measures compared to opinion polls. This research proposes a novel approach to measure the relative opinion of Twitter users towards public issues in order to accommodate more complex opinion structures and take advantage of the geography pertaining to the public issues. To ensure that this new measure is technically feasible, a modeling framework is developed including building a training dataset by adopting a state-of-the-art approach and devising a new deep learning method called Opinion-Oriented Word Embedding. With a case study of the tweets selected for the 2016 U.S. presidential election, we demonstrate the predictive superiority of our relative opinion approach and we show how it can aid visual analytics and support opinion predictions. Although the relative opinion measure is proved to be more robust compared to polling, our study also suggests that the former can advantageously complement the later in opinion prediction.


title: Microwave electrometry via electromagnetically induced absorption in cold Rydberg atoms

abstract: The atom-based traceable standard for microwave electrometry shows promising advantages by enabling stable and uniform measurement. Here we theoretically propose and then experimentally realize an alternative direct International System of Units (SI)-traceable and self-calibrated method for measuring a microwave electric field strength based on electromagnetically induced absorption (EIA) in cold Rydberg atoms. Comparing with the method of electromagnetically induced transparency, we show that the equivalence relation between microwave Rabi frequency and Autler-Townes splitting is more valid and is even more robust against the experimental parameters in the EIA's linear region. Furthermore, a narrower linewidth of cold Rydberg EIA enables us to realize a direct SI-traceable microwave-electric-field measurement as small as $\sim$100 $\mu\mathrm{\!V} \mathrm{cm}^{\!-\!1}$.

categories: quant-ph cond-mat.quant-gas physics.atom-ph physics.optics


title: Slippery flows of a Carbopol gel in a micro-channel

abstract: The ability to predict and/or control wall slip is a fundamental problem in the hydrodynamics of yield stress fluids, which is poorly understood to date and has important applications in bio- and micro-fluidic systems. Systematic measurements of steady flows of a simple yield stress fluid (Carbopol Ultrez 10) in a plane acrylic micro-channel are used to establish rigorous scaling laws for the wall velocity gradient and the slip velocity. By means of epi-fluorescent microscopy combined with a custom developed Digital Particle Image Velocimetry (DPIV) technique, times series of velocity fields were measured within a wide range of flow rates and three distinct flow regimes were identified: full plug, partial plug and fully yielded. Corresponding to each flow regime, wall velocity gradients and slip velocities were obtained by extrapolating the velocity profiles using a smoothing spline function. By combining the flow field measurements with the macro-rheological measurements, scaling laws for the wall velocity gradient and the slip velocity with the wall shear stress were identified, and compared with results from the literature. Detailed microscopic measurements of the velocity field enabled an assessment of the effectiveness of a chemical treatment of the channel walls meant to suppress wall slip proposed by Metivier and coworkers in Ref. [1].


title: Distributed Degenerate Band Edge Oscillator

abstract: We propose a new class of oscillators by engineering the dispersion of two-coupled periodic waveguides to exhibit a degenerate band edge (DBE). The DBE is an exceptional point of degeneracy (EPD) of order four, i.e., representing the coalescence of four eigenmodes of a waveguide system without loss and gain. We present a distributed DBE oscillator realized in periodic coupled transmission lines with a unique mode selection scheme that leads to a stable single-frequency oscillation, even in the presence of load variation. The DBE oscillator potentially leads to a boost of the efficiency and performance of RF sources, thanks to the unique features associated to the EPD concept. This class of oscillators is promising for improving discrete-distributed coherent sources and can be extended to radiating structures to achieve a new class of active integrated antenna arrays.


title: Enhanced carrier mobility in anisotropic 2D tetrahex-carbon through strain engineering

abstract: A recently predicted two dimensional (2D) carbon allotrope, tetrahex-carbon consisting of tetragonal and hexagonal rings, draws research interests due to its unique mechanical and electronic properties. Tetrahex-C shows ultrahigh strength, negative Poisson ratio, a direct band gap and high carrier mobility. In this work, we employ first-principles density-functional theory calculations to explore the directional dependence of electronic properties such as carrier effective mass and mobility in tetrahex-C. Tetrahex-C demonstrates strong anisotropicity in effective mass of charge carrier and therefore its mobility (electric conductance) exhibits a strong orientation preference. More interesting, we find that such unique anisotropic carrier effective mass and mobility can be controlled by simple uniaxial strain. The orientation dependence of effective mass can be dramatically rotated by 90 degrees through applying uniaxial tensile strain beyond ~ 7% (11%) in the armchair direction for the hole (electron). As a result, the intrinsic carrier mobility in tetrahex-C is significantly enhanced. The results are useful for potential electronic and mechanical applications in tetrahex-C.


title: Isometric study of Wasserstein spaces --- the real line

abstract: Recently Kloeckner described the structure of the isometry group of the quadratic Wasserstein space $\mathcal{W}_2\left(\mathbb{R}^n\right)$. It turned out that the case of the real line is exceptional in the sense that there exists an exotic isometry flow. Following this line of investigation, we compute $\mathrm{Isom}\left(\mathcal{W}_p(\mathbb{R})\right)$, the isometry group of the Wasserstein space $\mathcal{W}_p(\mathbb{R})$ for all $p \in [1, \infty)\setminus\{2\}$. We show that $\mathcal{W}_2(\mathbb{R})$ is also exceptional regarding the parameter $p$: $\mathcal{W}_p(\mathbb{R})$ is isometrically rigid if and only if $p\neq 2$. Regarding the underlying space, we prove that the exceptionality of $p=2$ disappears if we replace $\mathbb{R}$ by the compact interval $[0,1]$. Surprisingly, in that case, $\mathcal{W}_p\left([0,1]\right)$ is isometrically rigid if and only if $p\neq1$. Moreover, $\mathcal{W}_1\left([0,1]\right)$ admits isometries that split mass, and $\mathrm{Isom}\left(\mathcal{W}_1\left([0,1]\right)\right)$ cannot be embedded into $\mathrm{Isom}\left(\mathcal{W}_1(\mathbb{R})\right).$


title: Optimized spiking neurons classify images with high accuracy through temporal coding with two spikes

abstract: Spike-based neuromorphic hardware promises to reduce the energy consumption of image classification and other deep learning applications, particularly on mobile phones or other edge devices. However, direct training of deep spiking neural networks is difficult, and previous methods for converting trained artificial neural networks to spiking neurons were inefficient because the neurons had to emit too many spikes. We show that a substantially more efficient conversion arises when one optimizes the spiking neuron model for that purpose, so that it not only matters for information transmission how many spikes a neuron emits, but also when it emits those spikes. This advances the accuracy that can be achieved for image classification with spiking neurons, and the resulting networks need on average just two spikes per neuron for classifying an image. In addition, our new conversion method improves latency and throughput of the resulting spiking networks.


title: First search for dyons with the full MoEDAL trapping detector in 13 TeV pp collisions

abstract: The MoEDAL trapping detector, consists of approximately 800 kg of aluminium volumes. It was exposed during Run-2 of the LHC program to 6.46 fb^-1 of 13 TeV proton-proton collisions at the LHCb interaction point. Evidence for dyons (particles with electric and magnetic charge) captured in the trapping detector was sought by passing the aluminium volumes comprising the detector through a SQUID magnetometer. The presence of a trapped dyon would be signalled by a persistent current induced in the SQUID magnetometer. On the basis of a Drell-Yan production model, we exclude dyons with a magnetic charge ranging up to 5 Dirac charges, and an electric charge up to 200 times the fundamental electric charge for mass limits in the range 790 - 3130 GeV.


title: CMOS-Free Multilayer Perceptron Enabled by Four-Terminal MTJ Device

abstract: Neuromorphic computing promises revolutionary improvements over conventional systems for applications that process unstructured information. To fully realize this potential, neuromorphic systems should exploit the biomimetic behavior of emerging nanodevices. In particular, exceptional opportunities are provided by the non-volatility and analog capabilities of spintronic devices. While spintronic devices have previously been proposed that emulate neurons and synapses, complementary metal-oxide-semiconductor (CMOS) devices are required to implement multilayer spintronic perceptron crossbars. This work therefore proposes a new spintronic neuron that enables purely spintronic multilayer perceptrons, eliminating the need for CMOS circuitry and simplifying fabrication.

categories: cs.NE cond-mat.mes-hall cs.ET physics.app-ph


title: Supporting DNN Safety Analysis and Retraining through Heatmap-based Unsupervised Learning

abstract: Deep neural networks (DNNs) are increasingly important in safety-critical systems, for example in their perception layer to analyze images. Unfortunately, there is a lack of methods to ensure the functional safety of DNN-based components. We observe three major challenges with existing practices regarding DNNs in safety-critical systems: (1) scenarios that are underrepresented in the test set may lead to serious safety violation risks, but may, however, remain unnoticed; (2) characterizing such high-risk scenarios is critical for safety analysis; (3) retraining DNNs to address these risks is poorly supported when causes of violations are difficult to determine. To address these problems in the context of DNNs analyzing images, we propose HUDD, an approach that automatically supports the identification of root causes for DNN errors. HUDD identifies root causes by applying a clustering algorithm to heatmaps capturing the relevance of every DNN neuron on the DNN outcome. Also, HUDD retrains DNNs with images that are automatically selected based on their relatedness to the identified image clusters. We evaluated HUDD with DNNs from the automotive domain. HUDD was able to identify all the distinct root causes of DNN errors, thus supporting safety analysis. Also, our retraining approach has shown to be more effective at improving DNN accuracy than existing approaches.


title: Optimal Iterative Sketching with the Subsampled Randomized Hadamard Transform

abstract: Random projections or sketching are widely used in many algorithmic and learning contexts. Here we study the performance of iterative Hessian sketch for least-squares problems. By leveraging and extending recent results from random matrix theory on the limiting spectrum of matrices randomly projected with the subsampled randomized Hadamard transform, and truncated Haar matrices, we can study and compare the resulting algorithms to a level of precision that has not been possible before. Our technical contributions include a novel formula for the second moment of the inverse of projected matrices. We also find simple closed-form expressions for asymptotically optimal step-sizes and convergence rates. These show that the convergence rate for Haar and randomized Hadamard matrices are identical, and asymptotically improve upon Gaussian random projections. These techniques may be applied to other algorithms that employ randomized dimension reduction.


title: Designing GANs: A Likelihood Ratio Approach

abstract: We are interested in the design of generative networks. The training of these mathematical structures is mostly performed with the help of adversarial (min-max) optimization problems. We propose a simple methodology for constructing such problems assuring, at the same time, consistency of the corresponding solution. We give characteristic examples developed by our method, some of which can be recognized from other applications, and some are introduced here for the first time. We present a new metric, the likelihood ratio, that can be employed online to examine the convergence and stability during the training of different Generative Adversarial Networks (GANs). Finally, we compare various possibilities by applying them to well-known datasets using neural networks of different configurations and sizes.


title: To pipeline or not to pipeline, that is the question

abstract: In designing query processing primitives, a crucial design choice is the method for data transfer between two operators in a query plan. As we were considering this critical design mechanism for an in-memory database system that we are building, we quickly realized that (surprisingly) there isn't a clear definition of this concept. Papers are full or ad hoc use of terms like pipelining and blocking, but as these terms are not crisply defined, it is hard to fully understand the results attributed to these concepts. To address this limitation, we introduce a clear terminology for how to think about data transfer between operators in a query pipeline. We show that there isn't a clear definition of pipelining and blocking, and that there is a full spectrum of techniques based on a simple concept called unit-of-transfer. Next, we develop an analytical model for inter-operator communication, and highlight the key parameters that impact performance (for in-memory database settings). Armed with this model, we then apply it to the system we are designing and highlight the insights we gathered from this exercise. We find that the gap between pipelining and non-pipelining query execution, w.r.t. key factors such as performance and memory footprint is quite narrow, and thus system designers should likely rethink the notion of pipelining vs. blocking for in-memory database systems.


title: Deep Self-Supervised Representation Learning for Free-Hand Sketch

abstract: In this paper, we tackle for the first time, the problem of self-supervised representation learning for free-hand sketches. This importantly addresses a common problem faced by the sketch community -- that annotated supervisory data are difficult to obtain. This problem is very challenging in that sketches are highly abstract and subject to different drawing styles, making existing solutions tailored for photos unsuitable. Key for the success of our self-supervised learning paradigm lies with our sketch-specific designs: (i) we propose a set of pretext tasks specifically designed for sketches that mimic different drawing styles, and (ii) we further exploit the use of a textual convolution network (TCN) in a dual-branch architecture for sketch feature learning, as means to accommodate the sequential stroke nature of sketches. We demonstrate the superiority of our sketch-specific designs through two sketch-related applications (retrieval and recognition) on a million-scale sketch dataset, and show that the proposed approach outperforms the state-of-the-art unsupervised representation learning methods, and significantly narrows the performance gap between with supervised representation learning.


title: Parallel implicit-explicit general linear methods

abstract: High-order discretizations of partial differential equations (PDEs) necessitate high-order time integration schemes capable of handling both stiff and nonstiff operators in an efficient manner. Implicit-explicit (IMEX) integration based on general linear methods (GLMs) offers an attractive solution due to their high stage and method order, as well as excellent stability properties. The IMEX characteristic allows stiff terms to be treated implicitly and nonstiff terms to be efficiently integrated explicitly. This work develops two systematic approaches for the development of IMEX GLMs of arbitrary order with stages that can be solved in parallel. The first approach is based on diagonally implicit multistage integration methods (DIMSIMs) of types 3 and 4. The second is a parallel generalization of IMEX Euler and has the interesting feature that the linear stability is independent of the order of accuracy. Numerical experiments confirm the theoretical rates of convergence and reveal that the new schemes are more efficient than serial IMEX GLMs and IMEX Runge-Kutta methods.


title: Benjamini-Schramm convergence and spectrum of random hyperbolic surfaces of high genus

abstract: We study geometric and spectral properties of typical hyperbolic surfaces of high genus, excluding a set of small measure for the Weil-Petersson probability measure. We first prove Benjamini-Schramm convergence to the hyperbolic plane H as the genus g goes to infinity. An estimate for the number of eigenvalues in an interval [a,b] in terms of a, b and g is then proven using the Selberg trace formula. It implies the convergence of spectral measures to the spectral measure of H as g $\rightarrow$+$\infty$, and a uniform Weyl law as b $\rightarrow$+$\infty$. We deduce a bound on the number of small eigenvalues, and the multiplicity of any eigenvalue.


title: Boundary value problems in Euclidean space for Bosonic Laplacians

abstract: A bosonic Laplacian is a conformally invariant second order differential operator acting on smooth functions defined on domains in Euclidean space and taking values in higher order irreducible representations of the special orthogonal group. In this paper, we study boundary value problems involving bosonic Laplacians in the upper-half space and the unit ball. Poisson kernels in the upper-half space and the unit ball are constructed, which give us solutions to the Dirichlet problems with $L^p$ boundary data, $1 \leq p \leq \infty$. We also prove the uniqueness for solutions to the Dirichlet problems with continuous data for bosonic Laplacians and provide analogs of some properties of harmonic functions for null solutions of bosonic Laplacians, for instance, Cauchy's estimates, the mean-value property, Liouville's Theorem, etc.


title: Orbital period dependent spatial scatter of Galactic Black Hole Transients

abstract: There remain significant uncertainties in the origin and evolution of black holes in binary systems, in particular regarding their birth sites and the influence of natal kicks. These are long-standing issues, but their debate has been reinvigorated in the era of gravitational wave detections and the improving precision of astrometric measurements. Using recent and archival characterisation of Galactic black hole X-ray binaries (BHXBs), we report here an apparent anticorrelation between P{orb} (system orbital periods) and scatter in z (elevation above the Galactic plane). The absence of long period sources at high z is not an obvious observational bias, and two possible explanatory scenarios are qualitatively explored: (1) a disc origin for BHXBs followed by natal kicks producing the scatter in z, with only the tightest binaries preferentially surviving strong kicks; (2) a halo origin, with P{orb} shortening through dynamical interactions in globular clusters (GCs). For the latter case, we show a correspondence in z-scatter between BHXBs and the GCs with most compact core radii of <0.1pc. However, the known absence of outbursting BHXB transients within Galactic GCs remains puzzling in this case, in contrast to the multitude of known GC neutron star XRBs. These results provide an interesting observational constraint for any black hole binary evolutionary model to satisfy.


title: Scoring Graspability based on Grasp Regression for Better Grasp Prediction

abstract: Grasping objects is one of the most important abilities that a robot needs to master in order to interact with its environment. Current state-of-the-art methods rely on deep neural networks trained to jointly predict a graspability score together with a regression of an offset with respect to grasp reference parameters. However, these two predictions are performed independently, which can lead to a decrease in the actual graspability score when applying the predicted offset. Therefore, in this paper, we extend a state-of-the-art neural network with a scorer that evaluates the graspability of a given position, and introduce a novel loss function which correlates regression of grasp parameters with graspability score. We show that this novel architecture improves performance from 82.13% for a state-of-the-art grasp detection network to 85.74% on Jacquard dataset. When the learned model is transferred onto a real robot, the proposed method correlating graspability and grasp regression achieves a 92.4% rate compared to 88.1% for the baseline trained without the correlation.


title: Direct and real-time analysis in a plasma reactor using a compact FTICR/MS: degradation of acetone in nitrogen and by-products formation

abstract: Methods for reduction of Volatile Organic Compounds (VOCs) content in air depend on the 10 application considered. For low concentration and low flux, non-thermal plasma methods are often considered as efficient. However, the complex chemistry involved is still not well understood as there is a lack of datasets of byproducts formation. So as to overcome this issue, rapid analytical methods are needed. We present the coupling of a rapid chemical ionization mass spectrometer (CIMS) for the real-time analysis of the VOCs formed during a degradation experiment. The high resolution instrument used 15 allows for chemical ionization and direct quantification of non-targeted compounds. We present degradation experiments of acetone in a photo-triggered nitrogen plasma discharge. Two regimes were highlighted: efficient conversion at low concentrations (<100ppm) and moderate efficiency conversion at higher concentrations (>100ppm). Those two regimes were clearly delimited as the sum of two exponential curves occuring at respectively low and high concentrations. Many by-products were detected, in 20 particular HCN presented a significantly high yield. Nitrile compounds (acetonitrile, propionitrile,...) are formed as well. To a lower extent, ketene, acetaldehyde and formaldehyde are observed. The association of the high resolution mass spectrometer to the plasma reactor will allow further insights into the plasma chemistry and comparison to modelisation.


title: Finite-Sample Analysis of Stochastic Approximation Using Smooth Convex Envelopes

abstract: Stochastic Approximation (SA) is a popular approach for solving fixed-point equations where the information is corrupted by noise. In this paper, we consider an SA involving a contraction mapping with respect to an arbitrary norm, and show its finite-sample error bounds while using different stepsizes. The idea is to construct a smooth Lyapunov function using the generalized Moreau envelope, and show that the iterates of SA have negative drift with respect to that Lyapunov function. Our result is applicable in Reinforcement Learning (RL). In particular, we use it to establish the first-known convergence rate of the V-trace algorithm for off-policy TD-learning. Moreover, we also use it to study TD-learning in the on-policy setting, and recover the existing state-of-the-art results for $Q$-learning. Importantly, our construction results in only a logarithmic dependence of the convergence bound on the size of the state-space.


title: Absence of logarithmic divergence of the entanglement entropies at the phase transitions of a 2D classical hard rod model

abstract: Entanglement entropy is a powerful tool to detect continuous, discontinuous and even topological phase transitions in quantum as well as classical systems. In this work, von Neumann and Renyi entanglement entropies are studied numerically for classical lattice models in a square geometry. A cut is made from the center of the square to the midpoint of one of its edges, say the right edge. The entanglement entropies measure the entanglement between the left and right halves of the system. As in the strip geometry, von Neumann and Renyi entanglement entropies diverge logarithmically at the transition point while they display a jump for first-order phase transitions. The analysis is extended to a classical model of non-overlapping finite hard rods deposited on a square lattice for which Monte Carlo simulations have shown that, when the hard rods span over 7 or more lattice sites, a nematic phase appears in the phase diagram between two disordered phases. A new Corner Transfer Matrix Renormalization Group algorithm (CTMRG) is introduced to study this model. No logarithmic divergence of entanglement entropies is observed at the phase transitions in the CTMRG calculation discussed here. We therefore infer that the transitions neither can belong to the Ising universality class, as previously assumed in the literature, nor be discontinuous.


title: Torch-Struct: Deep Structured Prediction Library

abstract: The literature on structured prediction for NLP describes a rich collection of distributions and algorithms over sequences, segmentations, alignments, and trees; however, these algorithms are difficult to utilize in deep learning frameworks. We introduce Torch-Struct, a library for structured prediction designed to take advantage of and integrate with vectorized, auto-differentiation based frameworks. Torch-Struct includes a broad collection of probabilistic structures accessed through a simple and flexible distribution-based API that connects to any deep learning model. The library utilizes batched, vectorized operations and exploits auto-differentiation to produce readable, fast, and testable code. Internally, we also include a number of general-purpose optimizations to provide cross-algorithm efficiency. Experiments show significant performance gains over fast baselines and case-studies demonstrate the benefits of the library. Torch-Struct is available at https://github.com/harvardnlp/pytorch-struct.


title: Observability and unique continuation of the adjoint of a linearized compressible fluid-structure model in a 2d channel

abstract: Our objective is to study the observability and unique continuation property of the adjoint of a linearized compressible fluid structure interaction model in a 2D channel. Concerning the structure we will consider a damped Euler-Bernoulli beam located on a portion of the boundary. In the present article we establish an observability inequality for the adjoint of the linearized fluid structure interaction problem under consideration which in principle is equivalent with the null controllability of the linearized system. As a corollary of the derived observability inequality we also obtain a unique continuation property for the adjoint problem.


title: A Code for Unscented Kalman Filtering on Manifolds (UKF-M)

abstract: The present paper introduces a novel methodology for Unscented Kalman Filtering (UKF) on manifolds that extends previous work by the authors on UKF on Lie groups. Beyond filtering performance, the main interests of the approach are its versatility, as the method applies to numerous state estimation problems, and its simplicity of implementation for practitioners not being necessarily familiar with manifolds and Lie groups. We have developed the method on two independent open-source Python and Matlab frameworks we call UKF-M, for quickly implementing and testing the approach. The online repositories contain tutorials, documentation, and various relevant robotics examples that the user can readily reproduce and then adapt, for fast prototyping and benchmarking. The code is available at https://github.com/CAOR-MINES-ParisTech/ukfm.


title: Optimal l-one Rank One Matrix Decompositions

abstract: In this paper we consider the decomposition of positive semidefinite matrices as a sum of rank one matrices. We introduce and investigate the properties of various measures of optimality of such decompositions. For some classes of positive semidefinite matrices we give explicitly these optimal decompositions. These classes include diagonally dominant matrices and certain of their generalizations, $2\times 2$, and a class of $3\times 3$ matrices.


title: A Frame Decomposition of the Atmospheric Tomography Operator

abstract: We consider the problem of atmospheric tomography, as it appears for example in adaptive optics systems for extremely large telescopes. We derive a frame decomposition, i.e., a decomposition in terms of a frame, of the underlying atmospheric tomography operator, extending the singular-value-type decomposition results of Neubauer and Ramlau (2017) by allowing a mixture of both natural and laser guide stars, as well as arbitrary aperture shapes. Based on both analytical considerations as well as numerical illustrations, we provide insight into the properties of the derived frame decomposition and its building blocks.


title: FindBounce: package for multi-field bounce actions

abstract: We are launching FindBounce, a Mathematica package for the evaluation of the Euclidean bounce action that enters the decay rate of metastable states in quantum and thermal field theories. It is based on the idea of polygonal bounces, which is a semi-analytical approach to solving the bounce equation by discretizing the potential into piecewise linear segments. This allows for a fast and robust evaluation of arbitrary potentials with specified precision and any number of scalar fields. Time cost grows linearly with the number of fields and/or the number of segments. Computation with 20 fields takes $\sim 2$ seconds with $0.5\%$ accuracy of the action. The FindBounce function is simple to use with the native Mathematica look and feel, it is easy to install, and comes with detailed documentation and physical examples, such as the calculation of the nucleation temperature. We also provide timing benchmarks with comparisons to existing tools, where applicable.


title: Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS$_2$ nanopores

abstract: Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based in low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based in nanoporous materials have been suggested in recent years. One of the materials suggested is MoS$_2$, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct Force Fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were build based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through two nanopores sizes were used -- one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter, and two rigid water models from the TIP4P family -- the TIP4P/2005 and TIP4P/$\epsilon$. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only the water model can influence the outcomes, but the ion model plays a crucial role.


title: Adversarial-based neural networks for affect estimations in the wild

abstract: There is a growing interest in affective computing research nowadays given its crucial role in bridging humans with computers. This progress has been recently accelerated due to the emergence of bigger data. One recent advance in this field is the use of adversarial learning to improve model learning through augmented samples. However, the use of latent features, which is feasible through adversarial learning, is not largely explored, yet. This technique may also improve the performance of affective models, as analogously demonstrated in related fields, such as computer vision. To expand this analysis, in this work, we explore the use of latent features through our proposed adversarial-based networks for valence and arousal recognition in the wild. Specifically, our models operate by aggregating several modalities to our discriminator, which is further conditioned to the extracted latent features by the generator. Our experiments on the recently released SEWA dataset suggest the progressive improvements of our results. Finally, we show our competitive results on the Affective Behavior Analysis in-the-Wild (ABAW) challenge dataset


title: Ambient backscatters-friendly 5G networks: creating hot spots for tags and good spots for readers

abstract: In this paper, we present an ambient backscatters-friendly 5G network that creates locations with large power (hot spots) for tags and good reception locations (quiet spots or coherent spots) for readers. The massive multiple input multiple output (M-MIMO) antenna and beamforming capability of the 5G network is used as follows. In a first step, a training device (separate from tags and readers) is used to send pilots and train the 5G network to perform focusing and/or nulling onto marked locations. In a second step, tags and readers are positioned onto the marked locations. The robustness of M-MIMO beamforming to slight changes in the environment is exploited. Our initial simulation results, in a multipath propagation channel environment, show that creating a hot spot on a tag improves the tag-reader range however with a low probability of detection. Creating a hot spot on a tag and a good reception spot on a reader at the same time, improves the tag-reader range with 99% probability of detection. The study also shows that beamforming does not degrade the performance of a legacy 5G communication.


title: Diffusion bridges for stochastic Hamiltonian systems and shape evolutions

abstract: Stochastically evolving geometric systems are studied in shape analysis and computational anatomy for modelling random evolutions of human organ shapes. The notion of geodesic paths between shapes is central to shape analysis and has a natural generalisation as diffusion bridges in a stochastic setting. Simulation of such bridges is key to solve inference and registration problems in shape analysis. We demonstrate how to apply state-of-the-art diffusion bridge simulation methods to recently introduced stochastic shape deformation models thereby substantially expanding the applicability of such models. We exemplify these methods by estimating template shapes from observed shape configurations while simultaneously learning model parameters.


title: Operations on the Hochschild Bicomplex of a Diagram of Algebras

abstract: A diagram of algebras is a functor valued in a category of associative algebras. I construct an operad acting on the Hochschild bicomplex of a diagram of algebras. Using this operad, I give a direct proof that the Hochschild cohomology of a diagram of algebras is a Gerstenhaber algebra. I also show that the total complex is an $L_\infty$-algebra. The same results are true for the reduced and asimplicial subcomplexes and asimplicial cohomology. This structure governs deformations of diagrams of algebras through the Maurer-Cartan equation.


title: Hyperfine structure constants on the relativistic coupled cluster level with associated uncertainties

abstract: Accurate predictions of hyperfine structure (HFS) constants are important in many areas of chemistry and physics, from the determination of nuclear electric and magnetic moments to benchmarking of new theoretical methods. We present a detailed investigation of the performance of the relativistic coupled cluster method for calculating HFS constants withing the finite-field scheme. The two selected test systems are $^{133}$Cs and $^{137}$BaF. Special attention has been paid to construct a theoretical uncertainty estimate based on investigations on basis set, electron correlation and relativistic effects. The largest contribution to the uncertainty estimate comes from higher order correlation contributions. Our conservative uncertainty estimate for the calculated HFS constants is $\sim$ 5.5\%, while the actual deviation of our results from experimental values was $<1$\% in all cases.


title: Equal sums of two cubes of binary quadratic forms

abstract: We give a complete description of all solutions to the equation $f_1^3 + f_2^3 = f_3^3 + f_4^3$ for quadratic forms $f_j \in \mathbb C[x,y]$ and show how Ramanujan's example can be extended to three equal sums of pairs of cubes. We also give a complete census in counting the number of ways a sextic $p \in \mathbb C[x,y]$ can be written as a sum of two cubes. The extreme example is $p(x,y) = xy(x^4-y^4)$, which has six such representations.


title: Deformation-induced charge redistribution in Ceria thin film at room temperature

abstract: Tuning electronic properties through strain engineering of metal oxides is an important step toward understanding electrochemical and catalytic reactions in energy storage and conversion devices. Traditionally, strain engineering studies focused on movement of oxygen ions at high temperatures (500 {\deg}C and above), complicating electrical properties by introducing mixed electronic and ionic conductivity. In this study, we demonstrate room temperature charge redistribution in a CeO2 thin film as a result of phase transformation in a localized region by mechanical deformation. Mechanical indentation of the CeO2 thin film at room temperature results in irreversible deformation. Conductive-tip atomic force microscopy (C-AFM) analysis shows increased current passing through the locally deformed area of the CeO2 thin film. Electron energy loss spectroscopy (EELS) analysis equipped with transmission electron microscopy (TEM) suggests that the increase in current contrast in the deformed region arises from an increased concentration of Ce3+ ions. We herein discuss the fundamental reason behind the increased amount of Ce3+ ions in the deformed area, based on the atomic scale computational works performed by molecular dynamics (MD) simulations and first-principles density functional theory (DFT) calculations. Plastic deformation induces a phase transformation of cubic fluorite CeO2 into a newly formed T-CeO2 structure. This phase transformation occurs mainly by oxygen ions moving closer to cerium ions to release mechanical energy absorbed in the CeO2 thin film, followed by charge redistribution from the initial CeO2 to the newly created T-CeO2 structure.


title: Non-linear electrodynamics generalization of the rotating BTZ black hole

abstract: We obtain two exact solutions of Einstein gravity coupled to nonlinear electrodynamics (NLED) in $(2+ 1)$-dimensional Anti-de Sitter (AdS) spacetime. The solutions are characterized by the mass $M$, angular momentum $J$, cosmological constant or (anti) de Sitter parameter $\Lambda$, and an electromagnetic parameter $Q$, that is related to an electric field in the first solution and to a magnetic charge for the second solution. Depending on the range of the parameters, the solutions admit a charged rotating asymptotically AdS black hole (BH) interpretation or a charged rotating asymptotically AdS traversable wormhole (WH). If the electromagnetic field is turned off, the stationary Ba\~nados-Teitelboim-Zanelli (BTZ) BH is recovered; in such a way that our BH-WH solutions are nonlinear charged generalizations of the stationary BTZ-BH. Moreover, in contrast to the BTZ metric, the derived AdS solutions are singular at certain radius $r_{s} \neq 0$, resembling the ring singularity of the Kerr-Newman spacetime; while if $\Lambda$ is positive the curvature invariants of the second solution are finite.


title: Congruences on the Partial Automorphism Monoid of a Free Group Action

abstract: We study congruences on the partial automorphism monoid of a finite rank free group action. We give a decomposition of a congruence on this monoid into a Rees congruence, a congruence on a Brandt semigroup and an idempotent separating congruence. The constituent parts are further described in terms of subgroups of direct and semidirect products of groups. We utilize this description to demonstrate how the number of congruences on the partial automorphism monoid depends on the group and on the rank of the action.


title: Effect of top-down connections in Hierarchical Sparse Coding

abstract: Hierarchical Sparse Coding (HSC) is a powerful model to efficiently represent multi-dimensional, structured data such as images. The simplest solution to solve this computationally hard problem is to decompose it into independent layer-wise subproblems. However, neuroscientific evidence would suggest inter-connecting these subproblems as in the Predictive Coding (PC) theory, which adds top-down connections between consecutive layers. In this study, a new model called 2-Layers Sparse Predictive Coding (2L-SPC) is introduced to assess the impact of this inter-layer feedback connection. In particular, the 2L-SPC is compared with a Hierarchical Lasso (Hi-La) network made out of a sequence of independent Lasso layers. The 2L-SPC and the 2-layers Hi-La networks are trained on 4 different databases and with different sparsity parameters on each layer. First, we show that the overall prediction error generated by 2L-SPC is lower thanks to the feedback mechanism as it transfers prediction error between layers. Second, we demonstrate that the inference stage of the 2L-SPC is faster to converge than for the Hi-La model. Third, we show that the 2L-SPC also accelerates the learning process. Finally, the qualitative analysis of both models dictionaries, supported by their activation probability, show that the 2L-SPC features are more generic and informative.


title: An In-Situ Movable Calibration Source for Cryogenic Particle Detectors

abstract: A prototype device capable of moving a radioactive calibration source to multiple positions was operated at millikelvin temperatures using a modified commercial stepper motor. It was developed as an in-situ calibration strategy for cryogenic dark matter detectors. Data taken by scanning a calibration source across multiple radial positions of a prototype dark matter detector demonstrated its functionality. Construction, heat load, and operation of the device are discussed, as is the effect of the motor on the detector operation. A sample dataset taken over multiple positions of a SuperCDMS detector is presented as an example of the utility of such a device.


title: Surface roughness in finite element meshes

abstract: We present a practical approach for constructing meshes of general rough surfaces with given autocorrelation functions based on the unstructured meshes of nominally smooth surfaces. The approach builds on a well-known method to construct correlated random numbers from white noise using a decomposition of the autocorrelation matrix. We discuss important details arising in practical applications to the physicalmodeling of surface roughness and provide a software implementation to enable use of the approach with a broad range of numerical methods in various fields of science and engineering.


title: Inferring population statistics of receptor neurons sensitivities and firing-rates from general functional requirements

abstract: On the basis of the evident ability of neuronal olfactory systems to evaluate the intensity of an odorous stimulus and at the same time also recognise the identity of the odorant over a large range of concentrations, a few biologically-realistic hypotheses on some of the underlying neural processes are made. In particular, it is assumed that the receptor neurons mean firing-rate scale monotonically with odorant intensity, and that the receptor sensitivities range widely across odorants and receptor neurons hence leading to highly distributed representations of the stimuli. The mathematical implementation of the phenomenological postulates allows for inferring explicit functional relationships between some measurable quantities. It results that both the dependence of the mean firing-rate on odorant concentration and the statistical distribution of receptor sensitivity across the neuronal population are power-laws, whose respective exponents are in an arithmetic, testable relationship. In order to test quantitatively the prediction of power-law dependence of population mean firing-rate on odorant concentration, a probabilistic model is created to extract information from data available in the experimental literature. The values of the free parameters of the model are estimated by an info-geometric Bayesian maximum-likelihood inference which keeps into account the prior distribution of the parameters. The eventual goodness of fit is quantified by means of a distribution-independent test. [CONTINUES]


title: A duality between non-compact semisimple symmetric pairs and commutative compact semisimple symmetric triads and its general theory

abstract: The present paper investigates a natural generalization of the duality between Riemannian symmetric pairs of compact type and those of non-compact type \`a la \'E. Cartan. The main result of this paper is to construct an explicit description of a one-to-one correspondence between non-compact pseudo-Riemannian semisimple symmetric pairs and commutative compact semisimple symmetric triads, which is called the duality theorem. Further, we develop a general theory of the duality theorem.


title: Modular Simulation Framework for Process Variation Analysis of MRAM-based Deep Belief Networks

abstract: Magnetic Random-Access Memory (MRAM) based p-bit neuromorphic computing devices are garnering increasing interest as a means to compactly and efficiently realize machine learning operations in Restricted Boltzmann Machines (RBMs). When embedded within an RBM resistive crossbar array, the p-bit based neuron realizes a tunable sigmoidal activation function. Since the stochasticity of activation is dependent on the energy barrier of the MRAM device, it is essential to assess the impact of process variation on the voltage-dependent behavior of the sigmoid function. Other influential performance factors arise from varying energy barriers on power consumption requiring a simulation environment to facilitate the multi-objective optimization of device and network parameters. Herein, transportable Python scripts are developed to analyze the output variation under changes in device dimensions on the accuracy of machine learning applications. Evaluation with RBM circuits using the MNIST dataset reveal impacts and limits for processing variation of device fabrication in terms of the resulting energy vs. accuracy tradeoffs, and the resulting simulation framework is available via a Creative Commons license.


title: CMB and BBN constraints on evaporating primordial black holes revisited

abstract: We derive new CMB anisotropy power spectrum and BBN constraints for evaporating primordial black holes by explicitly solving the electromagnetic particle cascades of emitted particles and the deposition of this emitted energy to the background baryon-photon plasma. We show that the CMB anisotropies can provide stronger constraints compared to BBN and CMB spectral distortions on black holes with masses as small as $M_{\rm BH}=1.1\times 10^{13}$g, a slightly smaller mass than what has been considered in literature until now. We also show that, with more up-to-date data on abundances of deuterium and helium-3, BBN constraints are strengthened significantly.


title: EGO-CH: Dataset and Fundamental Tasks for Visitors BehavioralUnderstanding using Egocentric Vision

abstract: Equipping visitors of a cultural site with a wearable device allows to easily collect information about their preferences which can be exploited to improve the fruition of cultural goods with augmented reality. Moreover, egocentric video can be processed using computer vision and machine learning to enable an automated analysis of visitors' behavior. The inferred information can be used both online to assist the visitor and offline to support the manager of the site. Despite the positive impact such technologies can have in cultural heritage, the topic is currently understudied due to the limited number of public datasets suitable to study the considered problems. To address this issue, in this paper we propose EGOcentric-Cultural Heritage (EGO-CH), the first dataset of egocentric videos for visitors' behavior understanding in cultural sites. The dataset has been collected in two cultural sites and includes more than $27$ hours of video acquired by $70$ subjects, with labels for $26$ environments and over $200$ different Points of Interest. A large subset of the dataset, consisting of $60$ videos, is associated with surveys filled out by real visitors. To encourage research on the topic, we propose $4$ challenging tasks (room-based localization, point of interest/object recognition, object retrieval and survey prediction) useful to understand visitors' behavior and report baseline results on the dataset.


title: Double Hurwitz numbers: polynomiality, topological recursion and intersection theory

abstract: Double Hurwitz numbers enumerate branched covers of $\mathbb{CP}^1$ with prescribed ramification over two points and simple ramification elsewhere. In contrast to the single case, their underlying geometry is not well understood. In previous work by the second- and third-named authors, the double Hurwitz numbers were conjectured to satisfy a polynomiality structure and to be governed by the topological recursion, analogous to existing results concerning single Hurwitz numbers. In this paper, we resolve these conjectures by a careful analysis of the semi-infinite wedge representation for double Hurwitz numbers, by pushing further methods previously used for other Hurwitz problems. We deduce a preliminary version of an ELSV-like formula for double Hurwitz numbers, by deforming the Johnson-Pandharipande-Tseng formula for orbifold Hurwitz numbers and using properties of the topological recursion under variation of spectral curves. In the course of this analysis, we unveil certain vanishing properties of the Chiodo classes.


title: Fragmentation Coagulation Based Mixed Membership Stochastic Blockmodel

abstract: The Mixed-Membership Stochastic Blockmodel~(MMSB) is proposed as one of the state-of-the-art Bayesian relational methods suitable for learning the complex hidden structure underlying the network data. However, the current formulation of MMSB suffers from the following two issues: (1), the prior information~(e.g. entities' community structural information) can not be well embedded in the modelling; (2), community evolution can not be well described in the literature. Therefore, we propose a non-parametric fragmentation coagulation based Mixed Membership Stochastic Blockmodel (fcMMSB). Our model performs entity-based clustering to capture the community information for entities and linkage-based clustering to derive the group information for links simultaneously. Besides, the proposed model infers the network structure and models community evolution, manifested by appearances and disappearances of communities, using the discrete fragmentation coagulation process (DFCP). By integrating the community structure with the group compatibility matrix we derive a generalized version of MMSB. An efficient Gibbs sampling scheme with Polya Gamma (PG) approach is implemented for posterior inference. We validate our model on synthetic and real world data.


title: Sparse Magnetometer-free Inertial Motion Tracking -- A Condition for Observability in Double Hinge Joint Systems

abstract: Inertial measurement units are commonly used in a growing number of application fields to track or capture motions of kinematic chains, such as human limbs, exoskeletons or robotic actuators. A major challenge is the presence of magnetic disturbances that result in unreliable magnetometer readings. Recent research revealed that this problem can be overcome by exploitation of kinematic constraints. While typically each segment of the kinematic chain is equipped with an IMU, a novel approach called sparse inertial motion tracking aims at infering the complete motion states from measurements of a reduced set of sensors. In the present contribution, we combine the magnetometer-free and the sparse approach for real-time motion tracking of double-hinge joint systems with non-parallel joint axes. Analyzing the observability of the system, we find a condition which assures that the relative orientations between all segments are uniquely determined by a kinematic constraint, which contains only the gyroscope readings. Furthermore, we propose a moving-horizon estimator and validate it in a simulation study of three movements with different degrees of excitation. The results of this study confirm all theoretical conjectures and demonstrate that magnetometer-free sparse inertial real-time motion tracking is feasible under precise and simple excitation conditions.


title: Computational Synthesis of Substrates by Crystal Cleavage

abstract: The discovery of novel substrate materials has been dominated by trial and error, opening the opportunity for a systematic search. To identify stable crystal surfaces, we generate bonding networks for materials from the Materials Project database with one to five atoms in the primitive unit cell. For three-dimensional crystals in this set, we systematically break up to three bonds in the bonding network of the primitive cell. Successful cleavage reduces the bonding network to two periodic dimensions, creating a layer of the cleaved crystal. We identify 4,693 unique cleavage surfaces across 2,133 bulk crystals, 4,626 of which have a maximum Miller index of 1. To characterize the likelihood of cleavage and the thermodynamic stability of the cleaved surfaces, we create monolayers of these surfaces and calculate the work of cleavage and the partially-relaxed surface energy using density functional theory to discover 3,991 potential substrates, 2,307 of which do not contain f-valence electrons and 2,183 of which are derived from a bulk precursor with an entry in the Inorganic Crystal Structure Database. Following, we identify distinct trends in the work of cleavage of these layers and relate them to metallic and covalent/ionic bonding of the three-dimensional precursor. We also assembled a database of commercially available substrates and show that the database of predicted substrates significantly enhances the diversity and range of the distribution of electronic properties and lattice parameters, providing opportunities for the epitaxial growth of many materials. We illustrate the potential impact of the substrate database by identifying several new epitaxial substrates for the transparent conductor BaSnO3, which exhibit low cleavage energies and result in strains an order of magnitude lower than currently used substrates.


title: Siamese Neural Networks for EEG-based Brain-computer Interfaces

abstract: Motivated by the inconceivable capability of the human brain in simultaneously processing multi-modal signals and its real-time feedback to the outer world events, there has been a surge of interest in establishing a communication bridge between the human brain and a computer, which are referred to as Brain-computer Interfaces (BCI). To this aim, monitoring the electrical activity of brain through Electroencephalogram (EEG) has emerged as the prime choice for BCI systems. To discover the underlying and specific features of brain signals for different mental tasks, a considerable number of research works are developed based on statistical and data-driven techniques. However, a major bottleneck in the development of practical and commercial BCI systems is their limited performance when the number of mental tasks for classification is increased. In this work, we propose a new EEG processing and feature extraction paradigm based on Siamese neural networks, which can be conveniently merged and scaled up for multi-class problems. The idea of Siamese networks is to train a double-input neural network based on a contrastive loss-function, which provides the capability of verifying if two input EEG trials are from the same class or not. In this work, a Siamese architecture, which is developed based on Convolutional Neural Networks (CNN) and provides a binary output on the similarity of two inputs, is combined with OVR and OVO techniques to scale up for multi-class problems. The efficacy of this architecture is evaluated on a 4-class Motor Imagery (MI) dataset from BCI Competition IV-2a and the results suggest a promising performance compared to its counterparts.


title: Dynamical phase transitions for the activity biased Ising model in a magnetic field

abstract: We consider large deviations of the dynamical activity -- defined as the total number of configuration changes within a time interval -- for mean-field and one-dimensional Ising models, in the presence of a magnetic field. We identify several dynamical phase transitions that appear as singularities in the scaled cumulant generating function of the activity. In particular, we find low-activity ferromagnetic states and a novel high-activity phase, with associated first- and second-order phase transitions. The high-activity phase has a negative susceptibility to the magnetic field. In the mean-field case, we analyse the dynamical phase coexistence that occurs on first-order transition lines, including the optimal-control forces that reproduce the relevant large deviations. In the one-dimensional model, we use exact diagonalisation and cloning methods to perform finite-size scaling of the first-order phase transition at non-zero magnetic field.


title: The CARMENES search for exoplanets around M dwarfs. Dynamical characterization of the multiple planet system GJ 1148 and prospects of habitable exomoons around GJ 1148 b

abstract: Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.


title: The infrared view of dust and molecules around V4334 Sgr (Sakurai's Object): a 20-year retrospective

abstract: We present an analysis of the evolution of circumstellar dust and molecules in the environment of the very late thermal pulse object V4334 Sgr (Sakurai's Object) over a $\sim20$-year period, drawing on ground-, airborne- and space-based infrared photometry and spectroscopy. The dust emission, which started in 1997, resembles a blackbody that cooled from $\sim1200$K in 1998 August to $\sim180$K in 2016 July. The dust mass, assuming amorphous carbon, was $\sim5\times10^{-10}$M$_\odot$ in 1998 August, and we estimate that the total dust mass was $\sim2\times10^{-5}$M$_\odot$ by $\sim2016$. The appearance of a near infrared excess in 2008 suggests a new episode of (or renewed) mass loss began then. We infer lower limits on the bolometric luminosity of the embedded star from that of the dust shell, which rose to $\sim16000$L$_\odot$ before declining to $\sim3000$L$_\odot$. There is evidence for weak 6-7$\mu$m absorption, which we attribute to hydrogenated amorphous carbon formed in material ejected by Sakurai's Object during a mass ejection phase that preceded the 1997 event. We detect small hydrocarbon and other molecules in the spectra, and trace the column densities in hydrogen cyanide (HCN) and acetylene (C$_2$H$_2$). We use the former to determine the $^{12}$C/$^{13}$C ratio to be $6.4\pm0.7$, 14 times smaller than the Solar System value.


title: Infrared properties of Planetary Nebulae with [WR] and $wels$ central stars

abstract: We report the IR properties of planetary nebulae with WR type and $wels$ central stars known to date and compare them with the IR properties of a sample of PNe with H-rich central stars. We use near-, mid- and far-IR photometric data from archives to derive the IR properties of PNe. We have constructed IR colour-colour diagrams of PNe using measurements from $2MASS$, $IRAS$, $WISE$ and $Akari$ bands. [WR] PNe have a larger near-IR emission from the hot dust component and also show a tendency for stronger 12$\mu$m emission as compared to the other two groups. Cool AGB dust properties of all PNe are found to be similar. We derived the dust colour temperatures, dust masses, dust-to-gas mass ratios, IR luminosities and IR excess of PNe for these three groups. [WR] PNe and $wels$-PNe tend to have larger mean values for dust mass when compared to the third group. The average dust-to-gas mass ratio is found to be similar for the three groups of PNe. While there is a strong correlation of dust temperature and IR luminosity with the age for the three groups of PNe, the dust mass, dust-to-gas mass ratios and IR excess are found to be non-varying as the PNe evolve. [WR] PNe and $wels$-PNe show very similar distribution of excitation classes and also show similar distribution with Galactic latitude.


title: Towards Explainable Bit Error Tolerance of Resistive RAM-Based Binarized Neural Networks

abstract: Non-volatile memory, such as resistive RAM (RRAM), is an emerging energy-efficient storage, especially for low-power machine learning models on the edge. It is reported, however, that the bit error rate of RRAMs can be up to 3.3% in the ultra low-power setting, which might be crucial for many use cases. Binary neural networks (BNNs), a resource efficient variant of neural networks (NNs), can tolerate a certain percentage of errors without a loss in accuracy and demand lower resources in computation and storage. The bit error tolerance (BET) in BNNs can be achieved by flipping the weight signs during training, as proposed by Hirtzlin et al., but their method has a significant drawback, especially for fully connected neural networks (FCNN): The FCNNs overfit to the error rate used in training, which leads to low accuracy under lower error rates. In addition, the underlying principles of BET are not investigated. In this work, we improve the training for BET of BNNs and aim to explain this property. We propose straight-through gradient approximation to improve the weight-sign-flip training, by which BNNs adapt less to the bit error rates. To explain the achieved robustness, we define a metric that aims to measure BET without fault injection. We evaluate the metric and find that it correlates with accuracy over error rate for all FCNNs tested. Finally, we explore the influence of a novel regularizer that optimizes with respect to this metric, with the aim of providing a configurable trade-off in accuracy and BET.


title: Is It Safe Living in the Vicinity of Cellular Towers? Analysis of Long-Term Human EMF Exposure at Population Scale

abstract: We focus on the ElectroMagnetic Field (EMF) exposure safety for people living in the vicinity of cellular towers. To this aim, we analyze a large dataset of long-term EMF measurements collected over almost 20 years in more than 2000 measurement points spread over an Italian region. We evaluate the relationship between EMF exposure and the following factors: (i) distance from the closest installation(s), (ii) type of EMF sources in the vicinity, (iii) Base Station (BS) technology, and (iv) EMF regulation updates. Overall, the exposure levels from BSs in the vicinity are below the Italian EMF limits, thus ensuring safety for the population. Moreover, BSs represent the lowest exposure compared to Radio/TV repeaters and other EMF sources. However, the BS EMF exposure in proximity to users exhibits an increasing trend over the last years, which is likely due to the pervasive deployment of multiple technologies and to the EMF regulation updates. As a side consideration, if the EMF levels continue to increase with the current trends, the EMF exposure in proximity to BSs will saturate to the maximum EMF limit by the next 20 years at a distance of 30 meters from the closest BS.


title: L6DNet: Light 6 DoF Network for Robust and Precise Object Pose Estimation with Small Datasets

abstract: Estimating the 3D pose of an object is a challenging task that can be considered within augmented reality or robotic applications. In this paper, we propose a novel approach to perform 6 DoF object pose estimation from a single RGB-D image. We adopt a hybrid pipeline in two stages: data-driven and geometric respectively. The data-driven step consists of a classification CNN to estimate the object 2D location in the image from local patches, followed by a regression CNN trained to predict the 3D location of a set of keypoints in the camera coordinate system. To extract the pose information, the geometric step consists in aligning the 3D points in the camera coordinate system with the corresponding 3D points in world coordinate system by minimizing a registration error, thus computing the pose. Our experiments on the standard dataset LineMod show that our approach is more robust and accurate than state-of-the-art methods. The approach is also validated to achieve a 6 DoF positioning task by visual servoing.


title: A sign pattern with non-zero elements on the diagonal whose minimal rank realizations are not diagonalizable over the complex numbers

abstract: The rank of the $9\times 9$ matrix $$ \left( \begin{array}{cccc|c|cccc} 1&1&0&0&1&0&0&0&0\\ 1&1&0&0&0&0&0&0&0\\ 0&0&1&1&1&0&0&0&0\\ 0&0&1&1&0&0&0&0&0\\\hline 0&0&0&0&1&0&1&0&1\\\hline 0&0&0&0&0&1&1&0&0\\ 0&0&0&0&0&1&1&0&0\\ 0&0&0&0&0&0&0&1&1\\ 0&0&0&0&0&0&0&1&1 \end{array} \right) $$ is $6$. If we replace the ones by arbitrary non-zero numbers, we get a matrix $B$ with $\operatorname{rank} B\geqslant6$, and if $\operatorname{rank} B=6$, the $6\times 6$ principal minors of $B$ vanish.


title: An asymptotic vanishing theorem for the cohomology of thickenings

abstract: Let $X$ be a closed equidimensional local complete intersection subscheme of a smooth projective scheme $Y$ over a field, and let $X_t$ denote the $t$-th thickening of $X$ in $Y$. Fix an ample line bundle $\mathcal{O}_Y(1)$ on $Y$. We prove the following asymptotic formulation of the Kodaira vanishing theorem: there exists an integer $c$, such that for all integers $t \geqslant 1$, the cohomology group $H^k(X_t,\mathcal{O}_{X_t}(j))$ vanishes for $k < \dim X$ and $j < -ct$. Note that there are no restrictions on the characteristic of the field, or on the singular locus of $X$. We also construct examples illustrating that a linear bound is indeed the best possible, and that the constant $c$ is unbounded, even in a fixed dimension.


title: The relative isoperimetric inequality for minimal submanifolds with free boundary in the Euclidean space

abstract: In this paper, we mainly consider the relative isoperimetric inequalities for minimal submanifolds with free boundary. We first generalize ideas of restricted normal cones introduced by Choe-Ghomi-Ritor\'e in \cite{CGR06} and obtain an optimal area estimate for generalized restricted normal cones. This area estimate, together with the ABP method of Cabr\'e in \cite{Cabre2008}, provides a new proof of the relative isoperimetric inequality obtained by Choe-Ghomi-Ritor\'e in \cite{CGR07}. Furthermore, we use this estimate and the idea of Brendle in his recent work \cite{Brendle2019} to obtain a relative isoperimetric inequality for minimal submanifolds with free boundary on a convex support surface in $\mathbb{R}^{n+m}$, which is optimal and gives an affirmative answer to an open problem proposed by Choe in \cite{Choe2005}, Open Problem 12.6, when the codimension $m\leq 2$.


title: Complexity Guarantees for Polyak Steps with Momentum

abstract: In smooth strongly convex optimization, knowledge of the strong convexity parameter is critical for obtaining simple methods with accelerated rates. In this work, we study a class of methods, based on Polyak steps, where this knowledge is substituted by that of the optimal value, $f_*$. We first show slightly improved convergence bounds than previously known for the classical case of simple gradient descent with Polyak steps, we then derive an accelerated gradient method with Polyak steps and momentum, along with convergence guarantees.


title: Geometry of Bend: Singular Lines and Defects in Twist-Bend Nematics

abstract: We describe the geometry of bend distortions in twist-bend nematic liquid crystals in terms of their fundamental degeneracies, which we call $\beta$ lines. These represent a new class of line-like topological defect. We use them to construct and characterise novel structures, including grain boundary and focal conic smectic-like defects, Skyrmions, and knotted merons. We analyse their local geometry and global structure, showing that their intersection with any surface is twice the Skyrmion number. Finally, we demonstrate how arbitrary knots and links can be created and describe them in terms of merons, giving a new geometric perspective on the fractionalisation of Skyrmions.


title: A power Schur complement Low-Rank correction preconditioner for general sparse linear systems

abstract: An effective power based parallel preconditioner is proposed for general large sparse linear systems. The preconditioner combines a power series expansion method with some low-rank correction techniques, where the Sherman-Morrison-Woodbury formula is utilized. A matrix splitting of the Schur complement is proposed to expand the power series. The number of terms used in the power series expansion can control the approximation accuracy of the preconditioner to the inverse of the Schur complement. To construct the preconditioner, graph partitioning is invoked to reorder the original coefficient matrix, leading to a special block two-by-two matrix whose two off-diagonal submatrices are block diagonal. Variables corresponding to interface variables are obtained by solving a linear system with the coeffcient matrix being the Schur complement. For the variables related to the interior variables, one only needs to solve a block diagonal linear system. This can be performed efficiently in parallel. Various numerical examples are provided to illustrate that the efficiency of the proposed preconditioner.


title: BeCAPTCHA: Detecting Human Behavior in Smartphone Interaction using Multiple Inbuilt Sensors

abstract: We introduce a novel multimodal mobile database called HuMIdb (Human Mobile Interaction database) that comprises 14 mobile sensors acquired from 600 users. The heterogeneous flow of data generated during the interaction with the smartphones can be used to model human behavior when interacting with the technology. Based on this new dataset, we explore the capacity of smartphone sensors to improve bot detection. We propose a CAPTCHA method based on the analysis of the information obtained during a single drag and drop task. We evaluate the method generating fake samples synthesized with Generative Adversarial Neural Networks and handcrafted methods. Our results suggest the potential of mobile sensors to characterize the human behavior and develop a new generation of CAPTCHAs.


title: On the signs of Fourier coefficients of Hilbert cusp forms

abstract: We prove that given any $\epsilon > 0$ and a primitive adelic Hilbert cusp form $f$ of weight $k=(k_1,k_2,...,k_n) \in (2 \mathbb{Z})^n$ and full level, there exists an integral ideal $\mathfrak{m}$ with $N(\mathfrak{m}) \ll_{\epsilon} Q_{f}^{9/20+ \epsilon} $ such that the $\mathfrak{m}$-th Fourier coefficient of $C_{f} (\mathfrak{m})$ of $f$ is negative. Here $n$ is the degree of the associated number field, $N(\mathfrak{m})$ is the norm of integral ideal $\mathfrak{m}$ and $Q_{f}$ is the analytic conductor of $f$. In the case of arbitrary weights, we show that there is an integral ideal $\mathfrak{m}$ with $N(\mathfrak{m}) \ll_{\epsilon} Q_{f}^{1/2 + \epsilon}$ such that $C_{f}(\mathfrak{m}) <0$. We also prove that when $k=(k_1,k_2,...,k_n) \in (2 \mathbb{Z})^n$, asymptotically half of the Fourier coefficients are positive while half are negative.


title: Non-linear regression models for behavioral and neural data analysis

abstract: Regression models are popular tools in empirical sciences to infer the influence of a set of variables onto a dependent variable given an experimental dataset. In neuroscience and cognitive psychology, Generalized Linear Models (GLMs) -including linear regression, logistic regression, and Poisson GLM- is the regression model of choice to study the factors that drive participant's choices, reaction times and neural activations. These methods are however limited as they only capture linear contributions of each regressors. Here, we introduce an extension of GLMs called Generalized Unrestricted Models (GUMs), which allows to infer a much richer set of contributions of the regressors to the dependent variable, including possible interactions between the regressors. In a GUM, each regressor is passed through a linear or nonlinear function, and the contribution of the different resulting transformed regressors can be summed or multiplied to generate a predictor for the dependent variable. We propose a Bayesian treatment of these models in which we endow functions with Gaussian Process priors, and we present two methods to compute a posterior over the functions given a dataset: the Laplace method and a sparse variational approach, which scales better for large dataset. For each method, we assess the quality of the model estimation and we detail how the hyperparameters (defining for example the expected smoothness of the function) can be fitted. Finally, we illustrate the power of the method on a behavioral dataset where subjects reported the average perceived orientation of a series of gratings. The method allows to recover the mapping of the grating angle onto perceptual evidence for each subject, as well as the impact of the grating based on its position. Overall, GUMs provides a very rich and flexible framework to run nonlinear regression analysis in neuroscience, psychology, and beyond.


title: Repeated patterns in proper colourings

abstract: For a fixed graph $H$, what is the smallest number of colours $C$ such that there is a proper edge-colouring of the complete graph $K_n$ with $C$ colours containing no two vertex-disjoint colour-isomorphic copies, or repeats, of $H$? We study this function and its generalisation to more than two copies using a variety of combinatorial, probabilistic and algebraic techniques. For example, we show that for any tree $T$ there exists a constant $c$ such that any proper edge-colouring of $K_n$ with at most $c n^2$ colours contains two repeats of $T$, while there are colourings with at most $c' n^{3/2}$ colours for some absolute constant $c'$ containing no three repeats of any tree with at least two edges. We also show that for any graph $H$ containing a cycle there exist $k$ and $c$ such that there is a proper edge-colouring of $K_n$ with at most $c n$ colours containing no $k$ repeats of $H$, while, for a tree $T$ with $m$ edges, a colouring with $o(n^{(m+1)/m})$ colours contains $\omega(1)$ repeats of $T$.


title: A Neural-embedded Choice Model: TasteNet-MNL Modeling Taste Heterogeneity with Flexibility and Interpretability

abstract: Discrete choice models (DCMs) require a priori knowledge of the utility functions, especially how tastes vary across individuals. Utility misspecification may lead to biased estimates, inaccurate interpretations and limited predictability. In this paper, we utilize a neural network to learn taste representation. Our formulation consists of two modules: a neural network (TasteNet) that learns taste parameters (e.g., time coefficient) as flexible functions of individual characteristics; and a multinomial logit (MNL) model with utility functions defined with expert knowledge. Taste parameters learned by the neural network are fed into the choice model and link the two modules. Our approach extends the L-MNL model (Sifringer et al., 2020) by allowing the neural network to learn the interactions between individual characteristics and alternative attributes. Moreover, we formalize and strengthen the interpretability condition - requiring realistic estimates of behavior indicators (e.g., value-of-time, elasticity) at the disaggregated level, which is crucial for a model to be suitable for scenario analysis and policy decisions. Through a unique network architecture and parameter transformation, we incorporate prior knowledge and guide the neural network to output realistic behavior indicators at the disaggregated level. We show that TasteNet-MNL reaches the ground-truth model's predictability and recovers the nonlinear taste functions on synthetic data. Its estimated value-of-time and choice elasticities at the individual level are close to the ground truth. On a publicly available Swissmetro dataset, TasteNet-MNL outperforms benchmarking MNLs and Mixed Logit model's predictability. It learns a broader spectrum of taste variations within the population and suggests a higher average value-of-time.


title: Repr\'esentations de r\'eflexion de groupes de Coxeter Deuxi\`eme partie: outils pour des exemples

abstract: This part is made of three sections. In the first section we study the family of polynomials whose roots are 4cos2 k{\pi}, (n \geqslant 3,1 \leqslant k < \frac{n}{2}). We obtain n2 in this manner a family of orthogonal polynomials. This will permit us to study in details all the examples which follow. In the second section, we give technical formulae in order noto repeat calculations. In the third section,we give applications, first when the field K is a sub-field of R (presentations of W(H3) and W(H4)) then in the complex case (study of the complex reflection group G(p, p, n), G24 and G27).


title: Within-sample variability-invariant loss for robust speaker recognition under noisy environments

abstract: Despite the significant improvements in speaker recognition enabled by deep neural networks, unsatisfactory performance persists under noisy environments. In this paper, we train the speaker embedding network to learn the "clean" embedding of the noisy utterance. Specifically, the network is trained with the original speaker identification loss with an auxiliary within-sample variability-invariant loss. This auxiliary variability-invariant loss is used to learn the same embedding among the clean utterance and its noisy copies and prevents the network from encoding the undesired noises or variabilities into the speaker representation. Furthermore, we investigate the data preparation strategy for generating clean and noisy utterance pairs on-the-fly. The strategy generates different noisy copies for the same clean utterance at each training step, helping the speaker embedding network generalize better under noisy environments. Experiments on VoxCeleb1 indicate that the proposed training framework improves the performance of the speaker verification system in both clean and noisy conditions.


title: Extremes of the 2d scale-inhomogeneous discrete Gaussian free field: Extremal process in the weakly correlated regime

abstract: We prove convergence of the full extremal process of the two-dimensional scale-inhomogeneous discrete Gaussian free field in the weak correlation regime. The scale-inhomogeneous discrete Gaussian free field is obtained from the 2d discrete Gaussian free field by modifying the variance through a function $\mathcal{I}:[0,1]\rightarrow [0,1]$. The limiting process is a cluster Cox process. The random intensity of the Cox process depends on the $\mathcal{I}^\prime(0)$ through a random measure $Y$ and on the $\mathcal{I}^\prime(1)$ through a constant $\beta$. We describe the cluster process, which only depends on $\mathcal{I}^\prime(1)$, as points of a standard 2d discrete Gaussian free field conditioned to be unusually high.


title: Search for the $H \rightarrow \text{W}^+ \text{W}^-$ process at the LHeC Experiment

abstract: We consider the decay of the Higgs boson to W$^{+}$ W$^{-}$ at a proposed Large Hadron Electron Collider and determine the likelihood of detecting a signal for the Higgs mass from its decay product W jets by imposing cuts to select candidate jet pairs and optimizing the value of the angular separation $\Delta R$. It was found that at the LHeC experiment (CM energy $\sqrt{s}=1.3$ TeV and luminosity of 100 fb$^{-1}$ per year), the highest efficiency is obtained with $\Delta R = 0.4$, along with a selection scheme of $|\Delta \eta| <1, 10<m<85$ GeV, $p_T$ of jets 1 and 2 between $10 - 20$ GeV and $p_T$ of jets 3 and 4 $>10$ GeV: this led to an efficiency between $7.1 - 7.5 \%$ for finding the invariant 4-jet mass in a mass region $<140$ GeV. Under signal-to-background comparison, the signal showed a $3.8 \sigma$ excess compared to the charged current W$^{-}$ background.


title: Beurling numbers whose number of prime factors lies in a specified residue class

abstract: We find asymptotics for $S_{K,c}(x)$, the number of positive integers below $x$ whose number of prime factors is $c \; \mathrm{mod}\; K$. We study this question in the context of Beurling integers.


title: A Model of Electron Transport in Dense Plasmas Spanning Temperature Regimes

abstract: We present a new model of electron transport in warm and hot dense plasmas which combines the quantum Landau-Fokker-Planck equation with the concept of mean-force scattering. We obtain electrical and thermal conductivities across several orders of magnitude in temperature, from warm dense matter conditions to hot, nondegenerate plasma conditions, including the challenging crossover regime between the two. The small-angle approximation characteristic of Fokker-Planck collision theories is mitigated to good effect by the construction of accurate effective Coulomb logarithms based on mean-force scattering, which allows the theory to remain accurate even at low temperatures, as compared with high-fidelity quantum simulation results. Electron-electron collisions are treated on equal footing as electron-ion collisions. Their accurate treatment is found to be essential for hydrogen, and is expected to be important to other low-Z elements. We find that electron-electron scattering remains influential to the value of the thermal conductivity down to temperatures somewhat below the Fermi energy. The accuracy of the theory seems to falter only for the behavior of the thermal conductivity at very low temperatures due to a subtle interplay between the Pauli exclusion principle and the small-angle approximation as they pertain to electron-electron scattering. Even there, the model is in fair agreement with ab initio simulations.


title: 3D propagation of relativistic solar protons through interplanetary space

abstract: Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons through the Interplanetary Magnetic Field (IMF) using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. The latter are influenced by the IMF polarity. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the polarity and position of the observer, and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, in many cases resulting in a roll-over in the fluence spectrum.


title: Disorder Detection with Costly Observations

abstract: We study the Wiener disorder detection problem where each observation is associated with a positive cost. In this setting, a strategy is a pair consisting of a sequence of observation times and a stopping time corresponding to the declaration of disorder. We characterize the minimal cost of the disorder problem with costly observations as the unique fix-point of a certain jump operator, and we determine the optimal strategy.


title: Stochastic reaction networks in dynamic compartment populations

abstract: Compartmentalization of biochemical processes underlies all biological systems, from the organelle to the tissue scale. Theoretical models to study the interplay between noisy reaction dynamics and compartmentalization are sparse, and typically very challenging to analyze computationally. Recent studies have made progress towards addressing this problem in the context of concrete biological systems but general approaches remain lacking. In this work we propose a mathematical framework based on counting processes that allows us to study compartment populations with arbitrary interactions and internal biochemistry. We provide an efficient description of the population dynamics in terms of differential equations which capture moments of the population and their variability. We demonstrate the relevance of our approach using several case studies inspired by biological systems at different scales.


title: Evaluation of Energy- and Capacity-Market Revenues from Lithium-ion Battery Systems for Offshore Wind Using Advanced Battery Models

abstract: Revenue potential from offshore wind and energy storage systems for a Long Island node in the New York ISO (NYISO) is examined using advanced lithium-ion battery representations. These advanced mixed-integer-linear battery models account for the dynamic performance, as well as the degradation behavior of the batteries, which are usually not accounted for in power systems models. Multiple hybrid offshore wind and battery system designs are investigated to examine the impact of locating the battery offshore versus locating it onshore. For the examined systems, we explore different battery usable state-of-charge (SOC) windows, and corresponding dispatch of the battery to maximize energy- and capacity-market revenues. The impacts of variability of offshore wind output along with energy- and capacity-market prices are evaluated using publicly available data from 2010 to 2013. Locating the battery onshore resulted in higher revenues. For 2013, results highlight that without accurate battery representations, models can overestimate battery revenues by up to 155%, resulting primarily from degradation-related costs. Using advanced algorithms, net revenue can be increased by 29%. Results also indicate that wider useable SOC windows could lead to higher net revenues from the energy market, due to higher arbitrage opportunities that compensate for any additional degradation-tied costs in higher DODs. The added value of a MWh of energy storage varies from $2 to $3.5 per MWh of wind energy, which leads to a breakeven cost range of $50-$95 per kWh for the battery systems studied. As such, energy- and capacity-market revenues were found to be insufficient in recovering the investment costs of current battery systems for the applications considered in this analysis.


title: Hyperk\"ahler cones and instantons on quaternionic K\"ahler manifolds

abstract: We present a novel approach to the study of Yang-Mills instantons on quaternionic K\"ahler manifolds, based on an extension of the harmonic space method of constructing instantons on hyperk\"ahler manifolds. Our results establish a bijection between local equivalence classes of instantons on quaternionic K\"ahler manifolds M and equivalence classes of certain holomorphic maps on an appropriate SL_2(C)-bundle over the Swann bundle of M.


title: Stop Tracking Me Bro! Differential Tracking Of User Demographics On Hyper-partisan Websites

abstract: Websites with hyper-partisan, left or right-leaning focus offer content that is typically biased towards the expectations of their target audience. Such content often polarizes users, who are repeatedly primed to specific (extreme) content, usually reflecting hard party lines on political and socio-economic topics. Though this polarization has been extensively studied with respect to content, it is still unknown how it associates with the online tracking experienced by browsing users, especially when they exhibit certain demographic characteristics. For example, it is unclear how such websites enable the ad-ecosystem to track users based on their gender or age. In this paper, we take a first step to shed light and measure such potential differences in tracking imposed on users when visiting specific party-line's websites. For this, we design and deploy a methodology to systematically probe such websites and measure differences in user tracking. This methodology allows us to create user personas with specific attributes like gender and age and automate their browsing behavior in a consistent and repeatable manner. Thus, we systematically study how personas are being tracked by these websites and their third parties, especially if they exhibit particular demographic properties. Overall, we test 9 personas on 556 hyper-partisan websites and find that right-leaning websites tend to track users more intensely than left-leaning, depending on user demographics, using both cookies and cookie synchronization methods and leading to more costly delivered ads.