title: Current-mediated synchronization of a pair of beating non-identical flagella

abstract: The basic phenomenology of experimentally observed synchronization (i.e., a stochastic phase locking) of identical, beating flagella of a biflagellate alga is known to be captured well by a minimal model describing the dynamics of coupled, limit-cycle, noisy oscillators (known as the noisy Kuramoto model). As demonstrated experimentally, the amplitudes of the noise terms therein, which stem from fluctuations of the rotary motors, depend on the flagella length. Here we address the conceptually important question which kind of synchrony occurs if the two flagella have different lengths such that the noises acting on each of them have different amplitudes. On the basis of a minimal model, too, we show that a different kind of synchrony emerges, and here it is mediated by a current carrying, steady-state; it manifests itself via correlated "drifts" of phases. We quantify such a synchronization mechanism in terms of appropriate order parameters $Q$ and $Q_{\cal S}$ - for an ensemble of trajectories and for a single realization of noises of duration ${\cal S}$, respectively. Via numerical simulations we show that both approaches become identical for long observation times ${\cal S}$. This reveals an ergodic behavior and implies that a single-realization order parameter $Q_{\cal S}$ is suitable for experimental analysis for which ensemble averaging is not always possible.


title: Bouligand-Levenberg-Marquardt iteration for a non-smooth ill-posed inverse problem

abstract: In this paper, we consider a modified Levenberg--Marquardt method for solving an ill-posed inverse problem where the forward mapping is not G\^ateaux differentiable. By relaxing the standard assumptions for the classical smooth setting, we derive asymptotic stability estimates that are then used to prove the convergence of the proposed method. This method can be applied to an inverse source problem for a non-smooth semilinear elliptic PDE where a Bouligand subdifferential can be used in place of the non-existing Fr\'echet derivative, and we show that the corresponding Bouligand-Levenberg-Marquardt iteration is an iterative regularization scheme. Numerical examples illustrate the advantage over the corresponding Bouligand-Landweber iteration.


title: A tutorial on group effective connectivity analysis, part 1: first level analysis with DCM for fMRI

abstract: Dynamic Causal Modelling (DCM) is the predominant method for inferring effective connectivity from neuroimaging data. In the 15 years since its introduction, the neural models and statistical routines in DCM have developed in parallel, driven by the needs of researchers in cognitive and clinical neuroscience. In this tutorial, we step through an exemplar fMRI analysis in detail, reviewing the current implementation of DCM and demonstrating recent developments in group-level connectivity analysis. In the first part of the tutorial (current paper), we focus on issues specific to DCM for fMRI, unpacking the relevant theory and highlighting practical considerations. In particular, we clarify the assumptions (i.e., priors) used in DCM for fMRI and how to interpret the model parameters. This tutorial is accompanied by all the necessary data and instructions to reproduce the analyses using the SPM software. In the second part (in a companion paper), we move from subject-level to group-level modelling using the Parametric Empirical Bayes framework, and illustrate how to test for commonalities and differences in effective connectivity across subjects, based on imaging data from any modality.


title: On the Suppression and Distortion of Non-Equilibrium Fluctuations by Transpiration

abstract: A fluid in a non-equilibrium state exhibits long-ranged correlations of its hydrodynamic fluctuations. In this article, we examine the effect of a transpiration interface on these correlations -- specifically, we consider a dilute gas in a domain bisected by the interface. The system is held in a non-equilibrium steady state by using isothermal walls to impose a temperature gradient. The gas is simulated using both direct simulation Monte Carlo (DSMC) and fluctuating hydrodynamics (FHD). For the FHD simulations two models are developed for the interface based on master equation and Langevin approaches. For appropriate simulation parameters, good agreement is observed between DSMC and FHD results with the latter showing a significant advantage in computational speed. For each approach we quantify the effects of transpiration on long-ranged correlations in the hydrodynamic variables.


title: Fossil field decay due to nonlinear tides in massive binaries

abstract: Surface magnetic fields have been detected in 5 to 10% of isolated massive stars, hosting outer radiative envelopes. They are often thought to have a fossil origin, resulting from the stellar formation phase. Yet, magnetic massive stars are scarcer in (close) short-period binaries, as reported by the BinaMIcS (Binarity and Magnetic Interaction in various classes of Stars) collaboration. Different physical conditions in the molecular clouds giving birth to isolated stars and binaries are commonly invoked. In addition, we propose that the observed lower magnetic incidence in close binaries may be due to nonlinear tides. Indeed, close binaries are probably prone to tidal instability, a fluid instability growing upon the equilibrium tidal flow via nonlinear effects. Yet, stratified effects have hitherto been largely overlooked. We theoretically and numerically investigate tidal instability in rapidly rotating, stably stratified fluids permeated by magnetic fields. We use the short-wavelength stability method to propose a comprehensive (local) theory of tidal instability at the linear onset, discussing damping effects. Then, we propose a mixing-length theory for the mixing generated by tidal instability in the nonlinear regime. We successfully assess our theoretical predictions against proof-of-concept, direct numerical simulations. Finally, we compare our predictions with the observations of short-period, double-lined spectroscopic binary systems. Using new analytical results, cross-validated by a direct integration of the stability equations, we show that tidal instability can be generated by nonlinear couplings of inertia-gravity waves with the equilibrium tidal flow in short-period massive binaries, even against the Joule diffusion. In the nonlinear regime, a fossil magnetic field can be dissipated by the turbulent magnetic diffusion induced by the saturated tidal flows. We predict that the turbulent Joule diffusion of fossil fields would occur in a few million years for several short-period massive binaries. Therefore, turbulent tidal flows could explain the observed dearth of some short-period magnetic binaries.


title: A simple master Theorem for discrete divide and conquer recurrences

abstract: The aim of this note is to provide a Master Theorem for some discrete divide and conquer recurrences: $$X_{n}=a_n+\sum_{j=1}^m b_j X_{\lfloor p_j n \rfloor},$$ where the $p_i$'s belong to $(0,1)$. The main novelty of this work is there is no assumption of regularity or monotonicity for $(a_n)$. Then, this result can be applied to various sequences of random variables $(a_n)_{n\ge 0}$, for example such that $\sup_{n\ge 1}\mathbb{E}(|a_n|)<+\infty$.


title: Trion formation and unconventional superconductivity in a three-dimensional model with short-range attraction

abstract: A three-fermion problem in a three-dimensional lattice with anisotropic hopping is solved by discretizing the Schroedinger equation in momentum space. Interparticle interaction comprises on-site Hubbard repulsion and in-plane nearest-neighbor attraction. By comparing the energy of three-fermion bound clusters (trions) with the energy of one pair plus one free particle, a trion formation threshold is accurately determined, and the region of pair stability is mapped out. It is found that the "close-packed" density of fermion pairs is highest in a strongly anisotropic model. It is also argued that pair superconductivity with the highest critical temperature is always close to trion formation, which makes the system prone to phase separation and local charge ordering.


title: Functoriality of Moduli Spaces of Global $\mathbb G$-Shtukas

abstract: Moduli spaces of global $\mathbb G$-shtukas play a crucial role in the Langlands program for function fields. We analyze their functoriality properties following a change of the curve and a change of the group scheme $\mathbb G$ under various aspects. In particular, we prove two finiteness results which are of interest in the study of stratifications of these moduli spaces and which potentially allow the formulation of an analog of the Andr\'{e}-Oort conjecture for global $\mathbb G$-shtukas.


title: Multivariate Alexander quandles, II. The involutory medial quandle of a link (corrected)

abstract: Joyce showed that for a classical knot $K$, the involutory medial quandle $\text{IMQ}(K)$ is isomorphic to the core quandle of the homology group $H_1(X_2)$, where $X_2$ is the cyclic double cover of $\mathbb S ^3$, branched over $K$. It follows that $|\text{IMQ}(K)| = | \det K |$. In the present paper, the extension of Joyce's result to classical links is discussed. Among other things, we show that for a classical link $L$ of $\mu \geq 2$ components, the order of the involutory medial quandle is bounded as follows: \[ \frac{\mu | \det L |}{2} \geq |\text{IMQ}(L)| \geq \frac{ \mu | \det L |} {2^{\mu -1}}. \] In particular, $\text{IMQ}(L)$ is infinite if and only if $\det L =0$. We also show that in general, $\text{IMQ}(L)$ is a strictly stronger invariant than $H_1(X_2)$. That is, if $L$ and $L'$ are links with $\text{IMQ}(L) \cong \text{IMQ}(L')$, then $H_1(X_2) \cong H_1(X'_2)$; but it is possible to have $H_1(X_2) \cong H_1(X'_2)$ and $\text{IMQ}(L) \not \cong \text{IMQ}(L')$. In fact, it is possible to have $X_2 \cong X'_2$ and $\text{IMQ}(L) \not \cong \text{IMQ}(L')$.


title: A tutorial on group effective connectivity analysis, part 2: second level analysis with PEB

abstract: This tutorial provides a worked example of using Dynamic Causal Modelling (DCM) and Parametric Empirical Bayes (PEB) to characterise inter-subject variability in neural circuitry (effective connectivity). This involves specifying a hierarchical model with two or more levels. At the first level, state space models (DCMs) are used to infer the effective connectivity that best explains a subject's neuroimaging timeseries (e.g. fMRI, MEG, EEG). Subject-specific connectivity parameters are then taken to the group level, where they are modelled using a General Linear Model (GLM) that partitions between-subject variability into designed effects and additive random effects. The ensuing (Bayesian) hierarchical model conveys both the estimated connection strengths and their uncertainty (i.e., posterior covariance) from the subject to the group level; enabling hypotheses to be tested about the commonalities and differences across subjects. This approach can also finesse parameter estimation at the subject level, by using the group-level parameters as empirical priors. We walk through this approach in detail, using data from a published fMRI experiment that characterised individual differences in hemispheric lateralization in a semantic processing task. The preliminary subject specific DCM analysis is covered in detail in a companion paper. This tutorial is accompanied by the example dataset and step-by-step instructions to reproduce the analyses.


title: Maximum Likelihood Estimation of Sparse Networks with Missing Observations

