title: Achieving Secure and Differentially Private Computations in Multiparty Settings

abstract: Sharing and working on sensitive data in distributed settings from healthcare to finance is a major challenge due to security and privacy concerns. Secure multiparty computation (SMC) is a viable panacea for this, allowing distributed parties to make computations while the parties learn nothing about their data, but the final result. Although SMC is instrumental in such distributed settings, it does not provide any guarantees not to leak any information about individuals to adversaries. Differential privacy (DP) can be utilized to address this; however, achieving SMC with DP is not a trivial task, either. In this paper, we propose a novel Secure Multiparty Distributed Differentially Private (SM-DDP) protocol to achieve secure and private computations in a multiparty environment. Specifically, with our protocol, we simultaneously achieve SMC and DP in distributed settings focusing on linear regression on horizontally distributed data. That is, parties do not see each others' data and further, can not infer information about individuals from the final constructed statistical model. Any statistical model function that allows independent calculation of local statistics can be computed through our protocol. The protocol implements homomorphic encryption for SMC and functional mechanism for DP to achieve the desired security and privacy guarantees. In this work, we first introduce the theoretical foundation for the SM-DDP protocol and then evaluate its efficacy and performance on two different datasets. Our results show that one can achieve individual-level privacy through the proposed protocol with distributed DP, which is independently applied by each party in a distributed fashion. Moreover, our results also show that the SM-DDP protocol incurs minimal computational overhead, is scalable, and provides security and privacy guarantees.


title: Solutions of Nonlinear Polyharmonic Equation with Periodic Potential

abstract: Quasi-periodic solutions of a nonlinear periodic polyharmonic equation in $\R^n$, $n>1$, are studied. It is proven that there is an extensive "non-resonant" set ${\mathcal G}\subset \R^n$ such that for every $\vec k\in \mathcal G$ there is a solution asymptotically close to a plane wave $Ae^{i\langle{ \vec{k}, \vec{x} }\rangle}$ as $|\vec k|\to \infty $.


title: Blockchain Consensus Protocols in the Wild

abstract: A blockchain is a distributed ledger for recording transactions, maintained by many nodes without central authority through a distributed cryptographic protocol. All nodes validate the information to be appended to the blockchain, and a consensus protocol ensures that the nodes agree on a unique order in which entries are appended. Consensus protocols for tolerating Byzantine faults have received renewed attention because they also address blockchain systems. This work discusses the process of assessing and gaining confidence in the resilience of a consensus protocols exposed to faults and adversarial nodes. We advocate to follow the established practice in cryptography and computer security, relying on public reviews, detailed models, and formal proofs; the designers of several practical systems appear to be unaware of this. Moreover, we review the consensus protocols in some prominent permissioned blockchain platforms with respect to their fault models and resilience against attacks. The protocol comparison covers Hyperledger Fabric, Tendermint, Symbiont, R3~Corda, Iroha, Kadena, Chain, Quorum, MultiChain, Sawtooth Lake, Ripple, Stellar, and IOTA.


title: Maximizing the mean subtree order

abstract: This article focuses on properties and structures of trees with maximum mean subtree order in a given family; such trees are called optimal in the family. Our main goal is to describe the structure of optimal trees in $\mathcal{T}_n$ and $\mathcal{C}_n$, the families of all trees and caterpillars, respectively, of order $n$. We begin by establishing a powerful tool called the Gluing Lemma, which is used to prove several of our main results. In particular, we show that if $T$ is an optimal tree in $\mathcal{T}_n$ or $\mathcal{C}_n$ for $n\geq 4$, then every leaf of $T$ is adjacent to a vertex of degree at least $3$. We also use the Gluing Lemma to answer an open question of Jamison, and to provide a conceptually simple proof of Jamison's result that the path $P_n$ has minimum mean subtree order among all trees of order $n$. We prove that if $T$ is optimal in $\mathcal{T}_n$, then the number of leaves in $T$ is $\mathrm{O}(\log_2 n)$, and that if $T$ is optimal in $\mathcal{C}_n$, then the number of leaves in $T$ is $\mathrm{\Theta}(\log_2 n)$. Along the way, we describe the asymptotic structure of optimal trees in several narrower families of trees.


title: Calibrated Fairness in Bandits

abstract: We study fairness within the stochastic, \emph{multi-armed bandit} (MAB) decision making framework. We adapt the fairness framework of "treating similar individuals similarly" to this setting. Here, an `individual' corresponds to an arm and two arms are `similar' if they have a similar quality distribution. First, we adopt a {\em smoothness constraint} that if two arms have a similar quality distribution then the probability of selecting each arm should be similar. In addition, we define the {\em fairness regret}, which corresponds to the degree to which an algorithm is not calibrated, where perfect calibration requires that the probability of selecting an arm is equal to the probability with which the arm has the best quality realization. We show that a variation on Thompson sampling satisfies smooth fairness for total variation distance, and give an $\tilde{O}((kT)^{2/3})$ bound on fairness regret. This complements prior work, which protects an on-average better arm from being less favored. We also explain how to extend our algorithm to the dueling bandit setting.


title: Extended Hubbard model for mesoscopic transport in donor arrays in silicon

abstract: Arrays of dopants in silicon are promising platforms for the quantum simulation of the Fermi-Hubbard model. We show that the simplest model with only on-site interaction is insufficient to describe the physics of an array of phosphorous donors in silicon due to the strong intersite interaction in the system. We also study the resonant tunneling transport in the array at low temperature as a mean of probing the features of the Hubbard physics, such as the Hubbard bands and the Mott gap. Two mechanisms of localization which suppresses transport in the array are investigated: The first arises from the electron-ion core attraction and is significant at low filling; the second is due to the sharp oscillation in the tunnel coupling caused by the intervalley interference of the donor electron's wavefunction. This disorder in the tunnel coupling leads to a steep exponential decay of conductance with channel length in one-dimensional arrays, but its effect is less prominent in two-dimensional ones. Hence, it is possible to observe resonant tunneling transport in a relatively large array in two dimensions.


title: Changing Views on Curves and Surfaces

abstract: Visual events in computer vision are studied from the perspective of algebraic geometry. Given a sufficiently general curve or surface in 3-space, we consider the image or contour curve that arises by projecting from a viewpoint. Qualitative changes in that curve occur when the viewpoint crosses the visual event surface. We examine the components of this ruled surface, and observe that these coincide with the iterated singular loci of the coisotropic hypersurfaces associated with the original curve or surface. We derive formulas, due to Salmon and Petitjean, for the degrees of these surfaces, and show how to compute exact representations for all visual event surfaces using algebraic methods.


title: New Cameron-Liebler line classes with parameter $\frac{q^2+1}{2}$

abstract: New families of Cameron-Liebler line classes of ${\rm PG}(3,q)$, $q\ge 7$ odd, with parameter $(q^2+1)/2$ are constructed.


title: Internet of Things: Survey on Security and Privacy

abstract: The Internet of Things (IoT) is intended for ubiquitous connectivity among different entities or "things". While its purpose is to provide effective and efficient solutions, security of the devices and network is a challenging issue. The number of devices connected along with the ad-hoc nature of the system further exacerbates the situation. Therefore, security and privacy has emerged as a significant challenge for the IoT. In this paper,we aim to provide a thorough survey related to the privacy and security challenges of the IoT. This document addresses these challenges from the perspective of technologies and architecture used. This work focuses also in IoT intrinsic vulnerabilities as well as the security challenges of various layers based on the security principles of data confidentiality, integrity and availability. This survey analyzes articles published for the IoT at the time and relates it to the security conjuncture of the field and its projection to the future.


title: Project Makespan Estimation: Computational Load of Interval and Point Estimates

abstract: The estimation of project completion time is to be repeated several times in the project planning phase to reach the optimal tradeoff between time, cost, and quality. Estimation procedures provide either an interval or a point estimate. The computational load of several estimation procedures is reviewed. A multiple polynomial regression model is provided for major interval estimation procedures and shows that the accuracy in the probability model for activities is the most influential factor. The computational time does not appear to be an impeding factor, though it is larger for MonteCarlo simulation, so that the computational time can be traded off in search of a simpler estimation procedure.


title: Approximate degeneracy of $J=1$ spatial correlators in high temperature QCD

abstract: We study spatial isovector meson correlators in $N_f=2$ QCD with dynamical domain-wall fermions on $32^3\times 8$ lattices at temperatures $T=220-380$ MeV. We measure the correlators of spin-one ($J=1$) operators including vector, axial-vector, tensor and axial-tensor. Restoration of chiral $U(1)_A$ and $SU(2)_L \times SU(2)_R$ symmetries of QCD implies degeneracies in vector--axial-vector ($SU(2)_L \times SU(2)_R$) and tensor--axial-tensor ($U(1)_A$) pairs, which are indeed observed at temperatures above $T_c$. Moreover, we observe an approximate degeneracy of all $J=1$ correlators with increasing temperature. This approximate degeneracy suggests emergent $SU(2)_{CS}$ and $SU(4)$ symmeries at high temperatures, that mix left- and right-handed quarks.


title: A contemporary look at Hermann Hankel's 1861 pioneering work on Lagrangian fluid dynamics

