title: Scalable Robust Kidney Exchange

abstract: In barter exchanges, participants directly trade their endowed goods in a constrained economic setting without money. Transactions in barter exchanges are often facilitated via a central clearinghouse that must match participants even in the face of uncertainty---over participants, existence and quality of potential trades, and so on. Leveraging robust combinatorial optimization techniques, we address uncertainty in kidney exchange, a real-world barter market where patients swap (in)compatible paired donors. We provide two scalable robust methods to handle two distinct types of uncertainty in kidney exchange---over the quality and the existence of a potential match. The latter case directly addresses a weakness in all stochastic-optimization-based methods to the kidney exchange clearing problem, which all necessarily require explicit estimates of the probability of a transaction existing---a still-unsolved problem in this nascent market. We also propose a novel, scalable kidney exchange formulation that eliminates the need for an exponential-time constraint generation process in competing formulations, maintains provable optimality, and serves as a subsolver for our robust approach. For each type of uncertainty we demonstrate the benefits of robustness on real data from a large, fielded kidney exchange in the United States. We conclude by drawing parallels between robustness and notions of fairness in the kidney exchange setting.


title: Reducing Autocorrelation Times in Lattice Simulations with Generative Adversarial Networks

abstract: Short autocorrelation times are essential for a reliable error assessment in Monte Carlo simulations of lattice systems. In many interesting scenarios, the decay of autocorrelations in the Markov chain is prohibitively slow. Generative samplers can provide statistically independent field configurations, thereby potentially ameliorating these issues. In this work, the applicability of neural samplers to this problem is investigated. Specifically, we work with a generative adversarial network (GAN). We propose to address difficulties regarding its statistical exactness through the implementation of an overrelaxation step, by searching the latent space of the trained generator network. This procedure can be incorporated into a standard Monte Carlo algorithm, which then permits a sensible assessment of ergodicity and balance based on consistency checks. Numerical results for real, scalar $\phi^4$-theory in two dimensions are presented. We achieve a significant reduction of autocorrelations while accurately reproducing the correct statistics. We discuss possible improvements to the approach as well as potential solutions to persisting issues.


title: MOVPE growth of nitrogen- and aluminum-polar AlN on 4H-SiC

abstract: We present a comprehensive study on metal-organic vapor phase epitaxy growth of N-polar and Al-polar AlN on 4H-SiC with 4 deg miscut using constant growth parameters. At a high temperature of 1165 degC, N-polar AlN layers had high crystalline quality whereas the Al-polar AlN surfaces had a high density of etch pits. For N-polar AlN, the V/III ratio below 1000 forms hexagonal hillocks, while the V/III ratio over 1000 yields step bunching without the hillocks. 1-um-thick N-polar AlN layer grown in optimal conditions exhibited FWHMs of 307, 330 and 337 arcsec for (002), (102) and (201) reflections, respectively.


title: Learning Dense Stereo Matching for Digital Surface Models from Satellite Imagery

abstract: Digital Surface Model generation from satellite imagery is a difficult task that has been largely overlooked by the deep learning community. Stereo reconstruction techniques developed for terrestrial systems including self driving cars do not translate well to satellite imagery where image pairs vary considerably. In this work we present neural network tailored for Digital Surface Model generation, a ground truthing and training scheme which maximizes available hardware, and we present a comparison to existing methods. The resulting models are smooth, preserve boundaries, and enable further processing. This represents one of the first attempts at leveraging deep learning in this domain.


title: Iterative Filtering as a direct method for the decomposition of non-stationary signals

abstract: The Iterative Filtering method is a technique developed recently for the decomposition and analysis of non-stationary and non-linear signals. In this work we propose two alternative formulations of the original algorithm which allows to transform the Iterative Filtering method into a direct technique, making the algorithm closer to an online algorithm. We present a few numerical examples to show the effectiveness of the proposed approaches.


title: Iterative Classroom Teaching

abstract: We consider the machine teaching problem in a classroom-like setting wherein the teacher has to deliver the same examples to a diverse group of students. Their diversity stems from differences in their initial internal states as well as their learning rates. We prove that a teacher with full knowledge about the learning dynamics of the students can teach a target concept to the entire classroom using O(min{d,N} log(1/eps)) examples, where d is the ambient dimension of the problem, N is the number of learners, and eps is the accuracy parameter. We show the robustness of our teaching strategy when the teacher has limited knowledge of the learners' internal dynamics as provided by a noisy oracle. Further, we study the trade-off between the learners' workload and the teacher's cost in teaching the target concept. Our experiments validate our theoretical results and suggest that appropriately partitioning the classroom into homogenous groups provides a balance between these two objectives.


title: Integrating Security in Resource-Constrained Cyber-Physical Systems

abstract: Defense mechanisms against network-level attacks are commonly based on the use of cryptographic techniques, such as message authentication codes that provide data integrity guarantees. However, such mechanisms require significant resources, which prevents their continuous use in resource-constrained cyber-physical systems. Recently, it was shown how physical properties of plants can be exploited to relax these requirements for systems where sensor measurements and actuator commands are transmitted over a compromised network; specifically, intermittent use of data authentication, can still provide Quality-of-Control (QoC) guarantees even in the presence of false-data injection attacks. Consequently, in this work we focus on integrating security into existing systems, in order to protect against these attacks. We introduce a design-time methodology that incorporates requirements for QoC in the presence of attacks into end-to-end timing constraints for real-time control transactions, which include data acquisition and authentication, communication, and control. This allows us to formulate a mixed integer linear programming-based method for synthesis of schedulable task and message parameters (i.e., deadlines and offsets) that maintain timing requirements of deployed controllers, while adding a sufficient level of protection against attacks; specifically, this method provides suitable intermittent authentication policies that ensure the desired QoC levels under attack. To additionally reduce the security-related bandwidth overhead, we propose the use of cumulative message authentication. Furthermore, we introduce a method for opportunistic use of remaining resources to further improve the overall QoC guarantees while ensuring system schedulability. Finally, we demonstrate applicability of our methodology on synthetic automotive systems as well as an automotive case-study.


title: Uncovering the Social Interaction in Swarm Intelligence with Network Science

abstract: Swarm intelligence is the collective behavior emerging in systems with locally interacting components. Because of their self-organization capabilities, swarm-based systems show essential properties for handling real-world problems such as robustness, scalability, and flexibility. Yet, we do not know why swarm-based algorithms work well and neither we can compare the different approaches in the literature. The lack of a common framework capable of characterizing these several swarm-based algorithms, transcending their particularities, has led to a stream of publications inspired by different aspects of nature without a systematic comparison over existing approaches. Here, we address this gap by introducing a network-based framework---the interaction network---to examine computational swarm-based systems via the optics of the social dynamics of such interaction network; a clear example of network science being applied to bring further clarity to a complicated field within artificial intelligence. We discuss the social interactions of four well-known swarm-based algorithms and provide an in-depth case study of the Particle Swarm Optimization. The interaction network enables researchers to study swarm algorithms as systems, removing the algorithm particularities from the analyses while focusing on the structure of the social interactions.


title: A note on quasi-local energy

abstract: A definition of quasi-local energy in a gravitational field based upon its embedding into flat space is discussed. The outcome is not satisfactory from many points of view.


title: Towards Connecting Experiences during Collocated Events through Data Mining and Visualization

abstract: Themed collocated events, such as conferences, workshops, and seminars, invite people with related life experiences to connect with each other. In this era when people record lives through the Internet, individual experiences exist in different forms of digital contents. People share digital life records during collocated events, such as sharing blogs they wrote, Twitter posts they forwarded, and books they have read. However, connecting experiences during collocated events are challenging. Not only one is blind to the large contents of others, identifying related experiential items depends on how well experiences are retrieved. The collection of personal contents from all participants forms a valuable group repository, from which connections between experiences can be mined. Visualizing same or related experiences inspire conversations and support social exchange. Common topics in group content also help participants generate new perspectives about the collocated group. Advances in machine learning and data visualization provide automated approaches to process large data and enable interactions with data repositories. This position paper promotes the idea of event mining: how to utilize state-of-the-art data processing and visualization techniques to design event mining systems for connecting experiences during collocated activities. We discuss empirical and constructive problems in this design space, and our preliminary study of deploying a tabletop-based system, BlogCloud, which supports experience re-visitation and exchange with machine-learning and data visualization.


title: Adaptive Semantic Segmentation with a Strategic Curriculum of Proxy Labels

abstract: Training deep networks for semantic segmentation requires annotation of large amounts of data, which can be time-consuming and expensive. Unfortunately, these trained networks still generalize poorly when tested in domains not consistent with the training data. In this paper, we show that by carefully presenting a mixture of labeled source domain and proxy-labeled target domain data to a network, we can achieve state-of-the-art unsupervised domain adaptation results. With our design, the network progressively learns features specific to the target domain using annotation from only the source domain. We generate proxy labels for the target domain using the network's own predictions. Our architecture then allows selective mining of easy samples from this set of proxy labels, and hard samples from the annotated source domain. We conduct a series of experiments with the GTA5, Cityscapes and BDD100k datasets on synthetic-to-real domain adaptation and geographic domain adaptation, showing the advantages of our method over baselines and existing approaches.


title: Predicative well-ordering

abstract: Confusion over the predicativist conception of well-ordering pervades the literature and is responsible for widespread fundamental misconceptions about the nature of predicative reasoning. This short note aims to explain the principal fallacy, first noted in [N. Weaver, Predicativity beyond Gamma_0, arXiv:math/0509244], and some of its consequences.