abstract: Estimating the matrix of connections probabilities is one of the key questions when studying sparse networks. In this work, we consider networks generated under the sparse graphon model and the in-homogeneous random graph model with missing observations. Using the Stochastic Block Model as a parametric proxy, we bound the risk of the maximum likelihood estimator of network connections probabilities , and show that it is minimax optimal. When risk is measured in Frobenius norm, no estimator running in polynomial time has been shown to attain the minimax optimal rate of convergence for this problem. Thus, maximum likelihood estimation is of particular interest as computationally efficient approximations to it have been proposed in the literature and are often used in practice.


title: Vanishing relaxation time limit of the Jordan--Moore--Gibson--Thompson wave equation with Neumann and absorbing boundary conditions

abstract: We study the Jordan--Moore--Gibson--Thompson (JMGT) equation, a third order in time wave equation that models nonlinear sound propagation, in the practically relevant setting of Neumann and absorbing boundary conditions. In the analysis, we pay special attention to dependencies on the coefficient $\tau$ of the third order time derivative that plays the physical role of relaxation time. We establish local in time well-posedness and derive energy bounds that can be made independent of $\tau$ under appropriate conditions. This fact allows us to pass to the limit $\tau\to0$ and recover solutions of a classical model in nonlinear acoustics, the Westervelt equation, as singular limits of solutions to the JMGT equation.


title: Necessary and Sufficient Conditions for Passivity of Velocity-Sourced Impedance Control of Series Elastic Actuators

abstract: Series Elastic Actuation (SEA) has become prevalent in applications involving physical human-robot interaction as it provides considerable advantages over traditional stiff actuators in terms of stability robustness and fidelity of force control. Several impedance control architectures have been proposed for SEA. Among these alternatives, the cascaded controller with an inner-most velocity loop, an intermediate torque loop and an outer-most impedance loop is particularly favoured for its simplicity, robustness, and performance. In this paper, we derive the \emph{necessary and sufficient conditions} to ensure the passivity of this cascade-controller architecture for rendering two most common virtual impedance models. Based on the newly established passivity conditions, we provide non-conservative design guidelines to haptically display a null impedance and a pure spring while ensuring the passivity of interaction. We also demonstrate the importance of including physical damping in the actuator model during derivation of passivity conditions, when integral controllers are utilized. In particular, we show the adversary effect of physical damping on system passivity.


title: Energy Constraints for Climate Models from Hydrologic Partitioning

abstract: It is well known that evaporative cooling of Earth's surface water reduces the amount of radiation that goes into sensible heat, namely the portion of radiation that produces higher temperatures. However, a rigorous use of long-term hydrologic measurements and the related theories of hydrologic partitioning have not yet been fully exploited to quantify these effects on climate. Here, we show that the Budyko's curve, a well-known and efficient framework for water balance estimation, can be effectively utilized to partition the surface energy fluxes by expressing the long-term evaporative fraction as a function of the dryness index. The combination of this energy partitioning method with hydrological observations allows us to estimate the surface energy components at watershed and continental scales. Analyzing climate model outputs through this new lens reveals energy biases due to inaccurate parameterization of hydrological and atmospheric processes, offering insight into model parameterization and providing useful information for improving climate projections.


title: Enhanced convergence rates and asymptotics for a dispersive Boussinesq-type system with large ill-prepared data

abstract: In this article we prove highly improved and flexible Strichartz-type estimates allowing us to generalize the asymptotics we obtained for a stratified and rotating incompressible Navier-Stokes system: for large (and less regular) initial data, we obtain global well-posedness, asymptotics (as the Rossby number $\epsilon$ goes to zero) and convergence rates as a power of the small parameter $\epsilon$. Our approach is lead by the special structure of the limit system: the 3D quasi-geostrophic system.


title: Growth of layered Lu$_2$Fe$_3$O$_7$ and Lu$_3$Fe$_4$O$_{10} $ single crystals exhibiting long-range charge order via the optical floating-zone method

abstract: We report the controlled growth of single crystals of intercalated layered Lu$ _{1+n} $Fe$ _{2+n} $O$ _{4+3n-\delta} $ ($ n $=1,2) with different oxygen stoichiometries ${\delta}$. For the first time crystals sufficiently stoichiometric to exhibit superstructure reflections in X-ray diffraction attributable to charge ordering were obtained. The estimated correlation lengths tend to be smaller than for not intercalated LuFe$ _2 $O$ _4$. For Lu$ _2 $Fe$ _3 $O$ _7 $, two different superstructures were observed, one an incommensurate zigzag pattern similar to previous observations by electron diffraction, the other an apparently commensurate pattern with ($\frac{1}{3}\frac{1}{3}0$) propagation. Implications for the possible charge order in the bilayers are discussed. Magnetization measurements suggest reduced magnetic correlations and the absence of an antiferromagnetic phase.


title: Keck HIRES Spectroscopy of SkyMapper Commissioning Survey Candidate Extremely Metal-Poor Stars

abstract: We present results from the analysis of high-resolution spectra obtained with the Keck HIRES spectrograph for a sample of 17 candidate extremely metal-poor (EMP) stars originally selected from commissioning data obtained with the SkyMapper telescope. Fourteen of the stars have not been observed previously at high dispersion. Three have [Fe/H]<=-3.0 while the remainder, with two more metal-rich exceptions, have -3.0<=[Fe/H]<=-2.0 dex. Apart from Fe, we also derive abundances for the elements C, N, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, and Zn, and for n-capture elements Sr, Ba, and Eu. None of the current sample of stars is found to be carbon-rich. In general our chemical abundances follow previous trends found in the literature, although we note that two of the most metal-poor stars show very low [Ba/Fe] (~-1.7) coupled with low [Sr/Ba] (~-0.3). Such stars are relatively rare in the Galactic halo. One further star, and possibly two others, meet the criteria for classification as a r-I star. This study, together with that of Jacobson et al. (2015), completes the outcomes of the SkyMapper commissioning data survey for EMP stars.


title: Helical Twisting Number and Braiding Linkage Number of Solar Coronal Loops

abstract: Coronal loops in active regions are often characterized by quasi-circular and helically twisted (sigmoidal) geometries, which are consistent with dipolar potential field models in the former case, and with nonlinear force-free field models with vertical currents in the latter case. Alternatively, Parker-type nanoflare models of the solar corona hypothesize that a braiding mechanism operates between unresolved loop strands, which is a more complex topological model. In this study we use the vertical-current approximation of a nonpotential magnetic field solution (that fulfills the divergence-free and force-free conditions) to characterize the number of helical turns $N_{twist}$ in twisted coronal loops. We measure the helical twist in 15 active regions observed with AIA and HMI/SDO and find a mean nonpotentiality angle (between the potential and nonpotential field directions) of $\mu_{NP} = 15^\circ \pm 3^\circ$. The resulting mean rotational twist angle is $\varphi = 49^\circ \pm 11^\circ$, which corresponds to $N_{twist}=\varphi/360^\circ = 0.14\pm0.03$ turns with respect to the untwisted potential field, with an absolute upper limit of $N_{twist} \lapprox 0.5$, which is far below the kink instability limit of $|N_{twist}| \gapprox 1$. The number of twist turns $N_{twist}$ corresponds to the Gauss linkage number $N_{link}$ in braiding topologies. We conclude that any braided topology (with $|N_{link}| \ge 1$) cannot explain the observed stability of loops in a force-free corona, nor the observed low twist number. Parker-type nanoflaring can thus occur in non-forcefree environments only, such as in the chromosphere and transition region.


title: Bayesian data fusion for unmeasured confounding

abstract: Bayesian causal inference offers a principled approach to policy evaluation of proposed interventions on mediators or time-varying exposures. We outline a general approach to the estimation of causal quantities for settings with time-varying confounding, such as exposure-induced mediator-outcome confounders. We further extend this approach to propose two Bayesian data fusion (BDF) methods for unmeasured confounding. Using informative priors on quantities relating to the confounding bias parameters, our methods incorporate data from an external source where the confounder is measured in order to make inferences about causal estimands in the main study population. We present results from a simulation study comparing our data fusion methods to two common frequentist correction methods for unmeasured confounding bias in the mediation setting. We also demonstrate our method with an investigation of the role of stage at cancer diagnosis in contributing to Black-White colorectal cancer survival disparities.


title: Spherical centroid bodies

abstract: The spherical centroid body of a centrally-symmetric convex body in the Euclidean unit sphere is introduced. Two alternative definitions - one geometric, the other probabilistic in nature - are given and shown to lead to the same objects. The geometric approach is then used to establish a number of basic properties of spherical centroid bodies, while the probabilistic approach inspires the proof of a spherical analogue of the classical polar Busemann-Petty centroid inequality.


title: Upper-Confidence Bound for Channel Selection in LPWA Networks with Retransmissions

abstract: In this paper, we propose and evaluate different learning strategies based on Multi-Arm Bandit (MAB) algorithms. They allow Internet of Things (IoT) devices to improve their access to the network and their autonomy, while taking into account the impact of encountered radio collisions. For that end, several heuristics employing Upper-Confident Bound (UCB) algorithms are examined, to explore the contextual information provided by the number of retransmissions. Our results show that approaches based on UCB obtain a significant improvement in terms of successful transmission probabilities. Furthermore, it also reveals that a pure UCB channel access is as efficient as more sophisticated learning strategies.


title: Anharmonicity can enhance the performance of quantum refrigerators

abstract: We explore a thermodynamical effect of anharmonicity in quantum mechanical oscillators. We show that small quartic perturbations to the oscillator potential lead to an enhancement of performance of quantum refrigerators for both the Otto and Stirling cycles. A similar nonlinearity driven enhancement of performance is also observed for an analogous spin-qubit model of quantum refrigerators. We further demonstrate the robustness of improvement of the coefficient of performance versus the energy cost for creating anharmonicity. It is shown that the anharmonicity driven improvement in performance is a generic effect at the quantum level for the experimentally realizable Otto refrigerator.


title: Detection of spike-like structures near the front of type-II bursts