abstract: The present paper is a companion to the paper by Villone and Rampf (2017), titled "Hermann Hankel's On the general theory of motion of fluids, an essay including an English translation of the complete Preisschrift from 1861" together with connected documents. Here we give a critical assessment of Hankel's work, which covers many important aspects of fluid dynamics considered from a Lagrangian-coordinates point of view: variational formulation in the spirit of Hamilton for elastic (barotropic) fluids, transport (we would now say Lie transport) of vorticity, the Lagrangian significance of Clebsch variables, etc. Hankel's work is also put in the perspective of previous and future work. Hence, the action spans about two centuries: from Lagrange's 1760-1761 Turin paper on variational approaches to mechanics and fluid mechanics problems to Arnold's 1966 founding paper on the geometrical/variational formulation of incompressible flow. The 22-year old Hankel - who was to die 12 years later - emerges as a highly innovative master of mathematical fluid dynamics, fully deserving Riemann's assessment that his Preisschrift contains "all manner of good things."


title: Hermann Hankel's "On the general theory of motion of fluids", an essay including an English translation of the complete Preisschrift from 1861

abstract: The present is a companion paper to "A contemporary look at Hermann Hankel's 1861 pioneering work on Lagrangian fluid dynamics" by Frisch, Grimberg and Villone (2017). Here we present the English translation of the 1861 prize manuscript from G\"ottingen University "Zur allgemeinen Theorie der Bewegung der Fl\"ussigkeiten" (On the general theory of the motion of the fluids) of Hermann Hankel (1839-1873), which was originally submitted in Latin and then translated into German by the Author for publication. We also provide the English translation of two important reports on the manuscript, one written by Bernhard Riemann and the other by Wilhelm Eduard Weber, during the assessment process for the prize. Finally we give a short biography of Hermann Hankel with his complete bibliography.


title: Pointwise estimate for the Bergman kernel of the weighted Bergman spaces with exponential type weights

abstract: Let $AL^{2}_{\phi}(\mathbb{D})$ denote the closed subspace of $L^{2}(\mathbb{D},e^{-2\phi}d\lambda)$ consisting of holomorphic functions in the unit disc ${\mathbb D}$. For certain class of subharmonic funcions $\phi : {\mathbb D}\rightarrow{\mathbb D}$, we prove upper pointwise estimate for the Bergman kernel for $AL^{2}_{\phi}(\mathbb{D})$


title: Regimes of the Pomeron and its Intrinsic Entropy

abstract: We suggest that the perturbative and non-perturbative descriptions of the Pomeron can be viewed as complementary descriptions of different phases in the Pomeron phase diagram, with a phase boundary where the proper description of the produced systems are "string balls". Their intrinsic entropy is calculated and turned out to be the same, as the recently reported perturbative entanglement entropy. The distribution of large multiplicities stemming from the string balls is also wide, with its moments close to those reported for hadrons in $pp$ collisions at the LHC. At low-x, the quantum string is so entangled that sufficiently weak string self-attraction can cause it to turn to a string ball dual to a black hole. We suggest that low-x saturation occurs when the density of wee-strings reaches the Bekenstein bound, with a proton size that freezes with increasing rapidity. Some of these observations maybe checked at the future eIC.


title: Buildup of Speaking Skills in an Online Learning Community: A Network-Analytic Exploration

abstract: In this study, we explore peer-interaction effects in online networks on speaking skill development. In particular, we present an evidence for gradual buildup of skills in a small-group setting that has not been reported in the literature. We introduce a novel dataset of six online communities consisting of 158 participants focusing on improving their speaking skills. They video-record speeches for 5 prompts in 10 days and exchange comments and performance-ratings with their peers. We ask (i) whether the participants' ratings are affected by their interaction patterns with peers, and (ii) whether there is any gradual buildup of speaking skills in the communities towards homogeneity. To analyze the data, we employ tools from the emerging field of Graph Signal Processing (GSP). GSP enjoys a distinction from Social Network Analysis in that the latter is concerned primarily with the connection structures of graphs, while the former studies signals on top of graphs. We study the performance ratings of the participants as graph signals atop underlying interaction topologies. Total variation analysis of the graph signals show that the participants' rating differences decrease with time (slope=-0.04, p<0.01), while average ratings increase (slope=0.07, p<0.05)--thereby gradually building up the ratings towards community-wide homogeneity. We provide evidence for peer-influence through a prediction formulation. Our consensus-based prediction model outperforms baseline network-agnostic regression models by about 23% in predicting performance ratings. This, in turn, shows that participants' ratings are affected by their peers' ratings and the associated interaction patterns, corroborating previous findings. Then, we formulate a consensus-based diffusion model that captures these observations of peer-influence from our analyses.


title: Effect of density on the physical aging of pressure-densified polymethylmethacrylate

abstract: The rate of physical aging of glassy polymethylmethacrylate (PMMA), followed from the change in the secondary relaxation with aging, is found to be independent of the density, the latter controlled by the pressure during glass formation. Thus, the aging behavior of the secondary relaxation is the same whether the glass is more compacted or less dense than the corresponding equilibrium liquid. This equivalence in aging of glasses formed under different pressures indicates that local packing is the dominant variable governing the glassy dynamics. One consequence is that pressure densification yields a reduction in the glass transition temperature. The fact that pressure densification yields different glass structures is at odds with a model for non-associated materials having dynamic properties exhibited by PMMA, such as density scaling of the relaxation time and isochronal superposition of the relaxation dispersion.


title: Batch Informed Trees (BIT*): Informed Asymptotically Optimal Anytime Search

abstract: Path planning in robotics often requires finding high-quality solutions to continuously valued and/or high-dimensional problems. These problems are challenging and most planning algorithms instead solve simplified approximations. Popular approximations include graphs and random samples, as respectively used by informed graph-based searches and anytime sampling-based planners. Informed graph-based searches, such as A*, traditionally use heuristics to search a priori graphs in order of potential solution quality. This makes their search efficient but leaves their performance dependent on the chosen approximation. If its resolution is too low then they may not find a (suitable) solution but if it is too high then they may take a prohibitively long time to do so. Anytime sampling-based planners, such as RRT*, traditionally use random sampling to approximate the problem domain incrementally. This allows them to increase resolution until a suitable solution is found but makes their search dependent on the order of approximation. Arbitrary sequences of random samples approximate the problem domain in every direction simultaneously and but may be prohibitively inefficient at containing a solution. This paper unifies and extends these two approaches to develop Batch Informed Trees (BIT*), an informed, anytime sampling-based planner. BIT* solves continuous path planning problems efficiently by using sampling and heuristics to alternately approximate and search the problem domain. Its search is ordered by potential solution quality, as in A*, and its approximation improves indefinitely with additional computational time, as in RRT*. It is shown analytically to be almost-surely asymptotically optimal and experimentally to outperform existing sampling-based planners, especially on high-dimensional planning problems.


title: Fourth moment theorems on the Poisson space in any dimension

abstract: We extend to any dimension the quantitative fourth moment theorem on the Poisson setting, recently proved by C. D\"obler and G. Peccati (2017). In particular, by adapting the exchangeable pairs couplings construction introduced by I. Nourdin and G. Zheng (2017) to the Poisson framework, we prove our results under the weakest possible assumption of finite fourth moments. This yields a Peccati-Tudor type theorem, as well as an optimal improvement in the univariate case. Finally, a transfer principle "from-Poisson-to-Gaussian" is derived, which is closely related to the universality phenomenon for homogeneous multilinear sums.


title: An Interactive Tool for Natural Language Processing on Clinical Text

abstract: Natural Language Processing (NLP) systems often make use of machine learning techniques that are unfamiliar to end-users who are interested in analyzing clinical records. Although NLP has been widely used in extracting information from clinical text, current systems generally do not support model revision based on feedback from domain experts. We present a prototype tool that allows end users to visualize and review the outputs of an NLP system that extracts binary variables from clinical text. Our tool combines multiple visualizations to help the users understand these results and make any necessary corrections, thus forming a feedback loop and helping improve the accuracy of the NLP models. We have tested our prototype in a formative think-aloud user study with clinicians and researchers involved in colonoscopy research. Results from semi-structured interviews and a System Usability Scale (SUS) analysis show that the users are able to quickly start refining NLP models, despite having very little or no experience with machine learning. Observations from these sessions suggest revisions to the interface to better support review workflow and interpretation of results.


title: Trust-PCL: An Off-Policy Trust Region Method for Continuous Control

abstract: Trust region methods, such as TRPO, are often used to stabilize policy optimization algorithms in reinforcement learning (RL). While current trust region strategies are effective for continuous control, they typically require a prohibitively large amount of on-policy interaction with the environment. To address this problem, we propose an off-policy trust region method, Trust-PCL. The algorithm is the result of observing that the optimal policy and state values of a maximum reward objective with a relative-entropy regularizer satisfy a set of multi-step pathwise consistencies along any path. Thus, Trust-PCL is able to maintain optimization stability while exploiting off-policy data to improve sample efficiency. When evaluated on a number of continuous control tasks, Trust-PCL improves the solution quality and sample efficiency of TRPO.


title: Thermodynamic Formalism for Iterated Function Systems with Weights

abstract: This paper introduces an intrinsic theory of Thermodynamic Formalism for Iterated Functions Systems with general positive continuous weights (IFSw).We study the spectral properties of the Transfer and Markov operators and one of our first results is the proof of the existence of at least one eigenprobability for the Markov operator associated to a positive eigenvalue. Sufficient conditions are provided for this eingenvalue to be the spectral radius of the transfer operator and we also prove in this general setting that positive eigenfunctions of the transfer operator are always associated to its spectral radius. We introduce variational formulations for the topological entropy of holonomic measures and the topological pressure of IFSw's with weights given by a potential. A definition of equilibrium state is then natural and we prove its existence for any continuous potential. We show, in this setting, a uniqueness result for the equilibrium state requiring only the G\^ateaux differentiability of the pressure functional. We also recover the classical formula relating the powers of the transfer operator and the topological pressure and establish its uniform convergence. In the last section we present some examples and show that the results obtained can be viewed as a generalization of several classical results in Thermodynamic Formalism for ordinary dynamical systems.