title: Random Phase Approximation in Projected Oscillator Orbitals

abstract: The projected oscillator orbitals (pOOs) are localized virtual orbitals constructed by multiplying localized occupied orbitals by harmonics. Following a recent paper by Mussard and \'Angy\'an (Theor Chem Acc 134:1, 2015), further developments of projected oscillator orbitals are shown, notably the equations for pOOs of general order as well as their overlaps are derived. The performance of these localized virtual orbitals is demonstrated up to third order. It is found that a good fraction of the aug-cc-pVQZ RPA correlation energy is recovered by use of a smaller number of pOOs. This is especially true where considering only the long-range correlation energy, which is important for the description of London dispersion forces.


title: Merits and constraints of low-${\bm K^2}$ experimental data for the proton radius determination

abstract: The question of the determination of the proton charge radius $R_p$ from electron scattering data led to an unprecedented experimental effort for measurements of the electric form factor of the proton at low and very low momentum transfer in electron and muon elastic scattering. On the basis of basic properties of densities and fitting bias considerations, a procedure is developed in order to evaluate the impact of forthcoming data on $R_p$. Particularly, it is shown that a 0.1\% precision on these future cross section data is necessary to establish indisputably the $R_p$-value as determined from lepton scattering. The ProRad (Proton Radius) experiment at the PRAE (Platform for Research and Applications with Electrons) facility in Orsay is further discussed, especially the experimental method to meet this stringent constraint.


title: On simple spectral modules over Leavitt path algebras

abstract: In this work we describe all simple modules over Leavitt path algebras as induced modules from irreducible representations of the isotropy groups.


title: On the Erd\H{o}s Covering Problem: the density of the uncovered set

abstract: Since their introduction by Erd\H{o}s in 1950, covering systems (that is, finite collections of arithmetic progressions that cover the integers) have been extensively studied, and numerous questions and conjectures have been posed regarding the existence of covering systems with various properties. In particular, Erd\H{o}s asked if the moduli can be distinct and all arbitrarily large, Erd\H{o}s and Selfridge asked if the moduli can be distinct and all odd, and Schinzel conjectured that in any covering system there exists a pair of moduli, one of which divides the other. Another beautiful conjecture, proposed by Erd\H{o}s and Graham in 1980, states that if the moduli are distinct elements of the interval $[n,Cn]$, and $n$ is sufficiently large, then the density of integers uncovered by the union is bounded below by a constant (depending only on $C$). This conjecture was confirmed (in a strong form) by Filaseta, Ford, Konyagin, Pomerance and Yu in 2007, who moreover asked whether the same conclusion holds if the moduli are distinct and sufficiently large, and $\sum_{i=1}^k \frac{1}{d_i} < C$. Although this condition turns out not to be sufficiently strong to imply the desired conclusion, as the main result of this paper we will give an essentially best possible condition which is sufficient. Our method has a number of further applications. Most importantly, we prove the conjecture of Schinzel stated above, which was made in 1967. We moreover give an alternative (somewhat simpler) proof of a breakthrough result of Hough, who resolved Erd\H{o}s' minimum modulus problem, with an improved bound on the smallest difference. Finally, we make further progress on the problem of Erd\H{o}s and Selfridge.


title: Blazar VHE spectral alterations induced by photon-ALP oscillations

abstract: Prompted by the increasing interest of axion-like particles (ALPs) for very-high-energy (VHE) astrophysics, we have considered a full scenario for the propagation of a VHE photon/ALP beam emitted by a BL Lac and reaching us in the light of the most up-to-date astrophysical information and for energies up to above $100 \, \rm TeV$. During its trip, the beam -- generated in a small region of a BL Lac jet -- crosses a variety of magnetic structures in very different astronomical environments: the BL Lac jet, the host elliptical galaxy, the extragalactic space and the Milky Way. We have taken an effort to model all these magnetic fields in the most realistic fashion and using a new model developed by us concerning the extragalactic magnetic field. Assuming an intrinsic spectrum with a power law exponentially truncated at a fixed cut-off energy, we have evaluated the resulting observed spectra of Markarian 501, the extreme BL Lac 1ES 0229+200 and a similar source located at $z = 0.6$ up to above $100 \, \rm TeV$. We obtain interesting results: the model with photon-ALP oscillations possesses features (spectral energy oscillatory behaviour and photon excess above $20 \, \rm TeV$) which can be tested by $\gamma$-ray observatories like CTA, HAWC, GAMMA 400, LHAASO, TAIGA-HiSCORE and HERD. In addition, our ALP can be detected in dedicated laboratory experiments like the upgrade of ALPS II at DESY, the planned IAXO and STAX experiments, as well as with other techniques developed by Avignone and collaborators.


title: An End-to-end Approach to Semantic Segmentation with 3D CNN and Posterior-CRF in Medical Images

abstract: Fully-connected Conditional Random Field (CRF) is often used as post-processing to refine voxel classification results by encouraging spatial coherence. In this paper, we propose a new end-to-end training method called Posterior-CRF. In contrast with previous approaches which use the original image intensity in the CRF, our approach applies 3D, fully connected CRF to the posterior probabilities from a CNN and optimizes both CNN and CRF together. The experiments on white matter hyperintensities segmentation demonstrate that our method outperforms CNN, post-processing CRF and different end-to-end training CRF approaches.


title: A novel approach to the computation of one-loop three- and four-point functions. I -- The real mass case

abstract: This article is the first of a series of three presenting an alternative method to compute the one-loop scalar integrals. This novel method enjoys a couple of interesting features as compared with the method closely following 't Hooft and Veltman adopted previously. It directly proceeds in terms of the quantities driving algebraic reduction methods. It applies to the three-point functions and, in a similar way, to the four-point functions. It also extends to complex masses without much complication. Lastly, it extends to kinematics more general than the one of physical e.g.\ collider processes relevant at one loop. This last feature may be useful when considering the application of this method beyond one loop using generalised one-loop integrals as building blocks.


title: Quasi-Local Energy of a Charged Rotating Object Described by the Kerr-Newman Metric

abstract: The Brown-York quasi-local energy of a charged rotating black hole described by the Kerr-Newman metric and enclosed by a fixed-radius surface is computed. No further assumptions on the angular momentum or the radial coordinate in Boyer-Lindquist coordinates were made. The result can be expressed in terms of incomplete elliptic integrals and is used to analyze the self-energy of an electron which is assumed to be described by the Kerr-Newman metric. For this purpose the small sphere limit is investigated thoroughly comparing the analysis with known results. Evaluating the energy using the mass, angular momentum and charge of an electron a value in the order of the Planck energy is obtained in the small sphere limit as long as the Kerr-Newman metric acting as exterior solution can be used as a description of spacetime.


title: Dyson equations for correlators of Wilson loops

abstract: By considering a Gaussian truncation of ${\cal N}=4$ super Yang-Mills, we derive a set of Dyson equations that account for the ladder diagram contribution to connected correlators of circular Wilson loops. We consider different numbers of loops, with different relative orientations. We show that the Dyson equations admit a spectral representation in terms of eigenfunctions of a Schr\"odinger problem, whose classical limit describes the strong coupling limit of the ladder resummation. We also verify that in supersymmetric cases the exact solution to the Dyson equations reproduces known matrix model results.


title: Mixed moment of $GL(2)$ and $GL(3)$ $L$-functions

abstract: Let $ \mathfrak{f} $ run over the space $ H_{4k} $ of primitive cusp forms of level one and weight $ 4k $, $ k \in N $. We prove an explicit formula for the mixed moment of the Hecke $ L $-function $ L(\mathfrak{f}, 1/2) $ and the symmetric square $L$-function $ L(sym^2\mathfrak{f}, 1/2)$, relating it to the dual mixed moment of the double Dirichlet series and the Riemann zeta function weighted by the ${}_3F_{2}$ hypergeometric function. Analysing the corresponding special functions by the means of the Liouville-Green approximation followed by the saddle point method, we prove that the initial mixed moment is bounded by $\log^3k$.


title: Implicit Argument Prediction as Reading Comprehension

abstract: Implicit arguments, which cannot be detected solely through syntactic cues, make it harder to extract predicate-argument tuples. We present a new model for implicit argument prediction that draws on reading comprehension, casting the predicate-argument tuple with the missing argument as a query. We also draw on pointer networks and multi-hop computation. Our model shows good performance on an argument cloze task as well as on a nominal implicit argument prediction task.


title: Modular Architecture for StarCraft II with Deep Reinforcement Learning

abstract: We present a novel modular architecture for StarCraft II AI. The architecture splits responsibilities between multiple modules that each control one aspect of the game, such as build-order selection or tactics. A centralized scheduler reviews macros suggested by all modules and decides their order of execution. An updater keeps track of environment changes and instantiates macros into series of executable actions. Modules in this framework can be optimized independently or jointly via human design, planning, or reinforcement learning. We apply deep reinforcement learning techniques to training two out of six modules of a modular agent with self-play, achieving 94% or 87% win rates against the "Harder" (level 5) built-in Blizzard bot in Zerg vs. Zerg matches, with or without fog-of-war.


title: Topological strings, strips and quivers

