title: Accelerating DNN Training in Wireless Federated Edge Learning Systems

abstract: Training task in classical machine learning models, such as deep neural networks, is generally implemented at a remote cloud center for centralized learning, which is typically time-consuming and resource-hungry. It also incurs serious privacy issue and long communication latency since a large amount of data are transmitted to the centralized node. To overcome these shortcomings, we consider a newly-emerged framework, namely federated edge learning, to aggregate local learning updates at the network edge in lieu of users' raw data. Aiming at accelerating the training process, we first define a novel performance evaluation criterion, called learning efficiency. We then formulate a training acceleration optimization problem in the CPU scenario, where each user device is equipped with CPU. The closed-form expressions for joint batchsize selection and communication resource allocation are developed and some insightful results are highlighted. Further, we extend our learning framework to the GPU scenario. The optimal solution in this scenario is manifested to have the similar structure as that of the CPU scenario, recommending that our proposed algorithm is applicable in more general systems. Finally, extensive experiments validate the theoretical analysis and demonstrate that the proposed algorithm can reduce the training time and improve the learning accuracy simultaneously.


title: Extended Uncertainty Principle for Rindler and cosmological horizons

abstract: We find exact formulas for the Extended Uncertainty Principle (EUP) for the Rindler and Friedmann horizons and show that they can be expanded to obtain asymptotic forms known from the previous literature. We calculate the corrections to Hawking temperature and Bekenstein entropy of a black hole in the universe due to Rindler and Friedmann horizons. The effect of the EUP is similar to the canonical corrections of thermal fluctuations and so it rises the entropy signalling further loss of information.


title: Coulomb final state interaction in heavy ion collisions for Levy sources

abstract: Investigation of momentum space correlations of particles produced in high energy reactions requires taking final state interactions into account, a crucial point of any such analysis. Coulomb interaction between charged particles is the most important such effect. In small systems like those created in e+e- or p+p collisions, the so-called Gamow factor (valid for a point-like particle source) gives an acceptable description of the Coulomb interaction. However, in larger systems such as central or mid-central heavy ion collisions, more involved approaches are needed. In this paper we investigate the Coulomb final state interaction for Levy-type source functions that were recently shown to be of much interest for a refined description of the space-time picture of particle production in heavy-ion collisions.


title: An illustration of the risk of borrowing information via a shared likelihood

abstract: A concrete, stylized example illustrates that inferences may be degraded, rather than improved, by incorporating supplementary data via a joint likelihood. In the example, the likelihood is assumed to be correctly specified, as is the prior over the parameter of interest; all that is necessary for the joint modeling approach to suffer is misspecification of the prior over a nuisance parameter.


title: A Convolutional Cost-Sensitive Crack Localization Algorithm for Automated and Reliable RC Bridge Inspection

abstract: Bridges are an essential part of the transportation infrastructure and need to be monitored periodically. Visual inspections by dedicated teams have been one of the primary tools in structural health monitoring (SHM) of bridge structures. However, such conventional methods have certain shortcomings. Manual inspections may be challenging in harsh environments and are commonly biased in nature. In the last decade, camera-equipped unmanned aerial vehicles (UAVs) have been widely used for visual inspections; however, the task of automatically extracting useful information from raw images is still challenging. In this paper, a deep learning semantic segmentation framework is proposed to automatically localize surface cracks. Due to the high imbalance of crack and background classes in images, different strategies are investigated to improve performance and reliability. The trained models are tested on real-world crack images showing impressive robustness in terms of the metrics defined by the concepts of precision and recall. These techniques can be used in SHM of bridges to extract useful information from the unprocessed images taken from UAVs.


title: Network Pruning via Transformable Architecture Search

abstract: Network pruning reduces the computation costs of an over-parameterized network without performance damage. Prevailing pruning algorithms pre-define the width and depth of the pruned networks, and then transfer parameters from the unpruned network to pruned networks. To break the structure limitation of the pruned networks, we propose to apply neural architecture search to search directly for a network with flexible channel and layer sizes. The number of the channels/layers is learned by minimizing the loss of the pruned networks. The feature map of the pruned network is an aggregation of K feature map fragments (generated by K networks of different sizes), which are sampled based on the probability distribution.The loss can be back-propagated not only to the network weights, but also to the parameterized distribution to explicitly tune the size of the channels/layers. Specifically, we apply channel-wise interpolation to keep the feature map with different channel sizes aligned in the aggregation procedure. The maximum probability for the size in each distribution serves as the width and depth of the pruned network, whose parameters are learned by knowledge transfer, e.g., knowledge distillation, from the original networks. Experiments on CIFAR-10, CIFAR-100 and ImageNet demonstrate the effectiveness of our new perspective of network pruning compared to traditional network pruning algorithms. Various searching and knowledge transfer approaches are conducted to show the effectiveness of the two components. Code is at: https://github.com/D-X-Y/NAS-Projects.


title: Meta-GNN: On Few-shot Node Classification in Graph Meta-learning

abstract: Meta-learning has received a tremendous recent attention as a possible approach for mimicking human intelligence, i.e., acquiring new knowledge and skills with little or even no demonstration. Most of the existing meta-learning methods are proposed to tackle few-shot learning problems such as image and text, in rather Euclidean domain. However, there are very few works applying meta-learning to non-Euclidean domains, and the recently proposed graph neural networks (GNNs) models do not perform effectively on graph few-shot learning problems. Towards this, we propose a novel graph meta-learning framework -- Meta-GNN -- to tackle the few-shot node classification problem in graph meta-learning settings. It obtains the prior knowledge of classifiers by training on many similar few-shot learning tasks and then classifies the nodes from new classes with only few labeled samples. Additionally, Meta-GNN is a general model that can be straightforwardly incorporated into any existing state-of-the-art GNN. Our experiments conducted on three benchmark datasets demonstrate that our proposed approach not only improves the node classification performance by a large margin on few-shot learning problems in meta-learning paradigm, but also learns a more general and flexible model for task adaption.


title: Price of Dependence: Stochastic Submodular Maximization with Dependent Items

abstract: In this paper, we study the stochastic submodular maximization problem with dependent items subject to packing constraints such as matroid and knapsack constraints. The input of our problem is a finite set of items, and each item is in a particular state from a set of possible states. After picking an item, we are able to observe its state. We assume a monotone and submodular utility function over items and states, and our objective is to select a group of items adaptively so as to maximize the expected utility. Previous studies on stochastic submodular maximization often assume that items' states are independent, however, this assumption may not hold in general. This motivates us to study the stochastic submodular maximization problem with dependent items. We first introduce the concept of \emph{degree of independence} to capture the degree to which one item's state is dependent on others'. Then we propose a non-adaptive policy that approximates the optimal adaptive policy within a factor of $\alpha(1-e^{-\frac{\kappa}{2}+\frac{\kappa}{18m^2}}-\frac{\kappa+2}{3m\kappa})$ where the value of $\alpha$ is depending on the type of constraints, e.g., $\alpha=1$ for matroid constraint, $\kappa$ is the degree of independence, e.g., $\kappa=1$ for independent items, and $m$ is the number of items. We also analyze the adaptivity gap, i.e., the ratio of the values of best adaptive policy and best non-adaptive policy, of our problem with prefix-closed constraints.


title: Weak solutions obtained by the vortex method for the 2D Euler equations are Lagrangian and conserve the energy

abstract: We discuss the Lagrangian property and the conservation of the kinetic energy for solutions of the 2D incompressible Euler equations. Existence of Lagrangian solutions is known when the initial vorticity is in $L^p$ with $1\leq p\leq \infty$. Moreover, if $p\geq 3/2$ all weak solutions are conservative. In this work we prove that solutions obtained via the vortex method are Lagrangian, and that they are conservative if $p>1$.


title: Statistical Assertions for Validating Patterns and Finding Bugs in Quantum Programs

abstract: In support of the growing interest in quantum computing experimentation, programmers need new tools to write quantum algorithms as program code. Compared to debugging classical programs, debugging quantum programs is difficult because programmers have limited ability to probe the internal states of quantum programs; those states are difficult to interpret even when observations exist; and programmers do not yet have guidelines for what to check for when building quantum programs. In this work, we present quantum program assertions based on statistical tests on classical observations. These allow programmers to decide if a quantum program state matches its expected value in one of classical, superposition, or entangled types of states. We extend an existing quantum programming language with the ability to specify quantum assertions, which our tool then checks in a quantum program simulator. We use these assertions to debug three benchmark quantum programs in factoring, search, and chemistry. We share what types of bugs are possible, and lay out a strategy for using quantum programming patterns to place assertions and prevent bugs.


title: Random Norming Aids Analysis of Non-linear Regression Models with Sequential Informative Dose Selection

abstract: A two-stage adaptive optimal design is an attractive option for increasing the efficiency of clinical trials. In these designs, based on interim data, the locally optimal dose is chosen for further exploration, which induces dependencies between data from the two stages. When the maximum likelihood estimator (MLE) is used under nonlinear regression models with independent normal errors in a pilot study where the first stage sample size is fixed, and the second stage sample size is large, the Fisher information fails to normalize the estimator adequately asymptotically, because of dependencies. In this situation, we present three alternative random information measures and show that they provide better normalization of the MLE asymptotically. The performance of random information measures is investigated in simulation studies, and the results suggest that the observed information performs best when the sample size is small.


title: Site percolation and isoperimetric inequalities for plane graphs