title: Richardson's solutions in the real- and complex-energy spectrum

abstract: The constant pairing Hamiltonian holds exact solutions worked out by Richardson in the early Sixties. This exact solution of the pairing Hamiltonian regained interest at the end of the Nineties. The discret complex-energy states had been included in the Richardson's solutions by Hasegawa et al. [1]. In this contribution we reformulate the problem of determining the exact eigenenergies of the pairing Hamiltonian when the continuum is included through the single particle level density. The solutions with discret complex-energy states is recovered by analytic continuation of the equations to the complex energy plane. This formulation may be applied to loosely bound system where the correlations with the continuum-spectrum of energy is really important. Some details are given to show how the many-body eigenenergy emerges as sum of the pair-energies.

categories: math-ph cond-mat.supr-con math.MP nucl-th


title: The rank of Mazur's Eisenstein ideal

abstract: We use pseudodeformation theory to study Mazur's Eisenstein ideal. Given prime numbers $N$ and $p>3$, we study the Eisenstein part of the $p$-adic Hecke algebra for $\Gamma_0(N)$. We compute the rank of this Hecke algebra (and, more generally, its Newton polygon) in terms of Massey products in Galois cohomology, answering a question of Mazur and generalizing a result of Calegari-Emerton. We also also give new proofs of Merel's result on this rank and of Mazur's results on the structure of the Hecke algebra.


title: Adaptive user support in educational environments: A Bayesian Network approach

abstract: This paper is concerned with the design and implementation of an innovative user support system in the frame of an open educational environment. The environment adapted is ModelsCreator (MC), an educational system supporting learning through modelling activities. The pupils typical interaction with the system was modelled us-ing Bayesian Belief Networks (BBN). This model has been used in ModelsCreator to build an adaptive help system providing the most useful guidelines according to the current state of interaction. A brief description of the system and an overview of application of Bayesian techniques to educational systems is presented together with discussion about the process of building of the Bayesian Network derived from actual student interaction data. A preliminary evaluation of the developed prototype indicates that the proposed approach produces systems with promising performance.


title: Deformation conditions for pseudorepresentations

abstract: Given a property of representations satisfying a basic stability condition, Ramakrishna developed a variant of Mazur's Galois deformation theory for representations with that property. We introduce an axiomatic definition of pseudorepresentations with such a property. Among other things, we show that pseudorepresentations with a property enjoy a good deformation theory, generalizing Ramakrishna's theory to pseudorepresentations.


title: Pressure densification of a simple liquid

abstract: The magnitude of the high frequency, static dielectric permittivity is used to determine the density of tetramethyl tetraphenyl trisiloxane, a non-associated glass-forming liquid, as a function of temperature and pressure. We demonstrate that the properties in the glassy state are affected by the pressure applied to the liquid during vitrification. This behavior is normal for hydrogen-bonded liquids and polymers, but unanticipated by models of simple liquids.


title: Adaptive Modular Exponentiation Methods v.s. Python's Power Function

abstract: In this paper we use Python to implement two efficient modular exponentiation methods: the adaptive m-ary method and the adaptive sliding-window method of window size k, where both m's are adaptively chosen based on the length of exponent. We also conduct the benchmark for both methods. Evaluation results show that compared to the industry-standard efficient implementations of modular power function in CPython and Pypy, our algorithms can reduce 1-5% computing time for exponents with more than 3072 bits.


title: On Calculation of Bounds for Greedy Algorithms when Applied to Sensor Selection Problems

abstract: We consider the problem of studying the performance of greedy algorithm on sensor selection problem for stable linear systems with Kalman Filter. Specifically, the objective is to find the system parameters that affects the performance of greedy algorithms and conditions where greedy algorithm always produces optimal solutions. In this paper, we developed an upper bound for performance ratio of greedy algorithm, which is based on the work of Dr.Zhang \cite{Sundaram} and offers valuable insight into the system parameters that affects the performance of greedy algorithm. We also proposes a set of conditions where greedy algorithm will always produce the optimal solution. We then show in simulations how the system parameters mentioned by the performance ratio bound derived in this work affects the performance of greedy algorithm.


title: On the Distinct Periodicities of Sunspot Counts in Flaring and Non-flaring Active Regions

abstract: In a recent work, Kilcik et al. (2017), have detected the temporal and periodic behavior of sunspot counts (SSC) in flaring (i.e. C, M, or X class flares), and non-flaring active regions for the last two solar cycles, covering the period: 1996 - 2016. The main results obtained are: 1) The temporal behavior of monthly means of daily total SSC in flaring and non-flaring active regions are different and these differences are also varying from cycle to cycle; 2) The periodicities detected in SSC of flaring and non-flaring active regions are quite different and these variations are also different from one cycle to another; the highest detected period in the flaring active regions is 113 days, while there are much higher periodicities (327, 312, and 256 days) in non-flaring regions. The detection of typical different periodicities in flaring and non-flaring regions can suggests both important differences and physical interpretation in the magneto-hydrodynamic behavior of the Sun. For this reason in the present paper we show a further periodicity analysis of the sunspot counts in flaring and in non-flaring active regions using the same data source of that used by the above cited authors and applying a powerful wavelet analysis tool which is particularly useful to detect multiscale features of complex unsteady and unevenly sampled time series. In order to futher support the differences and similarities found in the time behavior of SSC in flaring and non-flaring regions, we also computed the behavior of the wavelet entropy, a proper time function which allow us to measure the degree of complexity in the dynamics of the related time series.


title: Entropy of Spatial Network Ensembles

abstract: We analyze complexity in spatial network ensembles through the lens of graph entropy. Mathematically, we model a spatial network as a soft random geometric graph, i.e., a graph with two sources of randomness, namely nodes located randomly in space and links formed independently between pairs of nodes with probability given by a specified function (the "pair connection function") of their mutual distance. We consider the general case where randomness arises in node positions as well as pairwise connections (i.e., for a given pair distance, the corresponding edge state is a random variable). Classical random geometric graph and exponential graph models can be recovered in certain limits. We derive a simple bound for the entropy of a spatial network ensemble and calculate the conditional entropy of an ensemble given the node location distribution for hard and soft (probabilistic) pair connection functions. Under this formalism, we derive the connection function that yields maximum entropy under general constraints. Finally, we apply our analytical framework to study two practical examples: ad hoc wireless networks and the US flight network. Through the study of these examples, we illustrate that both exhibit properties that are indicative of nearly maximally entropic ensembles.

categories: physics.soc-ph cond-mat.stat-mech cs.IT math.IT


title: Limits on Spectral Resolution Measurements by Quantum Probes

abstract: The limits of frequency resolution in nano-NMR experiments have been discussed extensively in recent years. It is believed that there is a crucial difference between the ability to resolve a few frequencies and the precision of estimating a single one. Whereas the efficiency of single frequency estimation gradually increases with the square root of the number of measurements, the ability to resolve two frequencies is limited by the specific timescale of the signal and cannot be compensated for by extra measurements. Here we show theoretically and demonstrate experimentally that the relationship between these quantities is more subtle and both are only limited by the Cram\'er-Rao bound of a single frequency estimation.


title: Structure and Topology of Band Structures in the 1651 Magnetic Space Groups

abstract: The properties of electrons in magnetically ordered crystals are of interest both from the viewpoint of realizing novel topological phases, such as magnetic Weyl semimetals, and from the applications perspective of creating energy-efficient memories. A systematic study of symmetry and topology in magnetic materials has been challenging given that there are 1651 magnetic space groups (MSGs). Here, by using an efficient representation of allowed band structures, we obtain a systematic description of several basic properties of free electrons in all MSGs in three dimensions as well as in the 528 magnetic layer groups relevant to two dimensional magnetic materials. We compute constraints on electron fillings and band connectivity compatible with insulating behavior. Also, by contrasting with atomic insulators, we identify band topology entailed by the symmetry transformation of bands, as determined by the MSG alone. We give an application of our results to identifying topological semimetals arising in periodic arrangements of hedgehog-like magnetic textures.


title: Individual Stellar Halos of Massive Galaxies Measured to 100 kpc at $0.3<z<0.5$ using Hyper Suprime-Cam

abstract: Massive galaxies display extended light profiles that can reach several hundreds of kilo parsecs. These stellar halos provide a fossil record of galaxy assembly histories. Using data that is both wide (~100 square degree) and deep (i>28.5 mag/arcsec^2 in i-band), we present a systematic study of the stellar halos of a sample of more than 3000 galaxies at 0.3 < z < 0.5 with $\log M_{\star}/M_{\odot} > 11.4$. Our study is based on high-quality (0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP), which enables us to individually estimate surface mass density profiles to 100 kpc without stacking. As in previous work, we find that more massive galaxies exhibit more extended outer profiles. When this extended light is not properly accounted for as a result of shallow imaging or inadequate profile modeling, the derived stellar mass function can be significantly underestimated at the highest masses. Across our sample, the ellipticity of outer light profiles increases substantially as we probe larger radii. We show for the first time that these ellipticity gradients steepen dramatically as a function of galaxy mass, but we detect no mass-dependence in outer color gradients. Our results support the two-phase formation scenario for massive galaxies in which outer envelopes are built up at late times from a series of merging events. We provide surface mass surface mass density profiles in a convenient tabulated format to facilitate comparisons with predictions from numerical simulations of galaxy formation.


title: Pulsar TeV Halos Explain the TeV Excess Observed by Milagro