abstract: We find a direct relation between quiver representation theory and open topological string theory on a class of toric Calabi-Yau manifolds without compact four-cycles, also referred to as strip geometries. We show that various quantities that characterize open topological string theory on these manifolds, such as partition functions, Gromov-Witten invariants, or open BPS invariants, can be expressed in terms of characteristics of the moduli space of representations of the corresponding quiver. This has various deep consequences; in particular, expressing open BPS invariants in terms of motivic Donaldson-Thomas invariants, immediately proves integrality of the former ones. Taking advantage of the relation to quivers we also derive explicit expressions for classical open BPS invariants for an arbitrary strip geometry, which lead to a large set of number theoretic integrality statements. Furthermore, for a specific framing, open topological string partition functions for strip geometries take form of generalized $q$-hypergeometric functions, which leads to a novel representation of these functions in terms of quantum dilogarithms and integral invariants. We also study quantum curves and A-polynomials associated to quivers, various limits thereof, and their specializations relevant for strip geometries. The relation between toric manifolds and quivers can be regarded as a generalization of the knots-quivers correspondence to more general Calabi-Yau geometries.


title: Efficient Numerical Algorithms based on Difference Potentials for Chemotaxis Systems in 3D

abstract: In this work, we propose efficient and accurate numerical algorithms based on Difference Potentials Method for numerical solution of chemotaxis systems and related models in 3D. The developed algorithms handle 3D irregular geometry with the use of only Cartesian meshes and employ Fast Poisson Solvers. In addition, to further enhance computational efficiency of the methods, we design a Difference-Potentials-based domain decomposition approach which allows mesh adaptivity and easy parallelization of the algorithm in space. Extensive numerical experiments are presented to illustrate the accuracy, efficiency and robustness of the developed numerical algorithms.


title: Iterated Integrals and Population Time Series Analysis

abstract: One of the core advantages topological methods for data analysis provide is that the language of (co)chains can be mapped onto the semantics of the data, providing a natural avenue for human understanding of the results. Here, we describe such a semantic structure on Chen's classical iterated integral cochain model for paths in Euclidean space. Specifically, in the context of population time series data, we observe that iterated integrals provide a model-free measure of pairwise influence that can be used for causality inference. Along the way, we survey recent results and applications, review the current standard methods for causality inference, and briefly provide our outlook on generalizations to go beyond time series data.


title: A Feature Complete SPIKE Banded Algorithm and Solver

abstract: New features and enhancements for the SPIKE banded solver are presented. Among all the SPIKE algorithm versions, we focus our attention on the recursive SPIKE technique which provides the best trade-off between generality and parallel efficiency, but was known for its lack of flexibility. Its application was essentially limited to power of two number of cores/processors. This limitation is successfully addressed in this paper. In addition, we present a new transpose solve option, a standard feature of most numerical solver libraries which has never been addressed by the SPIKE algorithm so far. A pivoting recursive SPIKE strategy is finally presented as an alternative to non-pivoting scheme for systems with large condition numbers. All these new enhancements participate to create a feature complete SPIKE algorithm and a new black-box SPIKE-OpenMP package that significantly outperforms the performance and scalability obtained with other state-of-the-art banded solvers.


title: Solar Sources of Interplanetary Magnetic Clouds Leading to Helicity Prediction

abstract: This study identifies the solar origins of magnetic clouds that are observed at 1 AU and predicts the helical handedness of these clouds from the solar surface magnetic fields. We started with the magnetic clouds listed by the Magnetic Field Investigation (MFI) team supporting NASA's WIND spacecraft in what is known as the MFI table and worked backwards in time to identify solar events that produced these clouds. Our methods utilize magnetograms from the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO) spacecraft so that we could only analyze MFI entries after the beginning of 2011. This start date and the end date of the MFI table gave us 37 cases to study. Of these we were able to associate only eight surface events with clouds detected by WIND at 1 AU. We developed a simple algorithm for predicting the cloud helicity which gave the correct handedness in all eight cases. The algorithm is based on the conceptual model that an ejected flux tube has two magnetic origination points at the positions of the strongest radial magnetic field regions of opposite polarity near the places where the ejected arches end at the solar surface. We were unable to find events for the remaining 29 cases: lack of a halo or partial halo CME in an appropriate time window, lack of magnetic and/or filament activity in the proper part of the solar disk, or the event was too far from disk center. The occurrence of a flare was not a requirement for making the identification but in fact flares, often weak, did occur for seven of the eight cases.


title: $R_{D^{(*)}}$ motivated $\mathcal{S}_1$ leptoquark scenarios: Impact of interference on the exclusion limits from LHC data

abstract: Motivated by the persistent anomalies in the semileptonic $B$-meson decays, we investigate the competency of LHC data to constrain the $R_{D^{(*)}}$-favoured parameter space in a charge $-1/3$ scalar leptoquark ($\mathcal S_1$) model. We consider some scenarios with one large free coupling to accommodate the $R_{D^{(*)}}$ anomalies. As a result, some of them dominantly yield nonresonant $\tau\tau$ and $\tau\nu$ events at the LHC through the $t$-channel $\mathcal S_1$ exchange. So far, no experiment has searched for leptoquarks using these signatures and the relevant resonant leptoquark searches are yet to put any strong exclusion limit on the parameter space. We recast the latest $\tau\tau$ and $\tau\nu$ resonance search data to obtain new exclusion limits. The nonresonant processes strongly interfere (destructively in our case) with the Standard Model background and play the determining role in setting the exclusion limits. To obtain precise limits, we include non-negligible effects coming from the subdominant (resonant) pair and inclusive single leptoquark productions systematically in our analysis. To deal with large destructive interference, we make use of the transverse mass distributions from the experiments in our statistical analysis. In addition, we also recast the relevant direct search results to obtain the most stringent collider bounds on these scenarios to date. These are independent bounds and are competitive to other known bounds. Finally, we indicate how one can adopt these bounds to a wide class of models with $\mathcal S_1$ that are proposed to accommodate the $R_{D^{(*)}}$ anomalies.


title: Real time Traffic Flow Parameters Prediction with Basic Safety Messages at Low Penetration of Connected Vehicles

abstract: The expected low market penetration of connected vehicles (CVs) in the near future could be a constraint in estimating traffic flow parameters, such as average travel speed of a roadway segment and average space headway between vehicles from the CV broadcasted data. This estimated traffic flow parameters from low penetration of connected vehicles become noisy compared to 100 percent penetration of CVs, and such noise reduces the real time prediction accuracy of a machine learning model, such as the accuracy of long short term memory (LSTM) model in terms of predicting traffic flow parameters. The accurate prediction of the parameters is important for future traffic condition assessment. To improve the prediction accuracy using noisy traffic flow parameters, which is constrained by limited CV market penetration and limited CV data, we developed a real time traffic data prediction model that combines LSTM with Kalman filter based Rauch Tung Striebel (RTS) noise reduction model. We conducted a case study using the Enhanced Next Generation Simulation (NGSIM) dataset, which contains vehicle trajectory data for every one tenth of a second, to evaluate the performance of this prediction model. Compared to a baseline LSTM model performance, for only 5 percent penetration of CVs, the analyses revealed that combined LSTM and RTS model reduced the mean absolute percentage error (MAPE) from 19 percent to 5 percent for speed prediction and from 27 percent to 9 percent for space-headway prediction. The statistical significance test with a 95 percent confidence interval confirmed no significant difference in predicted average speed and average space headway using this LSTM and RTS combination with only 5 percent CV penetration rate.


title: LAAIR: A Layered Architecture for Autonomous Interactive Robots

abstract: When developing general purpose robots, the overarching software architecture can greatly affect the ease of accomplishing various tasks. Initial efforts to create unified robot systems in the 1990s led to hybrid architectures, emphasizing a hierarchy in which deliberative plans direct the use of reactive skills. However, since that time there has been significant progress in the low-level skills available to robots, including manipulation and perception, making it newly feasible to accomplish many more tasks in real-world domains. There is thus renewed optimism that robots will be able to perform a wide array of tasks while maintaining responsiveness to human operators. However, the top layer in traditional hybrid architectures, designed to achieve long-term goals, can make it difficult to react quickly to human interactions during goal-driven execution. To mitigate this difficulty, we propose a novel architecture that supports such transitions by adding a top-level reactive module which has flexible access to both reactive skills and a deliberative control module. To validate this architecture, we present a case study of its application on a domestic service robot platform.


title: Description of two-electron atoms with correct cusp conditions

abstract: New sets of functions with arbitrary large finite cardinality are constructed for two-electron atoms. Functions from these sets exactly satisfy the Kato's cusp conditions. The new functions are special linear combinations of Hylleraas- and/or Kinoshita-type terms. Standard variational calculation, leading to matrix eigenvalue problem, can be carried out to calculate the energies of the system. There is no need for optimization with constraints to satisfy the cusp conditions. In the numerical examples the ground state energy of the He atom is considered.


title: The influence of the distribution of cosmic star formation at different metallicities on the properties of merging double compact objects

abstract: Binaries that merge within the local Universe originate from progenitor systems that formed at different times and in various environments. The efficiency of formation of double compact objects is highly sensitive to metallicity of the star formation. Therefore, to confront the theoretical predictions with observational limits resulting from gravitational waves observations one has to account for the formation and evolution of progenitor stars in the chemically evolving Universe. In particular, this requires knowledge of the distribution of cosmic star formation rate at different metallicities and times, probed by redshift (SFR(Z,z)). We investigate the effect of the assumed SFR(Z,z) on the properties of merging double compact objects, in particular on their merger rate densities. Using a set of binary evolution models from Chruslinska et al. (2018) we demonstrate that the reported tension between the merger rates of different types of double compact objects and current observational limits in some cases can be resolved if a SFR(Z,z) closer to that expected based on observations of local star-forming galaxies is used, without the need for changing the assumptions about the evolution of progenitor stars of different masses. This highlights the importance of finding tighter constraints on SFR(Z,z) and understanding the associated uncertainties.