abstract: We use isoperimetric inequalities combined with a new technique to prove upper bounds for the site percolation threshold of plane graphs with given minimum degree conditions. In the process we prove tight new isoperimetric bounds for certain classes of hyperbolic graphs. This establishes the vertex isoperimetric constant for all triangular and square hyperbolic lattices, answering a question of Lyons and Peres. We prove that plane graphs of minimum degree at least $7$ have site percolation threshold bounded away from $1/2$, which was conjectured by Benjamini and Schramm, and make progress on a conjecture of Angel, Benjamini and Horesh that the critical probability is at most $1/2$ for plane triangulations of minimum degree $6$. We prove additional bounds for stronger minimum degree conditions, and for graphs without triangular faces.


title: Two different analyses on derivation of PYTHAGORAS Theorem (569-479 BC): Discrete continuum states

abstract: We propose two different derivations of Pythagoras Theorem and apply the same to study discrete and continuum states.


title: A new algorithm that generates the image of the attractor of a generalized iterated function system

abstract: We provide a new algorithm (called the grid algorithm) designed to generate the image of the attractor of a generalized iterated function system on a finite dimensional space and we compare it with the deterministic algorithm regarding generalized iterated function systems presented by P. Jaros, L. Maslanka and F. Strobin in [Algorithms generating images of attractors of generalized iterated function systems, Numer. Algorithms, 73 (2016), 477-499].


title: Verification and Synthesis of Symmetric Uni-Rings for Leads-To Properties

abstract: This paper investigates the verification and synthesis of parameterized protocols that satisfy leadsto properties $R \leadsto Q$ on symmetric unidirectional rings (a.k.a. uni-rings) of deterministic and constant-space processes under no fairness and interleaving semantics, where $R$ and $Q$ are global state predicates. First, we show that verifying $R \leadsto Q$ for parameterized protocols on symmetric uni-rings is undecidable, even for deterministic and constant-space processes, and conjunctive state predicates. Then, we show that surprisingly synthesizing symmetric uni-ring protocols that satisfy $R \leadsto Q$ is actually decidable. We identify necessary and sufficient conditions for the decidability of synthesis based on which we devise a sound and complete polynomial-time algorithm that takes the predicates $R$ and $Q$, and automatically generates a parameterized protocol that satisfies $R \leadsto Q$ for unbounded (but finite) ring sizes. Moreover, we present some decidability results for cases where leadsto is required from multiple distinct $R$ predicates to different $Q$ predicates. To demonstrate the practicality of our synthesis method, we synthesize some parameterized protocols, including agreement and parity protocols.


title: Deep Drone Racing: From Simulation to Reality with Domain Randomization

abstract: Dynamically changing environments, unreliable state estimation, and operation under severe resource constraints are fundamental challenges that limit the deployment of small autonomous drones. We address these challenges in the context of autonomous, vision-based drone racing in dynamic environments. A racing drone must traverse a track with possibly moving gates at high speed. We enable this functionality by combining the performance of a state-of-the-art planning and control system with the perceptual awareness of a convolutional neural network (CNN). The resulting modular system is both platform- and domain-independent: it is trained in simulation and deployed on a physical quadrotor without any fine-tuning. The abundance of simulated data, generated via domain randomization, makes our system robust to changes of illumination and gate appearance. To the best of our knowledge, our approach is the first to demonstrate zero-shot sim-to-real transfer on the task of agile drone flight. We extensively test the precision and robustness of our system, both in simulation and on a physical platform, and show significant improvements over the state of the art.


title: Non-Abelian Statistics in one dimension: topological momentum spacings and SU(2) level $k$ fusion rules

abstract: We use a family of critical spin chain models discovered recently by one of us [M. Greiter, Mapping of Parent Hamiltonians, Springer, Berlin/Heidelberg 2011] to propose and elaborate that non-Abelian, SU(2) level $k=2S$ anyon statistics manifests itself in one dimension through topological selection rules for fractional shifts in the spacings of linear momenta, which yield an internal Hilbert space of, in the thermodynamic limit degenerate states. These shifts constitute the equivalent to the fractional shifts in the relative angular momenta of anyons in two dimensions. We derive the rules first for Ising anyons, and then generalize them to SU(2) level $k$ anyons. We establish a one-to-one correspondence between the topological choices for the momentum spacings and the fusion rules of spin \half spinons in the SU(2) level $k$ Wess--Zumino--Witten model, where the internal Hilbert space is spanned by the manifold of allowed fusion trees in the Bratelli diagrams. Finally, we show that the choices in the fusion trees may be interpreted as the choices between different domain walls between the $2S+1$ possible, degenerate dimer configurations of the spin $S$ chains at the multicritical point.


title: 2-factors with k cycles in Hamiltonian graphs

abstract: A well known generalisation of Dirac's theorem states that if a graph $G$ on $n\ge 4k$ vertices has minimum degree at least $n/2$ then $G$ contains a $2$-factor consisting of exactly $k$ cycles. This is easily seen to be tight in terms of the bound on the minimum degree. However, if one assumes in addition that $G$ is Hamiltonian it has been conjectured that the bound on the minimum degree may be relaxed. This was indeed shown to be true by S\'ark\"ozy. In subsequent papers, the minimum degree bound has been improved, most recently to $(2/5+\varepsilon)n$ by DeBiasio, Ferrara, and Morris. On the other hand no lower bounds close to this are known, and all papers on this topic ask whether the minimum degree needs to be linear. We answer this question, by showing that the required minimum degree for large Hamiltonian graphs to have a $2$-factor consisting of a fixed number of cycles is sublinear in $n.$


title: On modelling the emergence of logical thinking

abstract: Recent progress in machine learning techniques have revived interest in building artificial general intelligence using these particular tools. There has been a tremendous success in applying them for narrow intellectual tasks such as pattern recognition, natural language processing and playing Go. The latter application vastly outperforms the strongest human player in recent years. However, these tasks are formalized by people in such ways that it has become "easy" for automated recipes to find better solutions than humans do. In the sense of John Searle's Chinese Room Argument, the computer playing Go does not actually understand anything from the game. Thinking like a human mind requires to go beyond the curve fitting paradigm of current systems. There is a fundamental limit to what they can achieve currently as only very specific problem formalization can increase their performances in particular tasks. In this paper, we argue than one of the most important aspects of the human mind is its capacity for logical thinking, which gives rise to many intellectual expressions that differentiate us from animal brains. We propose to model the emergence of logical thinking based on Piaget's theory of cognitive development.


title: Edge states and corner modes in second-order topological phononic crystal plates

abstract: We realize an elastic second-order topological insulator hosting both one-dimensional gapped edge states and zero-dimensional in-gap corner modes in the double-sided pillared phononic crystal plates with square lattice. Changing the width of two neighbor pillars breaks the inversion symmetry and induces the band inversion to emulate the quantum spin Hall effect where the gapless edge states are obtained. Further breaking the space-symmetry at interface, the gapless edge states are gapped and inducing the edge topological transitions and then giving rise to the zero-dimensional in-gap corner modes. Our work offers a novel way for elastic wave trapping and robustly guiding.


title: Shear Thickening of Dense Suspensions: The Role of Friction

abstract: Shear thickening of particle suspensions is characterized by a transition between lubricated and frictional contacts between the particles. Using 3D numerical simulations, we study how the inter-particle friction coefficient influences the effective macroscopic friction coefficient and hence the microstructure and rheology of dense shear thickening suspensions. We propose expressions for effective friction coefficient in terms of distance to jamming for varying shear stresses and particle friction coefficient values. We find effective friction coefficient to be rather insensitive to interparticle friction, which is perhaps surprising but agrees with recent theory and experiments.


title: Real-Space Renormalization for disordered systems at the level of Large Deviations

abstract: The real-space renormalization procedures on hierarchical lattices have been much studied for many disordered systems in the past at the level of their typical fluctuations. In the present paper, the goal is to analyze instead the renormalization flows for the tails of probability distributions in order to extract the scalings of their large deviations and the tails behaviors of the corresponding rate functions. We focus on the renormalization rule for the ground-state energy of the Directed Polymer model in a random medium, and study the various renormalization flows that can emerge for the tails as a function of the tails of the initial condition.

categories: cond-mat.dis-nn cond-mat.stat-mech math-ph math.MP


title: Acceleration of SVRG and Katyusha X by Inexact Preconditioning

abstract: Empirical risk minimization is an important class of optimization problems with many popular machine learning applications, and stochastic variance reduction methods are popular choices for solving them. Among these methods, SVRG and Katyusha X (a Nesterov accelerated SVRG) achieve fast convergence without substantial memory requirement. In this paper, we propose to accelerate these two algorithms by \textit{inexact preconditioning}, the proposed methods employ \textit{fixed} preconditioners, although the subproblem in each epoch becomes harder, it suffices to apply \textit{fixed} number of simple subroutines to solve it inexactly, without losing the overall convergence. As a result, this inexact preconditioning strategy gives provably better iteration complexity and gradient complexity over SVRG and Katyusha X. We also allow each function in the finite sum to be nonconvex while the sum is strongly convex. In our numerical experiments, we observe an on average $8\times$ speedup on the number of iterations and $7\times$ speedup on runtime.


title: Digital Normativity: A challenge for human subjectivization and free will

abstract: Over the past decade, artificial intelligence has demonstrated its efficiency in many different applications and a huge number of algorithms have become central and ubiquitous in our life. Their growing interest is essentially based on their capability to synthesize and process large amounts of data, and to help humans making decisions in a world of increasing complexity. Yet, the effectiveness of algorithms in bringing more and more relevant recommendations to humans may start to compete with human-alone decisions based on values other than pure efficacy. Here, we examine this tension in light of the emergence of several forms of digital normativity, and analyze how this normative role of AI may influence the ability of humans to remain subject of their life. The advent of AI technology imposes a need to achieve a balance between concrete material progress and progress of the mind to avoid any form of servitude. It has become essential that an ethical reflection accompany the current developments of intelligent algorithms beyond the sole question of their social acceptability. Such reflection should be anchored where AI technologies are being developed as well as in educational programs where their implications can be explained.