abstract: Aims. We examine high time resolution dynamic spectra for fine structures in type II solar radio bursts Methods. We used data obtained with the (SAO) receiver of the Artemis-JLS (ARTEMIS-IV) solar radio spectrograph in the 450-270 MHz range at 10 ms cadence and identified more than 600 short, narrowband features. Their characteristics, such as instantaneous relative bandwidth and total duration were measured and compared with those of spikes embedded in type IV emissions. Results. Type II associated spikes occur mostly in chains inside or close to the slowly drifting type II emission. These spikes coexist with herringbone and pulsating structures. Their average duration is 96 ms and their average relative bandwidth 1.7%. These properties are not different from those of type IV embedded spikes. It is therefore possible that they are signatures of small-scale reconnection along the type II shock front.


title: Still a Pain in the Neck: Evaluating Text Representations on Lexical Composition

abstract: Building meaningful phrase representations is challenging because phrase meanings are not simply the sum of their constituent meanings. Lexical composition can shift the meanings of the constituent words and introduce implicit information. We tested a broad range of textual representations for their capacity to address these issues. We found that as expected, contextualized word representations perform better than static word embeddings, more so on detecting meaning shift than in recovering implicit information, in which their performance is still far from that of humans. Our evaluation suite, including 5 tasks related to lexical composition effects, can serve future research aiming to improve such representations.


title: Learning Factored Markov Decision Processes with Unawareness

abstract: Methods for learning and planning in sequential decision problems often assume the learner is aware of all possible states and actions in advance. This assumption is sometimes untenable. In this paper, we give a method to learn factored markov decision problems from both domain exploration and expert assistance, which guarantees convergence to near-optimal behaviour, even when the agent begins unaware of factors critical to success. Our experiments show our agent learns optimal behaviour on small and large problems, and that conserving information on discovering new possibilities results in faster convergence.


title: Singular stochastic integral operators

abstract: In this paper we introduce Calder\'on-Zygmund theory for singular stochastic integrals with operator-valued kernel. In particular, we prove $L^p$-extrapolation results under a H\"ormander condition on the kernel. Sparse domination and sharp weighted bounds are obtained under a Dini condition on the kernel, leading to a stochastic version of the solution to the $A_2$-conjecture. The results are applied to obtain $p$-independence and weighted bounds for stochastic maximal $L^p$-regularity both in the complex and real interpolation scale. As a consequence we obtain several new regularity results for the stochastic heat equation on $\mathbb{R}^d$ and smooth and angular domains.


title: Coexistence of ferromagnetic and stripe antiferromagnetic spin fluctuations in SrCo$_2$As$_2$

abstract: We use inelastic neutron scattering to study energy and wave vector dependence of spin fluctuations in SrCo$_2$As$_2$, derived from SrFe$_{2-x}$Co$_x$As$_2$ iron pnictide superconductors. Our data reveals the coexistence of antiferromagnetic (AF) and ferromagnetic (FM) spin fluctuations at wave vectors $\textbf{Q}_{\rm AF}$=(1,0) and $\textbf{Q}_{\rm FM}$=(0,0)/(2,0), respectively. By comparing neutron scattering results with those of dynamic mean field theory calculation and angle-resolved photoemission spectroscopy experiments, we conclude that both AF and FM spin fluctuations in SrCo$_2$As$_2$ are closely associated with a flat band of the $e_g$ orbitals near the Fermi level, different from the $t_{2g}$ orbitals in superconducting SrFe$_{2-x}$Co$_x$As$_2$. Therefore, Co-substitution in SrFe$_{2-x}$Co$_x$As$_2$ induces a $t_{2g}$ to $e_g$ orbital switching, and is responsible for FM spin fluctuations detrimental to the singlet pairing superconductivity.


title: On the radius of spatial analyticity for defocusing nonlinear Schr\"odinger equations

abstract: In this paper we study spatial analyticity of solutions to the defocusing nonlinear Schr\"odinger equations $iu_t + \Delta u = |u|^{p-1}u$, given initial data which is analytic with fixed radius. It is shown that the uniform radius of spatial analyticity of solutions at later time $t$ cannot decay faster than $1/|t|$ as $|t|\rightarrow\infty$. This extends the previous work of Tesfahun for the cubic case $p=3$ to the cases where $p$ is any odd integer greater than $3$.


title: Zoho at SemEval-2019 Task 9: Semi-supervised Domain Adaptation using Tri-training for Suggestion Mining

abstract: This paper describes our submission for the SemEval-2019 Suggestion Mining task. A simple Convolutional Neural Network (CNN) classifier with contextual word representations from a pre-trained language model was used for sentence classification. The model is trained using tri-training, a semi-supervised bootstrapping mechanism for labelling unseen data. Tri-training proved to be an effective technique to accommodate domain shift for cross-domain suggestion mining (Subtask B) where there is no hand labelled training data. For in-domain evaluation (Subtask A), we use the same technique to augment the training set. Our system ranks thirteenth in Subtask A with an $F_1$-score of 68.07 and third in Subtask B with an $F_1$-score of 81.94.


title: Markovian explorations of random planar maps are roundish

abstract: The infinite discrete stable Boltzmann maps are "heavy-tailed" generalisations of the well-known Uniform Infinite Planar Quadrangulation. Very efficient tools to study these objects are Markovian step-by-step explorations of the lattice called peeling processes. Such a process depends on an algorithm which selects at each step the next edge where the exploration continues. We prove here that, whatever this algorithm, a peeling process always reveals about the same portion of the map, thus growing roughly metric balls. Applied to well-designed algorithms, this easily enables us to compare distances in the map and in its dual, as well as to control the so-called pioneer points of the simple random walk, both on the map and on its dual.


title: Phenomenology of $^{63}$Cu nuclear relaxation in cuprate superconductors

abstract: Nuclear relaxation is an important thermodynamic probe of electronic excitations, in particular in conducting and superconducting systems. Here, an empirical phenomenology based on all available literature data for planar Cu in hole-doped cuprates is developed. It is found that most of the seemingly different relaxation rates among the systems are due to a temperature independent anisotropy that affects the mostly measured $1/T_{1\parallel}$, the rate with an external magnetic field along the crystal $c$-axis, while $1/T_{1\perp}$ is largely independent on doping and material above the critical temperature of superconductivity ($T_c$). This includes very strongly overdoped systems that show Fermi liquid behavior and obey the Korringa law. Below $T_c$ the relaxation rates are similar, as well, if plotted against the reduced temperature $T/T_c$. Thus, planar Cu nuclear relaxation is governed by a simple, dominant mechanism that couples the nuclei with varying anisotropy to a rather ubiquitous bath of electronic excitations that appear Fermi liquid-like irrespective of doping and family. In particular, there is no significant enhancement of the relaxation due to electronic spin fluctuations, different from earlier conclusions. Only the La$_{2-x}$Sr$_x$CuO$_4$ family appears to be an outlier as additional relaxation is present, however, the anisotropy remains temperature independent. Also systems with very low doping levels, for which there is a lack of data, may behave differently.


title: Identifying Galaxy Mergers in Observations and Simulations with Deep Learning

abstract: Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test biases in visually identified merger catalogues. A convolutional neural network architecture was developed and trained in two ways: one with observations from SDSS and one with simulated galaxies from EAGLE, processed to mimic the SDSS observations. The SDSS images were also classified by the simulation trained network and the EAGLE images classified by the observation trained network. The observationally trained network achieves an accuracy of 91.5% while the simulation trained network achieves 65.2% on the visually classified SDSS and physically classified EAGLE images respectively. Classifying the SDSS images with the simulation trained network was less successful, only achieving an accuracy of 64.6%, while classifying the EAGLE images with the observation network was very poor, achieving an accuracy of only 53.0% with preferential assignment to the non-merger classification. This suggests that most of the simulated mergers do not have conspicuous merger features and visually identified merger catalogues from observations are incomplete and biased towards certain merger types. The networks trained and tested with the same data perform the best, with observations performing better than simulations, a result of the observational sample being biased towards conspicuous mergers. Classifying SDSS observations with the simulation trained network has proven to work, providing tantalizing prospects for using simulation trained networks for galaxy identification in large surveys.


title: Mixed cohomology of Lie superalgebras

abstract: We investigate a new cohomology of Lie superalgebras, which may be compared to a de Rham cohomology of Lie supergroups involving both differential and integral forms. It is defined by a BRST complex of Lie superalgebra modules, which is formulated in terms of a Weyl superalgebra and incorporates inequivalent representations of the bosonic Weyl subalgebra. The new cohomology includes the standard Lie superalgebra cohomology as a special case. Examples of new cohomology groups are computed.


title: Proceedings to the 21st workshop "What comes beyond the standard models", Bled June 23--July 1, 2018, with links to the talks appearing in the proceedings, and/or in arXiv

abstract: The contribution contains the preface to the Proceedings to the 21st Workshop "What Comes Beyond the Standard Models", Bled, June 23 - July 1, 2018, published in Bled workshops in physics, Vol.19, No. 2, DMFA-Zalo\v{z}nistvo, Ljubljana, Dec. 2018, links to (most of) the published contributions and section (by M.Yu. Khlopov) on VIA at Bled 2018.


title: Optimal transient growth in an incompressible flow past a backward-slanted step

abstract: With the aim of providing a first step in the quest for a reduction of the aerodynamic drag on the rear-end of a car, we study the phenomena of separation and reattachment of an incompressible flow focusing on a specific aerodynamic geometry, namely a backward-slanted step at 25 degrees of inclination. The ensuing recirculation bubble provides the basis for an analytical and numerical investigation of streamwise-streak generation, lift-up effect, and turbulent-wake and Kelvin-Helmholtz instabilities. A linear stability analysis is performed, and an optimal control problem with a steady volumic forcing is tackled by means of variational formulation, adjoint method, penalization scheme and orthogonalization algorithm. Dealing with the transient growth of spanwise-periodic perturbations and inspired by the need of physically-realizable disturbances, we finally provide a procedure attaining a kinetic-energy maximal gain of the order of one million with respect to the power introduced by the external forcing.


title: Alternating Synthetic and Real Gradients for Neural Language Modeling