title: A Natural Language Interface with Relayed Acoustic Communications for Improved Command and Control of AUVs

abstract: Autonomous underwater vehicles (AUVs) are being tasked with increasingly complex missions. The acoustic communications required for AUVs are, by the nature of the medium, low bandwidth while adverse environmental conditions underwater often mean they are also intermittent. This has motivated development of highly autonomous systems, which can operate independently of their operators for considerable periods of time. These missions often involve multiple vehicles leading not only to challenges in communications but also in command and control (C2). Specifically operators face complexity in controlling multi-objective, multi-vehicle missions, whilst simultaneously facing uncertainty over the current status and safety of several remote high value assets. Additionally, it may be required to perform command and control of these complex missions in a remote control room. In this paper, we propose a combination of an intuitive, natural language operator interface combined with communications that use platforms from multiple domains to relay data over different mediums and transmission modes, improving command and control of collaborative and fully autonomous missions. In trials, we have demonstrated an integrated system combining working prototypes with established commercial C2 software that enables the use of a natural language interface to monitor an AUV survey mission in an on-shore command and control centre.


title: Biologically-plausible learning algorithms can scale to large datasets

abstract: The backpropagation (BP) algorithm is often thought to be biologically implausible in the brain. One of the main reasons is that BP requires symmetric weight matrices in the feedforward and feedback pathways. To address this "weight transport problem" (Grossberg, 1987), two more biologically plausible algorithms, proposed by Liao et al. (2016) and Lillicrap et al. (2016), relax BP's weight symmetry requirements and demonstrate comparable learning capabilities to that of BP on small datasets. However, a recent study by Bartunov et al. (2018) evaluate variants of target-propagation (TP) and feedback alignment (FA) on MINIST, CIFAR, and ImageNet datasets, and find that although many of the proposed algorithms perform well on MNIST and CIFAR, they perform significantly worse than BP on ImageNet. Here, we additionally evaluate the sign-symmetry algorithm (Liao et al., 2016), which differs from both BP and FA in that the feedback and feedforward weights share signs but not magnitudes. We examine the performance of sign-symmetry and feedback alignment on ImageNet and MS COCO datasets using different network architectures (ResNet-18 and AlexNet for ImageNet, RetinaNet for MS COCO). Surprisingly, networks trained with sign-symmetry can attain classification performance approaching that of BP-trained networks. These results complement the study by Bartunov et al. (2018), and establish a new benchmark for future biologically plausible learning algorithms on more difficult datasets and more complex architectures.


title: A Geometric Approach of Gradient Descent Algorithms in Linear Neural Networks

abstract: In this paper, we propose a geometric framework to analyze the convergence properties of gradient descent trajectories in the context of linear neural networks. We translate a well-known empirical observation of linear neural nets into a conjecture that we call the \emph{overfitting conjecture} which states that, for almost all training data and initial conditions, the trajectory of the corresponding gradient descent system converges to a global minimum. This would imply that the solution achieved by vanilla gradient descent algorithms is equivalent to that of the least-squares estimation, for linear neural networks of an arbitrary number of hidden layers. Built upon a key invariance property induced by the network structure, we first establish convergence of gradient descent trajectories to critical points of the square loss function in the case of linear networks of arbitrary depth. Our second result is the proof of the \emph{overfitting conjecture} in the case of single-hidden-layer linear networks with an argument based on the notion of normal hyperbolicity and under a generic property on the training data (i.e., holding for almost all training data).


title: An Axiomatic Study of Query Terms Order in Ad-hoc Retrieval

abstract: Classic retrieval methods use simple bag-of-word representations for queries and documents. This representation fails to capture the full semantic richness of queries and documents. More recent retrieval models have tried to overcome this deficiency by using approaches such as incorporating dependencies between query terms, using bi-gram representations of documents, proximity heuristics, and passage retrieval. While some of these previous works have implicitly accounted for term order, to the best of our knowledge, term order has not been the primary focus of any research. In this paper, we focus solely on the effect of term order in information retrieval. We will show that documents that have two query terms in the same order as in the query have a higher probability of being relevant than documents that have two query terms in the reverse order. Using the axiomatic framework for information retrieval, we introduce a constraint that retrieval models must adhere to in order to effectively utilize term order dependency among query terms. We modify existing retrieval models based on this constraint so that if the order of a pair of query terms is semantically important, a document that includes these query terms in the same order as the query should receive a higher score compared to a document that includes them in the reverse order. Our empirical evaluation using both TREC newswire and web corpora demonstrates that the modified retrieval models significantly outperform their original counterparts.


title: Dynamics and stationary configurations of heterogeneous foams

abstract: We consider the variational foam model, where the goal is to minimize the total surface area of a collection of bubbles subject to the constraint that the volume of each bubble is prescribed. We apply sharp interface methods to develop an efficient computational method for this problem. In addition to simulating time dynamics, we also report on stationary states of this flow for <22 bubbles in two dimensions and <18 bubbles in three dimensions. For small numbers of bubbles, we recover known analytical results, which we briefly discuss. In two dimensions, we also recover the previous numerical results of Cox et. al. (2003), computed using other methods. Particular attention is given to locally optimal foam configurations and heterogeneous foams, where the volumes of the bubbles are not equal. Configurational transitions are reported for the quasi-stationary flow where the volume of one of the bubbles is varied and, for each volume, the stationary state is computed. The results from these numerical experiments are described and accompanied by many figures and videos.


title: Intrinsic Geometric Vulnerability of High-Dimensional Artificial Intelligence

abstract: The success of modern Artificial Intelligence (AI) technologies depends critically on the ability to learn non-linear functional dependencies from large, high dimensional data sets. Despite recent high-profile successes, empirical evidence indicates that the high predictive performance is often paired with low robustness, making AI systems potentially vulnerable to adversarial attacks. In this report, we provide a simple intuitive argument suggesting that high performance and vulnerability are intrinsically coupled, and largely dependent on the geometry of typical, high-dimensional data sets. Our work highlights a major potential pitfall of modern AI systems, and suggests practical research directions to ameliorate the problem.


title: Probing Pseudoscalars with Pulsar Polarisation Data Sets

abstract: Recently a data set containing linear and circular polarisation information of a collection of six hundred pulsars has been released. The operative radio wavelength for the same was 21 cm. Pulsars radio emission process is modelled either with synchroton / superconducting self-Compton route or with curvature radiation route. These theories fall short of accounting for the circular polarisation observed, as they are predisposed towards producing, solely, linear polarisation. Here we invoke (pseudo)scalars and their interaction with photons mediated by colossal magnetic fields of pulsars, to account for the circular part of polarisation data. This enables us to estimate the pseudoscalar parameters such as its coupling to photons and its mass in conjunction as product. To obtain these values separately, we turn our attention to recent observation on 47 pulsars, whose absolute polarisation position angles have been made available. Except, a third of the latter set, the rest of it overlaps with the expansive former data set on polarisation type \& degree. This helps us figure out, both the pseudoscalar parameters individually, that we report here.


title: Scale-variant topological information for characterizing the structure of complex networks

abstract: The structure of real-world networks is usually difficult to characterize owing to the variation of topological scales, the nondyadic complex interactions, and the fluctuations in the network. We aim to address these problems by introducing a general framework using a method based on topological data analysis. By considering the diffusion process at a single specified timescale in a network, we map the network nodes to a finite set of points that contains the topological information of the network at a single scale. Subsequently, we study the shape of these point sets over variable timescales that provide scale-variant topological information, to understand the varying topological scales and the complex interactions in the network. We conduct experiments on synthetic and real-world data to demonstrate the effectiveness of the proposed framework in identifying network models, classifying real-world networks, and detecting transition points in time-evolving networks. Overall, our study presents a unified analysis that can be applied to more complex network structures, as in the case of multilayer and multiplex networks.


title: Magnetic-field controlled anomalous refraction in doped semiconductors

abstract: We predict here that a slab made of a doped semiconductor can exhibit anomalous refraction under the application of a static magnetic field. This anomalous refraction takes place in the far-infrared range and it occurs for any angle of incidence. We show that this effect is due to the fact that a doped semiconductor under a magnetic field can behave, to some extent, as a hyperbolic metamaterial. We also show that the occurrence of this anomalous refraction enables a semiconductor slab under a magnetic field to partially focus the electromagnetic radiation. The remarkable thing in our case is that we deal with naturally occurring materials and the anomalous refraction can be tuned at will with an external field.


title: Large-Scale Visual Active Learning with Deep Probabilistic Ensembles

abstract: Annotating the right data for training deep neural networks is an important challenge. Active learning using uncertainty estimates from Bayesian Neural Networks (BNNs) could provide an effective solution to this. Despite being theoretically principled, BNNs require approximations to be applied to large-scale problems, where both performance and uncertainty estimation are crucial. In this paper, we introduce Deep Probabilistic Ensembles (DPEs), a scalable technique that uses a regularized ensemble to approximate a deep BNN. We conduct a series of large-scale visual active learning experiments to evaluate DPEs on classification with the CIFAR-10, CIFAR-100 and ImageNet datasets, and semantic segmentation with the BDD100k dataset. Our models require significantly less training data to achieve competitive performances, and steadily improve upon strong active learning baselines as the annotation budget is increased.


title: Interplay between electromagnetically induced transparency(EIT), absorption (EIA) and Autler-Townes (AT) splitting in N-type atomic system: An experiment and theory