title: Consistent Dynamic Mode Decomposition

abstract: We propose a new method for computing Dynamic Mode Decomposition (DMD) evolution matrices, which we use to analyze dynamical systems. Unlike the majority of existing methods, our approach is based on a variational formulation consisting of data alignment penalty terms and constitutive orthogonality constraints. Our method does not make any assumptions on the structure of the data or their size, and thus it is applicable to a wide range of problems including non-linear scenarios or extremely small observation sets. In addition, our technique is robust to noise that is independent of the dynamics and it does not require input data to be sequential. Our key idea is to introduce a regularization term for the forward and backward dynamics. The obtained minimization problem is solved efficiently using the Alternating Method of Multipliers (ADMM) which requires two Sylvester equation solves per iteration. Our numerical scheme converges empirically and is similar to a provably convergent ADMM scheme. We compare our approach to various state-of-the-art methods on several benchmark dynamical systems.


title: Aligned SICs and embedded tight frames in even dimensions

abstract: Alignment is a geometric relation between pairs of Weyl-Heisenberg SICs, one in dimension $d$ and another in dimension $d(d-2)$, manifesting a well-founded conjecture about a number-theoretical connection between the SICs. In this paper, we prove that if $d$ is even, the SIC in dimension $d(d-2)$ of an aligned pair can be partitioned into $(d-2)^2$ tight $d^2$-frames of rank $d(d-1)/2$ and, alternatively, into $d^2$ tight $(d-2)^2$-frames of rank $(d-1)(d-2)/2$. The corresponding result for odd $d$ is already known, but the proof for odd $d$ relies on results which are not available for even $d$. We develop methods that allow us to overcome this issue. In addition, we provide a relatively detailed study of parity operators in the Clifford group, emphasizing differences in the theory of parity operators in even and odd dimensions and discussing consequences due to such differences. In a final section, we study implications of alignment for the symmetry of the SIC.


title: Rotating binary Bose-Einstein condensates and vortex clusters in quantum droplets

abstract: Quantum droplets may form out of a gaseous Bose-Einstein condensate, stabilized by quantum fluctuations beyond mean field. We show that multiple singly-quantized vortices may form in these droplets at moderate angular momenta in two dimensions. Droplets carrying these precursors of an Abrikosov lattice remain self-bound for certain timescales after switching off an initial harmonic confinement. Furthermore, we examine how these vortex-carrying droplets can be formed in a more pertubation-resistant setting, by starting from a rotating binary Bose-Einstein condensate and inducing a metastable persistent current via a non-monotonic trapping potential.


title: Photoionization of Molecular Endohedrals

abstract: We calculate the photoionization cross-section of a molecular endohedral. We limit ourselves to two-atomic molecules. The consideration is much more complex than for atomic endohedrals because the system even for almost spherical fullerenes has only cylindrical instead of spherical symmetry. On the other hand, molecular endohedral is more interesting since the interelectron interaction in molecules is relatively stronger than in similar atoms. We present here results of calculations of molecular hydrogen stuffed inside almost spherical fullerene. For comparison, we perform calculations also for atomic endohedral with Helium inside fullerene. The results are obtained both in the single-electron Hartree-Fock approximation and with account of multi-electron correlations in the frame of so-called random phase approximation with exchange. The presence of the fullerenes shell results in prominent oscillations in the endohedrals photoionization cross section. The role of interelectron correlations becomes clear by comparing HF and RPAE results for molecular and atomic endohedral on the one side with that for corresponding isolated molecule and atom on the other.


title: Near-infrared optical properties and proposed phase-change usefulness of transition metal disulfides

abstract: The development of photonic integrated circuits would benefit from a wider selection of materials that can strongly-control near-infrared (NIR) light. Transition metal dichalcogenides (TMDs) have been explored extensively for visible spectrum opto-electronics, but the NIR properties of these layered materials have been less-studied. The measurement of optical constants is the foremost step to qualify TMDs for use in NIR photonics. Here we measure the complex optical constants for select sulfide TMDs (bulk crystals of MoS2, TiS2 and ZrS2) via spectroscopic ellipsometry in the visible-to-NIR range. Through Mueller matrix measurements and generalized ellipsometry, we explicitly measure the direction of the ordinary optical axis. We support our measurements with density functional theory (DFT) calculations, which agree with our measurements and predict giant birefringence. We further propose that TMDs could find use as photonic phase-change materials, by designing alloys that are thermodynamically adjacent to phase boundaries between competing crystal structures, to realize martensitic (i.e. displacive, order-order) switching.


title: PAC structures in nutshell

abstract: An expository paper written down after RIMS Model Theory Workshop 2018. To appear in RIMS Kokyuroku.


title: Explicit a posteriori local error estimation for FEM solutions

abstract: For the finite element solution of Poisson's equation, a local a posteriori error estimation based on the Hypercircle method is proposed. Even for the solution of Poisson's equation without the $H^2$ regularity, this method can provide explicit local error estimation. The efficiency of the proposed method is demonstrated by numerical experiments for the boundary value problem of Poisson's equation defined on the 2D and 3D domains.


title: Cross-chain Deals and Adversarial Commerce

abstract: Modern distributed data management systems face a new challenge: how can autonomous, mutually-distrusting parties cooperate safely and effectively? Addressing this challenge brings up questions familiar from classical distributed systems: how to combine multiple steps into a single atomic action, how to recover from failures, and how to synchronize concurrent access to data. Nevertheless, each of these issues requires rethinking when participants are autonomous and potentially adversarial. We propose the notion of a *cross-chain deal*, a new way to structure complex distributed computations that manage assets in an adversarial setting. Deals are inspired by classical atomic transactions, but are necessarily different, in important ways, to accommodate the decentralized and untrusting nature of the exchange. We describe novel safety and liveness properties, along with two alternative protocols for implementing cross-chain deals in a system of independent blockchain ledgers. One protocol, based on synchronous communication, is fully decentralized, while the other, based on semi-synchronous communication, requires a globally shared ledger.


title: A consistent and comprehensive computational approach for general Fluid-Structure-Contact Interaction problems