abstract: Training recurrent neural networks (RNNs) with backpropagation through time (BPTT) has known drawbacks such as being difficult to capture longterm dependencies in sequences. Successful alternatives to BPTT have not yet been discovered. Recently, BP with synthetic gradients by a decoupled neural interface module has been proposed to replace BPTT for training RNNs. On the other hand, it has been shown that the representations learned with synthetic and real gradients are different though they are functionally identical. In this project, we explore ways of combining synthetic and real gradients with application to neural language modeling tasks. Empirically, we demonstrate the effectiveness of alternating training with synthetic and real gradients after periodic warm restarts on language modeling tasks.


title: Measurement of the 235U(n,f) cross section relative to the 6Li(n,t) and 10B(n,alpha) standards from thermal to 170 keV neutron energy range at n_TOF

abstract: The 235U(n,f) cross section was measured in a wide energy range at n_TOF relative to 6Li(n,t) and 10B(n,alpha), with high resolution and in a wide energy range, with a setup based on a stack of six samples and six silicon detectors placed in the neutron beam. This allowed us to make a direct comparison of the reaction yields under the same experimental conditions, and taking into account the forward/backward emission asymmetry. A hint of an anomaly in the 10{\div}30 keV neutron energy range had been previously observed in other experiments, indicating a cross section systematically lower by several percent relative to major evaluations. The present results indicate that the evaluated cross section in the 9{\div}18 keV neutron energy range is indeed overestimated, both in the recent updates of ENDF/B-VIII.0 and of the IAEA reference data. Furthermore, these new high-resolution data confirm the existence of resonance-like structures in the keV neutron energy region. The new, high accuracy results here reported may lead to a reduction of the uncertainty in the 1{\div}100 keV neutron energy region. Finally, the present data provide additional confidence on the recently re-evaluated cross section integral between 7.8 and 11 eV.


title: Bias implies low rank for quartic polynomials

abstract: We investigate the structure of polynomials of degree four in many variables over a fixed prime field $\mathbb{F}=\mathbb{F}_{p}$. In 2007, Green and Tao proved that if a polynomial $f:\mathbb{F}^{n}\rightarrow\mathbb{F}$ is poorly distributed, then it is a function of a few polynomials of smaller degree. In 2009, Haramaty and Shpilka found an effective bound for $f$ of degree four: If $bias\left(f\right)\geq\delta$, then the number of lower degree polynomials required is at most polynomial in $1/\delta$ and $f$ has a simple presentation as a sum of their products. We make a step towards showing that in fact the number of lower degree polynomials required is at most log-polynomial in $1/\delta$, with the same simple presentation of $f$. This result was a Master's thesis supervised by T. Ziegler at the Hebrew University of Jerusalem, submitted in October 2018. A log-polynomial bound for polynomials of arbitrary degree was recently proved independently by Milicevic and by Janzer.


title: Dimension-independent Sparse Fourier Transform

abstract: The Discrete Fourier Transform (DFT) is a fundamental computational primitive, and the fastest known algorithm for computing the DFT is the FFT (Fast Fourier Transform) algorithm. One remarkable feature of FFT is the fact that its runtime depends only on the size $N$ of the input vector, but not on the dimensionality of the input domain: FFT runs in time $O(N\log N)$ irrespective of whether the DFT in question is on $\mathbb{Z}_N$ or $\mathbb{Z}_n^d$ for some $d>1$, where $N=n^d$. The state of the art for Sparse FFT, i.e. the problem of computing the DFT of a signal that has at most $k$ nonzeros in Fourier domain, is very different: all current techniques for sublinear time computation of Sparse FFT incur an exponential dependence on the dimension $d$ in the runtime. In this paper we give the first algorithm that computes the DFT of a $k$-sparse signal in time $\text{poly}(k, \log N)$ in any dimension $d$, avoiding the curse of dimensionality inherent in all previously known techniques. Our main tool is a new class of filters that we refer to as adaptive aliasing filters: these filters allow isolating frequencies of a $k$-Fourier sparse signal using $O(k)$ samples in time domain and $O(k\log N)$ runtime per frequency, in any dimension $d$. We also investigate natural average case models of the input signal: (1) worst case support in Fourier domain with randomized coefficients and (2) random locations in Fourier domain with worst case coefficients. Our techniques lead to an $\widetilde O(k^2)$ time algorithm for the former and an $\widetilde O(k)$ time algorithm for the latter.


title: Synthesis of the Morphological Description of Cometary Dust at Comet 67P

abstract: Before Rosetta, the space missions Giotto and Stardust shaped our view on cometary dust, supported by plentiful data from Earth based observations and interplanetary dust particles collected in the Earth's atmosphere. The Rosetta mission at comet 67P/Churyumov-Gerasimenko was equipped with a multitude of instruments designed to study cometary dust. While an abundant amount of data was presented in several individual papers, many focused on a dedicated measurement or topic. Different instruments, methods, and data sources provide different measurement parameters and potentially introduce different biases. This can be an advantage if the complementary aspect of such a complex data set can be exploited. However, it also poses a challenge in the comparison of results in the first place. The aim of this work therefore is to summarise dust results from Rosetta and before. We establish a simple classification as a common framework for inter-comparison. This classification is based on a dust particle's structure, porosity, and strength as well as its size. Depending on the instrumentation, these are not direct measurement parameters but we chose them as they were the most reliable to derive our model. The proposed classification already proved helpful in the Rosetta dust community and we propose to take it into consideration also beyond. In this manner we hope to better identify synergies between different instruments and methods in the future.


title: Discovery of a nitrogen-enhanced mildly metal-poor binary system: Possible evidence for pollution from an extinct AGB Star

abstract: We report the serendipitous discovery of a nitrogen-rich, mildly metal-poor ([Fe/H]=-1.08) giant star in a single-lined spectroscopic binary system found in the SDSS-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) survey, Data Release 14 (DR14). Previous work has assumed that the two percent of halo giants with unusual elemental abundances have been evaporated from globular clusters, but other origins for their abundance signatures, including binary mass transfer, must also be explored. We present the results of an abundance re-analysis of the APOGEE-2 high-resolution near-infrared spectrum of 2M12451043+1217401 with the Brussels Automatic Stellar Parameter (BACCHUS) automated spectral analysis code, and re-derive manually the main element families, namely the light elements (C, N), elements (O, Mg, Si), iron-peak element (Fe), \textit{s}-process element (Ce), and the light odd-Z element (Al). Our analysis confirm the N-rich nature of 2M12451043+1217401, which has a [N/Fe] ratio of $+0.69$, and shows that the abundances of C and Al are slightly discrepant from that of a typical mildly metal-poor RGB star, but exhibit Mg, Si, O and \textit{s}-process abundances (Ce) of typical field stars. We also detect a particularly large variability in its radial velocity over the period of the APOGEE-2 observations, and the most likely orbit fit to the radial velocity data has a period of $730.89\pm106.86$ days, a velocity semi-amplitude of $9.92 \pm 0.14$ km s$^{-1}$, and an eccentricity of $\sim 0.1276 \pm0.1174$, which support the hypothesis of a binary companion, and that has probably been polluted by a now-extinct AGB star.


title: Alleviating the $H_0$ and $\sigma_8$ anomalies with a decaying dark matter model

abstract: The Hubble tension between the $\Lambda$CDM-model-dependent prediction of the current expansion rate $H_0$ using Planck data and direct, model-independent measurements in the local universe from the SH0ES collaboration disagree at $>3.5\sigma$. Moreover, there exists a milder $\sim 2\sigma$ tension between similar predictions for the amplitude $S_8$ of matter fluctuations and its measurement in the local universe. As explanations relying on unresolved systematics have not been found, theorists have been exploring explanations for these anomalies that modify the cosmological model, altering early-universe-based predictions for these parameters. However, new cosmological models that attempt to resolve one tension often worsen the other. In this paper, we investigate a decaying dark matter (DDM) model as a solution to both tensions simultaneously. Here, a fraction of dark matter density decays into dark radiation. The decay rate $\Gamma$ is proportional to the Hubble rate $H$ through the constant $\alpha_{\rm dr}$, the only additional parameter of this model. Then, this model deviates most from $\Lambda$CDM in the early universe, with $\alpha_{\rm dr}$ being positively correlated with $H_0$ and negatively with $S_8$. Hence, increasing $\alpha_{\rm dr}$ (and allowing dark matter to decay in this way) can then diminish both tensions simultaneously. When only considering Planck CMB data and the local SH0ES prior on $H_0$, $\sim 1$\% dark matter decays, decreasing the $S_8$ tension to $0.3\sigma$ and increasing the best-fit $H_0$ by $1.6$ km/s/Mpc. However, the addition of intermediate-redshift data (the JLA supernova dataset and baryon acoustic oscillation data) weakens the effectiveness of this model. Only $\sim 0.5$\% of the dark matter decays bringing the $S_8$ tension back up to $\sim 1.5 \sigma$ and the increase in the best-fit $H_0$ down to $0.4$ km/s/Mpc.


title: Space-time fractional stochastic partial differential equations with L\'evy Noise

abstract: We consider non-linear time-fractional stochastic heat type equation $$\frac{\partial^\beta u}{\partial t^\beta}+\nu(-\Delta)^{\alpha/2} u=I^{1-\beta}_t \bigg[\int_{\mathbb{R}^d}\sigma(u(t,x),h) \stackrel{\cdot}{\tilde N }(t,x,h)\bigg]$$ and $$\frac{\partial^\beta u}{\partial t^\beta}+\nu(-\Delta)^{\alpha/2} u=I^{1-\beta}_t \bigg[\int_{\mathbb{R}^d}\sigma(u(t,x),h) \stackrel{\cdot}{N }(t,x,h)\bigg]$$ in $(d+1)$ dimensions, where $\alpha\in (0,2]$ and $d<\min\{2,\beta^{-1}\}\alpha$, $\nu>0$, $\partial^\beta_t$ is the Caputo fractional derivative, $-(-\Delta)^{\alpha/2} $ is the generator of an isotropic stable process, $I^{1-\beta}_t$ is the fractional integral operator, ${N}(t,x)$ are Poisson random measure with $\tilde{N}(t,x)$ being the compensated Poisson random measure. $\sigma:{\mathbb{R}}\to{\mathbb{R}}$ is a Lipschitz continuous function. We prove existence and uniqueness of mild solutions to this equation. Our results extend the results in the case of parabolic stochastic partial differential equations obtained in "M. Foondun and D. Khoshnevisan. Intermittence and nonlinear parabolic stochastic partial differential equations. \emph{ Electron. J. Probab.} {\bf14} (2009), 548--568" and " J. B. Walsh. An Introduction to Stochastic Partial Differential Equations, \'Ecoled'\'et\'e de Probabilit\'es de Saint-Flour, XIV|1984, Lecture Notes in Math., vol. 1180, Springer, Berlin, (1986), 265--439". Under the linear growth of $\sigma$, we show that the solution of the time fractional stochastic partial differential equation follows an exponential growth with respect to the time. We also show the nonexistence of the random field solution of both stochastic partial differential equations when $\sigma$ grows faster than linear.