abstract: In this article we have shown that the atomic states can be engineered by tunning the coupling Rabi frequency for a system with $\mathcal{N}$-type configuration. Electromagnetically induced transparency (EIT), Electromagnetically induced absorption (EIA) and Autler-Townes (AT) splitting has been observed experimentally in a four level $\mathcal{N}$-type atomic vapor of $^{85}Rb$ atoms in the hyperfine levels of $D_2$ transition. It has been shown that the response of the atomic medium can be tunned from highly transparent to highly absorptive in our case. The evolution of the atomic states from the dark state |D> to the non-coupled state |NC> has been studied with the partial dressed state approach which makes the backbone of the modification of the atomic response. In addition, transient solutions in the time domain and steady state solution in the frequency domain has been studied. The population dynamics and the coherence contribution in each case has been analyzed by the time dependent solutions. The experimentally observed steady line-shape profiles has been supported by the steady state solution of optical-Bloch equations considering the Maxwell Boltzmann velocity distributions of the atoms. It has been observed that the crossover between the EIT and the AT splitting has been replaced by the interference contribution of the EIA in this $\mathcal{N}$-type system.


title: Labeling Bias in Galaxy Morphologies

abstract: We present a metric to quantify systematic labeling bias in galaxy morphology data sets stemming from the quality of the labeled data. This labeling bias is independent from labeling errors and requires knowledge about the intrinsic properties of the data with respect to the observed properties. We conduct a relative comparison of label bias for different low redshift galaxy morphology data sets. We show our metric is able to recover previous de-biasing procedures based on redshift as biasing parameter. By using the image resolution instead, we find biases that have not been addressed. We find that the morphologies based on supervised machine-learning trained over features such as colors, shape, and concentration show significantly less bias than morphologies based on expert or citizen-science classifiers. This result holds even when there is underlying bias present in the training sets used in the supervised machine learning process. We use catalog simulations to validate our bias metric, and show how to bin the multidimensional intrinsic and observed galaxy properties used in the bias quantification. Our approach is designed to work on any other labeled multidimensional data sets and the code is publicly available.


title: The ASCCR Frame for Learning Essential Collaboration Skills

abstract: Statistics and data science are especially collaborative disciplines that typically require practitioners to interact with many different people or groups. Consequently, interdisciplinary collaboration skills are part of the personal and professional skills essential for success as an applied statistician or data scientist. These skills are learnable and teachable, and learning and improving collaboration skills provides a way to enhance one's practice of statistics and data science. To help individuals learn these skills and organizations to teach them, we have developed a framework covering five essential components of statistical collaboration: Attitude, Structure, Content, Communication, and Relationship. We call this the ASCCR Frame. This framework can be incorporated into formal training programs in the classroom or on the job and can also be used by individuals through self-study. We show how this framework can be applied specifically to statisticians and data scientists to improve their collaboration skills and their interdisciplinary impact. We believe that the ASCCR Frame can help organize and stimulate research and teaching in interdisciplinary collaboration and call on individuals and organizations to begin generating evidence regarding its effectiveness.


title: Mechanism Design with Limited Commitment

abstract: We develop a tool akin to the revelation principle for dynamic mechanism-selection games in which the designer can only commit to short-term mechanisms. We identify a canonical class of mechanisms rich enough to replicate the outcomes of any equilibrium in a mechanism-selection game between an uninformed designer and a privately informed agent. A cornerstone of our methodology is the idea that a mechanism should encode not only the rules that determine the allocation, but also the information the designer obtains from the interaction with the agent. Therefore, how much the designer learns, which is the key tension in design with limited commitment, becomes an explicit part of the design. Our result simplifies the search for the designer-optimal outcome by reducing the agent's behavior to a series of participation, truthtelling, and Bayes' plausibility constraints the mechanisms must satisfy.


title: $J/\psi$ -nucleon scattering in $P_{c}^{+}$ pentaquarks channel

abstract: Two pentaquarks $P_{c}^{+}$ were discovered by LHCb collaboration as peaks in the proton-$J/\psi$ invariant mass. We perform the lattice QCD study of the scattering between $J/\psi$ meson and nucleon in the channels with $J^{P}=\frac{3}{2}^{+},\frac{3}{2}^{-}, \frac{5}{2}^{+}, \frac{5}{2}^{-}$, where $P_{c}^{+}$ was discovered. This is the first lattice simulation that reaches the energies $4.3-4.5~$GeV where pentaquarks reside. The higher partial waves $L>0$ are also explored for the first time. In this study we consider the single-channel approximation for scattering of $NJ/\psi$. Energies and eigenstates are extracted for the $NJ/\psi$ system at the zero total momentum for all six irreducible representations of the lattice irreducible representation. No significant energy shifts are observed. The number of eigenstates agrees with the number expected from non-interacting limit for scattering. This could possibly indicate that the $P_{c}$ resonances seen in experiment are a consequence of a coupling of the $NJ/\psi$ channel with other two-hadron channels.


title: Decidability in Robot Manipulation Planning

abstract: Consider the problem of planning collision-free motion of $n$ objects in the plane movable through contact with a robot that can autonomously translate in the plane and that can move a maximum of $m \leq n$ objects simultaneously. This represents the abstract formulation of a manipulation planning problem that is proven to be decidable in this paper. The tools used for proving decidability of this simplified manipulation planning problem are, in fact, general enough to handle the decidability problem for the wider class of systems characterized by a stratified configuration space. These include, for example, problems of legged and multi-contact locomotion, bi-manual manipulation. In addition, the described approach does not restrict the dynamics of the manipulation system to be considered.


title: Comment on the article "Anisotropies in the astrophysical gravitational-wave background: The impact of black hole distributions" by A.C. Jenkins et al. [arXiv:1810.13435]

abstract: We investigate the discrepancy pointed out by Jenkins et al. in Ref. [1] between the predictions of anisotropies of the astrophysical gravitational wave (GW) background, derived using different methods in Cusin et al. [2] and in Jenkins et al. [3]. We show that this discrepancy is not due to our treatment of galaxy clustering, contrary to the claim made in Ref. [1] and we show that our modeling of clustering gives results in very good agreement with observations. Furthermore we show that the power law spectrum used in Refs. [1] and [3] to describe galaxy clustering is incorrect on large scales and leads to a different scaling of the multipoles $C_\ell$. Moreover, we also explain that the analytic derivation of the gravitational wave background correlation function in Refs. [1] and [3] is mathematically ill-defined and predicts an amplitude of the angular power spectrum which depends on the (arbitrary) choice of a non-physical cut-off.


title: The low-energy TQFT of the generalized double semion model

abstract: The generalized double semion (GDS) model, introduced by Freedman and Hastings, is a lattice system similar to the toric code, with a gapped Hamiltonian whose definition depends on a triangulation of the ambient manifold $M$, but whose space of ground states does not depend on the triangulation, but only on the underlying manifold. In this paper, we use topological quantum field theory (TQFT) to investigate the low-energy limit of the GDS model. We define and study a functorial TQFT $Z_{\mathrm{GDS}}$ in every dimension $n$ such that for every closed $(n - 1)$-manifold $M$, $Z_{\mathrm{GDS}}(M)$ is isomorphic to the space of ground states of the GDS model on $M$; the isomorphism can be chosen to intertwine the actions of the mapping class group of $M$ that arise on both sides. Throughout this paper, we compare our constructions and results with their known analogues for the toric code.

categories: math-ph cond-mat.str-el hep-th math.AT math.MP


title: Thermodynamics at Strong Coupling on Anisotropic Lattices

abstract: Lattice QCD with staggered fermions at strong coupling has long been studied in a dual representation to circumvent the finite baryon density sign problem. Monte Carlo simulations at finite temperature and density require anisotropic lattices. Recent results that established the non-perturbative functional dependence between the bare anisotropy $\gamma$ and the physical anisotropy $a_s/a_t$ in the chiral limit are now extended to finite quark mass. We illustrate how the calibration of the anisotropy works and discuss the consequences of the anisotropy on thermodynamic observables. We also show first results on the energy density and pressure in the QCD phase diagram in the strong coupling regime.


title: Quantitative Robustness Analysis of Quantum Programs (Extended Version)

abstract: Quantum computation is a topic of significant recent interest, with practical advances coming from both research and industry. A major challenge in quantum programming is dealing with errors (quantum noise) during execution. Because quantum resources (e.g., qubits) are scarce, classical error correction techniques applied at the level of the architecture are currently cost-prohibitive. But while this reality means that quantum programs are almost certain to have errors, there as yet exists no principled means to reason about erroneous behavior. This paper attempts to fill this gap by developing a semantics for erroneous quantum while-programs, as well as a logic for reasoning about them. This logic permits proving a property we have identified, called $\epsilon$-robustness, which characterizes possible "distance" between an ideal program and an erroneous one. We have proved the logic sound, and showed its utility on several case studies, notably: (1) analyzing the robustness of noisy versions of the quantum Bernoulli factory (QBF) and quantum walk (QW); (2) demonstrating the (in)effectiveness of different error correction schemes on single-qubit errors; and (3) analyzing the robustness of a fault-tolerant version of QBF.


title: A version of Putinar's Positivstellensatz for cylinders