abstract: We present a consistent approach that allows to solve challenging general nonlinear fluid-structure-contact interaction (FSCI) problems. The underlying continuous formulation includes both "no-slip" fluid-structure interaction as well as frictionless contact between multiple elastic bodies. The respective interface conditions in normal and tangential orientation and especially the role of the fluid stress within the region of closed contact are discussed for the general problem of FSCI. To ensure continuity of the tangential constraints from no-slip to frictionless contact, a transition is enabled by using the general Navier condition with varying slip length. Moreover, the fluid stress in the contact zone is obtained by an extension approach as it plays a crucial role for the lift-off behavior of contacting bodies. With the given continuity of the spatially continuous formulation, continuity of the discrete problem (which is essential for the convergence of Newton's method) is reached naturally. As topological changes of the fluid domain are an inherent challenge in FSCI configurations, a non-interface fitted Cut Finite Element Method (CutFEM) is applied to discretize the fluid domain. All interface conditions, that is the `no-slip' FSI, the general Navier condition, and frictionless contact are incorporated using Nitsche based methods, thus retaining the continuity and consistency of the model. To account for the strong interaction between the fluid and solid discretization, the overall coupled discrete system is solved monolithically. Numerical examples of varying complexity are presented to corroborate the developments.


title: Kuroda's formula and arithmetic statistics

abstract: Kuroda's formula relates the class number of a multi-quadratic number field $K$ to the class numbers of its quadratic subfields $k_i$. A key component in this formula is the unit group index $Q(K) = [\mathcal{O}_{K}^{\times}: \prod_i\mathcal{O}_{k_i}^{\times}]$. We study how $Q(K)$ behaves on average in certain natural families of totally real biquadratic fields $K$ parametrized by prime numbers.


title: A Predictive Model for Steady-State Multiphase Pipe Flow: Machine Learning on Lab Data

abstract: Engineering simulators used for steady-state multiphase pipe flows are commonly utilized to predict pressure drop. Such simulators are typically based on either empirical correlations or first-principles mechanistic models. The simulators allow evaluating the pressure drop in multiphase pipe flow with acceptable accuracy. However, the only shortcoming of these correlations and mechanistic models is their applicability. In order to extend the applicability and the accuracy of the existing accessible methods, a method of pressure drop calculation in the pipeline is proposed. The method is based on well segmentation and calculation of the pressure gradient in each segment using three surrogate models based on Machine Learning algorithms trained on a representative lab data set from the open literature. The first model predicts the value of a liquid holdup in the segment, the second one determines the flow pattern, and the third one is used to estimate the pressure gradient. To build these models, several ML algorithms are trained such as Random Forest, Gradient Boosting Decision Trees, Support Vector Machine, and Artificial Neural Network, and their predictive abilities are cross-compared. The proposed method for pressure gradient calculation yields $R^2 = 0.95$ by using the Gradient Boosting algorithm as compared with $R^2 = 0.92$ in case of Mukherjee and Brill correlation and $R^2 = 0.91$ when a combination of Ansari and Xiao mechanistic models is utilized. The method for pressure drop prediction is also validated on three real field cases. Validation indicates that the proposed model yields the following coefficients of determination: $R^2 = 0.806, 0.815$ and 0.99 as compared with the highest values obtained by commonly used techniques: $R^2 = 0.82$ (Beggs and Brill correlation), $R^2 = 0.823$ (Mukherjee and Brill correlation) and $R^2 = 0.98$ (Beggs and Brill correlation).

categories: physics.data-an cs.LG physics.comp-ph physics.flu-dyn


title: Adversarially Robust Distillation

abstract: Knowledge distillation is effective for producing small, high-performance neural networks for classification, but these small networks are vulnerable to adversarial attacks. This paper studies how adversarial robustness transfers from teacher to student during knowledge distillation. We find that a large amount of robustness may be inherited by the student even when distilled on only clean images. Second, we introduce Adversarially Robust Distillation (ARD) for distilling robustness onto student networks. In addition to producing small models with high test accuracy like conventional distillation, ARD also passes the superior robustness of large networks onto the student. In our experiments, we find that ARD student models decisively outperform adversarially trained networks of identical architecture in terms of robust accuracy, surpassing state-of-the-art methods on standard robustness benchmarks. Finally, we adapt recent fast adversarial training methods to ARD for accelerated robust distillation.


title: Co-theory of sorted profinite groups for PAC structures

abstract: We achieve several results. First, we develop a variant of the theory of absolute Galois groups in the context of many sorted structures. Second, we provide a method for coding absolute Galois groups of structures, so they can be interpreted in some monster model with an additional predicate. Third, we prove the "Weak Independence Theorem" for PAC substructures of an ambient structure with nfcp and the property B(3). Fourth, we describe Kim-dividing in these PAC substructures and show several results related to the SOPn hierarchy. Fifth, we characterize the algebraic closure in PAC structures


title: How to factor 2048 bit RSA integers in 8 hours using 20 million noisy qubits

abstract: We significantly reduce the cost of factoring integers and computing discrete logarithms in finite fields on a quantum computer by combining techniques from Shor 1994, Griffiths-Niu 1996, Zalka 2006, Fowler 2012, Eker{\aa}-H{\aa}stad 2017, Eker{\aa} 2017, Eker{\aa} 2018, Gidney-Fowler 2019, Gidney 2019. We estimate the approximate cost of our construction using plausible physical assumptions for large-scale superconducting qubit platforms: a planar grid of qubits with nearest-neighbor connectivity, a characteristic physical gate error rate of $10^{-3}$, a surface code cycle time of 1 microsecond, and a reaction time of 10 microseconds. We account for factors that are normally ignored such as noise, the need to make repeated attempts, and the spacetime layout of the computation. When factoring 2048 bit RSA integers, our construction's spacetime volume is a hundredfold less than comparable estimates from earlier works (Van Meter et al. 2009, Jones et al. 2010, Fowler et al. 2012, Gheorghiu et al. 2019). In the abstract circuit model (which ignores overheads from distillation, routing, and error correction) our construction uses $3 n + 0.002 n \lg n$ logical qubits, $0.3 n^3 + 0.0005 n^3 \lg n$ Toffolis, and $500 n^2 + n^2 \lg n$ measurement depth to factor $n$-bit RSA integers. We quantify the cryptographic implications of our work, both for RSA and for schemes based on the DLP in finite fields.


title: Approximation schemes for the generalized extensible bin packing problem

abstract: We present a new generalization of the extensible bin packing with unequal bin sizes problem. In our generalization the cost of exceeding the bin size depends on the index of the bin and not only on the amount in which the size of the bin is exceeded. This generalization does not satisfy the assumptions on the cost function that were used to present the existing polynomial time approximation scheme (PTAS) for the extensible bin packing with unequal bin sizes problem. In this work, we show the existence of an efficient PTAS (EPTAS) for this new generalization and thus in particular we improve the earlier PTAS for the extensible bin packing with unequal bin sizes problem into an EPTAS. Our new scheme is based on using the shifting technique followed by a solution of polynomial number of $n$-fold programming instances. In addition, we present an asymptotic fully polynomial time approximation scheme (AFPTAS) for the related bin packing type variant of the problem.


title: Learning When-to-Treat Policies

abstract: Many applied decision-making problems have a dynamic component: The policymaker needs not only to choose whom to treat, but also when to start which treatment. For example, a medical doctor may choose between postponing treatment (watchful waiting) and prescribing one of several available treatments during the many visits from a patient. We develop an "advantage doubly robust" estimator for learning such dynamic treatment rules using observational data under the assumption of sequential ignorability. We prove welfare regret bounds that generalize results for doubly robust learning in the single-step setting, and show promising empirical performance in several different contexts. Our approach is practical for policy optimization, and does not need any structural (e.g., Markovian) assumptions.


title: Measurement of distributions sensitive to the underlying event in inclusive $Z$-boson production in pp collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

abstract: This paper presents measurements of charged-particle distributions sensitive to the properties of the underlying event in events containing a $Z$ boson decaying into a muon pair. The data were obtained using the ATLAS detector at the LHC in proton-proton collisions at a centre-of-mass energy of 13 TeV with an integrated luminosity of 3.2 fb$^{-1}$. Distributions of the charged-particle multiplicity and of the charged-particle transverse momentum are measured in regions of the azimuth defined relative to the $Z$ boson direction. The measured distributions are compared with the predictions of various Monte Carlo generators which implement different underlying-event models. The Monte Carlo model predictions qualitatively describe the data well, but with some significant discrepancies.


title: An embedded-hybridized discontinuous Galerkin method for the coupled Stokes-Darcy system

abstract: We introduce an embedded-hybridized discontinuous Galerkin (EDG-HDG) method for the coupled Stokes-Darcy system. This EDG-HDG method is a pointwise mass-conserving discretization resulting in a divergence-conforming velocity field on the whole domain. In the proposed scheme, coupling between the Stokes and Darcy domains is achieved naturally through the EDG-HDG facet variables. \emph{A priori} error analysis shows optimal convergence rates, and that the velocity error does not depend on the pressure. The error analysis is verified through numerical examples on unstructured grids for different orders of polynomial approximation.


title: A Perceptual Weighting Filter Loss for DNN Training in Speech Enhancement

abstract: Single-channel speech enhancement with deep neural networks (DNNs) has shown promising performance and is thus intensively being studied. In this paper, instead of applying the mean squared error (MSE) as the loss function during DNN training for speech enhancement, we design a perceptual weighting filter loss motivated by the weighting filter as it is employed in analysis-by-synthesis speech coding, e.g., in code-excited linear prediction (CELP). The experimental results show that the proposed simple loss function improves the speech enhancement performance compared to a reference DNN with MSE loss in terms of perceptual quality and noise attenuation. The proposed loss function can be advantageously applied to an existing DNN-based speech enhancement system, without modification of the DNN topology for speech enhancement. The source code for the proposed approach is made available.


title: Misspelling Oblivious Word Embeddings

abstract: In this paper we present a method to learn word embeddings that are resilient to misspellings. Existing word embeddings have limited applicability to malformed texts, which contain a non-negligible amount of out-of-vocabulary words. We propose a method combining FastText with subwords and a supervised task of learning misspelling patterns. In our method, misspellings of each word are embedded close to their correct variants. We train these embeddings on a new dataset we are releasing publicly. Finally, we experimentally show the advantages of this approach on both intrinsic and extrinsic NLP tasks using public test sets.


title: Transition between characters of classical groups, decomposition of Gelfand-Tsetlin patterns and last passage percolation

abstract: We study the combinatorial structure of the irreducible characters of the classical groups ${\rm GL}_n(\mathbb{C})$, ${\rm SO}_{2n+1}(\mathbb{C})$, ${\rm Sp}_{2n}(\mathbb{C})$, ${\rm SO}_{2n}(\mathbb{C})$ and the "non-classical" odd symplectic group ${\rm Sp}_{2n+1}(\mathbb{C})$, finding new connections to the probabilistic model of Last Passage Percolation (LPP). Perturbing the expressions of these characters as generating functions of Gelfand-Tsetlin patterns, we produce two families of symmetric polynomials that interpolate between characters of ${\rm Sp}_{2n}(\mathbb{C})$ and ${\rm SO}_{2n+1}(\mathbb{C})$ and between characters of ${\rm SO}_{2n}(\mathbb{C})$ and ${\rm SO}_{2n+1}(\mathbb{C})$. We identify the first family as a one-parameter specialization of Koornwinder polynomials, for which we thus provide a novel combinatorial structure; on the other hand, the second family appears to be new. We next develop a method of Gelfand-Tsetlin pattern decomposition to establish identities between all these polynomials that, in the case of irreducible characters, can be viewed as branching rules. Through these formulas we connect orthogonal and symplectic characters, and more generally the interpolating polynomials, to LPP models with various symmetries, thus going beyond the link with classical Schur polynomials originally found by Baik and Rains (Duke Math. J., 2001). Taking the scaling limit of the LPP models, we finally provide an explanation of why the Tracy-Widom GOE and GSE distributions from random matrix theory admit formulations in terms of both Fredholm determinants and Fredholm Pfaffians.


title: Inverse Scattering and Stability Estimates for The Biharmonic Operator

abstract: We study an inverse scattering problem of a perturbed biharmonic operator. we show that the high-frequency asymptotic of scattering amplitude of the biharmonic operator uniquely determine the curl A and {V -1/2 nablaA}. We also study the near-field scattering problem and show that the high-frequency asymptotic expansion up to a certain error in terms of frequency lambda recovers the same two above quantities with no additional information about A and V. We also prove stability estimates for curl A and {V -1/2 nablaA}.


title: Network Density of States

abstract: Spectral analysis connects graph structure to the eigenvalues and eigenvectors of associated matrices. Much of spectral graph theory descends directly from spectral geometry, the study of differentiable manifolds through the spectra of associated differential operators. But the translation from spectral geometry to spectral graph theory has largely focused on results involving only a few extreme eigenvalues and their associated eigenvalues. Unlike in geometry, the study of graphs through the overall distribution of eigenvalues - the spectral density - is largely limited to simple random graph models. The interior of the spectrum of real-world graphs remains largely unexplored, difficult to compute and to interpret. In this paper, we delve into the heart of spectral densities of real-world graphs. We borrow tools developed in condensed matter physics, and add novel adaptations to handle the spectral signatures of common graph motifs. The resulting methods are highly efficient, as we illustrate by computing spectral densities for graphs with over a billion edges on a single compute node. Beyond providing visually compelling fingerprints of graphs, we show how the estimation of spectral densities facilitates the computation of many common centrality measures, and use spectral densities to estimate meaningful information about graph structure that cannot be inferred from the extremal eigenpairs alone.


title: Smoothness and monotonicity of the excursion set density of planar Gaussian fields

abstract: Nazarov and Sodin have shown that the number of connected components of the nodal set of a planar Gaussian field in a ball of radius $R$, normalised by area, converges to a constant as $R\to \infty $. This has been generalised to excursion/level sets at arbitrary levels, implying the existence of functionals $c_{ES}(\ell )$ and $c_{LS}(\ell )$ that encode the density of excursion/level set components at the level $\ell $. We prove that these functionals are continuously differentiable for a wide class of fields. This follows from a more general result, which derives differentiability of the functionals from the decay of the probability of `four-arm events' for the field conditioned to have a saddle point at the origin. For some fields, including the important special cases of the Random Plane Wave and the Bargmann-Fock field, we also derive stochastic monotonicity of the conditioned field, which allows us to deduce regions on which $c_{ES}(\ell )$ and $c_{LS}(\ell )$ are monotone.


title: Design Dimensions for Software Certification: A Grounded Analysis

abstract: In many domains, software systems cannot be deployed until authorities judge them fit for use in an intended operating environment. Certification standards and processes have been devised and deployed to regulate operations of software systems and prevent their failures. However, practitioners are often unsatisfied with the efficiency and value proposition of certification efforts. In this study, we compare two certification standards, Common Criteria and DO-178C, and collect insights from literature and from interviews with subject-matter experts to identify design options relevant to the design of standards. The results of the comparison of certification efforts---leading to the identification of design dimensions that affect their quality---serve as a framework to guide the comparison, creation, and revision of certification standards and processes. This paper puts software engineering research in context and discusses key issues around process and quality assurance and includes observations from industry about relevant topics such as recertification, timely evaluations, but also technical discussions around model-driven approaches and formal methods. Our initial characterization of the design space of certification efforts can be used to inform technical discussions and to influence the directions of new or existing certification efforts. Practitioners, technical commissions, and government can directly benefit from our analytical framework.


title: An Efficient Approach for Super and Nested Term Indexing and Retrieval

abstract: This paper describes a new approach, called Terminological Bucket Indexing (TBI), for efficient indexing and retrieval of both nested and super terms using a single method. We propose a hybrid data structure for facilitating faster indexing building. An evaluation of our approach with respect to widely used existing approaches on several publicly available dataset is provided. Compared to Trie based approaches, TBI provides comparable performance on nested term retrieval and far superior performance on super term retrieval. Compared to traditional hash table, TBI needs 80\% less time for indexing.


title: A Spectral Generalization of Von Neumann Minimax Theorem

abstract: Given $n \times n$ real symmetric matrices $A_1, \dots, A_m$, the following {\it spectral minimax} property holds: $$\min_{X \in \mathbf{\Delta}_n} \max_{y \in S_m} \sum_{i=1}^m y_iA_i \bullet X=\max_{y \in S_m} \min_{X \in \mathbf{\Delta}_n} \sum_{i=1}^m y_iA_i \bullet X,$$ where $S_m$ is the simplex and $\mathbf{\Delta}_n$ the spectraplex. For diagonal $A_i$'s this reduces to the classic minimax.


title: GLEE: Geometric Laplacian Eigenmap Embedding

abstract: Graph embedding seeks to build a low-dimensional representation of a graph G. This low-dimensional representation is then used for various downstream tasks. One popular approach is Laplacian Eigenmaps, which constructs a graph embedding based on the spectral properties of the Laplacian matrix of G. The intuition behind it, and many other embedding techniques, is that the embedding of a graph must respect node similarity: similar nodes must have embeddings that are close to one another. Here, we dispose of this distance-minimization assumption. Instead, we use the Laplacian matrix to find an embedding with geometric properties instead of spectral ones, by leveraging the so-called simplex geometry of G. We introduce a new approach, Geometric Laplacian Eigenmap Embedding (or GLEE for short), and demonstrate that it outperforms various other techniques (including Laplacian Eigenmaps) in the tasks of graph reconstruction and link prediction.


title: Invariant manifolds of partially normally hyperbolic invariant manifolds in Banach spaces

abstract: We investigate the existence and regularity of (locally) invariant manifolds nearby an approximately invariant set satisfying certain (geometric) hyperbolicity with respect to an abstract `generalized' dynamical system in a Banach space; such hyperbolicity is between normal hyperbolicity and partial hyperbolicity which has being studied in the finite-dimension and in some concrete PDEs. The `generalized' dynamical system is allowed to be non-smooth, non-Lipschitz, or even `non-mapping', making it applicable to both well-posed and ill-posed differential equations. As an illustration, we apply our results to study the dynamics of the whiskered tori.


title: Convergence Analysis of (Statistical) Inverse Problems under Conditional Stability Estimates

abstract: Conditional stability estimates require additional regularization for obtaining stable approximate solutions if the validity area of such estimates is not completely known. In this context, we consider ill-posed nonlinear inverse problems in Hilbert scales satisfying conditional stability estimates characterized by general concave index functions. For that case, we exploit Tikhonov regularization and provide convergence and convergence rates of regularized solutions for both deterministic and stochastic noise. We further discuss a priori and a posteriori parameter choice rules and illustrate the validity of our assumptions in different model and real world situations.


title: Workflow Design Analysis for High Resolution Satellite Image Analysis

abstract: Ecological sciences are using imagery from a variety of sources to monitor and survey populations and ecosystems. Very High Resolution (VHR) satellite imagery provide an effective dataset for large scale surveys. Convolutional Neural Networks have successfully been employed to analyze such imagery and detect large animals. As the datasets increase in volume, O(TB), and number of images, O(1k), utilizing High Performance Computing (HPC) resources becomes necessary. In this paper, we investigate a task-parallel data-driven workflows design to support imagery analysis pipelines with heterogeneous tasks on HPC. We analyze the capabilities of each design when processing a dataset of 3,000 VHR satellite images for a total of 4~TB. We experimentally model the execution time of the tasks of the image processing pipeline. We perform experiments to characterize the resource utilization, total time to completion, and overheads of each design. Based on the model, overhead and utilization analysis, we show which design approach to is best suited in scientific pipelines with similar characteristics.


title: The Role of Asymmetries in Thermal Non-Equilibrium

abstract: Thermal non-equilibrium (TNE) is a fascinating situation that occurs in coronal magnetic flux tubes (loops) for which no solution to the steady-state fluid equations exists. The plasma is constantly evolving even though the heating that produces the hot temperatures does not. This is a promising explanation for isolated phenomena such as prominences, coronal rain, and long-period pulsating loops, but it may also have much broader relevance. As known for some time, TNE requires that the heating be both (quasi) steady and concentrated at low coronal altitudes. Recent studies indicate that asymmetries are also important, with large enough asymmetries in the heating and/or cross-sectional area resulting in steady flow rather than TNE. Using reasonable approximations, we have derived two formulae for quantifying the conditions necessary for TNE. As a rough rule of thumb, the ratio of apex to footpoint heating rates must be less than about 0.1, and asymmetries must be less than about a factor of 3. The precise values are case dependent. We have tested our formulae with 1D hydrodynamic loop simulations and find a very acceptable agreement. These results are important for developing physical insight about TNE and assessing how widespread it may be on the Sun.


title: Zero-shot Knowledge Transfer via Adversarial Belief Matching

abstract: Performing knowledge transfer from a large teacher network to a smaller student is a popular task in modern deep learning applications. However, due to growing dataset sizes and stricter privacy regulations, it is increasingly common not to have access to the data that was used to train the teacher. We propose a novel method which trains a student to match the predictions of its teacher without using any data or metadata. We achieve this by training an adversarial generator to search for images on which the student poorly matches the teacher, and then using them to train the student. Our resulting student closely approximates its teacher for simple datasets like SVHN, and on CIFAR10 we improve on the state-of-the-art for few-shot distillation (with 100 images per class), despite using no data. Finally, we also propose a metric to quantify the degree of belief matching between teacher and student in the vicinity of decision boundaries, and observe a significantly higher match between our zero-shot student and the teacher, than between a student distilled with real data and the teacher. Code available at: https://github.com/polo5/ZeroShotKnowledgeTransfer


title: Interface states in polariton topological insulators

abstract: We address linear and nonlinear topological interface states in polariton condensates excited at the interface of the honeycomb and Lieb arrays of microcavity pillars in the presence of spin-orbit coupling and Zeeman splitting in the external magnetic field. Such interface states appear only in total energy gaps of the composite structure when parameters of the honeycomb and Lieb arrays are selected such that some topological gaps in the spectrum of one of the arrays overlap with topological or nontopological gaps in the spectrum of the other array. This is in contrast to conventional edge states at the interface of periodic topological and uniform trivial insulators, whose behavior is determined exclusively by the spectrum of the topological medium. The number of emerging interface states is determined by the difference of the Chern numbers of the overlapping gaps. Illustrative examples with one or two coexisting unidirectional interface states are provided. The representative feature of the system is the possibility of wide tuning of the concentration of power of the interface states between two limiting cases when practically all power is concentrated either in the Lieb or the honeycomb array. Localization of the interface states and their penetration depth into arrays drastically vary with Bloch momentum or upon modification of the amplitude of the interface state in the nonlinear regime. We illustrate topological protection of the interface states manifested in the absence of backscattering on interface defects, and study their modulation instability in the nonlinear regime. The latter leads to formation of quasisolitons whose penetration into different arrays also depends on Bloch momentum. In addition, we discuss the impact of losses and coherent pump leading to bistability of the interface states.


title: Polynomial-time proofs that groups are hyperbolic

abstract: It is undecidable in general whether a given finitely presented group is word hyperbolic. We use the concept of pregroups, introduced by Stallings, to define a new class of van Kampen diagrams, which represent groups as quotients of virtually free groups. We then present a polynomial-time procedure which analyses these diagrams, and either returns an explicit linear Dehn function for the presentation, or returns fail, together with its reasons for failure. Furthermore, if our procedure succeeds we are often able to produce in polynomial time a word problem solver for the presentation that runs in linear time. Our algorithms have been implemented, and are often many orders of magnitude faster than KBMAG, the only comparable publicly available software.


title: Multi-Service Mobile Traffic Forecasting via Convolutional Long Short-Term Memories

abstract: Network slicing is increasingly used to partition network infrastructure between different mobile services. Precise service-wise mobile traffic forecasting becomes essential in this context, as mobile operators seek to pre-allocate resources to each slice in advance, to meet the distinct requirements of individual services. This paper attacks the problem of multi-service mobile traffic forecasting using a sequence-to-sequence (S2S) learning paradigm and convolutional long short-term memories (ConvLSTMs). The proposed architecture is designed so as to effectively extract complex spatiotemporal features of mobile network traffic and predict with high accuracy the future demands for individual services at city scale. We conduct experiments on a mobile traffic dataset collected in a large European metropolis, demonstrating that the proposed S2S-ConvLSTM can forecast the mobile traffic volume produced by tens of different services in advance of up to one hour, by just using measurements taken during the past hour. In particular, our solution achieves mean absolute errors (MAE) at antenna level that are below 13KBps, outperforming other deep learning approaches by up to 31.2%.


title: What is the Entropy of a Social Organization?

abstract: We quantify a social organization's potentiality, that is its ability to attain different configurations. The organization is represented as a network in which nodes correspond to individuals and (multi-)edges to their multiple interactions. Attainable configurations are treated as realizations from a network ensemble. To encode interaction preferences between individuals, we choose the generalized hypergeometric ensemble of random graphs, which is described by a closed-form probability distribution. From this distribution we calculate Shannon entropy as a measure of potentiality. This allows us to compare different organizations as well different stages in the development of a given organization. The feasibility of the approach is demonstrated using data from 3 empirical and 2 synthetic systems.


title: Speech2Face: Learning the Face Behind a Voice

abstract: How much can we infer about a person's looks from the way they speak? In this paper, we study the task of reconstructing a facial image of a person from a short audio recording of that person speaking. We design and train a deep neural network to perform this task using millions of natural Internet/YouTube videos of people speaking. During training, our model learns voice-face correlations that allow it to produce images that capture various physical attributes of the speakers such as age, gender and ethnicity. This is done in a self-supervised manner, by utilizing the natural co-occurrence of faces and speech in Internet videos, without the need to model attributes explicitly. We evaluate and numerically quantify how--and in what manner--our Speech2Face reconstructions, obtained directly from audio, resemble the true face images of the speakers.


title: Charmed hadron production in an improved quark coalescence model

abstract: We study the production of charmed hadrons $D^{0}$ and $\Lambda_c^+$ in relativistic heavy-ion collisions using an improved quark coalescence model. In particular, we extend the usual coalescence model by letting a produced hadron to have the same velocity as the center-of-mass velocity of coalesced constituent quarks during hadronization to take into account the effect of collective flow in produced quark-gluon plasma. This results in a shift of charmed resonances of higher masses to larger transverse momenta ($p_T^{}$). Requiring all charm quarks of very low $p_T^{}$ to be converted to hadrons via coalescence and letting charm quarks not undergoing coalescence to hadronize by independent fragmentation, we obtain a good description of the measured yield ratio $\Lambda_c^+/D^0$ as a function of $p_T^{}$ in $\text{Au} + \text{Au}$ collisions at $\sqrt{s_{NN}}^{}=200$~GeV by the STAR Collaboration at the Relativistic Heavy Ion Collider.


title: Determination of topological properties of thin samples by the van der Pauw method

abstract: We solve the problem of determining basic topological properties of flat samples by performing measurements on their outer edge. The global maximum of four probe resistances shows a characteristic behaviour, which is dependent on the genus (i.e., the number of holes) of the domain. An extension of the van der Pauw method on domains having zero, one, or two holes is presented and discussed. A possibility of measuring topological properties of condensed matter is demonstrated. Experimental results for triply connected domains are presented and explained by continuous symmetry breaking caused by the presence of two holes. The results are consistent with the topological theorem of Hurwitz on the number of automorphisms of Riemann surfaces.


title: Can fermionic dark matter mimic supermassive black holes?

abstract: We analyze the intriguing possibility to explain both dark mass components in a galaxy: the dark matter (DM) halo and the supermassive dark compact object lying at the center, by a unified approach in terms of a quasi-relaxed system of massive, neutral fermions in general relativity. The solutions to the mass distribution of such a model that fulfill realistic halo boundary conditions inferred from observations, develop a highly-density core supported by the fermion degeneracy pressure able to mimic massive black holes at the center of galaxies. Remarkably, these dense core-diluted halo configurations can explain the dynamics of the closest stars around Milky Way's center (SgrA*) all the way to the halo rotation curve, without spoiling the baryonic bulge-disk components, for a narrow particle mass range $mc^2 \sim 10$-$10^2$~keV.


title: A Smoothness Energy without Boundary Distortion for Curved Surfaces

abstract: Current quadratic smoothness energies for curved surfaces either exhibit distortions near the boundary due to zero Neumann boundary conditions, or they do not correctly account for intrinsic curvature, which leads to unnatural-looking behavior away from the boundary. This leads to an unfortunate trade-off: one can either have natural behavior in the interior, or a distortion-free result at the boundary, but not both. We introduce a generalized Hessian energy for curved surfaces, expressed in terms of the covariant one-form Dirichlet energy, the Gaussian curvature, and the exterior derivative. Energy minimizers solve the Laplace-Beltrami biharmonic equation, correctly accounting for intrinsic curvature, leading to natural-looking isolines. On the boundary, minimizers are as-linear-as-possible, which reduces the distortion of isolines at the boundary. We discretize the covariant one-form Dirichlet energy using Crouzeix-Raviart finite elements, arriving at a discrete formulation of the Hessian energy for applications on curved surfaces. We observe convergence of the discretization in our experiments.


title: Privacy-Preserving Obfuscation of Critical Infrastructure Networks

abstract: The paper studies how to release data about a critical infrastructure network (e.g., the power network or a transportation network) without disclosing sensitive information that can be exploited by malevolent agents, while preserving the realism of the network. It proposes a novel obfuscation mechanism that combines several privacy-preserving building blocks with a bi-level optimization model to significantly improve accuracy. The obfuscation is evaluated for both realism and privacy properties on real energy and transportation networks. Experimental results show the obfuscation mechanism substantially reduces the potential damage of an attack exploiting the released data to harm the real network.


title: A general constitutive model for dense, fine particle suspensions validated in many geometries

abstract: Fine particle suspensions (such as cornstarch mixed with water) exhibit dramatic changes in viscosity when sheared, producing fascinating behaviors that captivate children and rheologists alike. Recent examination of these mixtures in simple flow geometries suggests inter-granular repulsion is central to this effect --- for mixtures at rest or shearing slowly, repulsion prevents frictional contacts from forming between particles, whereas, when sheared more forcefully, granular stresses overcome the repulsion allowing particles to interact frictionally and form microscopic structures that resist flow. Previous constitutive studies of these mixtures have focused on particular cases, typically limited to two-dimensional, steady, simple shearing flows. In this work, we introduce a predictive and general, three-dimensional continuum model for this material, using mixture theory to couple the fluid and particle phases. Playing a central role in the model, we introduce a micro-structural state variable, whose evolution is deduced from small-scale physical arguments and checked with existing data. Our space- and time-dependent model is implemented numerically in a variety of unsteady, non-uniform flow configurations where it is shown to accurately capture a variety of key behaviors: (i) the continuous shear thickening (CST) and discontinuous shear thickening (DST) behavior observed in steady flows, (ii) the time-dependent propagation of `shear jamming fronts', (iii) the time-dependent propagation of `impact activated jamming fronts', and (iv) the non-Newtonian, `running on oobleck' effect wherein fast locomotors stay afloat while slow ones sink.


title: Bayesian Optimization with Approximate Set Kernels

abstract: We propose a practical Bayesian optimization method over sets, to minimize a black-box function that takes a set as a single input. Because set inputs are permutation-invariant, traditional Gaussian process-based Bayesian optimization strategies which assume vector inputs can fall short. To address this, we develop a Bayesian optimization method with \emph{set kernel} that is used to build surrogate functions. This kernel accumulates similarity over set elements to enforce permutation-invariance, but this comes at a greater computational cost. To reduce this burden, we propose two key components: (i) a more efficient approximate set kernel which is still positive-definite and is an unbiased estimator of the true set kernel with upper-bounded variance in terms of the number of subsamples, (ii) a constrained acquisition function optimization over sets, which uses symmetry of the feasible region that defines a set input. Finally, we present several numerical experiments which demonstrate that our method outperforms other methods.


title: Model-independent determination of $H_0$ and $\Omega_{K0}$ from strong lensing and type Ia supernovae

abstract: We present the first determination of the Hubble constant $H_0$ from strong lensing time delay data and type Ia supernova luminosity distances that is independent of the cosmological model. We also determine the spatial curvature model-independently. We assume that light propagation over long distances is described by the FLRW metric and geometrical optics holds, but make no assumption about the contents of the Universe or the theory of gravity on cosmological scales. We find $H_0=75.7^{+4.5}_{-4.4}$ km/s/Mpc and $\Omega_{K0}=0.12^{+0.27}_{-0.25}$. This is a 6\% determination of $H_0$. A weak prior from the cosmic microwave background on the distance to the last scattering surface improves this to $H_0=76.8^{+4.2}_{-3.8}$ km/s/Mpc and $\Omega_{K0}=0.18^{+0.25}_{-0.18}$. Assuming zero spatial curvature, we get $H_0=74.2^{+3.0}_{-2.9}$ km/s/Mpc, a precision of $4\%$. The measurements also provide a consistency test of the FLRW metric: we find no evidence against it.


title: On the Bourbaki's fixed point theorem and the axiom of choice

abstract: In this note we generalize the Moroianu's fixed point theorem. We propose a very elegant common proof of the Bourbaki's fixed point theorem and our result. We apply our result to give a very elegant proof of the fact that, in the Zermelo-Fraenkel system, the axiom of choice is equivalent to each of the following statements: H. Kneser's Lemma, Zorn's Lemma, Zermelo's Lemma


title: Hinode EIS line widths in the quiet corona up to 1.5 Rsun

abstract: We present an analysis of several Hinode EIS observations of coronal line widths in the quiet Sun, up to 1.5 Rsun radial distances. No significant variations are found, which indicates no damping of Alfv\'en waves in the quiescent corona. However, the uncertainties in estimating the instrumental width mean that a firm conclusion cannot be reached. We present a discussion of various EIS instrumental issues and suggest that the strongest lines, from Fe XII at 193.5 and 195.1 A, have anomalous instrumental widths. We show how line widths in EIS are uncertain when the signal is low, and that the instrumental variation along the slit is also uncertain. We also found an anomalous decrease (up to 40%) in the intensities of these lines in many off-limb and active region observations, and suggest that this is due to opacity effects. We find that the most reliable measurements are obtained from the weaker lines.


title: Puzzling Blue Dips in the Black Hole Candidate Swift J1357.2-0933, from ULTRACAM, SALT, ATCA, Swift and NuSTAR

abstract: We present rapid, multiwavelength photometry of the low-mass X-ray binary Swift J1357.2-0933 during its 2017 outburst. Using several sets of quasi-simultaneous ULTRACAM/NTT (optical), NuSTAR (X-ray), XRT/Swift (X-ray), SALT (optical) and ATCA (radio) observations taken during outburst decline, we confirm the frequent optical dipping that has previously been noted both in outburst and in quiescence. We also find: 1) that the dip frequency decreases as the outburst decays, similar to what was seen in the previous outburst, 2) that the dips produce a shape similar to that in binary systems with partial disc occultations, 3) that the source becomes significantly bluer during these dips, indicating an unusual geometry compared to other LMXB dippers, and 4) that dip superposition analysis confirms the lack of an X-ray response to the optical dips. These very unusual properties appear to be unique to Swift J1357.2-0933, and are likely the result of a high binary inclination, as inferred from features such as its very low outburst X-ray luminosity. From this analysis as well as X-ray/optical timing correlations, we suggest a model with multi-component emission/absorption features with differing colours. This could include the possible presence of a sporadically occulted jet base and a recessed disc. This source still hosts many puzzling features, with consequences for the very faint X-ray transients population.


title: ALMA Polarimetry of AT2018cow

abstract: We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model reported by Ho et al. (2019) and indicate that the spectral peaks ($\sim110$ GHz at the first epoch and $\sim67$ GHz at the second epoch) represent the synchrotron self-absorption frequency. The non-detection of linear polarization with $<$0.15% in the 230 GHz band at the phase when the effect of synchrotron self-absorption was quite small in the band may be explained by internal Faraday depolarization with high circumburst density and strong magnetic field. This result supports the stellar explosion scenario rather than the tidal disruption model. The maximum energy of accelerating particles at the shocks of AT2018cow-like objects is also discussed.


title: Positional Encoding by Robots with Non-Rigid Movements

abstract: Consider a set of autonomous computational entities, called \emph{robots}, operating inside a polygonal enclosure (possibly with holes), that have to perform some collaborative tasks. The boundary of the polygon obstructs both visibility and mobility of a robot. Since the polygon is initially unknown to the robots, the natural approach is to first explore and construct a map of the polygon. For this, the robots need an unlimited amount of persistent memory to store the snapshots taken from different points inside the polygon. However, it has been shown by Di Luna et al. [DISC 2017] that map construction can be done even by oblivious robots by employing a positional encoding strategy where a robot carefully positions itself inside the polygon to encode information in the binary representation of its distance from the closest polygon vertex. Of course, to execute this strategy, it is crucial for the robots to make accurate movements. In this paper, we address the question whether this technique can be implemented even when the movements of the robots are unpredictable in the sense that the robot can be stopped by the adversary during its movement before reaching its destination. However, there exists a constant $\delta > 0$, unknown to the robot, such that the robot can always reach its destination if it has to move by no more than $\delta$ amount. This model is known in literature as \emph{non-rigid} movement. We give a partial answer to the question in the affirmative by presenting a map construction algorithm for robots with non-rigid movement, but having $O(1)$ bits of persistent memory and ability to make circular moves.


title: Search for heavy neutral leptons in decays of $W$ bosons produced in 13 TeV $pp$ collisions using prompt and displaced signatures with the ATLAS detector

abstract: The problems of neutrino masses, matter-antimatter asymmetry, and dark matter could be successfully addressed by postulating right-handed neutrinos with Majorana masses below the electroweak scale. In this work, leptonic decays of $W$ bosons extracted from 32.9 fb$^{-1}$ to 36.1 fb$^{-1}$ of 13 TeV proton-proton collisions at the LHC are used to search for heavy neutral leptons (HNLs) that are produced through mixing with muon or electron neutrinos. The search is conducted using the ATLAS detector in both prompt and displaced leptonic decay signatures. The prompt signature requires three leptons produced at the interaction point (either $\mu\mu e$ or $e e\mu$) with a veto on same-flavour opposite-charge topologies. The displaced signature comprises a prompt muon from the $W$ boson decay and the requirement of a dilepton vertex (either $\mu\mu$ or $\mu e$) displaced in the transverse plane by 4-300 mm from the interaction point. The search sets constraints on the HNL mixing to muon and electron neutrinos for HNL masses in the range 4.5-50 GeV.


title: Multi-Sample Dropout for Accelerated Training and Better Generalization

abstract: Dropout is a simple but efficient regularization technique for achieving better generalization of deep neural networks (DNNs); hence it is widely used in tasks based on DNNs. During training, dropout randomly discards a portion of the neurons to avoid overfitting. This paper presents an enhanced dropout technique, which we call multi-sample dropout, for both accelerating training and improving generalization over the original dropout. The original dropout creates a randomly selected subset (called a dropout sample) from the input in each training iteration while the multi-sample dropout creates multiple dropout samples. The loss is calculated for each sample, and then the sample losses are averaged to obtain the final loss. This technique can be easily implemented by duplicating a part of the network after the dropout layer while sharing the weights among the duplicated fully connected layers. Experimental results using image classification tasks including ImageNet, CIFAR-10, and CIFAR-100 showed that multi-sample dropout accelerates training. Moreover, the networks trained using multi-sample dropout achieved lower error rates compared to networks trained with the original dropout. The additional computation cost due to the duplicated operations is not significant for deep convolutional networks because most of the computation time is consumed in the convolution layers before the dropout layer, which are not duplicated.


title: A qubit model for U(1) lattice gauge theory

abstract: A conceptually simple model for strongly interacting compact U(1) lattice gauge theory is expressed as operators acting on qubits. The number of independent gauge links is reduced to its minimum through the use of Gauss's law. The model can be implemented with any number of qubits per gauge link, and a choice as small as two is shown to be useful. Real-time propagation and real-time collisions are observed on lattices in two spatial dimensions. The extension to three spatial dimensions is also developed, and a first look at 3-dimensional real-time dynamics is presented.


title: Cross-verification of independent quantum devices

abstract: Quantum computers are on the brink of surpassing the capabilities of even the most powerful classical computers. This naturally raises the question of how one can trust the results of a quantum computer when they cannot be compared to classical simulation. Here we present a verification technique that exploits the principles of measurement-based quantum computation to link quantum circuits of different input size, depth, and structure. Our approach enables consistency checks of quantum computations within a device, as well as between independent devices. We showcase our protocol by applying it to five state-of-the-art quantum processors, based on four distinct physical architectures: nuclear magnetic resonance, superconducting circuits, trapped ions, and photonics, with up to 6 qubits and 200 distinct circuits.


title: Multi-relational Poincar\'e Graph Embeddings

abstract: Hyperbolic embeddings have recently gained attention in machine learning due to their ability to represent hierarchical data more accurately and succinctly than their Euclidean analogues. However, multi-relational knowledge graphs often exhibit multiple simultaneous hierarchies, which current hyperbolic models do not capture. To address this, we propose a model that embeds multi-relational graph data in the Poincar\'e ball model of hyperbolic space. Our Multi-Relational Poincar\'e model (MuRP) learns relation-specific parameters to transform entity embeddings by M\"obius matrix-vector multiplication and M\"obius addition. Experiments on the hierarchical WN18RR knowledge graph show that our Poincar\'e embeddings outperform their Euclidean counterpart and existing embedding methods on the link prediction task, particularly at lower dimensionality.


title: Scheme for Majorana Manipulation Using Magnetic Force Microscopy

abstract: We propose a scheme for the use of magnetic force microscopy to manipulate Majorana zero modes emergent in vortex cores of topological superconductors in the Fe(Se,Te) family. We calculate the pinning forces necessary to drag two vortices together and the resulting change in current and charge density of the composite fermion. A possible algorithm for measuring and altering Majorana pair parity is demonstrated.


title: Electricity Market Equilibrium under Information Asymmetry

abstract: We study a competitive electricity market equilibrium with two trading stages, day-ahead and real-time. The welfare of each market agent is exposed to uncertainty (here from renewable energy production), while agent information on the probability distribution of this uncertainty is not identical at the day-ahead stage. We show a high sensitivity of the equilibrium solution to the level of information asymmetry and demonstrate economic, operational, and computational value for the system stemming from potential information sharing.


title: Evaluating the Effects of Control Surfaces Failure on the GTM

abstract: Despite the advances in aircraft guidance and control systems technology, Loss of Control remains as the main cause of the fatal accidents of large transport aircraft. Loss of Control is defined as excursion beyond the allowable flight envelope and is often a consequence of upset condition if improper maneuver is implemented by the pilot. Hence, extensive research in recent years has focused on improving the current fault tolerant control systems and developing new strategies for loss of control prevention and recovery systems. However, success of such systems requires the perception of the damaged aircraft's dynamic behavior and performance, and understanding of its new flight envelope. This paper provides a comprehensive understanding of lateral control surfaces' failure effect on the NASA Generic Transport Model's maneuvering flight envelope; which is a set of attainable steady state maneuvers herein referred to as trim points. The study utilizes a massive database of the Generic Transport Model's high-fidelity maneuvering flight envelopes computed for the unimpaired case and wide ranges of aileron and rudder failure cases at different flight conditions. Flight envelope boundary is rigorously investigated and the key parameters confining the trim points at different boundary sections are identified. Trend analysis of the impaired flight envelopes and the corresponding limiting factors is performed which demonstrates the effect of various failure degrees on the remaining feasible trim points. Results of the post-failure analysis can be employed in emergency path planning and have potential uses in the development of aircraft resilient control and upset recovery systems.


title: Nature-Inspired Computational Model of Population Desegregation under Group Leaders Influence

abstract: This paper presents an agent-based model of population desegregation and provides a thorough analysis of the social behavior leading to it, namely, the contact hypothesis. Based on the parameters of frequency and intensity of influence of group leaders on the population, the proposed model is constituted by two layers: 1) a physical layer of the population that is influenced by and 2) a virtual layer of group leaders. The model of negotiation and survival of group leaders are governed by the nature-inspired evolutionary process of queen ants, also known as Foundress Dilemma. The motivation of using a virtual grouping concept (instead of taking a subset of population as the group leaders) is to stay focused on finding the conditions leading individuals in a society tolerating a significantly diversified (desegregated) neighborhood, rather than, indulging into complex details, which would be more relevant to studies targeting the evolution of societal group and leaders. A geographic information system-driven simulation is performed, which reveals that: 1) desegregation is directly proportional to the frequency of group leaders' contact with the population and 2) mostly, it remains ineffective with an increase in the intensity of group leaders' contact with the population. The mechanism of group selection (the conflict resolution model resolving the Foundress Dilemma) reveals an exciting result concerning negative influence of cooperative group leaders. Most of the time, desegregation decreases with increase in cooperative leaders (the leaders enforcing desegregation) when compared with fierce leaders (the leaders enforcing segregation).


title: Augmenting correlation structures in spatial data using deep generative models

abstract: State-of-the-art deep learning methods have shown a remarkable capacity to model complex data domains, but struggle with geospatial data. In this paper, we introduce SpaceGAN, a novel generative model for geospatial domains that learns neighbourhood structures through spatial conditioning. We propose to enhance spatial representation beyond mere spatial coordinates, by conditioning each data point on feature vectors of its spatial neighbours, thus allowing for a more flexible representation of the spatial structure. To overcome issues of training convergence, we employ a metric capturing the loss in local spatial autocorrelation between real and generated data as stopping criterion for SpaceGAN parametrization. This way, we ensure that the generator produces synthetic samples faithful to the spatial patterns observed in the input. SpaceGAN is successfully applied for data augmentation and outperforms compared to other methods of synthetic spatial data generation. Finally, we propose an ensemble learning framework for the geospatial domain, taking augmented SpaceGAN samples as training data for a set of ensemble learners. We empirically show the superiority of this approach over conventional ensemble learning approaches and rivaling spatial data augmentation methods, using synthetic and real-world prediction tasks. Our findings suggest that SpaceGAN can be used as a tool for (1) artificially inflating sparse geospatial data and (2) improving generalization of geospatial models.


title: Interpreting Adversarially Trained Convolutional Neural Networks

abstract: We attempt to interpret how adversarially trained convolutional neural networks (AT-CNNs) recognize objects. We design systematic approaches to interpret AT-CNNs in both qualitative and quantitative ways and compare them with normally trained models. Surprisingly, we find that adversarial training alleviates the texture bias of standard CNNs when trained on object recognition tasks, and helps CNNs learn a more shape-biased representation. We validate our hypothesis from two aspects. First, we compare the salience maps of AT-CNNs and standard CNNs on clean images and images under different transformations. The comparison could visually show that the prediction of the two types of CNNs is sensitive to dramatically different types of features. Second, to achieve quantitative verification, we construct additional test datasets that destroy either textures or shapes, such as style-transferred version of clean data, saturated images and patch-shuffled ones, and then evaluate the classification accuracy of AT-CNNs and normal CNNs on these datasets. Our findings shed some light on why AT-CNNs are more robust than those normally trained ones and contribute to a better understanding of adversarial training over CNNs from an interpretation perspective.


title: Deuterated methanol toward NGC 7538-IRS1

abstract: We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.


title: Conformally Soft Theorem In Gravity

abstract: A central feature of scattering amplitudes in gravity or gauge theory is the existence of a variety of energetically soft theorems which put constraints on the amplitudes. Celestial amplitudes which are obtained from momentum-space amplitudes by a Mellin transform over the external particle energies cannot obey the usual energetically soft theorems. Instead, the symmetries of the celestial sphere imply that the scattering of conformally soft particles whose conformal weights under the 4D Lorentz group SL(2,C) are taken to zero obey special relations. Such conformally soft theorems have recently been found for gauge theory. Here, I show conformally soft factorization of celestial amplitudes for gravity and identify it as the celestial analogue of Weinberg's soft graviton theorem.


title: Gaussbock: Fast parallel-iterative cosmological parameter estimation with Bayesian nonparametrics

abstract: We present and apply Gaussbock, a new embarrassingly parallel iterative algorithm for cosmological parameter estimation designed for an era of cheap parallel computing resources. Gaussbock uses Bayesian nonparametrics and truncated importance sampling to accurately draw samples from posterior distributions with an orders-of-magnitude speed-up in wall time over alternative methods. Contemporary problems in this area often suffer from both increased computational costs due to high-dimensional parameter spaces and consequent excessive time requirements, as well as the need for fine tuning of proposal distributions or sampling parameters. Gaussbock is designed specifically with these issues in mind. We explore and validate the performance and convergence of the algorithm on a fast approximation to the Dark Energy Survey Year 1 (DES Y1) posterior, finding reasonable scaling behavior with the number of parameters. We then test on the full DES Y1 posterior using large-scale supercomputing facilities, and recover reasonable agreement with previous chains, although the algorithm can underestimate the tails of poorly-constrained parameters. Additionally, we discuss and demonstrate how Gaussbock recovers complex posterior shapes very well at lower dimensions, but faces challenges to perform well on such distributions in higher dimensions. In addition, we provide the community with a user-friendly software tool for accelerated cosmological parameter estimation based on the methodology described in this paper.


title: Spanning Trees in Graphs of High Minimum Degree with a Universal Vertex I: An Asymptotic Result

abstract: In this paper and a companion paper, we prove that, if $m$ is sufficiently large, every graph on $m+1$ vertices that has a universal vertex and minimum degree at least $\lfloor \frac{2m}{3} \rfloor$ contains each tree $T$ with $m$ edges as a subgraph. Our result confirms, for large $m$, an important special case of a recent conjecture by Havet, Reed, Stein, and Wood. The present paper already contains an approximate version of the result.


title: MIFAL: fully automated Multiple-Image Finder ALgorithm for strong-lens modelling -- proof of concept

abstract: We outline a simple procedure designed for \emph{automatically} finding sets of multiple images in strong lensing (SL) clusters. We show that by combining (a) an arc-finding (or source extracting) program, (b) photometric redshift measurements, and (c) a preliminary light-traces-mass lens model, multiple-image systems can be identified in a fully automated (`blind') manner. The presented procedure yields an assessment of the likelihood of each arc to belong to one of the multiple-image systems, as well as the preferred redshift for the different systems. These could be then used to automatically constrain and refine the initial lens model for an accurate mass distribution. We apply this procedure to \emph{Cluster Lensing And Supernova with Hubble} observations of three galaxy clusters, MACS J0329.6-0211, MACS J1720.2+3536, and MACS J1931.8-2635, comparing the results to published SL analyses where multiple images were verified by eye on a particular basis. In the first cluster all originally identified systems are recovered by the automated procedure, and in the second and third clusters about half are recovered. Other known systems are not picked up, in part due to a crude choice of parameters, ambiguous photometric redshifts, or inaccuracy of the initial lens model. On top of real systems recovered, some false images are also mistakenly identified by the procedure, depending on the thresholds used. While further improvements to the procedure and a more thorough scrutinisation of its performance are warranted, the work constitutes another important step toward fully automatising SL analyses for studying mass distributions of large cluster samples.