title: Capturing sets of ordinals by normal ultrapowers

abstract: We investigate the extent to which ultrapowers by normal measures on $\kappa$ can be correct about powersets $\mathcal{P}(\lambda)$ for $\lambda>\kappa$. We consider two versions of this questions, the capturing property $\mathrm{CP}(\kappa,\lambda)$ and the local capturing property $\mathrm{LCP}(\kappa,\lambda)$. $\mathrm{CP}(\kappa,\lambda)$ holds if there is an ultrapower by a normal measure on $\kappa$ which correctly computes $\mathcal{P}(\lambda)$. $\mathrm{LCP}(\kappa,\lambda)$ is a weakening of $\mathrm{CP}(\kappa,\lambda)$ which holds if every subset of $\lambda$ is contained in some ultrapower by a normal measure on $\kappa$. After examining the basic properties of these two notions, we identify the exact consistency strength of $\mathrm{LCP}(\kappa,\kappa^+)$. Building on results of Cummings, who determined the exact consistency strength of $\mathrm{CP}(\kappa,\kappa^+)$, and using a forcing due to Apter and Shelah, we show that $\mathrm{CP}(\kappa,\lambda)$ can hold at the least measurable cardinal.


title: Forward model for quantitative pulse-echo speed-of-sound imaging

abstract: Computed ultrasound tomography in echo mode (CUTE) allows determining the spatial distribution of speed-of-sound (SoS) inside tissue using handheld pulse-echo ultrasound (US). This technique is based on measuring the changing phase of beamformed echoes obtained under varying transmit (Tx) and/or receive (Rx) steering angles. The SoS is reconstructed by inverting a forward model describing how the spatial distribution of SoS is related to the spatial distribution of the echo phase shift. CUTE holds promise as a novel diagnostic modality that complements conventional US in a single, real-time handheld system. Here we demonstrate that, in order to obtain robust quantitative results, the forward model must contain two features that were not taken into account so far: a) the phase shift must be detected between pairs of Tx and Rx angles that are centred around a set of common mid-angles, and b) it must account for an additional phase shift induced by the error of the reconstructed position of echoes. In a phantom study mimicking liver imaging, this new model leads to a substantially improved quantitative SoS reconstruction compared to the model that has been used so far. The importance of the new model as a prerequisite for an accurate diagnosis is corroborated in preliminary volunteer results.


title: Efficient Video Classification Using Fewer Frames

abstract: Recently,there has been a lot of interest in building compact models for video classification which have a small memory footprint (<1 GB). While these models are compact, they typically operate by repeated application of a small weight matrix to all the frames in a video. E.g. recurrent neural network based methods compute a hidden state for every frame of the video using a recurrent weight matrix. Similarly, cluster-and-aggregate based methods such as NetVLAD, have a learnable clustering matrix which is used to assign soft-clusters to every frame in the video. Since these models look at every frame in the video, the number of floating point operations (FLOPs) is still large even though the memory footprint is small. We focus on building compute-efficient video classification models which process fewer frames and hence have less number of FLOPs. Similar to memory efficient models, we use the idea of distillation albeit in a different setting. Specifically, in our case, a compute-heavy teacher which looks at all the frames in the video is used to train a compute-efficient student which looks at only a small fraction of frames in the video. This is in contrast to a typical memory efficient Teacher-Student setting, wherein both the teacher and the student look at all the frames in the video but the student has fewer parameters. Our work thus complements the research on memory efficient video classification. We do an extensive evaluation with three types of models for video classification,viz.(i) recurrent models (ii) cluster-and-aggregate models and (iii) memory-efficient cluster-and-aggregate models and show that in each of these cases, a see-it-all teacher can be used to train a compute efficient see-very-little student. We show that the proposed student network can reduce the inference time by 30% and the number of FLOPs by approximately 90% with a negligible drop in the performance.


title: All minimal Cantor systems are slow

abstract: We show that every (invertible, or noninvertible) minimal Cantor system embeds in $\mathbb{R}$ with vanishing derivative everywhere. We also study relations between local shrinking and periodic points.


title: Order of contact and ruled submanifolds

abstract: We prove a generalization of the Monge-Cayley-Salmon theorem on osculation and ruled submanifolds using elementary geometric measure theory.


title: Some facts about the life and the scientific work of the Belgian mathematician Paul Mansion (1844-1919) Centenary Paul Mansion (working paper, 2019)

abstract: In this article, we are interested in the life and scientific work of the Belgian mathematician Paul Mansion. The year 2019 marks the centenary of his passing. We bring some new insights into Paul Mansion's work thanks to his scientific correspondence at the Royal Library in Brussels.


title: Provable Guarantees for Gradient-Based Meta-Learning

abstract: We study the problem of meta-learning through the lens of online convex optimization, developing a meta-algorithm bridging the gap between popular gradient-based meta-learning and classical regularization-based multi-task transfer methods. Our method is the first to simultaneously satisfy good sample efficiency guarantees in the convex setting, with generalization bounds that improve with task-similarity, while also being computationally scalable to modern deep learning architectures and the many-task setting. Despite its simplicity, the algorithm matches, up to a constant factor, a lower bound on the performance of any such parameter-transfer method under natural task similarity assumptions. We use experiments in both convex and deep learning settings to verify and demonstrate the applicability of our theory.


title: S-PRAC: Fast Partial Packet Recovery with Network Coding in Very Noisy Wireless Channels

abstract: Well-known error detection and correction solutions in wireless communications are slow or incur high transmission overhead. Recently, notable solutions like PRAC and DAPRAC, implementing partial packet recovery with network coding, could address these problems. However, they perform slowly when there are many errors. We propose S-PRAC, a fast scheme for partial packet recovery, particularly designed for very noisy wireless channels. S-PRAC improves on DAPRAC. It divides each packet into segments consisting of a fixed number of small RLNC encoded symbols and then attaches a CRC code to each segment and one to each coded packet. Extensive simulations show that S-PRAC can detect and correct errors quickly. It also outperforms DAPRAC significantly when the number of errors is high.


title: Unifying Ensemble Methods for Q-learning via Social Choice Theory

abstract: Ensemble methods have been widely applied in Reinforcement Learning (RL) in order to enhance stability, increase convergence speed, and improve exploration. These methods typically work by employing an aggregation mechanism over actions of different RL algorithms. We show that a variety of these methods can be unified by drawing parallels from committee voting rules in Social Choice Theory. We map the problem of designing an action aggregation mechanism in an ensemble method to a voting problem which, under different voting rules, yield popular ensemble-based RL algorithms like Majority Voting Q-learning or Bootstrapped Q-learning. Our unification framework, in turn, allows us to design new ensemble-RL algorithms with better performance. For instance, we map two diversity-centered committee voting rules, namely Single Non-Transferable Voting Rule and Chamberlin-Courant Rule, into new RL algorithms that demonstrate excellent exploratory behavior in our experiments.


title: Deep Learning-based Concept Detection in vitrivr at the Video Browser Showdown 2019 - Final Notes

abstract: This paper presents an after-the-fact summary of the participation of the vitrivr system to the 2019 Video Browser Showdown. Analogously to last year's report, the focus of this paper lies on additions made since the original publication and the system's performance during the competition.


title: Probabilistic Parity Shaping for Linear Codes

abstract: Linear layered probabilistic shaping (LLPS) is proposed, an architecture for linear codes to efficiently encode to shaped code words. In the previously proposed probabilistic amplitude shaping (PAS) architecture, a distribution matcher (DM) maps information bits to shaped bits, which are then systematically encoded by appending uniformly distributed parity bits. LLPS extends PAS by probabilistic parity shaping (PPS), which uses a syndrome DM to calculate shaped parity bits. LLPS enables the transmission with any desired distribution using linear codes, furthermore, by LLPS, a given linear code with rate $R_\text{fec}$ can be operated at any rate $R\leq R_\text{fec}$ by changing the distribution. LLPS is used with an LDPC code for dirty paper coding against an interfering BPSK signal, improving the energy efficiency by 0.8 dB.


title: F10-SGD: Fast Training of Elastic-net Linear Models for Text Classification and Named-entity Recognition

abstract: Voice-assistants text classification and named-entity recognition (NER) models are trained on millions of example utterances. Because of the large datasets, long training time is one of the bottlenecks for releasing improved models. In this work, we develop F10-SGD, a fast optimizer for text classification and NER elastic-net linear models. On internal datasets, F10-SGD provides 4x reduction in training time compared to the OWL-QN optimizer without loss of accuracy or increase in model size. Furthermore, we incorporate biased sampling that prioritizes harder examples towards the end of the training. As a result, in addition to faster training, we were able to obtain statistically significant accuracy improvements for NER. On public datasets, F10-SGD obtains 22% faster training time compared to FastText for text classification. And, 4x reduction in training time compared to CRFSuite OWL-QN for NER.


title: Quantum Anomalous Hall Effect through Canted Antiferromagnetism