abstract: We prove that, under some additional assumption, Putinar's Positivstellensatz holds on cylinders of type $S \times {\mathbb R}$ with $S = \{x \in {\mathbb R}^n | g_1(x) \ge 0, ..., g_s(x) \ge 0\}$ such that the quadratic module generated by $g_1, ..., g_s$ in ${\mathbb R}[X_1, ..., X_n]$ is archimedean, and we provide a degree bound for the representation of a polynomial $f \in {\mathbb R}[X_1, ..., X_n, Y]$ which is positive on $S \times {\mathbb R}$ as an explicit element of the quadratic module generated by $g_1, ..., g_s$ in ${\mathbb R}[X_1, ..., X_n, Y]$. We also include an example to show that an additional assumption is necessary for Putinar's Positivstellensatz to hold on cylinders of this type.


title: Identities And Relations On The Hermite-based Tangent Polynomials

abstract: In this note, we introduce and investigate the Hermite-based Tangent numbers and polynomials, Hermite-based modifieed degenerate- Tangent polynomials, poly-Tangent polynomials. We give some identities and relations for these polynomials.


title: Constrained Information Design

abstract: We provide tools to analyze information design problems subject to constraints. We do so by extending the insight in Le Treust and Tomala (2019) to the case of multiple inequality and equality constraints. Namely, that an information design problem subject to constraints can be represented as an unconstrained information design problem with a additional states, one for each constraint. Thus, without loss of generality, optimal solutions induce as many posteriors as the number of states and constraints. We provide results that refine this upper bound. Furthermore, we provide conditions under which there is no duality gap in constrained information design, thus validating a Lagrangian approach. We illustrate our results with applications to mechanism design with limited commitment (Doval and Skreta, 2022a) and persuasion of a privately informed receiver (Kolotilin et al., 2017).


title: Cellular Structure of Wreath Product Algebras

abstract: We apply the method of iterated inflations to show that the wreath product of a cellular algebra with a symmetric group is cellular, and obtain descriptions of the cell and simple modules together with a semisimplicity condition for such wreath products.


title: Social media cluster dynamics create resilient global hate highways

abstract: Online social media allows individuals to cluster around common interests - including hate. We show that tight-knit social clusters interlink to form resilient 'global hate highways' that bridge independent social network platforms, countries, languages and ideologies, and can quickly self-repair and rewire. We provide a mathematical theory that reveals a hidden resilience in the global axis of hate; explains a likely ineffectiveness of current control methods; and offers improvements. Our results reveal new science for networks-of-networks driven by bipartite dynamics, and should apply more broadly to illicit networks.


title: Nonlinear Dimension Reduction via Outer Bi-Lipschitz Extensions

abstract: We introduce and study the notion of an outer bi-Lipschitz extension of a map between Euclidean spaces. The notion is a natural analogue of the notion of a Lipschitz extension of a Lipschitz map. We show that for every map $f$ there exists an outer bi-Lipschitz extension $f'$ whose distortion is greater than that of $f$ by at most a constant factor. This result can be seen as a counterpart of the classic Kirszbraun theorem for outer bi-Lipschitz extensions. We also study outer bi-Lipschitz extensions of near-isometric maps and show upper and lower bounds for them. Then, we present applications of our results to prioritized and terminal dimension reduction problems. * We prove a prioritized variant of the Johnson-Lindenstrauss lemma: given a set of points $X\subset \mathbb{R}^d$ of size $N$ and a permutation ("priority ranking") of $X$, there exists an embedding $f$ of $X$ into $\mathbb{R}^{O(\log N)}$ with distortion $O(\log \log N)$ such that the point of rank $j$ has only $O(\log^{3 + \varepsilon} j)$ non-zero coordinates - more specifically, all but the first $O(\log^{3+\varepsilon} j)$ coordinates are equal to $0$; the distortion of $f$ restricted to the first $j$ points (according to the ranking) is at most $O(\log\log j)$. The result makes a progress towards answering an open question by Elkin, Filtser, and Neiman about prioritized dimension reductions. * We prove that given a set $X$ of $N$ points in $\mathbb{R}^d$, there exists a terminal dimension reduction embedding of $\mathbb{R}^d$ into $\mathbb{R}^{d'}$, where $d' = O\left(\frac{\log N}{\varepsilon^4}\right)$, which preserves distances $\|x-y\|$ between points $x\in X$ and $y \in \mathbb{R}^{d}$, up to a multiplicative factor of $1 \pm \varepsilon$. This improves a recent result by Elkin, Filtser, and Neiman. The dimension reductions that we obtain are nonlinear, and this nonlinearity is necessary.


title: An O^*(2.619^k) algorithm for 4-path vertex cover

abstract: In the 4-path vertex cover problem, the input is an undirected graph $G$ and an integer $k$. The goal is to decide whether there is a set of vertices $S$ of size at most $k$ such that every path with 4 vertices in $G$ contains at least one vertex of $S$. In this paper we give a parameterized algorithm for 4-path vertex cover whose time complexity is $O^*(2.619^k)$.


title: The Borg Cube Simulation: Cosmological Hydrodynamics with CRK-SPH

abstract: A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity and the other with non-radiative hydrodynamics. Each contains 2x2304^3 cold dark matter plus baryon particles in an 800 Mpc/h box. We compare statistics to previous non-radiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of non-thermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10^15 (10^13) Msun/h. We compute angular power spectra for the thermal and kinematic Sunyaev-Zel'dovich effects and find good agreement with the low-l Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH with future work including subresolution gas treatments.


title: Nonlocal black hole evaporation and quantum metric fluctuations via inhomogeneous vacuum density

abstract: Inhomogeneity of the actual value of the vacuum energy density is considered in a black hole background. We examine the back-reaction of a Schwarzschild black hole to the highly inhomogeneous vacuum density, and argue the fluctuations lead to deviations from general relativity in the near-horizon region. In particular, we found that vacuum fluctuations onto the horizon trigger adiabatic release of quantum information, while vacuum fluctuations in the vicinity of the horizon produce potentially observable metric fluctuations of order the Schwarzschild radius. Consequently, we propose a form of strong nonviolent nonlocality in which we simultaneously get nonlocal release of quantum information, and observable metric fluctuations.


title: The ordinal generated by an ordinal grammar is computable

abstract: A prefix grammar is a context-free grammar whose nonterminals generate prefix-free languages. A prefix grammar $G$ is an ordinal grammar if the language $L(G)$ is well-ordered with respect to the lexicographic ordering. It is known that from a finite system of parametric fixed point equations one can construct an ordinal grammar $G$ such that the lexicographic order of $G$ is isomorphic with the least solution of the system, if this solution is well-ordered. In this paper we show that given an ordinal grammar, one can compute (the Cantor normal form of) the order type of the lexicographic order of its language, yielding that least solutions of fixed point equation systems defining algebraic ordinals are effectively computable (and thus, their isomorphism problem is also decidable).


title: Updated Results on the Galaxy-Halo Connection from Satellite Kinematics in SDSS

abstract: We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.


title: Holographic Complexity for disentangled states

abstract: In this paper we consider the maximal volume and the action, which are conjectured to be gravity duals of the complexity, in the black hole geometries with end of the world branes. These geometries are duals of boundary states in CFTs which have small real space entanglement. When we raise the black hole temperature while keeping the cutoff radius, black hole horizons or end of the world branes come in contact with the cutoff surface. In this limit, holographic entanglement entropy reduces to 0. We studied the behavior of the volume and the action. We found that the volume reduces to 0 in this limit. The behavior of the action depends on their regularization. We study the implication of these results to the reference state of the holographic complexity both in the complexity = volume or the complexity = action conjectures.


title: Cross-comparative analysis of evacuation behavior after earthquakes using mobile phone data

abstract: Despite the importance of predicting evacuation mobility dynamics after large scale disasters for effective first response and disaster relief, our general understanding of evacuation behavior remains limited because of the lack of empirical evidence on the evacuation movement of individuals across multiple disaster instances. Here we investigate the GPS trajectories of a total of more than 1 million anonymized mobile phone users whose positions are tracked for a period of 2 months before and after four of the major earthquakes that occurred in Japan. Through a cross comparative analysis between the four disaster instances, we find that in contrast with the assumed complexity of evacuation decision making mechanisms in crisis situations, the individuals' evacuation probability is strongly dependent on the seismic intensity that they experience. In fact, we show that the evacuation probabilities in all earthquakes collapse into a similar pattern, with a critical threshold at around seismic intensity 5.5. This indicates that despite the diversity in the earthquakes profiles and urban characteristics, evacuation behavior is similarly dependent on seismic intensity. Moreover, we found that probability density functions of the distances that individuals evacuate are not dependent on seismic intensities that individuals experience. These insights from empirical analysis on evacuation from multiple earthquake instances using large scale mobility data contributes to a deeper understanding of how people react to earthquakes, and can potentially assist decision makers to simulate and predict the number of evacuees in urban areas with little computational time and cost, by using population density information and seismic intensity which can be observed instantaneously after the shock.


title: Stabilization and Destabilization of Semiflexible Polymers Induced by Spherical Confinement

abstract: Spherical confinement can act either stabilizing or destabilizing on the collapsed state of a semi-flexible polymer. General free-energy arguments suggest that the order of the unconstrained collapse transition is the distinguishing factor: First order implies stabilization, second order causes destabilization. We confirm this conjecture by Monte Carlo simulations of a coarse-grained model for semiflexible polymers whose chain stiffness tunes the transition order. The resulting physical picture is potentially relevant for other systems under strong confinement such as proteins and DNA.


title: Measuring the Effects of Data Parallelism on Neural Network Training