abstract: Milagro observations have found bright, diffuse TeV emission concentrated along the galactic plane of the Milky Way. The intensity and spectrum of this emission is difficult to explain with current models where gamma-ray production is dominated by hadronic mechanisms, and has been named the "TeV excess". We show that TeV emission from pulsars naturally explains this excess. In particular, recent observations have detected "TeV halos" surrounding pulsars that are either nearby or particularly luminous. Here, we show that the full population of Milky Way pulsars will produce diffuse TeV emission concentrated along the Milky Way plane. The total gamma-ray flux from TeV halos is expected to exceed the hadronic gamma-ray flux at energies above ~500 GeV. Moreover, the spectrum and intensity of TeV halo emission naturally matches the TeV excess. If this scenario is common to all galaxies, it will decrease the contribution of star-forming galaxies to the IceCube neutrino flux. Finally, we show that upcoming HAWC observations will resolve a significant fraction of the TeV excess into individual TeV halos, conclusively confirming, or ruling out, this model.


title: Joint Bayesian Estimation of Quasar Continua and the Lyman-Alpha Forest Flux Probability Distribution Function

abstract: We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.


title: Constraints on the mass-richness relation from the abundance and weak lensing of SDSS clusters

abstract: We constrain the scaling relation between optical richness ($\lambda$) and halo mass ($M$) for a sample of SDSS redMaPPer galaxy clusters within the context of the {\it Planck} cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, $P(\ln \lambda| M)$. To model the abundance and the stacked lensing profiles, we use an emulator specifically built to interpolate the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, which is calibrated based on a suite of high-resolution $N$-body simulations. We apply our method to 8,312 SDSS redMaPPer clusters with $20\le \lambda \le 100$ and $0.10\le z_{\lambda}\le0.33$, and show that the log-normal distribution model for $P(\lambda|M)$, with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, $\left\langle \ln{\lambda}\right\rangle(M)=A+B\ln(M/M_{\rm pivot})$, with $A=3.207^{+0.044}_{-0.046}$ and $B=0.993^{+0.041}_{-0.055}$ (68\%~CL), where the pivot mass scale $M_{\rm pivot}=3\times 10^{14} h^{-1}M_\odot$, and the scatter $\sigma_{\mathrm{\ln\lambda}|M}=\sigma_0+q\ln(M/M_{\rm pivot})$ with $\sigma_0=0.456^{+0.047}_{-0.039}$ and $q=-0.169^{+0.035}_{-0.026}$. We find that a large scatter in halo masses is required at the lowest richness bins ($20\le \lambda \lesssim 30$) in order to reproduce the measurements. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of intrinsically low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of $30\le \lambda \le 100$.


title: Microlensing Makes Lensed Quasar Time Delays Significantly Time Variable

abstract: The time delays of gravitationally lensed quasars are generally believed to be unique numbers whose measurement is limited only by the quality of the light curves and the models for the contaminating contribution of gravitational microlensing to the light curves. This belief is incorrect -- gravitational microlensing also produces changes in the actual time delays on the ~day(s) light-crossing time scale of the emission region. This is due to a combination of the inclination of the disk relative to the line of sight and the differential magnification of the temperature fluctuations producing the variability. We demonstrate this both mathematically and with direct calculations using microlensing magnification patterns. Measuring these delay fluctuations can provide a physical scale for microlensing observations, removing the need for priors on either the microlens masses or the component velocities. That time delays in lensed quasars are themselves time variable likely explains why repeated delay measurements of individual lensed quasars appear to vary by more than their estimated uncertainties. This effect is also an important new systematic problem for attempts to use time delays in lensed quasars for cosmology or to detect substructures (satellites) in lens galaxies.


title: Simulating Neutron Star Mergers as r-process Sources in Ultra Faint Dwarf Galaxies

abstract: To explain the high observed abundances of r-process elements in local ultra- faint dwarf (UFD) galaxies, we perform cosmological zoom simulations that include r-process production from neutron star mergers (NSMs). We model star-formation stochastically and simulate two different halos with total masses $\approx 10^8 M_{\odot}$ at z = 6. We find that the final distribution of [Eu/H] vs. [Fe/H] is relatively insensitive to the energy by which the r-process material is ejected into the interstellar medium, but strongly sensitive to the environment in which the NSM event occurs. In one halo the NSM event takes place at the center of the stellar distribution, leading to high-levels of r-process enrichment such as seen in a local UFD, Reticulum II (Ret II). In a second halo, the NSM event takes place outside of the densest part of the galaxy, leading to a more extended r-process distribution. The subsequent star formation occurs in an interstellar medium with shallow levels of r-process enrichment which results in stars with low levels of [Eu/H] compared to Ret II stars even when the maximum possible r-process mass is assumed to be ejected. This suggests that the natal kicks of neutron stars may also play an important role in determining the r-process abundances in UFD galaxies, a topic that warrants further theoretical investigation.


title: Distortions of the Cosmic Microwave Background through cooling lines during the epoch of Reionization

abstract: By using N-body hydrodynamical cosmological simulations in which the chemistry of major metals and molecules is consistently solved for, we study the interaction of metallic fine-structure lines with the CMB. Our analysis shows that the collisional induced emissions in the OI 145 $\mu$m and CII 158 $\mu$m lines during reionization introduce a distortion of the CMB spectrum at low frequencies ($\nu < 300$ GHz) with amplitudes up to $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-8}$-$10^{-7}$, i.e., at the $\sim 0.1$ percent level of FIRAS upper limits. Shorter wavelength fine-structure transitions (OI 63 $\mu$m, FeII 26 $\mu$m, and SiII 35 $\mu$m) typically sample the reionization epoch at higher observing frequencies ($\nu > 400$ GHz). This corresponds to the Wien tail of the CMB spectrum and the distortion level induced by those lines may be as high as $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-4}$. The angular anisotropy produced by these lines should be more relevant at higher frequencies: while practically negligible at $\nu=145 $GHz, signatures from CII 158 $\mu$m and OI 145 $\mu$m should amount to 1%-5% of the anisotropy power measured at $l \sim 5000$ and $\nu=220 $GHz by the ACT and SPT collaborations (after assuming $\Delta \nu_{\rm obs}/\nu_{\rm obs}\simeq 0.005$ for the line observations). Our simulations show that anisotropy maps from different lines (e.g., OI 145 $\mu$m and CII 158 $\mu$m) at the same redshift show a very high degree ($>0.8$) of spatial correlation, allowing for the use of observations at different frequencies to unveil the same snapshot of the reionization epoch. Finally, our simulations demonstrate that line-emission anisotropies extracted in narrow frequency/redshift shells are practically uncorrelated in frequency space, thus enabling standard methods for removal of foregrounds that vary smoothly in frequency, just as in HI 21 cm studies.


title: Small-N collisional dynamics III: The battle for the realm of not-so-small-N

abstract: In this paper, the third in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We study the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and we focus on interactions that produce direct collisions between any two stars. Our long-term goal is to construct an equation that gives the probability of a given collision event occurring over the course of the interaction, as a function of the total encounter energy and angular momentum as well as the numbers and properties of the particles. In previous papers, we varied the number of interacting particles and the distribution of particle radii, for all equal mass particles. Here, we focus on the effects of different combinations of particle masses. We develop an analytic formalism for calculating the time-scales for different collision scenarios to occur. Our analytic time-scales reproduce the simulated time-scales when gravitational focusing is included. We present a method for calculating the relative rates for different types of collisions to occur, assuming two different limits for the particle orbits; radial and tangential. These limits yield relative collision probabilities that bracket the probabilities we obtain directly from numerical scattering experiments, and are designed to reveal important information about the (time-averaged) trajectories of the particles as a function of the interaction parameters. Finally, we present a Collision Rate Diagram (CRD), which directly compares the predictions of our analytic rates to the simulations and quantifies the quality of the agreement. The CRD will facilitate refining our analytic collision rates in future work, as we expand in to the remaining parameter space.


title: Observational properties of massive black hole binary progenitors

abstract: The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ~36Msun and ~29Msun. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (PoWR), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. We provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.


title: Do galaxies that leak ionizing photons have extreme outflows?

abstract: To reionize the early universe, high-energy photons must escape the galaxies that produce them. It has been suggested that stellar feedback drives galactic outflows out of star-forming regions, creating low density channels through which ionizing photons escape into the inter-galactic medium. We compare the galactic outflow properties of confirmed Lyman continuum (LyC) leaking galaxies to a control sample of nearby star-forming galaxies to explore whether the outflows from leakers are extreme as compared to the control sample. We use data from the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure the equivalent widths and velocities of Si II and Si III absorption lines, tracing neutral and ionized galactic outflows. We find that the Si II and Si III equivalent widths of the LyC leakers reside on the low-end of the trend established by the control sample. The leakers' velocities are not statistically different than the control sample, but their absorption line profiles have a different asymmetry: their central velocities are closer to their maximum velocities. The outflow kinematics and equivalent widths are consistent with the scaling relations between outflow properties and host galaxy properties -- most notably metallicity -- defined by the control sample. Additionally, we use the Ly\alpha\ profiles to show that the Si II equivalent width scales with the Ly\alpha\ peak velocity separation. We determine that the low equivalent widths of the leakers are likely driven by low metallicities and low H I column densities, consistent with a density-bounded ionization region, although we cannot rule out significant variations in covering fraction. While we do not find that the LyC leakers have extreme outflow velocities, the low maximum-to-central velocity ratios demonstrate the importance of the acceleration and density profiles for LyC and Ly\alpha\ escape. [abridged]


title: Collapsar Gamma-Ray Bursts: how the luminosity function dictates the duration distribution