abstract: Most canted antiferromagnets are also anomalous Hall antiferromagnets (AHE AFMs), i.e. they have an anomalous Hall response and other responses with the same symmetry requirements. We suggest that AHE AFMs are promising materials as hosts for high-temperature quantum anomalous Hall effects. By considering models of two-dimensional (001) perovskite layers with strong spin-orbit coupling that isolates an effective total angular momentum $\tilde{j}= \frac{1}{2}$ subspace of the $t_{2g}$ manifold, we propose a strategy to engineering quantum anomalous Hall antiferromagnets.


title: Lorentzian Geodesic Flows and Interpolation between Hypersurfaces in Euclidean Spaces

abstract: We consider geodesic flows between hypersurfaces in $\R^n$. However, rather than consider using geodesics in $\R^n$, which are straight lines, we consider an induced flow using geodesics between the tangent spaces of the hypersurfaces viewed as affine hyperplanes. For naturality, we want the geodesic flow to be invariant under rigid transformations and homotheties. Consequently, we do not use the dual projective space, as the geodesic flow in this space is not preserved under translations. Instead we give an alternate approach using a Lorentzian space, which is semi-Riemannian with a metric of index $1$. For this space for points corresponding to affine hyperplanes in $\R^n$, we give a formula for the geodesic between two such points. As a consequence, we show the geodesic flow is preserved by rigid transformations and homotheties of $\R^n$. Furthermore, we give a criterion that a vector field in a smoothly varying family of hyperplanes along a curve yields a Lorentzian parallel vector field for the corresponding curve in the Lorentzian space. As a result this provides a method to extend an orthogonal frame in one affine hyperplane to a smoothly "Lorentzian varying" family of orthogonal frames in a family of affine hyperplanes along a smooth curve, as well as a interpolating between two such frames with a smooth " minimally Lorentzian varying" family of orthogonal frames. We further give sufficient conditions that the Lorentzian flow from a hypersurface is nonsingular and that the resulting corresponding flow in $\R^n$ is nonsingular. This is illustrated for surfaces in $\R^3$.


title: Learning nonlinear level sets for dimensionality reduction in function approximation

abstract: We developed a Nonlinear Level-set Learning (NLL) method for dimensionality reduction in high-dimensional function approximation with small data. This work is motivated by a variety of design tasks in real-world engineering applications, where practitioners would replace their computationally intensive physical models (e.g., high-resolution fluid simulators) with fast-to-evaluate predictive machine learning models, so as to accelerate the engineering design processes. There are two major challenges in constructing such predictive models: (a) high-dimensional inputs (e.g., many independent design parameters) and (b) small training data, generated by running extremely time-consuming simulations. Thus, reducing the input dimension is critical to alleviate the over-fitting issue caused by data insufficiency. Existing methods, including sliced inverse regression and active subspace approaches, reduce the input dimension by learning a linear coordinate transformation; our main contribution is to extend the transformation approach to a nonlinear regime. Specifically, we exploit reversible networks (RevNets) to learn nonlinear level sets of a high-dimensional function and parameterize its level sets in low-dimensional spaces. A new loss function was designed to utilize samples of the target functions' gradient to encourage the transformed function to be sensitive to only a few transformed coordinates. The NLL approach is demonstrated by applying it to three 2D functions and two 20D functions for showing the improved approximation accuracy with the use of nonlinear transformation, as well as to an 8D composite material design problem for optimizing the buckling-resistance performance of composite shells of rocket inter-stages.


title: Possible Realization and Protection of Valley-Polarized Quantum Hall Effect in Mn/WS2