abstract: Recent hardware developments have dramatically increased the scale of data parallelism available for neural network training. Among the simplest ways to harness next-generation hardware is to increase the batch size in standard mini-batch neural network training algorithms. In this work, we aim to experimentally characterize the effects of increasing the batch size on training time, as measured by the number of steps necessary to reach a goal out-of-sample error. We study how this relationship varies with the training algorithm, model, and data set, and find extremely large variation between workloads. Along the way, we show that disagreements in the literature on how batch size affects model quality can largely be explained by differences in metaparameter tuning and compute budgets at different batch sizes. We find no evidence that larger batch sizes degrade out-of-sample performance. Finally, we discuss the implications of our results on efforts to train neural networks much faster in the future. Our experimental data is publicly available as a database of 71,638,836 loss measurements taken over the course of training for 168,160 individual models across 35 workloads.


title: Deep BV: A Fully Automated System for Brain Ventricle Localization and Segmentation in 3D Ultrasound Images of Embryonic Mice

abstract: Volumetric analysis of brain ventricle (BV) structure is a key tool in the study of central nervous system development in embryonic mice. High-frequency ultrasound (HFU) is the only non-invasive, real-time modality available for rapid volumetric imaging of embryos in utero. However, manual segmentation of the BV from HFU volumes is tedious, time-consuming, and requires specialized expertise. In this paper, we propose a novel deep learning based BV segmentation system for whole-body HFU images of mouse embryos. Our fully automated system consists of two modules: localization and segmentation. It first applies a volumetric convolutional neural network on a 3D sliding window over the entire volume to identify a 3D bounding box containing the entire BV. It then employs a fully convolutional network to segment the detected bounding box into BV and background. The system achieves a Dice Similarity Coefficient (DSC) of 0.8956 for BV segmentation on an unseen 111 HFU volume test set surpassing the previous state-of-the-art method (DSC of 0.7119) by a margin of 25%.


title: Unraveling the structure of treelike networks from first-passage times of lazy random walkers

abstract: We study the problem of random search in finite networks with a tree topology, where it is expected that the distribution of the first-passage time F(t) decays exponentially. We show that the slope of the exponential tail is independent of the initial conditions of entering the tree in general, and scales exponentially or as a power law with the extent of the tree L, depending on the tendency p to jump toward the target node. It is unfeasible to uniquely determine L and p from measuring the tail slope or the mean first-passage time (MFPT) of an ordinary diffusion along the tree. To unravel the structure, we consider lazy random walkers that take steps with probability m when jumping on the nodes and return with probability q from the leaves. By deriving an exact analytical expression for the MFPT of the intermittent random walk, we verify that the structural information of the tree can be uniquely extracted by measuring the MFPT for two randomly chosen types of tracer particles with distinct experimental parameters m and q. We also address the applicability of our approach in the presence of disorder in the structure of the tree or statistical uncertainty in the experimental parameters.


title: Calculated magnetic exchange interactions in Dirac magnon material Cu3TeO6

abstract: Recently topological aspects of magnon band structure have attracted much interest, and especially, the Dirac magnons in Cu3TeO6 have been observed experimentally. In this work, we calculate the magnetic exchange interactions J's using the first-principles linear-response approach and find that these J's are short-range and negligible for the Cu-Cu atomic pair apart by longer than 7 Angstrom. Moreover there are only 5 sizable magnetic exchange interactions, and according to their signs and strengths, modest magnetic frustration is expected. Based on the obtained magnetic exchange couplings, we successfully reproduce the experimental spin-wave dispersions. The calculated neutron scattering cross section also agrees very well with the experiments. We also calculate Dzyaloshinskii-Moriya interactions (DMIs) and estimate the canting angle (about 1.3{\deg}) of the magnetic non-collinearity based on the competition between DMIs and J's, which is consistent with the experiment. The small canting angle agrees with that the current experiments cannot distinguish the DMI induced nodal line from a Dirac point in the spin-wave spectrum. Finally we analytically prove that the "sum rule" conjectured in [Nat. Phys. 14, 1011 (2018)] holds but only up to the 11th nearest neighbour.


title: Federated Learning for Mobile Keyboard Prediction

abstract: We train a recurrent neural network language model using a distributed, on-device learning framework called federated learning for the purpose of next-word prediction in a virtual keyboard for smartphones. Server-based training using stochastic gradient descent is compared with training on client devices using the Federated Averaging algorithm. The federated algorithm, which enables training on a higher-quality dataset for this use case, is shown to achieve better prediction recall. This work demonstrates the feasibility and benefit of training language models on client devices without exporting sensitive user data to servers. The federated learning environment gives users greater control over the use of their data and simplifies the task of incorporating privacy by default with distributed training and aggregation across a population of client devices.


title: Towards Lefschetz thimbles regularization of heavy-dense QCD

abstract: At finite density, lattice simulations are hindered by the well-known sign problem: for finite chemical potentials, the QCD action becomes complex and the Boltzmann weight $e^{-S}$ cannot be interpreted as a probability distribution to determine expectation values by Monte Carlo techniques. Different workarounds have been devised to study the QCD phase diagram, but their application is mostly limited to the region of small chemical potentials. The Lefschetz thimbles method takes a new approach in which one complexifies the theory and deforms the integration paths. By integrating over Lefschetz thimbles, the imaginary part of the action is kept constant and can be factored out, while $e^{-Re(S)}$ can be interpreted as a probability measure. The method has been applied in recent years to more or less difficult problems. Here we report preliminary results on Lefschetz thimbles regularization of heavy-dense QCD. While still simple, this is a very interesting problem. It is a first look at thimbles for QCD, although in a simplified, effective version. From an algorithmic point of view, it is a nice ground to test effectiveness of techniques we developed for multi thimbles simulations.


title: Bisimilarity of Open Terms in Stream GSOS

abstract: Stream GSOS is a specification format for operations and calculi on infinite sequences. The notion of bisimilarity provides a canonical proof technique for equivalence of closed terms in such specifications. In this paper, we focus on open terms, which may contain variables, and which are equivalent whenever they denote the same stream for every possible instantiation of the variables. Our main contribution is to capture equivalence of open terms as bisimilarity on certain Mealy machines, providing a concrete proof technique. Moreover, we introduce an enhancement of this technique, called bisimulation up-to substitutions, and show how to combine it with other up-to techniques to obtain a powerful method for proving equivalence of open terms.


title: Compact numerical solutions to the two-dimensional repulsive Hubbard model obtained via non-unitary similarity transformations

abstract: Similarity transformation of the Hubbard Hamiltonian using a Gutzwiller correlator leads to a non-Hermitian effective Hamiltonian, which can be expressed exactly in momentum-space representation, and contains three-body interactions. We apply this methodology to study the two-dimensional Hubbard model with repulsive interactions near half-filling in the intermediate interaction strength regime ($U/t=4$). We show that at optimal or near optimal strength of the Gutzwiller correlator, the similarity transformed Hamiltonian has extremely compact right eigenvectors, which can be sampled to high accuracy using the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) method, and its initiator approximation. Near-optimal correlators can be obtained using a simple projective equation, thus obviating the need for a numerical optimisation of the correlator. The FCIQMC method, as a projective technique, is well-suited for such non-Hermitian problems, and its stochastic nature can handle the 3-body interactions exactly without undue increase in computational cost. The highly compact nature of the right eigenvectors means that the initiator approximation in FCIQMC is not severe, and that large lattices can be simulated, well beyond the reach of the method applied to the original Hubbard Hamiltonian. Results are provided in lattice sizes upto 50 sites and compared to auxiliary-field QMC. New benchmark results are provided in the off half-filling regime, with no severe sign-problem being encountered. In addition, we show that methodology can be used to calculate excited states of the Hubbard model and lay the groundwork for the calculation of observables other than the energy.


title: Homeopathic Dark Matter, or how diluted heavy substances produce high energy cosmic rays

abstract: We point out that current and planned telescopes have the potential of probing annihilating Dark Matter (DM) with a mass of O(100) TeV and beyond. As a target for such searches, we propose models where DM annihilates into lighter mediators, themselves decaying into Standard Model (SM) particles. These models allow to reliably compute the energy spectra of the SM final states, and to naturally evade the unitarity bound on the DM mass. Indeed, long-lived mediators may cause an early matter-dominated phase in the evolution of the Universe and, upon decaying, dilute the density of preexisting relics thus allowing for very large DM masses. We compute this dilution in detail and provide results in a ready-to-use form. Considering for concreteness a model of dark U(1) DM, we then study both dilution and the signals at various high energy telescopes observing gamma rays, neutrinos and charged cosmic rays. This study enriches the physics case of these experiments, and opens a new observational window on heavy new physics sectors.


title: Symplectic cohomology rings of affine varieties in the topological limit

abstract: We construct a multiplicative spectral sequence converging to the symplectic cohomology ring of any affine variety $X$, with first page built out of topological invariants associated to strata of any fixed normal crossings compactification $(M,\mathbf{D})$ of $X$. We exhibit a broad class of pairs $(M,\mathbf{D})$ (characterized by the absence of relative holomorphic spheres or vanishing of certain relative GW invariants) for which the spectral sequence degenerates, and a broad subclass of pairs (similarly characterized) for which the ring structure on symplectic cohomology can also be described topologically. Sample applications include: (a) a complete topological description of the symplectic cohomology ring of the complement, in any projective $M$, of the union of sufficiently many generic ample divisors whose homology classes span a rank one subspace, (b) complete additive and partial multiplicative computations of degree zero symplectic cohomology rings of many log Calabi-Yau varieties, and (c) a proof in many cases that symplectic cohomology is finitely generated as a ring. A key technical ingredient in our results is a logarithmic version of the PSS morphism, introduced in our earlier work [GP1].