abstract: Jets in long-duration $\gamma$-ray bursts (GRBs) have to drill through the collapsing star in order to break out of it and produce the $\gamma$-ray signal while the central engine is still active. If the breakout time is shorter for more powerful engines, then the jet-collapsar interaction acts as a filter of less luminous jets. We show that the observed broken power-law GRB luminosity function is a natural outcome of this process. For a theoretically motivated breakout time that scales with jet luminosity as $L^{-\chi}$ with $\chi\sim 1/3-1/2$, we show that the shape of the $\gamma$-ray duration distribution can be uniquely determined by the GRB luminosity function and matches the observed one. This analysis has also interesting implications about the supernova-central engine connection. We show that not only successful jets can deposit sufficient energy in the stellar envelope to power the GRB-associated supernovae, but also failed jets may operate in all Type Ib/c supernovae.


title: Counting states and the Hadron Resonance Gas: Does X(3872) count?

abstract: We analyze how the renowned X(3872), a weakly bound state right below the $D \bar D^*$ threshold, should effectively be included in a hadronic representation of the QCD partition function. This can be decided by analyzing the $D \bar D^*$ scattering phase-shifts in the $J^{PC}=1^{++}$ channel and their contribution to the level density in the continuum from which the abundance in a hot medium can be determined. We show that in a purely molecular picture the bound state contribution cancels the continuum providing a vanishing occupation number density at finite temperature and the $X(3872)$ does not count below the Quark-Gluon Plasma crossover happening at $T \sim 150$MeV. In contrast, within a coupled-channels approach, for a non vanishing $c \bar c$ content the cancellation does not occur due to the onset of the $X(3940)$ which effectively counts as an elementary particle for temperatures above $T \gtrsim 250$MeV. Thus, a direct inclusion of the $X(3872)$ in the Hadron Resonance Gas is not justified. We also estimate the role of this cancellation in X(3872) production in heavy-ion collision experiments in terms of the corresponding $p_T$ distribution due to a finite energy resolution.


title: Luminous and Variable Stars in M31 and M33 V. The Upper HR Diagram

abstract: We present HR Diagrams for the massive star populations in M31 and M33 including several different types of emission-line stars: the confirmed Luminous Blue Variables (LBVs), candidate LBVs, B[e] supergiants and the warm hypergiants. We estimate their apparent temperatures and luminosities for comparison with their respective massive star populations and to evaluate the possible relationships of these different classes of evolved, massive stars, and their evolutionary state. Several of the LBV candidates lie near the LBV/S Dor instability strip which supports their classification. Most of the B[e] supergiants, however, are less luminous than the LBVs. Many are very dusty with the infrared flux contributing one-third or more to their total flux. They are also relatively isolated from other luminous OB stars. Overall, their spatial distribution suggests a more evolved state. Some may be post-RSGs like the warm hypergiants, and there may be more than one path to becoming a B[e] star. There are sufficient differences in the spectra, luminosities, spatial distribution, and the presence or lack of dust between the LBVs and B[e] supergiants to conclude that one group does not evolve into the other.


title: Higher-order Relation Schema Induction using Tensor Factorization with Back-off and Aggregation

abstract: Relation Schema Induction (RSI) is the problem of identifying type signatures of arguments of relations from unlabeled text. Most of the previous work in this area have focused only on binary RSI, i.e., inducing only the subject and object type signatures per relation. However, in practice, many relations are high-order, i.e., they have more than two arguments and inducing type signatures of all arguments is necessary. For example, in the sports domain, inducing a schema win(WinningPlayer, OpponentPlayer, Tournament, Location) is more informative than inducing just win(WinningPlayer, OpponentPlayer). We refer to this problem as Higher-order Relation Schema Induction (HRSI). In this paper, we propose Tensor Factorization with Back-off and Aggregation (TFBA), a novel framework for the HRSI problem. To the best of our knowledge, this is the first attempt at inducing higher-order relation schemata from unlabeled text. Using the experimental analysis on three real world datasets, we show how TFBA helps in dealing with sparsity and induce higher order schemata.


title: All $(4,0)$: Sigma Models with $(4,0)$ Off-Shell Supersymmetry

abstract: Off-shell $(4,0)$ supermultiplets in 2-dimensions are formulated. These are used to construct sigma models whose target spaces are vector bundles over manifolds that are hyperk\"ahler with torsion. The off-shell supersymmetry implies that the complex structures are simultaneously integrable and allows us to write actions using extended superspace and projective superspace, giving an explicit construction of the target space geometries.


title: Sonic metamaterials: reflection on the role of topology on dispersion surface morphology

abstract: Investigating dispersion surface morphology of sonic metamaterials is crucial in providing information on related phenomena as inertial coupling, acoustic transparency, polarisation, and absorption. In the present study, we look into frequency surface morphology of two-dimensional metamaterials of K3,3 and K6 topologies. The elastic structures under consideration consist of the same substratum lattice points and form a pair of sublattices with hexagonal symmetry. We show that, through introducing universal localised mass-in-mass phononic microstructures at lattice points, six single optical frequency-surfaces can be formed with required properties including negative group velocity. Splitting the frequency-surfaces is based on the classical analog of the quantum phenomenon of "energy-level repulsion", which can be achieved only through internal anisotropy of the nodes and allows us to obtain different frequency band gaps.


title: A General Perspective of Fe-Mn-Al-C Steels

abstract: During the last years, the scientific and industrial community has focused on the astonishing properties of Fe-Mn-Al-C steels. These high advanced steels allow high-density reductions about ~15% lighter than conventional steels, high corrosion resistance, high strength (ultimate tensile strength (UTS) ~1 Gpa) and at the same time ductilities above 60%. The increase of the tensile or yield strength and the ductility at the same time is almost a special feature of this kind of new steels, which makes them so interesting for many applications such as in the automotive, armor and mining industry. The control of these properties depends on a complex relationship between the chemical composition of the steel, the test temperature, the external loads and the processing parameters of the steel. This review has been conceived to tried to elucidate these complex relations and gather the most important aspects of Fe-Mn-Al-C steels developed so far.


title: A Visual Narrative Path from Switching to Resuming a Requirements Engineering Task

abstract: Requirements Engineering (RE) is closely tied to other development activities and is at the heart and foundation of every software development process. This makes RE the most data and communication-intensive activity compared to other development tasks. The highly demanding communication makes task switching and interruptions inevitable in RE activities. While task switching often allows us to perform tasks effectively, it imposes a cognitive load and can be detrimental to the primary task, particularly in complex tasks as the ones typical for RE activities. Visualization mechanisms enhanced with analytical methods and interaction techniques help software developers obtain a better cognitive understanding of the complexity of RE decisions, leading to timelier and higher quality decisions. In this paper, we propose to apply interactive visual analytics techniques for managing requirements decisions from various perspectives, including stakeholders communication, RE task switching, and interruptions. We propose a new layered visualization framework that supports the analytical reasoning process of task switching. This framework consists of both data analysis and visualization layers. The visual layers offer interactive knowledge visualization components for managing task interruption decisions at different stages of an interruption (i.e. before, during, and after). The analytical layers provide narrative knowledge about the consequences of task switching decisions and help requirements engineers to recall their reasoning process and decisions upon resuming a task. Moreover, we surveyed 53 software developers to test our visual prototype and to explore more required features for the visual and analytical layers of our framework.


title: Zero-Shot Deep Domain Adaptation

abstract: Domain adaptation is an important tool to transfer knowledge about a task (e.g. classification) learned in a source domain to a second, or target domain. Current approaches assume that task-relevant target-domain data is available during training. We demonstrate how to perform domain adaptation when no such task-relevant target-domain data is available. To tackle this issue, we propose zero-shot deep domain adaptation (ZDDA), which uses privileged information from task-irrelevant dual-domain pairs. ZDDA learns a source-domain representation which is not only tailored for the task of interest but also close to the target-domain representation. Therefore, the source-domain task of interest solution (e.g. a classifier for classification tasks) which is jointly trained with the source-domain representation can be applicable to both the source and target representations. Using the MNIST, Fashion-MNIST, NIST, EMNIST, and SUN RGB-D datasets, we show that ZDDA can perform domain adaptation in classification tasks without access to task-relevant target-domain training data. We also extend ZDDA to perform sensor fusion in the SUN RGB-D scene classification task by simulating task-relevant target-domain representations with task-relevant source-domain data. To the best of our knowledge, ZDDA is the first domain adaptation and sensor fusion method which requires no task-relevant target-domain data. The underlying principle is not particular to computer vision data, but should be extensible to other domains.


title: Facilitated exclusion process

abstract: We study the Facilitated TASEP, an interacting particle system on the one dimensional integer lattice. We prove that starting from step initial condition, the position of the rightmost particle has Tracy Widom GSE statistics on a cube root time scale, while the statistics in the bulk of the rarefaction fan are GUE. This uses a mapping with last-passage percolation in a half-quadrant which is exactly solvable through Pfaffian Schur processes. Our results further probe the question of how first particles fluctuate for exclusion processes with downward jump discontinuities in their limiting density profiles. Through the Facilitated TASEP and a previously studied MADM exclusion process we deduce that cube-root time fluctuations seem to be a common feature of such systems. However, the statistics which arise are shown to be model dependent (here they are GSE, whereas for the MADM exclusion process they are GUE). We also discuss a two-dimensional crossover between GUE, GOE and GSE distribution by studying the multipoint distribution of the first particles when the rate of the first one varies. In terms of half-space last passage percolation, this corresponds to last passage times close to the boundary when the size of the boundary weights is simultaneously scaled close to the critical point.

categories: math.PR cond-mat.stat-mech math-ph math.MP