abstract: By using the first-principles calculations and model analyses, we found that the combination of defected tungsten disulfide monolayer and sparse manganese adsorption may give a KK` valley spin splitting up to 210 meV. This system also has a tunable magnetic anisotropy energy, a clean band gap, and an appropriate band alignment, with the Fermi level sitting right above the top of valence bands at the K-valleys. Therefore, it can be used for the realization of the valley-polarized anomalous Hall effect and for the exploration of other valley related physics without using optical methods. A protective environment can be formed by covering it with a hexagonal BN layer, without much disturbance to the benign properties of Mn/WS2.


title: Register Games

abstract: The complexity of parity games is a long standing open problem that saw a major breakthrough in 2017 when two quasi-polynomial algorithms were published. This article presents a third, independent approach to solving parity games in quasi-polynomial time, based on the notion of register game, a parameterised variant of a parity game. The analysis of register games leads to a quasi-polynomial algorithm for parity games, a polynomial algorithm for restricted classes of parity games and a novel measure of complexity, the register index, which aims to capture the combined complexity of the priority assignement and the underlying game graph. We further present a translation of alternating parity word automata into alternating weak automata with only a quasi-polynomial increase in size, based on register games; this improves on the previous exponential translation. We also use register games to investigate the parity index hierarchy: while for words the index hierarchy of alternating parity automata collapses to the weak level, and for trees it is strict, for structures between trees and words, it collapses logarithmically, in the sense that any parity tree automaton of size n is equivalent, on these particular classes of structures, to an automaton with a number of priorities logarithmic in n.


title: Linear Time Visualization and Search in Big Data using Pixellated Factor Space Mapping

abstract: It is demonstrated how linear computational time and storage efficient approaches can be adopted when analyzing very large data sets. More importantly, interpretation is aided and furthermore, basic processing is easily supported. Such basic processing can be the use of supplementary, i.e. contextual, elements, or particular associations. Furthermore pixellated grid cell contents can be utilized as a basic form of imposed clustering. For a given resolution level, here related to an associated m-adic ($m$ here is a non-prime integer) or p-adic ($p$ is prime) number system encoding, such pixellated mapping results in partitioning. The association of a range of m-adic and p-adic representations leads naturally to an imposed hierarchical clustering, with partition levels corresponding to the m-adic-based and p-adic-based representations and displays. In these clustering embedding and imposed cluster structures, some analytical visualization and search applications are described


title: The loop-stable homotopy category of algebras

abstract: Let l be a commutative ring with unit. Garkusha constructed a functor from the category of l-algebras into a triangulated category D, that is a universal excisive and homotopy invariant homology theory. Later on, he provided different descriptions of D, as an application of his motivic homotopy theory of algebras. Using these, it can be shown that D is triangulated equivalent to a category, denote it by k, whose objects are pairs (A,m) with A an l-algebra and m an integer, and whose Hom-sets can be described in terms of homotopy classes of morphisms. All these computations, however, require a heavy machinery of homotopy theory. In this paper, we give a more explicit construction of the triangulated category k and prove its universal property, avoiding the homotopy-theoretic methods and using instead the ones developed by Corti\~nas-Thom for defining kk-theory. Moreover, we give a new description of the composition law in k, mimicking the one in the suspension-stable homotopy category of bornological algebras defined by Cuntz-Meyer-Rosenberg. We also prove that the triangulated structure in k can be defined using either extension or mapping path triangles.


title: From explanation to synthesis: Compositional program induction for learning from demonstration

abstract: Hybrid systems are a compact and natural mechanism with which to address problems in robotics. This work introduces an approach to learning hybrid systems from demonstrations, with an emphasis on extracting models that are explicitly verifiable and easily interpreted by robot operators. We fit a sequence of controllers using sequential importance sampling under a generative switching proportional controller task model. Here, we parameterise controllers using a proportional gain and a visually verifiable joint angle goal. Inference under this model is challenging, but we address this by introducing an attribution prior extracted from a neural end-to-end visuomotor control model. Given the sequence of controllers comprising a task, we simplify the trace using grammar parsing strategies, taking advantage of the sequence compositionality, before grounding the controllers by training perception networks to predict goals given images. Using this approach, we are successfully able to induce a program for a visuomotor reaching task involving loops and conditionals from a single demonstration and a neural end-to-end model. In addition, we are able to discover the program used for a tower building task. We argue that computer program-like control systems are more interpretable than alternative end-to-end learning approaches, and that hybrid systems inherently allow for better generalisation across task configurations.


title: Regularity Normalization: Neuroscience-Inspired Unsupervised Attention across Neural Network Layers

abstract: Inspired by the adaptation phenomenon of neuronal firing, we propose the regularity normalization (RN) as an unsupervised attention mechanism (UAM) which computes the statistical regularity in the implicit space of neural networks under the Minimum Description Length (MDL) principle. Treating the neural network optimization process as a partially observable model selection problem, the regularity normalization constrains the implicit space by a normalization factor, the universal code length. We compute this universal code incrementally across neural network layers and demonstrate the flexibility to include data priors such as top-down attention and other oracle information. Empirically, our approach outperforms existing normalization methods in tackling limited, imbalanced and non-stationary input distribution in image classification, classic control, procedurally-generated reinforcement learning, generative modeling, handwriting generation and question answering tasks with various neural network architectures. Lastly, the unsupervised attention mechanisms is a useful probing tool for neural networks by tracking the dependency and critical learning stages across layers and recurrent time steps of deep networks.

categories: cs.LG cs.CV cs.IT math.IT q-bio.NC stat.ML


title: Online adaptive basis refinement and compression for reduced-order models via vector-space sieving

abstract: In many applications, projection-based reduced-order models (ROMs) have demonstrated the ability to provide rapid approximate solutions to high-fidelity full-order models (FOMs). However, there is no a priori assurance that these approximate solutions are accurate; their accuracy depends on the ability of the low-dimensional trial basis to represent the FOM solution. As a result, ROMs can generate inaccurate approximate solutions, e.g., when the FOM solution at the online prediction point is not well represented by training data used to construct the trial basis. To address this fundamental deficiency of standard model-reduction approaches, this work proposes a novel online-adaptive mechanism for efficiently enriching the trial basis in a manner that ensures convergence of the ROM to the FOM, yet does not incur any FOM solves. The mechanism is based on the previously proposed adaptive $h$-refinement method for ROMs [12], but improves upon this work in two crucial ways. First, the proposed method enables basis refinement with respect to any orthogonal basis (not just the Kronecker basis), thereby generalizing the refinement mechanism and enabling it to be tailored to the physics characterizing the problem at hand. Second, the proposed method provides a fast online algorithm for periodically compressing the enriched basis via an efficient proper orthogonal decomposition (POD) method, which does not incur any operations that scale with the FOM dimension. These two features allow the proposed method to serve as (1) a failsafe mechanism for ROMs, as the method enables the ROM to satisfy any prescribed error tolerance online (even in the case of inadequate training), and (2) an efficient online basis-adaptation mechanism, as the combination of basis enrichment and compression enables the basis to adapt online while controlling its dimension.


title: Robust Decision Trees Against Adversarial Examples

abstract: Although adversarial examples and model robustness have been extensively studied in the context of linear models and neural networks, research on this issue in tree-based models and how to make tree-based models robust against adversarial examples is still limited. In this paper, we show that tree based models are also vulnerable to adversarial examples and develop a novel algorithm to learn robust trees. At its core, our method aims to optimize the performance under the worst-case perturbation of input features, which leads to a max-min saddle point problem. Incorporating this saddle point objective into the decision tree building procedure is non-trivial due to the discrete nature of trees --- a naive approach to finding the best split according to this saddle point objective will take exponential time. To make our approach practical and scalable, we propose efficient tree building algorithms by approximating the inner minimizer in this saddle point problem, and present efficient implementations for classical information gain based trees as well as state-of-the-art tree boosting models such as XGBoost. Experimental results on real world datasets demonstrate that the proposed algorithms can substantially improve the robustness of tree-based models against adversarial examples.


title: Maximum Wiener index of unicyclic graphs with given bipartition

abstract: The \emph{Wiener index} is a widely studied topological index of graphs. One of the main problems in the area is to determine which graphs of given properties attain the extremal values of Wiener index. In this paper we resolve an open problem posed by Du in 2012. To this end we determine the unicyclic bipartite graphs with given size of parts having the maximum Wiener index. This completes the previous research in which the minimum case was solved.


title: Non-Uniform Robot Densities in Vibration Driven Swarms Using Phase Separation Theory

abstract: In robot swarms operating under highly restrictive sensing and communication constraints, individuals may need to use direct physical proximity to facilitate information exchange. However, in certain task-related scenarios, this requirement might conflict with the need for robots to spread out in the environment, e.g., for distributed sensing or surveillance applications. This paper demonstrates how a swarm of minimally-equipped robots can form high-density robot aggregates which coexist with lower robot densities in the domain. We envision a scenario where a swarm of vibration-driven robots---which sit atop bristles and achieve directed motion by vibrating them---move somewhat randomly in an environment while colliding with each other. Theoretical techniques from the study of far-from-equilibrium collectives and statistical mechanics clarify the mechanisms underlying the formation of these high and low density regions. Specifically, we capitalize on a transformation that connects the collective properties of a system of self-propelled particles with that of a well-studied molecular fluid system, thereby inheriting the rich theory of equilibrium thermodynamics. This connection is a formal one and is a relatively recent result in studies of motility induced phase separation; it is previously unexplored in the context of robotics. Real robot experiments as well as simulations illustrate how inter-robot collisions can precipitate the formation of non-uniform robot densities in a closed and bounded region.


title: Semi-inclusive Deep-Inelastic Scattering, Parton Distributions and Fragmentation Functions at a Future Electron-Ion Collider

abstract: We present a quantitative assessment of the impact a future Electron-Ion Collider would have in the determination of parton distribution functions in the proton and parton-to-hadron fragmentation functions through semi-inclusive deep-inelastic electron-proton scattering data. Specifically, we estimate the kinematic configurations for which the forthcoming data are expected to have the most significant impact in the precision of these distributions, computing the respective correlation and sensitivity coefficients. Implementing the reweighting technique over sets of simulated data with realistic uncertainties for two different center-of-mass energies, we analyse the resulting upgraded sets of PDFs and FFs, which have significantly reduced uncertainties.


title: Local Function Complexity for Active Learning via Mixture of Gaussian Processes

abstract: Inhomogeneities in real-world data, e.g., due to changes in the observation noise level or variations in the structural complexity of the source function, pose a unique set of challenges for statistical inference. Accounting for them can greatly improve predictive power when physical resources or computation time is limited. In this paper, we draw on recent theoretical results on the estimation of local function complexity (LFC), derived from the domain of local polynomial smoothing (LPS), to establish a notion of local structural complexity, which is used to develop a model-agnostic active learning (AL) framework. Due to its reliance on pointwise estimates, the LPS model class is not robust and scalable concerning large input space dimensions that typically come along with real-world problems. Here, we derive and estimate the Gaussian process regression (GPR)-based analog of the LPS-based LFC and use it as a substitute in the above framework to make it robust and scalable. We assess the effectiveness of our LFC estimate in an AL application on a prototypical low-dimensional synthetic dataset, before taking on the challenging real-world task of reconstructing a quantum chemical force field for a small organic molecule and demonstrating state-of-the-art performance with a significantly reduced training demand.


title: Quartic graphs which are Bakry-\'Emery curvature sharp

abstract: We give a classification of all connected quartic graphs which are (infinity) curvature sharp in all vertices with respect to Bakry-\'Emery curvature. The result is based on a computer classification by F. Gurr and L. Watson May and a combinatorial case by case investigation.


title: Improving Missing Data Imputation with Deep Generative Models

abstract: Datasets with missing values are very common on industry applications, and they can have a negative impact on machine learning models. Recent studies introduced solutions to the problem of imputing missing values based on deep generative models. Previous experiments with Generative Adversarial Networks and Variational Autoencoders showed interesting results in this domain, but it is not clear which method is preferable for different use cases. The goal of this work is twofold: we present a comparison between missing data imputation solutions based on deep generative models, and we propose improvements over those methodologies. We run our experiments using known real life datasets with different characteristics, removing values at random and reconstructing them with several imputation techniques. Our results show that the presence or absence of categorical variables can alter the selection of the best model, and that some models are more stable than others after similar runs with different random number generator seeds.


title: Bridging the Gap: Attending to Discontinuity in Identification of Multiword Expressions

abstract: We introduce a new method to tag Multiword Expressions (MWEs) using a linguistically interpretable language-independent deep learning architecture. We specifically target discontinuity, an under-explored aspect that poses a significant challenge to computational treatment of MWEs. Two neural architectures are explored: Graph Convolutional Network (GCN) and multi-head self-attention. GCN leverages dependency parse information, and self-attention attends to long-range relations. We finally propose a combined model that integrates complementary information from both through a gating mechanism. The experiments on a standard multilingual dataset for verbal MWEs show that our model outperforms the baselines not only in the case of discontinuous MWEs but also in overall F-score.


title: Wide band spectroscopic response of monocrystallines to low dose neutron and gamma radiation

abstract: We identify a number of crystalline structures with promising characteristics to serve as a detection medium for a novel Dark Matter (DM) detector with a low threshold energy. A detector of this kind can be specifically useful in application requiring the detection of nuclear recoils, such as in direct detection of low mass DM, coherent neutrino scattering and neutrons. We describe a broad band, high sensitivity optical setup designed and constructed for the purpose of this search and future investigations of specific crystals. We report on the fluorescent signals produced from exposure to low doses of neutrons and $\gamma$ rays and find potential targets in Quartz, Sapphire, LiF, CaF$_{2}$ and BaF$_{2}$. These crystals and specific signals will be the subject of further study to establish the various traits relevant for a full scale DM detector. In this paper we identify the most interesting signals that will be promoted to significantly more detailed studies, including their production mechanism.


title: Localization of a Model With U(1) Kinetic Gauge Mixing

abstract: In this paper we study the localization of a model with kinetic gauge mixing on a thin membrane. The model we discuss is a theory for millicharged particles, proposed by Holdom in 1985, and that now is object of study in the LHC. We propose a geometrical coupling between the gauge fields, the Ricci scalar and the Ricci tensor. We show that it is possible to localize such a model by regarding specific values for the coupling constants. We find the solutions for the two gauge fields and discuss the localization for scalar fields that appears naturally in the process. We show that not necessarily the gauge and scalar fields are localized at the same time.


title: From radio-quiet to radio-silent: low luminosity Seyfert radio cores

abstract: A strong effort has been devoted to understand the physical origin of radio emission from low-luminosity AGN (LLAGN), but a comprehensive picture is still missing. We used high-resolution ($\le$1 arcsec), multi-frequency (1.5, 5.5, 9 and 14 GHz) NSF's Karl G. Jansky Very Large Array (VLA) observations to characterise the state of the nuclear region of ten Seyfert nuclei, which are the faintest members of a complete, distance-limited sample of 28 sources. With the sensitivity and resolution guaranteed by the VLA-A configuration, we measured radio emission for six sources (NGC3185, NGC3941, NGC4477, NGC4639, NGC4698 and NGC4725), while for the remaining four (NGC0676, NGC1058, NGC2685 and NGC3486) we put upper limits at tens uJy/beam level, below the previous 0.12 mJy/beam level of Ho&Ulvestad (2001), corresponding to luminosities down to L$\le10^{19}$ W/Hz at 1.5 GHz for the highest RMS observation. Two sources, NGC4639 and NGC4698, exhibit spectral slopes compatible with inverted spectra ($\alpha\le$0, $S_{\nu}\,\propto\,{\nu}^{-\alpha}$), hint for radio emission from an optically-thick core, while NGC4477 exhibits a steep (+0.52$\pm$0.09) slope. The detected sources are mainly compact on scales $\le$ arcseconds, predominantly unresolved, except NGC3185 and NGC3941, in which the resolved radio emission could be associated to star-formation processes. A significant X-ray - radio luminosities correlation is extended down to very low luminosities, with slope consistent with inefficient accretion, expected at such low Eddington ratios. Such sources will be one of the dominant Square Kilometre Array (SKA) population, allowing a deeper understanding of the physics underlying such faint AGN.


title: Customizing Object Detectors for Indoor Robots

abstract: Object detection models based on convolutional neural networks (CNNs) demonstrate impressive performance when trained on large-scale labeled datasets. While a generic object detector trained on such a dataset performs adequately in applications where the input data is similar to user photographs, the detector performs poorly on small objects, particularly ones with limited training data or imaged from uncommon viewpoints. Also, a specific room will have many objects that are missed by standard object detectors, frustrating a robot that continually operates in the same indoor environment. This paper describes a system for rapidly creating customized object detectors. Data is collected from a quadcopter that is teleoperated with an interactive interface. Once an object is selected, the quadcopter autonomously photographs the object from multiple viewpoints to %create training data that is used by DUNet (Dense Upscaled Net), collect data to train DUNet (Dense Upscaled Network), our proposed model for learning customized object detectors from scratch given limited data. Our experiments compare the performance of learning models from scratch with DUNet vs.\ fine tuning existing state of the art object detectors, both on our indoor robotics domain and on standard datasets.


title: On Carmichael and polygonal numbers, Bernoulli polynomials, and sums of base-$p$ digits

abstract: We give a new characterization of the set $\mathcal{C}$ of Carmichael numbers in the context of $p$-adic theory, independently of the classical results of Korselt and Carmichael. The characterization originates from a surprising link to the denominators of the Bernoulli polynomials via the sum-of-base-$p$-digits function. More precisely, we show that such a denominator obeys a triple-product identity, where one factor is connected with a $p$-adically defined subset $\mathcal{S}$ of the squarefree integers that contains $\mathcal{C}$. This leads to the definition of a new subset $\mathcal{C}'$ of $\mathcal{C}$, called the "primary Carmichael numbers". Subsequently, we establish that every Carmichael number equals an explicitly determined polygonal number. Finally, the set $\mathcal{S}$ is covered by modular subsets $\mathcal{S}_d$ ($d \geq 1$) that are related to the Kn\"odel numbers, where $\mathcal{C} = \mathcal{S}_1$ is a special case.


title: Improved Fault-Tolerant Quantum Simulation of Condensed-Phase Correlated Electrons via Trotterization

abstract: Recent work has deployed linear combinations of unitaries techniques to reduce the cost of fault-tolerant quantum simulations of correlated electron models. Here, we show that one can sometimes improve upon those results with optimized implementations of Trotter-Suzuki-based product formulas. We show that low-order Trotter methods perform surprisingly well when used with phase estimation to compute relative precision quantities (e.g. energies per unit cell), as is often the goal for condensed-phase systems. In this context, simulations of the Hubbard and plane-wave electronic structure models with $N < 10^5$ fermionic modes can be performed with roughly $O(1)$ and $O(N^2)$ T complexities. We perform numerics revealing tradeoffs between the error and gate complexity of a Trotter step; e.g., we show that split-operator techniques have less Trotter error than popular alternatives. By compiling to surface code fault-tolerant gates and assuming error rates of one part per thousand, we show that one can error-correct quantum simulations of interesting, classically intractable instances with a few hundred thousand physical qubits.


title: The Best Defense Is a Good Offense: Adversarial Attacks to Avoid Modulation Detection

abstract: We consider a communication scenario, in which an intruder tries to determine the modulation scheme of the intercepted signal. Our aim is to minimize the accuracy of the intruder, while guaranteeing that the intended receiver can still recover the underlying message with the highest reliability. This is achieved by perturbing channel input symbols at the encoder, similarly to adversarial attacks against classifiers in machine learning. In image classification, the perturbation is limited to be imperceptible to a human observer, while in our case the perturbation is constrained so that the message can still be reliably decoded by the legitimate receiver, which is oblivious to the perturbation. Simulation results demonstrate the viability of our approach to make wireless communication secure against state-of-the-art intruders (using deep learning or decision trees) with minimal sacrifice in the communication performance. On the other hand, we also demonstrate that using diverse training data and curriculum learning can significantly boost the accuracy of the intruder.


title: High-dimensional Bayesian optimization using low-dimensional feature spaces

abstract: Bayesian optimization (BO) is a powerful approach for seeking the global optimum of expensive black-box functions and has proven successful for fine tuning hyper-parameters of machine learning models. However, BO is practically limited to optimizing 10--20 parameters. To scale BO to high dimensions, we usually make structural assumptions on the decomposition of the objective and\slash or exploit the intrinsic lower dimensionality of the problem, e.g. by using linear projections. We could achieve a higher compression rate with nonlinear projections, but learning these nonlinear embeddings typically requires much data. This contradicts the BO objective of a relatively small evaluation budget. To address this challenge, we propose to learn a low-dimensional feature space jointly with (a) the response surface and (b) a reconstruction mapping. Our approach allows for optimization of BO's acquisition function in the lower-dimensional subspace, which significantly simplifies the optimization problem. We reconstruct the original parameter space from the lower-dimensional subspace for evaluating the black-box function. For meaningful exploration, we solve a constrained optimization problem.


title: Real-time City-scale Ridesharing via Linear Assignment Problems

abstract: In this paper, we propose a novel, computational efficient, dynamic ridesharing algorithm. The beneficial computational properties of the algorithm arise from casting the ridesharing problem as a linear assignment problem between fleet vehicles and customer trip requests within a federated optimization architecture. The resulting algorithm is up to four times faster than the state-of-the-art, even if it is implemented on a less dedicated hardware, and achieves similar service quality. Current literature showcases the ability of state-of-the-art ridesharing algorithms to tackle very large fleets and customer requests in almost near real-time, but the benefits of ridesharing seem limited to centralized systems. Our algorithm suggests that this does not need to be the case. The algorithm that we propose is fully distributable among multiple ridesharing companies. By leveraging two datasets, the New York city taxi dataset and the Melbourne Metropolitan Area dataset, we show that with our algorithm, real-time ridesharing offers clear benefits with respect to more traditional taxi fleets in terms of level of service, even if one considers partial adoption of the system. In fact, e.g., the quality of the solutions obtained in the state-of-the-art works that tackle the whole customer set of the New York city taxi dataset is achieved, even if one considers only a proportion of the fleet size and customer requests. This could make real-time urban-scale ridesharing very attractive to small enterprises and city authorities alike. However, in some cases, e.g., in multi-company scenarios where companies have predefined market shares, we show that the number of vehicles needed to achieve a comparable performance to the monopolistic setting increases, and this raises concerns on the possible negative effects of multi-company ridesharing.


title: On Constrained Open-World Probabilistic Databases

abstract: Increasing amounts of available data have led to a heightened need for representing large-scale probabilistic knowledge bases. One approach is to use a probabilistic database, a model with strong assumptions that allow for efficiently answering many interesting queries. Recent work on open-world probabilistic databases strengthens the semantics of these probabilistic databases by discarding the assumption that any information not present in the data must be false. While intuitive, these semantics are not sufficiently precise to give reasonable answers to queries. We propose overcoming these issues by using constraints to restrict this open world. We provide an algorithm for one class of queries, and establish a basic hardness result for another. Finally, we propose an efficient and tight approximation for a large class of queries.


title: Programmable linear quantum networks with a multimode fibre

abstract: Reconfigurable quantum circuits are fundamental building blocks for the implementation of scalable quantum technologies. Their implementation has been pursued in linear optics through the engineering of sophisticated interferometers. While such optical networks have been successful in demonstrating the control of small-scale quantum circuits, scaling up to larger dimensions poses significant challenges. Here, we demonstrate a potentially scalable route towards reconfigurable optical networks based on the use of a multimode fibre and advanced wavefront-shaping techniques. We program networks involving spatial and polarisation modes of the fibre and experimentally validate the accuracy and robustness of our approach using two-photon quantum states. In particular, we illustrate the reconfigurability of our platform by emulating a tunable coherent absorption experiment. By demonstrating reliable reprogrammable linear transformations, with the prospect to scale, our results highlight the potential of complex media driven by wavefront shaping for quantum information processing.


title: Increasing the representation accuracy of quantum simulations of chemistry without extra quantum resources

abstract: Proposals for near-term experiments in quantum chemistry on quantum computers leverage the ability to target a subset of degrees of freedom containing the essential quantum behavior, sometimes called the active space. This approximation allows one to treat more difficult problems using fewer qubits and lower gate depths than would otherwise be possible. However, while this approximation captures many important qualitative features, it may leave the results wanting in terms of absolute accuracy (basis error) of the representation. In traditional approaches, increasing this accuracy requires increasing the number of qubits and an appropriate increase in circuit depth as well. Here we introduce a technique requiring no additional qubits or circuit depth that is able to remove much of this approximation in favor of additional measurements. The technique is constructed and analyzed theoretically, and some numerical proof of concept calculations are shown. As an example, we show how to achieve the accuracy of a 20 qubit representation using only 4 qubits and a modest number of additional measurements for a simple hydrogen molecule. We close with an outlook on the impact this technique may have on both near-term and fault-tolerant quantum simulations.


title: Analyzing the Perceived Severity of Cybersecurity Threats Reported on Social Media

abstract: Breaking cybersecurity events are shared across a range of websites, including security blogs (FireEye, Kaspersky, etc.), in addition to social media platforms such as Facebook and Twitter. In this paper, we investigate methods to analyze the severity of cybersecurity threats based on the language that is used to describe them online. A corpus of 6,000 tweets describing software vulnerabilities is annotated with authors' opinions toward their severity. We show that our corpus supports the development of automatic classifiers with high precision for this task. Furthermore, we demonstrate the value of analyzing users' opinions about the severity of threats reported online as an early indicator of important software vulnerabilities. We present a simple, yet effective method for linking software vulnerabilities reported in tweets to Common Vulnerabilities and Exposures (CVEs) in the National Vulnerability Database (NVD). Using our predicted severity scores, we show that it is possible to achieve a Precision@50 of 0.86 when forecasting high severity vulnerabilities, significantly outperforming a baseline that is based on tweet volume. Finally we show how reports of severe vulnerabilities online are predictive of real-world exploits.


title: An efficient protocol of quantum walk in circuit QED

abstract: Implementation of discrete-time quantum walk (DTQW) with superconducting qubits is difficult since on-chip superconducting qubits cannot hop between lattice sites. We propose an efficient protocol for the implementation of DTQW in circuit quantum electrodynamics (QED), in which only $N+1$ qutrits and $N$ assistant cavities are needed for an $N$-step DTQW. The operation of each DTQW step is very quick because only resonant processes are adopted. The numerical simulations show that high-similarity DTQW with the number of step up to $20$ is feasible with present-day circuit QED technique. This protocol can help to study properties and applications of large-step DTQW in experiments, which is important for the development of quantum computation and quantum simulation in circuit QED.