title: On the Filamentation Instability in Degenerate Relativistic Plasmas

abstract: The filamentation instability of the electromagnetic (EM) beam in an underdense plasma with high level of degeneracy is examined by means of the momentum equation, continuity equation and Maxwell's equations. It has been demonstrated that the instability develops for weakly as well as strongly relativistic degenerate plasma and arbitrary strong amplitude of EM beams.


title: Evolution of Hubble wedges in episodic protostellar outflows

abstract: Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin of these features, we have developed a new episodic outflow model for the SPH code {\sc gandalf}, which mimics the accretion and ejection behaviour of FU Ori type stars. We apply this model to simulations of star formation, invoking two types of initial conditions: spherically symmetric cores in solid-body rotation with $\rho\propto r^{-2}$, and spherically symmetric turbulent cores with density proportional to the density of a Bonnor-Ebert sphere. For a wide range of model parameters, we find that episodic outflows lead to self-regulation of the ejected mass and momentum, and we achieve acceptable results, even with relatively low resolution. Using this model, we find that recently ejected outflow bullets produce a `Hubble wedge' in the position-velocity relation. However, once such a bullet hits the leading shock front, it decelerates and aligns with older bullets to form a `Hubble-law'. Bullets can be identified as bumps in the mass-velocity relation, which can be fit with a power-law, $dM/d\upsilon_{_{\rm RAD}}\propto\upsilon_{_{\rm RAD}}^{-1.5}$.


title: Communication and Information Theory of Single Action Potential Signals in Plants

abstract: Many plants, such as Mimosa pudica (the sensitive plant), employ electrochemical signals known as action potentials (APs) for rapid intercellular communication. In this paper, we consider a reaction diffusion model of individual AP signals to analyze APs from a communication and information theoretic perspective. We use concepts from molecular communication to explain the underlying process of information transfer in a plant for a single AP pulse that is shared with one or more receiver cells. We also use the chemical Langevin equation to accommodate the deterministic as well as stochastic component of the system. Finally we present an information theoretic analysis of single action potentials, obtaining achievable information rates for these signals. We show that, in general, the presence of an AP signal can increase the mutual information and information propagation speed among neighboring cells with receivers in different settings.


title: The Transition Function of $G_2$ over $S^6$

abstract: We obtain explicit formulas for the trivialization functions of the ${\rm SU}(3)$ principal bundle $G_2 \to S^6$ over two affine charts. We also calculate the explicit transition function of this fibration over the equator of the six-sphere. In this way we obtain a new proof of the known fact that this fibration corresponds to a generator of $\pi_{5}({\rm SU}(3))$.


title: New limits on Anomalous Spin-Spin Interactions

abstract: We report the results of a new search for long range spin-dependent interactions using a Rb -$^{21}$Ne atomic comagnetometer and a rotatable electron spin source based on a SmCo$_{5}$ magnet with an iron flux return. By looking for signal correlations with the orientation of the spin source we set new constrains on the product of the pseudoscalar electron and neutron couplings $g^e_p g^n_p/\hbar c<1.7\times10^{-14}$ and on the product of their axial couplings $g^e_A g^n_A/\hbar c<5\times10^{-42}$ to a new particle with a mass of less than about $1~\mu$eV. Our measurements improve by about 2 orders of magnitude previous constraints on such spin-dependent interactions.


title: BRST invariant RG flows

abstract: A mass parameter for the gauge bosons in gauge-fixed four-dimensional Yang-Mills theory can be accommodated in a local and manifestly BRST-invariant action. The construction is based on the Faddeev-Popov method involving a nonlinear gauge-fixing and a background Nakanishi-Lautrup field. When applied to momentum-dependent masslike deformations, this formalism leads to a full regularization of the theory which explicitly preserves BRST symmetry. We deduce a functional renormalization group equation for the one-particle-irreducible effective action, which has a one-loop form. The master equation is compatible with it, i.e. BRST symmetry is preserved along the flow, and it has a standard regulator-independent Zinn-Justin form. As a first application, we compute the leading-order gluon wave-function renormalization.


title: Meter- to Millimeter Emission from Cool Stellar Systems: Latest Results, Synergies Across the Spectrum, and Outlook for the Next Decade

abstract: Radio observations of cool stellar systems provide unique information on their magnetic fields, high-energy processes, and chemistry. Buoyed by powerful new instruments (e.g. ALMA, JVLA, LOFAR), advances in related fields (e.g., the Gaia astrometric revolution), and above all a renewed interest in the relevant stellar astrophysics, stellar radio astronomy is experiencing a renaissance. In this splinter session, participants took stock of the present state of stellar radio astronomy to chart a course for the field's future.


title: GradiVeQ: Vector Quantization for Bandwidth-Efficient Gradient Aggregation in Distributed CNN Training

abstract: Data parallelism can boost the training speed of convolutional neural networks (CNN), but could suffer from significant communication costs caused by gradient aggregation. To alleviate this problem, several scalar quantization techniques have been developed to compress the gradients. But these techniques could perform poorly when used together with decentralized aggregation protocols like ring all-reduce (RAR), mainly due to their inability to directly aggregate compressed gradients. In this paper, we empirically demonstrate the strong linear correlations between CNN gradients, and propose a gradient vector quantization technique, named GradiVeQ, to exploit these correlations through principal component analysis (PCA) for substantial gradient dimension reduction. GradiVeQ enables direct aggregation of compressed gradients, hence allows us to build a distributed learning system that parallelizes GradiVeQ gradient compression and RAR communications. Extensive experiments on popular CNNs demonstrate that applying GradiVeQ slashes the wall-clock gradient aggregation time of the original RAR by more than 5X without noticeable accuracy loss, and reduces the end-to-end training time by almost 50%. The results also show that GradiVeQ is compatible with scalar quantization techniques such as QSGD (Quantized SGD), and achieves a much higher speed-up gain under the same compression ratio.


title: Demonstrating Advantages of Neuromorphic Computation: A Pilot Study

abstract: Neuromorphic devices represent an attempt to mimic aspects of the brain's architecture and dynamics with the aim of replicating its hallmark functional capabilities in terms of computational power, robust learning and energy efficiency. We employ a single-chip prototype of the BrainScaleS 2 neuromorphic system to implement a proof-of-concept demonstration of reward-modulated spike-timing-dependent plasticity in a spiking network that learns to play the Pong video game by smooth pursuit. This system combines an electronic mixed-signal substrate for emulating neuron and synapse dynamics with an embedded digital processor for on-chip learning, which in this work also serves to simulate the virtual environment and learning agent. The analog emulation of neuronal membrane dynamics enables a 1000-fold acceleration with respect to biological real-time, with the entire chip operating on a power budget of 57mW. Compared to an equivalent simulation using state-of-the-art software, the on-chip emulation is at least one order of magnitude faster and three orders of magnitude more energy-efficient. We demonstrate how on-chip learning can mitigate the effects of fixed-pattern noise, which is unavoidable in analog substrates, while making use of temporal variability for action exploration. Learning compensates imperfections of the physical substrate, as manifested in neuronal parameter variability, by adapting synaptic weights to match respective excitability of individual neurons.


title: Pipe-SGD: A Decentralized Pipelined SGD Framework for Distributed Deep Net Training

abstract: Distributed training of deep nets is an important technique to address some of the present day computing challenges like memory consumption and computational demands. Classical distributed approaches, synchronous or asynchronous, are based on the parameter server architecture, i.e., worker nodes compute gradients which are communicated to the parameter server while updated parameters are returned. Recently, distributed training with AllReduce operations gained popularity as well. While many of those operations seem appealing, little is reported about wall-clock training time improvements. In this paper, we carefully analyze the AllReduce based setup, propose timing models which include network latency, bandwidth, cluster size and compute time, and demonstrate that a pipelined training with a width of two combines the best of both synchronous and asynchronous training. Specifically, for a setup consisting of a four-node GPU cluster we show wall-clock time training improvements of up to 5.4x compared to conventional approaches.


title: Diffusion Kurtosis Imaging maps neural damage in the EAE model of multiple sclerosis

abstract: Diffusion kurtosis imaging (DKI), is an imaging modality that yields novel disease biomarkers and in combination with nervous tissue modeling, provides access to microstructural parameters. Recently, DKI and subsequent estimation of microstructural model parameters has been used for assessment of tissue changes in neurodegenerative diseases and their animal models. In this study, mouse spinal cords from the experimental autoimmune encephalomyelitis (EAE) model of multiple sclerosis (MS) were investigated for the first time using DKI in combination with biophysical modeling to study the relationship between microstructural metrics and degree of animal dysfunction. Thirteen spinal cords were extracted from animals of variable disability and scanned in a high-field MRI scanner along with five control specimen. Diffusion weighted data were acquired together with high resolution T2* images. Diffusion data were fit to estimate diffusion and kurtosis tensors and white matter modeling parameters, which were all used for subsequent statistical analysis using a linear mixed effects model. T2* images were used to delineate focal demyelination/inflammation. Our results unveil a strong relationship between disability and measured microstructural parameters in normal appearing white matter and gray matter. The changes we found in biophysical modeling parameters and in particular in extra-axonal axial diffusivity were clearly different from previous studies employing other animal models of MS. In conclusion, our data suggest that DKI and microstructural modeling can provide a unique contrast capable of detecting EAE-specific changes correlating with clinical disability. These findings could close the gap between MRI findings and clinical presentation in patients and deepen our understanding of EAE and the MS mechanisms.