title: Non-vanishing of certain cyclotomic multiple harmonic sums and application to the non-vanishing of certain $p$-adic cyclotomic multiple zeta values

abstract: We define and apply a method to study the non-vanishing of $p$-adic cyclotomic multiple zeta values. We prove the non-vanishing of certain cyclotomic multiple harmonic sums, and, via a formula proved in another paper, which expresses a cyclotomic multiple harmonic sums as an infinite sum of products of $p$-adic cyclotomic multiple zeta values, this implies the non-vanishing of certain $p$-adic cyclotomic multiple zeta values.


title: Biharmonic Immersion in Cartan Hadamard manifold

abstract: If $(N^{m+p},h)$ is a Cartan-Hadamard manifold such that $Ric(h)\geq -G(r_{N}(x))$ where $G(0)\geq 1, G^{'}\geq 0$ and $G^{-1/2}\not\in L^{1}(+\infty)$ then every proper biharmonic isometric immersion $\phi : M^{m}\rightarrow(N^{m+p},h)$ is a harmonic map.


title: Diffusion Convolutional Recurrent Neural Network: Data-Driven Traffic Forecasting

abstract: Spatiotemporal forecasting has various applications in neuroscience, climate and transportation domain. Traffic forecasting is one canonical example of such learning task. The task is challenging due to (1) complex spatial dependency on road networks, (2) non-linear temporal dynamics with changing road conditions and (3) inherent difficulty of long-term forecasting. To address these challenges, we propose to model the traffic flow as a diffusion process on a directed graph and introduce Diffusion Convolutional Recurrent Neural Network (DCRNN), a deep learning framework for traffic forecasting that incorporates both spatial and temporal dependency in the traffic flow. Specifically, DCRNN captures the spatial dependency using bidirectional random walks on the graph, and the temporal dependency using the encoder-decoder architecture with scheduled sampling. We evaluate the framework on two real-world large scale road network traffic datasets and observe consistent improvement of 12% - 15% over state-of-the-art baselines.


title: Let's hear it from RETTA: A Requirements Elicitation Tool for TrAffic management systems

abstract: The area of Traffic Management (TM) is characterized by uncertainty, complexity, and imprecision. The complexity of software systems in the TM domain which contributes to a more challenging Requirements Engineering (RE) job mainly stems from the diversity of stakeholders and complexity of requirements elicitation in this domain. This work brings an interactive solution for exploring functional and non-functional requirements of software-reliant systems in the area of traffic management. We prototyped the RETTA tool which leverages the wisdom of the crowd and combines it with machine learning approaches such as Natural Language Processing and Naive Bayes to help with the requirements elicitation and classification task in the TM domain. This bridges the gap among stakeholders from both areas of software development and transportation engineering. The RETTA prototype is mainly designed for requirements engineers and software developers in the area of TM and can be used on Android-based devices.


title: Transient Mass Loss Analysis of Solar Observations using Stellar Methods

abstract: Low frequency dynamic spectra of radio bursts from nearby stars offer the best chance to directly detect the stellar signature of transient mass loss on low mass stars. Crosley et al. (2016) proposes a multi-wavelength methodology to determine coronal mass ejection parameters, such as Coronal Mass Ejection (CME) speed, mass, and kinetic energy. We test the validity and accuracy of the results derived from the methodology by using Geostationary Operational Environmental Satellite X-ray observations and Bruny Island Radio Spectrometer radio observations. These are analogous observations to those which would be found in the stellar studies. Derived results from these observations are compared to direct white light measurements of the Large Angle and Spectrometric Coronagraph. We find that, when a pre-event temperature can be determined, that the accuracy of CME speeds are within a few hundred km/s, and are reliable when specific criteria has been met. CME mass and kinetic energies are only useful in determining approximate order of magnitude measurements when considering the large errors associated to them. These results will be directly applicable to interpretation of any detected stellar events and derivation of stellar CME properties.


title: Conformal metrics with prescribed fractional scalar curvature on conformal infinities with positive fractional Yamabe constants

abstract: Let $(X, g^+)$ be an asymptotically hyperbolic manifold and $(M, [\hat{h}])$ its conformal infinity. Our primary aim in this paper is to introduce the prescribed fractional scalar curvature problem on $M$ and provide solutions under various geometric conditions on $X$ and $M$. We also obtain the existence results for the fractional Yamabe problem in the endpoint case, e.g., $n = 3$, $\gamma = 1/2$ and $M$ is non-umbilic, etc. Every solution we find turns out to be smooth on $M$.


title: A tale of stars and cliques

abstract: We show that for an infinitely many natural numbers $k$ there are $k$-uniform hypergraphs which admit a `rescaling phenomenon' as described in [9]. More precisely, let $\mathcal{A}(k,I, n)$ denote the class of $k$-graphs on $n$ vertices in which the sizes of all pairwise intersections of edges belong to a set $I$. We show that if $k=rt^2$ for some $r\ge 1$ and $t\ge 2$, and~$I$ is chosen in some special way, the densest graphs in $\mathcal{A}(rt^2,I, n)$ are either dominated by stars of large degree, or basically, they are `$t$-thick' $rt^2$-graphs in which vertices are partitioned into groups of $t$ vertices each and every edge is a union of $tr$ such groups. It is easy to see that, unlike in stars, the maximum degree of $t$-thick graphs is of a lower order than the number of its edges. Thus, if we study the graphs from $\mathcal{A}(rt^2,I, n)$ with a prescribed number of edges $m$ which minimize the maximum degree, around the value of $m$ which is the number of edges of the largest $t$-thick graph, a rapid, discontinuous phase transition can be observed. Interestingly, these two types of $k$-graphs determine the structure of all hypergraphs in $\mathcal{A}(rt^2,I, n)$. Namely, we show that each such hypergraph can be decomposed into a $t$-thick graph $H_T$, a special collection $H_S$ of stars, and a sparse `left-over' graph $H_R$.


title: Lattice paths with catastrophes

abstract: In queuing theory, it is usual to have some models with a "reset" of the queue. In terms of lattice paths, it is like having the possibility of jumping from any altitude to zero. These objects have the interesting feature that they do not have the same intuitive probabilistic behaviour as classical Dyck paths (the typical properties of which are strongly related to Brownian motion theory), and this article quantifies some relations between these two types of paths. We give a bijection with some other lattice paths and a link with a continued fraction expansion. Furthermore, we prove several formulae for related combinatorial structures conjectured in the On-Line Encyclopedia of Integer Sequences. Thanks to the kernel method and via analytic combinatorics, we provide the enumeration and limit laws of these "lattice paths with catastrophes" for any finite set of jumps. We end with an algorithm to generate such lattice paths uniformly at random.


title: End-to-End Learning of Semantic Grasping

abstract: We consider the task of semantic robotic grasping, in which a robot picks up an object of a user-specified class using only monocular images. Inspired by the two-stream hypothesis of visual reasoning, we present a semantic grasping framework that learns object detection, classification, and grasp planning in an end-to-end fashion. A "ventral stream" recognizes object class while a "dorsal stream" simultaneously interprets the geometric relationships necessary to execute successful grasps. We leverage the autonomous data collection capabilities of robots to obtain a large self-supervised dataset for training the dorsal stream, and use semi-supervised label propagation to train the ventral stream with only a modest amount of human supervision. We experimentally show that our approach improves upon grasping systems whose components are not learned end-to-end, including a baseline method that uses bounding box detection. Furthermore, we show that jointly training our model with auxiliary data consisting of non-semantic grasping data, as well as semantically labeled images without grasp actions, has the potential to substantially improve semantic grasping performance.


title: Application of the Kronecker product to simple spin systems

abstract: We show that the well known Kronecker product is a suitable tool for the construction of matrix representations of widely used spin Hamiltonians. In this way we avoid the explicit use of basis sets for the construction of the matrix elements. As illustrative examples we discuss two isotropic models and an anisotropic one.


title: Elliptic equations involving the $1$-Laplacian and a subcritical source term

abstract: In this paper we deal with a Dirichlet problem for an elliptic equation involving the $1$-Laplacian operator and a source term. We prove that, when the growth of the source is subcritical, there exist two bounded nontrivial solutions to our problem. Moreover, a Pohozaev type identity is proved, which holds even when the growth is supercritical. We also show explicit examples of our results.


title: Root data with group actions

abstract: Suppose $k$ is a field, $G$ is a connected reductive algebraic $k$-group, $T$ is a maximal $k$-torus in $G$, and $\Gamma$ is a finite group that acts on $(G,T)$. From the above, one obtains a root datum $\Psi$ on which $\text{Gal}(k)\times\Gamma$ acts. Provided that $\Gamma$ preserves a positive system in $\Psi$, not necessarily invariant under $\text{Gal}(k)$, we construct an inverse to this process. That is, given a root datum on which $\text{Gal}(k)\times\Gamma$ acts appropriately, we show how to construct a pair $(G,T)$, on which $\Gamma$ acts as above. Although the pair $(G,T)$ and the action of $\Gamma$ are canonical only up to an equivalence relation, we construct a particular pair for which $G$ is $k$-quasisplit and $\Gamma$ fixes a $\text{Gal}(k)$-stable pinning of $G$. Using these choices, we can define a notion of taking "$\Gamma$-fixed points" at the level of equivalence classes, and this process is compatible with a general "restriction" process for root data with $\Gamma$-action.


title: Improving superconducting resonators in magnetic fields by reduced field-focussing and engineered flux screening

abstract: We experimentally investigate superconducting coplanar waveguide resonators in external magnetic fields and present two strategies to reduce field-induced dissipation channels and resonance frequency shifts. One of our approaches is to significantly reduce the superconducting ground-plane areas, which leads to reduced magnetic field-focussing and thus to lower effective magnetic fields inside the waveguide cavity. By this measure, the field-induced losses can be reduced by more than one order of magnitude in mT out-of-plane magnetic fields. When these resonators are additionally coupled inductively instead of capacitively to the microwave feedlines, an intrinsic closed superconducting loop is effectively shielding the heart of the resonator from magnetic fields by means of flux conservation. In total, we achieve a reduction of the field-induced resonance frequency shift by up to two orders of magnitude. We combine systematic parameter variations on the experimental side with numerical magnetic field calculations to explain the effects of our approaches and to support our conclusions. The presented results are relevant for all areas, where high-performance superconducting resonators need to be operated in magnetic fields, e.g. for quantum hybrid devices with superconducting circuits or electron spin resonance detectors based on coplanar waveguide cavities.


title: Simultaneous Monitoring of X-ray and Radio Variability in Sagittarius A*

abstract: Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on sub-mm and radio variability hints at a lag between X-ray flares and their candidate sub-millimeter or radio counterparts, with the long wavelength data lagging the X-ray. However, there is only one published time lag between an X-ray flare and a possible radio counterpart. Here we report 9 contemporaneous X-ray and radio observations of Sgr A*. We detect significant radio variability peaking $\gtrsim$176 minutes after the brightest X-ray flare ever detected from Sgr A*. We also report other potentially associated X-ray and radio variability, with the radio peaks appearing $\lesssim$80 minutes after these weaker X-ray flares. Taken at face value, these results suggest that stronger X-ray flares lead to longer time lags in the radio. However, we also test the possibility that the variability at X-ray and radio wavelengths is not temporally correlated. We cross-correlate data from mismatched X-ray and radio epochs and obtain comparable correlations to the matched data. Hence, we find no overall statistical evidence that X-ray flares and radio variability are correlated, underscoring a need for more simultaneous, long duration X-ray--radio monitoring of Sgr A*.


title: Euler vs. Lagrange: The role of coordinates in practical Evans-function computations

abstract: The Evans function has become a standard tool in the mathematical study of nonlinear wave stability. In particular, computation of its zero set gives a convenient numerical method for determining the point spectrum of the associated linear operator (and thus the spectral stability of the wave in question). We report on an unexpected complication that frustrates this computation for viscous shock profiles in gas dynamics. Although this phenomenon---related to the choice of Eulerian or Lagrangian coordinate system used to describe the gas---is present already in the one-dimensional setting, its implications are especially important in the multidimensional case where no computationally viable Lagrangian description of the gas is readily available. We introduce new "pseudo-Lagrangian" coordinates that allow us to overcome this difficulty, and we illustrate the utility of these coordinates in the setting of isentropic gas dynamics in two space dimensions.


title: High-Performance FPGA Implementation of Equivariant Adaptive Separation via Independence Algorithm for Independent Component Analysis

abstract: Independent Component Analysis (ICA) is a dimensionality reduction technique that can boost efficiency of machine learning models that deal with probability density functions, e.g. Bayesian neural networks. Algorithms that implement adaptive ICA converge slower than their nonadaptive counterparts, however, they are capable of tracking changes in underlying distributions of input features. This intrinsically slow convergence of adaptive methods combined with existing hardware implementations that operate at very low clock frequencies necessitate fundamental improvements in both algorithm and hardware design. This paper presents an algorithm that allows efficient hardware implementation of ICA. Compared to previous work, our FPGA implementation of adaptive ICA improves clock frequency by at least one order of magnitude and throughput by at least two orders of magnitude. Our proposed algorithm is not limited to ICA and can be used in various machine learning problems that use stochastic gradient descent optimization.


title: Information Landscape and Flux, Mutual Information Rate Decomposition and Entropy Production

abstract: We explore the dynamics of information systems. We show that the driving force for information dynamics is determined by both the information landscape and information flux which determines the equilibrium time reversible and the nonequilibrium time-irreversible behaviours of the system respectively. We further demonstrate that the mutual information rate between the two subsystems can be decomposed into the time-reversible and time-irreversible parts respectively, analogous to the information landscape-flux decomposition for dynamics. Finally, we uncover the intimate relation between the nonequilibrium thermodynamics in terms of the entropy production rates and the time-irreversible part of the mutual information rate. We demonstrate the above features by the dynamics of a bivariate Markov chain.


title: MPG - A Framework for Reasoning on 6 DOF Pose Uncertainty

abstract: Reasoning about the pose, i.e. position and orientation of objects is one of the cornerstones of robotic manipulation under uncertainty. In a number of joint research projects our group is developing a robotic perception system that perceives and models an unprepared kitchen scenario with many objects. Since no single sensor or measurement provides sufficient information, a technique is needed to fuse a number of uncertain estimates of the pose, i.e. estimates with a widely stretched probability density function ($pdf$). The most frequently used approaches to describe the $pdfs$ are sample based description and multivariate normal (Gaussian) distributions. Sample based descriptions in 6D can describe basically any type of $pdfs$, but they require a large number of samples and there are no analytic formulae to fuse several $pdfs$. For Gaussian distributions these formulae exist, but the Gaussian distributions are unimodal and don't model widely spread distributions well. In this paper we present a framework for probabilistic modeling of 6D poses that combines the expressive power of the sample based description with the conciseness and algorithmic power of the Gaussian models. As parameterization of the 6D poses we select the dual quaternions, i.e. any pose is represented by two quaternions. The orientation part of a pose is described by a unit quaternion. The translation part is described by a purely imaginary quaternion. A basic probability density function over the poses is constructed by selecting a tangent point on the 3D sphere representing unit quaternions and taking the Cartesian set product of the tangent space with the 3D space of translations. In this 6D Euclidean space a 6D Gaussian distribution is defined. Projecting this Gaussian back to the unit sphere and renormalizing induces a distribution over 6D poses, called a Projected Gaussian.


title: The Densities of Planets in Multiple Stellar Systems

abstract: We analyze the effect of companion stars on the bulk density of 29 planets orbiting 15 stars in the Kepler field. These stars have at least one stellar companion within 2", and the planets have measured masses and radii, allowing an estimate of their bulk density. The transit dilution by the companion star requires the planet radii to be revised upward, even if the planet orbits the primary star; as a consequence, the planetary bulk density decreases. We find that, if planets orbited a faint companion star, they would be more volatile-rich, and in several cases their densities would become unrealistically low, requiring large, inflated atmospheres or unusually large mass fractions in a H/He envelope. In addition, for planets detected in radial velocity data, the primary star has to be the host. We can exclude 14 planets from orbiting the companion star; the remaining 15 planets in seven planetary systems could orbit either the primary or the secondary star, and for five of these planets the decrease in density would be substantial even if they orbited the primary, since the companion is of almost equal brightness as the primary. Substantial follow-up work is required in order to accurately determine the radii of transiting planets. Of particular interest are small, rocky planets that may be habitable; a lower mean density might imply a more volatile-rich composition. Reliable radii, masses, and thus bulk densities will allow us to identify which small planets are truly Earth-like.


title: A causal framework for explaining the predictions of black-box sequence-to-sequence models

abstract: We interpret the predictions of any black-box structured input-structured output model around a specific input-output pair. Our method returns an "explanation" consisting of groups of input-output tokens that are causally related. These dependencies are inferred by querying the black-box model with perturbed inputs, generating a graph over tokens from the responses, and solving a partitioning problem to select the most relevant components. We focus the general approach on sequence-to-sequence problems, adopting a variational autoencoder to yield meaningful input perturbations. We test our method across several NLP sequence generation tasks.


title: Narrow-line Laser Cooling by Adiabatic Transfer

abstract: We propose and demonstrate a novel laser cooling mechanism applicable to particles with narrow-linewidth optical transitions. By sweeping the frequency of counter-propagating laser beams in a sawtooth manner, we cause adiabatic transfer back and forth between the ground state and a long-lived optically excited state. The time-ordering of these adiabatic transfers is determined by Doppler shifts, which ensures that the associated photon recoils are in the opposite direction to the particle's motion. This ultimately leads to a robust cooling mechanism capable of exerting large forces via a weak transition and with reduced reliance on spontaneous emission. We present a simple intuitive model for the resulting frictional force, and directly demonstrate its efficacy for increasing the total phase-space density of an atomic ensemble. We rely on both simulation and experimental studies using the 7.5~kHz linewidth $^1$S$_0$ to $^3$P$_1$ transition in $^{88}$Sr. The reduced reliance on spontaneous emission may allow this adiabatic sweep method to be a useful tool for cooling particles that lack closed cycling transitions, such as molecules.


title: Simple Classification using Binary Data

abstract: Binary, or one-bit, representations of data arise naturally in many applications, and are appealing in both hardware implementations and algorithm design. In this work, we study the problem of data classification from binary data and propose a framework with low computation and resource costs. We illustrate the utility of the proposed approach through stylized and realistic numerical experiments, and provide a theoretical analysis for a simple case. We hope that our framework and analysis will serve as a foundation for studying similar types of approaches.


title: Solitary waves in the Nonlinear Dirac Equation

abstract: In the present work, we consider the existence, stability, and dynamics of solitary waves in the nonlinear Dirac equation. We start by introducing the Soler model of self-interacting spinors, and discuss its localized waveforms in one, two, and three spatial dimensions and the equations they satisfy. We present the associated explicit solutions in one dimension and numerically obtain their analogues in higher dimensions. The stability is subsequently discussed from a theoretical perspective and then complemented with numerical computations. Finally, the dynamics of the solutions is explored and compared to its non-relativistic analogue, which is the nonlinear Schr{\"o}dinger equation. A few special topics are also explored, including the discrete variant of the nonlinear Dirac equation and its solitary wave properties, as well as the PT-symmetric variant of the model.


title: Efficient simulation of DC-DC switch-mode power converters by multirate partial differential equations

abstract: In this paper, Multirate Partial Differential Equations (MPDEs) are used for the efficient simulation of problems with 2-level pulsed excitations as they often occur in power electronics, e.g., DC-DC switch-mode converters. The differential equations describing the problem are reformulated as MPDEs which are solved by a Galerkin approach and time discretization. For the solution expansion two types of basis functions are proposed, namely classical Finite Element (FE) nodal functions and the recently introduced excitation-specific pulse width modulation (PWM) basis functions. The new method is applied to the example of a buck converter. Convergence, accuracy of the solution and computational efficiency of the method are numerically analyzed.


title: Emergence of dynamical dark energy from polynomial $f(R)$ theory in Palatini formalism

abstract: We consider FRW cosmology in $f(R)= R+ \gamma R^2+\delta R^3$ modified framework. The Palatini approach reduces its dynamics to the simple generalization of Friedmann equation. Thus we study the dynamics in two-dimensional phase space with some details. After reformulation of the model in the Einstein frame, it reduces to the FRW cosmological model with a homogeneous scalar field and vanishing kinetic energy term. This potential determines the running cosmological constant term as a function of the Ricci scalar. As a result we obtain the emergent dark energy parametrization from the covariant theory. We study also singularities of the model and demonstrate that in the Einstein frame some undesirable singularities disappear.


title: Elements of $C^*$-algebras Attaining Their Norm in a Finite-Dimensional Representation

abstract: We characterize the class of RFD $C^*$-algebras as those containing a dense subset of elements that attain their norm under a finite-dimensional representation. We show further that this subset is the whole space precisely when every irreducible representation of the $C^*$-algebra is finite-dimensional, which is equivalent to the $C^*$-algebra having no simple infinite-dimensional AF subquotient. We apply techniques from this proof to show the existence of elements in more general classes of $C^*$-algebras whose norms in finite-dimensional representations fit certain prescribed properties.


title: Mufasa:The strength and evolution of galaxy conformity in various tracers

abstract: We investigate galaxy conformity using the Mufasa cosmological hydrodynamical simulation. We show a bimodal distribution in galaxy colour with radius, albeit with too many low-mass quenched satellite galaxies compared to observations. Mufasa produces conformity in observed properties such as colour, sSFR, and HI content; i.e neighbouring galaxies have similar properties. We see analogous trends in other properties such as in environment, stellar age, H$_2$ content, and metallicity. We introduce quantifying conformity using $S(R)$, measuring the relative difference in upper and lower quartile properties of the neighbours. We show that low-mass and non-quenched haloes have weak conformity ($S(R)\leq 0.5$) extending to large projected radii $R$ in all properties, while high-mass and quenched haloes have strong conformity ($S(R)\sim 1$) that diminishes rapidly with $R$ and disappears at $R\geq 1$ Mpc. $S(R)$ is strongest for environment in low-mass haloes, and sSFR (or colour) in high-mass haloes, and is dominated by one-halo conformity with the exception of HI in small haloes. Metallicity shows a curious anti-conformity in massive haloes. Tracking the evolution of conformity for $z=0$ galaxies back in time shows that conformity broadly emerges as a late-time ($z\leq 1$) phenomenon. However, for fixed halo mass bins, conformity is fairly constant with redshift out to $z\geq 2$. These trends are consistent with the idea that strong conformity only emerges once haloes grow above Mufasa's quenching mass scale of $\sim 10^{12}M_\odot$. A quantitative measure of conformity in various properties, along with its evolution, thus represents a new and stringent test of the impact of quenching on environment within current galaxy formation models.


title: Robust Doubly Protected Estimators for Quantiles with Missing Data

abstract: Doubly protected estimators are widely used for estimating the population mean of an outcome Y from a sample where the response is missing in some individuals. To compensate for the missing responses, a vector X of covariates is observed at each individual, and the missing mechanism is assumed to be independent of the response, conditioned on X (missing at random). In recent years, many authors have moved from the mean to the median, and more generally, doubly protected estimators of the quantiles have been proposed, assuming a parametric regression model for the relationship between X and Y and a parametric form for the propensity score. In this work, we present doubly protected estimators for the quantiles that are also robust, in the sense that they are resistant to the presence of outliers in the sample. We also flexibilize the model for the relationship between X and Y . Thus we present robust doubly protected estimators for the quantiles of the response in the presence of missing observations, postulating a semiparametric regression model for the relationship between the response and the covariates and a parametric model for the propensity score.


title: Significance of Disk Failure Prediction in Datacenters

abstract: Modern datacenters assemble a very large number of disk drives under a single roof. Even if economic and technical factors where to make individual drives more reliable (which is not at all clear, given the commoditization of the technology), their sheer numbers combined with their ever increasing utilization in a well-balanced design makes achieving storage reliability a major challenge. In this paper, we assess the challenge of storage system reliability in the modern datacenter, and demonstrate how good disk failure prediction models can significantly improve the reliability of such systems.


title: A tensor network study of the complete ground state phase diagram of the spin-1 bilinear-biquadratic Heisenberg model on the square lattice

abstract: Using infinite projected entangled pair states, we study the ground state phase diagram of the spin-1 bilinear-biquadratic Heisenberg model on the square lattice directly in the thermodynamic limit. We find an unexpected partially nematic partially magnetic phase in between the antiferroquadrupolar and ferromagnetic regions. Furthermore, we describe all observed phases and discuss the nature of the phase transitions involved.


title: Convergence and normal continuity analysis of non-stationary subdivision schemes near extraordinary vertices and faces

abstract: Convergence and normal continuity analysis of a bivariate non-stationary (level-dependent) subdivision scheme for 2-manifold meshes with arbitrary topology is still an open issue. Exploiting ideas from the theory of asymptotically equivalent subdivision schemes, in this paper we derive new sufficient conditions for establishing convergence and normal continuity of any rotationally symmetric, non-stationary, subdivision scheme near an extraordinary vertex/face.


title: Renormalized Reduced Order Models with Memory for Long Time Prediction

abstract: We examine the challenging problem of constructing reduced models for the long time prediction of systems where there is no timescale separation between the resolved and unresolved variables. In previous work we focused on the case where there was only transfer of activity (e.g. energy, mass) from the resolved to the unresolved variables. Here we investigate the much more difficult case where there is two-way transfer of activity between the resolved and unresolved variables. Like in the case of activity drain out of the resolved variables, even if one starts with an exact formalism, like the Mori-Zwanzig (MZ) formalism, the constructed reduced models can become unstable. We show how to remedy this situation by using dynamic information from the full system to renormalize the MZ reduced models. In addition to being stabilized, the renormalized models can be accurate for very long times. We use the Korteweg-de Vries equation to illustrate the approach. The coefficients of the renormalized models exhibit rich structure, including algebraic time dependence and incomplete similarity.


title: Periodically Driven Array of Single Rydberg Atoms

abstract: An array of single Rydberg atoms driven by a frequency modulated light field is studied. The periodic modulation effectively modifies the Rabi coupling, leading to unprecedented dynamics in the presence of Rydberg-Rydberg interactions. They include state dependent population trapping, the Rydberg blockade for small and anti-blockades at large interaction strengths. Interestingly, the Schrieffer-wolf transformation reveals a fundamental process in Rydberg gases, correlated Rabi oscillations, arising from the long-range interactions, provides an alternative depiction for Rydberg blockade and it exhibits a nontrivial behaviour in the presence of periodic modulation. The dynamical localization of a many body configuration in a driven Rydberg-lattice is discussed.


title: A 3D DLM/FD method for simulating the motion of spheres in a bounded shear flow of Oldroyd-B fluids

abstract: We present a novel distributed Lagrange multiplier/fictitious domain (DLM/FD) method for simulating fluid-particle interaction in Oldroyd-B fluids under creeping conditions. The results concerning two ball interaction in a three dimensional (3D) bounded shear flow are obtained for Weissenberg numbers up to 1 . The pass and return trajectories of the two ball mass centers are similar to those in a Newtonian fluid; but they lose the symmetry due to the effect of elastic force arising from viscoelastic fluids. A tumbling chain of two balls (a dipole) may occur, depending on the value of the Weissenberg number and the initial vertical displacement of the ball mass center to the middle plane between two walls.


title: Non-Gaussian Limit Theorem for Non-Linear Langevin Equations Driven by L\'evy Noise

abstract: In this paper, we study the small noise behaviour of solutions of a non-linear second order Langevin equation $\ddot x^\varepsilon_t +|\dot x^\varepsilon_t|^\beta=\dot Z^\varepsilon_{\varepsilon t}$, $\beta\in\mathbb R$, driven by symmetric non-Gaussian L\'evy processes $Z^\varepsilon$. This equation describes the dynamics of a one-degree-of-freedom mechanical system subject to non-linear friction and noisy vibrations. For a compound Poisson noise, the process $x^\varepsilon$ on the macroscopic time scale $t/\varepsilon$ has a natural interpretation as a non-linear filter which responds to each single jump of the driving process. We prove that a system driven by a general symmetric L\'evy noise exhibits essentially the same asymptotic behaviour under the principal condition $\alpha+2\beta<4$, where $\alpha\in [0,2]$ is the ``uniform'' Blumenthal--Getoor index of the family $\{Z^\varepsilon\}_{\varepsilon>0}$.