-
Vortex confinement through an unquantized magnetic flux
Authors:
Geunyong Kim,
Jinyoung Yun,
Jinho Yang,
Ilkyu Yang,
Dirk Wulferding,
Roman Movshovich,
Gil Young Cho,
Ki-Seok Kim,
Garam Hahn,
Jeehoon Kim
Abstract:
Geometrically confined superconductors often experience a breakdown in the quantization of magnetic flux owing to the incomplete screening of the supercurrent against the field penetration. In this study, we report that the confinement of a magnetic field occurs regardless of the dimensionality of the system, extending even to 1D linear potential systems. By utilizing a vector-field magnetic force…
▽ More
Geometrically confined superconductors often experience a breakdown in the quantization of magnetic flux owing to the incomplete screening of the supercurrent against the field penetration. In this study, we report that the confinement of a magnetic field occurs regardless of the dimensionality of the system, extending even to 1D linear potential systems. By utilizing a vector-field magnetic force microscope, we successfully create a vortex-antivortex pair connected by a 1D unquantized magnetic flux in ultra-thin superconducting films. Through an investigation of the manipulation and thermal behavior of the vortex pair, we uncover a long-range interaction mediated by the unquantized magnetic flux. These findings suggest a universal phenomenon of unquantized magnetic flux formation, independent of the geometry of the system. Our results present an experimental route for probing the impact of confinement on superconducting properties and order parameters in unconventional superconductors characterized by extremely low dimensionality.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Comparing Quantum Annealing and Spiking Neuromorphic Computing for Sampling Binary Sparse Coding QUBO Problems
Authors:
Kyle Henke,
Elijah Pelofske,
Garrett Kenyon,
Georg Hahn
Abstract:
We consider the problem of computing a sparse binary representation of an image. To be precise, given an image and an overcomplete, non-orthonormal basis, we aim to find a sparse binary vector indicating the minimal set of basis vectors that when added together best reconstruct the given input. We formulate this problem with an $L_2$ loss on the reconstruction error, and an $L_0$ (or, equivalently…
▽ More
We consider the problem of computing a sparse binary representation of an image. To be precise, given an image and an overcomplete, non-orthonormal basis, we aim to find a sparse binary vector indicating the minimal set of basis vectors that when added together best reconstruct the given input. We formulate this problem with an $L_2$ loss on the reconstruction error, and an $L_0$ (or, equivalently, an $L_1$) loss on the binary vector enforcing sparsity. This yields a quadratic binary optimization problem (QUBO), whose optimal solution(s) in general is NP-hard to find. The method of unsupervised and unnormalized dictionary feature learning for a desired sparsity level to best match the data is presented. Next, we solve the sparse representation QUBO by implementing it both on a D-Wave quantum annealer with Pegasus chip connectivity via minor embedding, as well as on the Intel Loihi 2 spiking neuromorphic processor. On the quantum annealer, we sample from the sparse representation QUBO using parallel quantum annealing combined with quantum evolution Monte Carlo, also known as iterated reverse annealing. On Loihi 2, we use a stochastic winner take all network of neurons. The solutions are benchmarked against simulated annealing, a classical heuristic, and the optimal solutions are computed using CPLEX. Iterated reverse quantum annealing performs similarly to simulated annealing, although simulated annealing is always able to sample the optimal solution whereas quantum annealing was not always able to. The Loihi 2 solutions that are sampled are on average more sparse than the solutions from any of the other methods. Loihi 2 outperforms a D-Wave quantum annealer standard linear-schedule anneal, while iterated reverse quantum annealing performs much better than both unmodified linear-schedule quantum annealing and iterated warm starting on Loihi 2.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
New Beam Dynamics Code for Cyclotron Analysis
Authors:
G-H. Kim,
H-J. Cho,
B-H. Oh,
G-R. Hahn,
M. Chung,
S. Park,
S. Shin
Abstract:
This paper describes the beam dynamic simulation with transfer matrix method for cyclotron. Starting from a description on the equation of motion in the cyclotron, lattice functions were determined from transfer matrix method and the solutions for the 2nd-order nonlinear Hamiltonian were introduced and used in phase space particle tracking. Based on the description of beam dynamics in the cyclotro…
▽ More
This paper describes the beam dynamic simulation with transfer matrix method for cyclotron. Starting from a description on the equation of motion in the cyclotron, lattice functions were determined from transfer matrix method and the solutions for the 2nd-order nonlinear Hamiltonian were introduced and used in phase space particle tracking. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design.
△ Less
Submitted 19 January, 2024;
originally announced January 2024.
-
Penalized Principal Component Analysis using Nesterov Smoothing
Authors:
Rebecca M. Hurwitz,
Georg Hahn
Abstract:
Principal components computed via PCA (principal component analysis) are traditionally used to reduce dimensionality in genomic data or to correct for population stratification. In this paper, we explore the penalized eigenvalue problem (PEP) which reformulates the computation of the first eigenvector as an optimization problem and adds an L1 penalty constraint. The contribution of our article is…
▽ More
Principal components computed via PCA (principal component analysis) are traditionally used to reduce dimensionality in genomic data or to correct for population stratification. In this paper, we explore the penalized eigenvalue problem (PEP) which reformulates the computation of the first eigenvector as an optimization problem and adds an L1 penalty constraint. The contribution of our article is threefold. First, we extend PEP by applying Nesterov smoothing to the original LASSO-type L1 penalty. This allows one to compute analytical gradients which enable faster and more efficient minimization of the objective function associated with the optimization problem. Second, we demonstrate how higher order eigenvectors can be calculated with PEP using established results from singular value decomposition (SVD). Third, using data from the 1000 Genome Project dataset, we empirically demonstrate that our proposed smoothed PEP allows one to increase numerical stability and obtain meaningful eigenvectors. We further investigate the utility of the penalized eigenvector approach over traditional PCA.
△ Less
Submitted 24 September, 2023;
originally announced September 2023.
-
Posiform Planting: Generating QUBO Instances for Benchmarking
Authors:
Georg Hahn,
Elijah Pelofske,
Hristo N. Djidjev
Abstract:
We are interested in benchmarking both quantum annealing and classical algorithms for minimizing Quadratic Unconstrained Binary Optimization (QUBO) problems. Such problems are NP-hard in general, implying that the exact minima of randomly generated instances are hard to find and thus typically unknown. While brute forcing smaller instances is possible, such instances are typically not interesting…
▽ More
We are interested in benchmarking both quantum annealing and classical algorithms for minimizing Quadratic Unconstrained Binary Optimization (QUBO) problems. Such problems are NP-hard in general, implying that the exact minima of randomly generated instances are hard to find and thus typically unknown. While brute forcing smaller instances is possible, such instances are typically not interesting due to being too easy for both quantum and classical algorithms. In this contribution, we propose a novel method, called posiform planting, for generating random QUBO instances of arbitrary size with known optimal solutions, and use those instances to benchmark the sampling quality of four D-Wave quantum annealers utilizing different interconnection structures (Chimera, Pegasus, and Zephyr hardware graphs) as well as the simulated annealing algorithm. Posiform planting differs from many existing methods in two key ways. It ensures the uniqueness of the planted optimal solution, thus avoiding groundstate degeneracy, and it enables the generation of QUBOs that are tailored to a given hardware connectivity structure, provided that the connectivity is not too sparse. Posiform planted QUBOs are a type of 2-SAT boolean satisfiability combinatorial optimization problems. Our experiments demonstrate the capability of the D-Wave quantum annealers to sample the optimal planted solution of combinatorial optimization problems with up to $5627$ qubits.
△ Less
Submitted 28 October, 2023; v1 submitted 10 August, 2023;
originally announced August 2023.
-
A Decade of Near-Infrared Variability in NGC4388: Insights into the AGN Structure
Authors:
Luis G. Dahmer-Hahn,
Alberto Rodríguez-Ardila,
Marina Bianchin,
Rogemar A. Riffel,
Rogério Riffel,
Thaisa Storchi-Bergmann,
Lei Hao
Abstract:
Variability studies have proven to be a powerful diagnostic tool for understanding the physics and properties of of Active Galactic Nuclei (AGNs). They provide insights into the spatial and temporal distribution of the emitting regions, the structure and dynamics of the accretion disk, and the properties of the central black hole. Here, we have analysed the K-band spectral variability of the Seyfe…
▽ More
Variability studies have proven to be a powerful diagnostic tool for understanding the physics and properties of of Active Galactic Nuclei (AGNs). They provide insights into the spatial and temporal distribution of the emitting regions, the structure and dynamics of the accretion disk, and the properties of the central black hole. Here, we have analysed the K-band spectral variability of the Seyfert 1.9/2 galaxy NGC4388 spanning five epochs over a period of ten years. We have performed spectral synthesis of the nuclear region and found that the contribution of warm dust (T~800K) declined by 88% during these 10 years. In the same period, the [CaVIII] coronal line decreased 61%, whereas BrG emission declined 35%. For the HeI and H2, we did not detect any significant variation beyond their uncertainties. Based on the time span of these changes, we estimate that the region where the warm dust is produced is smaller than 0.6pc, which suggests that this spectral feature comes from the innermost part of the region sampled, directly from the AGN torus. On the other hand, the bulk of [CaVIII] is produced in the inner ~2pc and the nuclear BrG region is more extended, spanning a region larger than 3pc. Lastly, HeI and H2 are even more external, with most of the emission probably being produced in the host galaxy rather than in the AGN. This is the first spectroscopic variability study in the NIR for an AGN where the central source is not directly visible.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Sampling binary sparse coding QUBO models using a spiking neuromorphic processor
Authors:
Kyle Henke,
Elijah Pelofske,
Georg Hahn,
Garrett T. Kenyon
Abstract:
We consider the problem of computing a sparse binary representation of an image. To be precise, given an image and an overcomplete, non-orthonormal basis, we aim to find a sparse binary vector indicating the minimal set of basis vectors that when added together best reconstruct the given input. We formulate this problem with an $L_2$ loss on the reconstruction error, and an $L_0$ (or, equivalently…
▽ More
We consider the problem of computing a sparse binary representation of an image. To be precise, given an image and an overcomplete, non-orthonormal basis, we aim to find a sparse binary vector indicating the minimal set of basis vectors that when added together best reconstruct the given input. We formulate this problem with an $L_2$ loss on the reconstruction error, and an $L_0$ (or, equivalently, an $L_1$) loss on the binary vector enforcing sparsity. This yields a so-called Quadratic Unconstrained Binary Optimization (QUBO) problem, whose solution is generally NP-hard to find. The contribution of this work is twofold. First, the method of unsupervised and unnormalized dictionary feature learning for a desired sparsity level to best match the data is presented. Second, the binary sparse coding problem is then solved on the Loihi 1 neuromorphic chip by the use of stochastic networks of neurons to traverse the non-convex energy landscape. The solutions are benchmarked against the classical heuristic simulated annealing. We demonstrate neuromorphic computing is suitable for sampling low energy solutions of binary sparse coding QUBO models, and although Loihi 1 is capable of sampling very sparse solutions of the QUBO models, there needs to be improvement in the implementation in order to be competitive with simulated annealing.
△ Less
Submitted 2 August, 2023; v1 submitted 2 June, 2023;
originally announced June 2023.
-
Initial State Encoding via Reverse Quantum Annealing and h-gain Features
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
Quantum annealing is a specialized type of quantum computation that aims to use quantum fluctuations in order to obtain global minimum solutions of combinatorial optimization problems. Programmable D-Wave quantum annealers are available as cloud computing resources which allow users low level access to quantum annealing control features. In this paper, we are interested in improving the quality of…
▽ More
Quantum annealing is a specialized type of quantum computation that aims to use quantum fluctuations in order to obtain global minimum solutions of combinatorial optimization problems. Programmable D-Wave quantum annealers are available as cloud computing resources which allow users low level access to quantum annealing control features. In this paper, we are interested in improving the quality of the solutions returned by a quantum annealer by encoding an initial state into the annealing process. We explore two D-Wave features allowing one to encode such an initial state: the reverse annealing and the h-gain features. Reverse annealing (RA) aims to refine a known solution following an anneal path starting with a classical state representing a good solution, going backwards to a point where a transverse field is present, and then finishing the annealing process with a forward anneal. The h-gain (HG) feature allows one to put a time-dependent weighting scheme on linear ($h$) biases of the Hamiltonian, and we demonstrate that this feature likewise can be used to bias the annealing to start from an initial state. We also consider a hybrid method consisting of a backward phase resembling RA, and a forward phase using the HG initial state encoding. Importantly, we investigate the idea of iteratively applying RA and HG to a problem, with the goal of monotonically improving on an initial state that is not optimal. The HG encoding technique is evaluated on a variety of input problems including the edge-weighted Maximum Cut problem and the vertex-weighted Maximum Clique problem, demonstrating that the HG technique is a viable alternative to RA for some problems. We also investigate how the iterative procedures perform for both RA and HG initial state encoding on random whole-chip spin glasses with the native hardware connectivity of the D-Wave Chimera and Pegasus chips.
△ Less
Submitted 24 September, 2023; v1 submitted 23 March, 2023;
originally announced March 2023.
-
The James Webb Space Telescope Mission: Optical Telescope Element Design, Development, and Performance
Authors:
Michael W. McElwain,
Lee D. Feinberg,
Marshall D. Perrin,
Mark Clampin,
C. Matt Mountain,
Matthew D. Lallo,
Charles-Philippe Lajoie,
Randy A. Kimble,
Charles W. Bowers,
Christopher C. Stark,
D. Scott Acton,
Ken Aiello,
Charles Atkinson,
Beth Barinek,
Allison Barto,
Scott Basinger,
Tracy Beck,
Matthew D. Bergkoetter,
Marcel Bluth,
Rene A. Boucarut,
Gregory R. Brady,
Keira J. Brooks,
Bob Brown,
John Byard,
Larkin Carey
, et al. (104 additional authors not shown)
Abstract:
The James Webb Space Telescope (JWST) is a large, infrared space telescope that has recently started its science program which will enable breakthroughs in astrophysics and planetary science. Notably, JWST will provide the very first observations of the earliest luminous objects in the Universe and start a new era of exoplanet atmospheric characterization. This transformative science is enabled by…
▽ More
The James Webb Space Telescope (JWST) is a large, infrared space telescope that has recently started its science program which will enable breakthroughs in astrophysics and planetary science. Notably, JWST will provide the very first observations of the earliest luminous objects in the Universe and start a new era of exoplanet atmospheric characterization. This transformative science is enabled by a 6.6 m telescope that is passively cooled with a 5-layer sunshield. The primary mirror is comprised of 18 controllable, low areal density hexagonal segments, that were aligned and phased relative to each other in orbit using innovative image-based wavefront sensing and control algorithms. This revolutionary telescope took more than two decades to develop with a widely distributed team across engineering disciplines. We present an overview of the telescope requirements, architecture, development, superb on-orbit performance, and lessons learned. JWST successfully demonstrates a segmented aperture space telescope and establishes a path to building even larger space telescopes.
△ Less
Submitted 4 January, 2023;
originally announced January 2023.
-
Noise Dynamics of Quantum Annealers: Estimating the Effective Noise Using Idle Qubits
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealing is a type of analog computation that aims to use quantum mechanical fluctuations in search of optimal solutions of QUBO (quadratic unconstrained binary optimization) or, equivalently, Ising problems. Since NP-hard problems can in general be mapped to Ising and QUBO formulations, the quantum annealing paradigm has the potential to help solve various NP-hard problems. Current quant…
▽ More
Quantum annealing is a type of analog computation that aims to use quantum mechanical fluctuations in search of optimal solutions of QUBO (quadratic unconstrained binary optimization) or, equivalently, Ising problems. Since NP-hard problems can in general be mapped to Ising and QUBO formulations, the quantum annealing paradigm has the potential to help solve various NP-hard problems. Current quantum annealers, such as those manufactured by D-Wave Systems, Inc., have various practical limitations including the size (number of qubits) of the problem that can be solved, the qubit connectivity, and error due to the environment or system calibration, which can reduce the quality of the solutions. Typically, for an arbitrary problem instance, the corresponding QUBO (or Ising) structure will not natively embed onto the available qubit architecture on the quantum chip. Thus, in these cases, a minor embedding of the problem structure onto the device is necessary. However, minor embeddings on these devices do not always make use of the full sparse chip hardware graph, and a large portion of the available qubits stay unused during quantum annealing. In this work, we embed a disjoint random QUBO on the unused parts of the chip alongside the QUBO to be solved, which acts as an indicator of the solution quality of the device over time. Using experiments on three different D-Wave quantum annealers, we demonstrate that (i) long term trends in solution quality exist on the D-Wave device, and (ii) the unused qubits can be used to measure the current level of noise of the quantum system.
△ Less
Submitted 17 February, 2023; v1 submitted 12 September, 2022;
originally announced September 2022.
-
Solving Larger Maximum Clique Problems Using Parallel Quantum Annealing
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealing has the potential to find low energy solutions of NP-hard problems that can be expressed as quadratic unconstrained binary optimization problems. However, the hardware of the quantum annealer manufactured by D-Wave Systems, which we consider in this work, is sparsely connected and moderately sized (on the order of thousands of qubits), thus necessitating a minor-embedding of a lo…
▽ More
Quantum annealing has the potential to find low energy solutions of NP-hard problems that can be expressed as quadratic unconstrained binary optimization problems. However, the hardware of the quantum annealer manufactured by D-Wave Systems, which we consider in this work, is sparsely connected and moderately sized (on the order of thousands of qubits), thus necessitating a minor-embedding of a logical problem onto the physical qubit hardware. The combination of relatively small hardware sizes and the necessity of a minor-embedding can mean that solving large optimization problems is not possible on current quantum annealers. In this research, we show that a hybrid approach combining parallel quantum annealing with graph decomposition allows one to solve larger optimization problem accurately. We apply the approach on the Maximum Clique problem on graphs with up to 120 nodes and 6395 edges.
△ Less
Submitted 14 March, 2023; v1 submitted 24 May, 2022;
originally announced May 2022.
-
Improving the Security of the IEEE 802.15.6 Standard for Medical BANs
Authors:
Muhammad Ali Siddiqi,
Georg Hahn,
Said Hamdioui,
Wouter A. Serdijn,
Christos Strydis
Abstract:
A Medical Body Area Network (MBAN) is an ensemble of collaborating, potentially heterogeneous, medical devices located inside, on the surface of or around the human body with the objective of tackling one or multiple medical conditions of the MBAN host. These devices -- which are a special category of Wireless Body Area Networks (WBANs) -- collect, process and transfer medical data outside of the…
▽ More
A Medical Body Area Network (MBAN) is an ensemble of collaborating, potentially heterogeneous, medical devices located inside, on the surface of or around the human body with the objective of tackling one or multiple medical conditions of the MBAN host. These devices -- which are a special category of Wireless Body Area Networks (WBANs) -- collect, process and transfer medical data outside of the network, while in some cases they also administer medical treatment autonomously. Since communication is so pivotal to their operation, the newfangled IEEE 802.15.6 standard is aimed at the communication aspects of WBANs. It places a set of physical and communication constraints while it also includes association/disassociation protocols and security services that WBAN applications need to comply with. However, the security specifications put forward by the standard can be easily shown to be insufficient when considering realistic MBAN use cases and need further enhancements. The present work addresses these shortcomings by, first, providing a structured analysis of the IEEE 802.15.6 security features and, afterwards, proposing comprehensive and tangible recommendations on improving the standard's security.
△ Less
Submitted 27 June, 2022; v1 submitted 17 January, 2022;
originally announced January 2022.
-
Stellar populations in local AGNs: evidence for enhanced star formation in the inner 100pc
Authors:
L. G. Dahmer-Hahn,
R. Riffel,
A. Rodríguez-Ardila,
R. A. Riffel,
T. Storchi-Bergmann,
M. Marinello,
R. I. Davies,
L. Burtscher,
D. Ruschel-Dutra,
D. J. Rosario
Abstract:
In modern models and simulations of galactic evolution, the star formation in massive galaxies is regulated by an ad hoc active galactic nuclei (AGN) feedback process. However, the physics and the extension of such effects on the star formation history of galaxies is matter of vivid debate. In order to shed some light in the AGN effects over the star formation, we analyzed the inner 500$\times$500…
▽ More
In modern models and simulations of galactic evolution, the star formation in massive galaxies is regulated by an ad hoc active galactic nuclei (AGN) feedback process. However, the physics and the extension of such effects on the star formation history of galaxies is matter of vivid debate. In order to shed some light in the AGN effects over the star formation, we analyzed the inner 500$\times$500pc of a sample of 14 Seyfert galaxies using GMOS and MUSE integral field spectroscopy. We fitted the continuum spectra in order to derive stellar age, metallicity, velocity and velocity dispersion maps in each source. After stacking our sample and averaging their properties, we found that the contribution of young SP, as well as that of AGN featureless continuum both peak at the nucleus. The fraction of intermediate-age SPs is smaller in the nucleus if compared to outer regions, and the contribution of old SPs vary very little within our field of view (FoV). We also found no variation of velocity dispersion or metallicity within our FoV. Lastly, we detected an increase in the dust reddening towards the center of the galaxies. These results lead us to conclude that AGN phenomenon is usually related to a recent star formation episode in the circumnuclear region.
△ Less
Submitted 18 November, 2021;
originally announced November 2021.
-
Parallel Quantum Annealing
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealers of D-Wave Systems, Inc., offer an efficient way to compute high quality solutions of NP-hard problems. This is done by mapping a problem onto the physical qubits of the quantum chip, from which a solution is obtained after quantum annealing. However, since the connectivity of the physical qubits on the chip is limited, a minor embedding of the problem structure onto the chip is r…
▽ More
Quantum annealers of D-Wave Systems, Inc., offer an efficient way to compute high quality solutions of NP-hard problems. This is done by mapping a problem onto the physical qubits of the quantum chip, from which a solution is obtained after quantum annealing. However, since the connectivity of the physical qubits on the chip is limited, a minor embedding of the problem structure onto the chip is required. In this process, and especially for smaller problems, many qubits will stay unused. We propose a novel method, called parallel quantum annealing, to make better use of available qubits, wherein either the same or several independent problems are solved in the same annealing cycle of a quantum annealer, assuming enough physical qubits are available to embed more than one problem. Although the individual solution quality may be slightly decreased when solving several problems in parallel (as opposed to solving each problem separately), we demonstrate that our method may give dramatic speed-ups in terms of Time-to-Solution (TTS) for solving instances of the Maximum Clique problem when compared to solving each problem sequentially on the quantum annealer. Additionally, we show that solving a single Maximum Clique problem using parallel quantum annealing reduces the TTS significantly.
△ Less
Submitted 28 November, 2022; v1 submitted 10 November, 2021;
originally announced November 2021.
-
Quantum Annealing Algorithms for Boolean Tensor Networks
Authors:
Elijah Pelofske,
Georg Hahn,
Daniel O'Malley,
Hristo N. Djidjev,
Boian S. Alexandrov
Abstract:
Quantum annealers manufactured by D-Wave Systems, Inc., are computational devices capable of finding high-quality solutions of NP-hard problems. In this contribution, we explore the potential and effectiveness of such quantum annealers for computing Boolean tensor networks. Tensors offer a natural way to model high-dimensional data commonplace in many scientific fields, and representing a binary t…
▽ More
Quantum annealers manufactured by D-Wave Systems, Inc., are computational devices capable of finding high-quality solutions of NP-hard problems. In this contribution, we explore the potential and effectiveness of such quantum annealers for computing Boolean tensor networks. Tensors offer a natural way to model high-dimensional data commonplace in many scientific fields, and representing a binary tensor as a Boolean tensor network is the task of expressing a tensor containing categorical (i.e., {0, 1}) values as a product of low dimensional binary tensors. A Boolean tensor network is computed by Boolean tensor decomposition, and it is usually not exact. The aim of such decomposition is to minimize the given distance measure between the high-dimensional input tensor and the product of lower-dimensional (usually three-dimensional) tensors and matrices representing the tensor network. In this paper, we introduce and analyze three general algorithms for Boolean tensor networks: Tucker, Tensor Train, and Hierarchical Tucker networks. The computation of a Boolean tensor network is reduced to a sequence of Boolean matrix factorizations, which we show can be expressed as a quadratic unconstrained binary optimization problem suitable for solving on a quantum annealer. By using a novel method we introduce called \textit{parallel quantum annealing}, we demonstrate that tensor with up to millions of elements can be decomposed efficiently using a DWave 2000Q quantum annealer.
△ Less
Submitted 27 March, 2022; v1 submitted 28 July, 2021;
originally announced July 2021.
-
Using machine learning for quantum annealing accuracy prediction
Authors:
Aaron Barbosa,
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealers, such as the device built by D-Wave Systems, Inc., offer a way to compute solutions of NP-hard problems that can be expressed in Ising or QUBO (quadratic unconstrained binary optimization) form. Although such solutions are typically of very high quality, problem instances are usually not solved to optimality due to imperfections of the current generations quantum annealers. In th…
▽ More
Quantum annealers, such as the device built by D-Wave Systems, Inc., offer a way to compute solutions of NP-hard problems that can be expressed in Ising or QUBO (quadratic unconstrained binary optimization) form. Although such solutions are typically of very high quality, problem instances are usually not solved to optimality due to imperfections of the current generations quantum annealers. In this contribution, we aim to understand some of the factors contributing to the hardness of a problem instance, and to use machine learning models to predict the accuracy of the D-Wave 2000Q annealer for solving specific problems. We focus on the Maximum Clique problem, a classic NP-hard problem with important applications in network analysis, bioinformatics, and computational chemistry. By training a machine learning classification model on basic problem characteristics such as the number of edges in the graph, or annealing parameters such as D-Wave's chain strength, we are able to rank certain features in the order of their contribution to the solution hardness, and present a simple decision tree which allows to predict whether a problem will be solvable to optimality with the D-Wave 2000Q. We extend these results by training a machine learning regression model that predicts the clique size found by D-Wave.
△ Less
Submitted 31 May, 2021;
originally announced June 2021.
-
Witness electron beam injection using an active plasma lens for a proton beam-driven plasma wakefield accelerator
Authors:
S. -Y. Kim,
K. Moon,
M. Chung,
K. N. Sjobak,
E. Adli,
S. Doebert,
M. Dayyani,
E. S. Yoon,
I. Nam,
G. Hahn
Abstract:
An active plasma lens focuses the beam in both the horizontal and vertical planes simultaneously using a magnetic field generated by a discharge current through the plasma. A beam size of 5--10 $μ$m can be achieved within a short distance using a focusing gradient on the order of 100 T/m. The active plasma lens is therefore an attractive element for plasma wakefield acceleration, because an ultra-…
▽ More
An active plasma lens focuses the beam in both the horizontal and vertical planes simultaneously using a magnetic field generated by a discharge current through the plasma. A beam size of 5--10 $μ$m can be achieved within a short distance using a focusing gradient on the order of 100 T/m. The active plasma lens is therefore an attractive element for plasma wakefield acceleration, because an ultra-small size of the witness electron beam is required for injection into the plasma wakefield to minimize emittance growth and to enhance the capturing efficiency. When the drive beam and witness electron beam co-propagate through the active plasma lens, interactions between the drive and witness beams, and the plasma must be considered. In this paper, through particle-in-cell simulations, we discuss the possibility of using an active plasma lens for the final focusing of the electron beam for the AWAKE RUN 2 experiments. It is confirmed that the amplitude of the plasma wakefield excited by proton bunches remains the same even after propagation through the active plasma lens. The emittance of the witness electron beam increases rapidly in the plasma density ramp regions of the lens. Nevertheless, when the witness electron beam has a charge of 100 pC, emittance of 10 mm mrad, and bunch length of 60 $μ$m, its emittance growth is not significant along the active plasma lens. For small emittance, such as 2 mm mrad, the emittance growth is found to be strongly dependent on the RMS beam size, plasma density, and multiple Coulomb scattering.
△ Less
Submitted 10 December, 2021; v1 submitted 20 April, 2021;
originally announced April 2021.
-
Boolean Hierarchical Tucker Networks on Quantum Annealers
Authors:
Elijah Pelofske,
Georg Hahn,
Daniel O'Malley,
Hristo N. Djidjev,
Boian S. Alexandrov
Abstract:
Quantum annealing is an emerging technology with the potential to solve some of the computational challenges that remain unresolved as we approach an era beyond Moore's Law. In this work, we investigate the capabilities of the quantum annealers of D-Wave Systems, Inc., for computing a certain type of Boolean tensor decomposition called Boolean Hierarchical Tucker Network (BHTN). Boolean tensor dec…
▽ More
Quantum annealing is an emerging technology with the potential to solve some of the computational challenges that remain unresolved as we approach an era beyond Moore's Law. In this work, we investigate the capabilities of the quantum annealers of D-Wave Systems, Inc., for computing a certain type of Boolean tensor decomposition called Boolean Hierarchical Tucker Network (BHTN). Boolean tensor decomposition problems ask for finding a decomposition of a high-dimensional tensor with categorical, [true, false], values, as a product of smaller Boolean core tensors. As the BHTN decompositions are usually not exact, we aim to approximate an input high-dimensional tensor by a product of lower-dimensional tensors such that the difference between both is minimized in some norm. We show that BHTN can be calculated as a sequence of optimization problems suitable for the D-Wave 2000Q quantum annealer. Although current technology is still fairly restricted in the problems they can address, we show that a complex problem such as BHTN can be solved efficiently and accurately.
△ Less
Submitted 12 March, 2021;
originally announced March 2021.
-
Reducing quantum annealing biases for solving the graph partitioning problem
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealers offer an efficient way to compute high quality solutions of NP-hard problems when expressed in a QUBO (quadratic unconstrained binary optimization) or an Ising form. This is done by mapping a problem onto the physical qubits and couplers of the quantum chip, from which a solution is read after a process called quantum annealing. However, this process is subject to multiple source…
▽ More
Quantum annealers offer an efficient way to compute high quality solutions of NP-hard problems when expressed in a QUBO (quadratic unconstrained binary optimization) or an Ising form. This is done by mapping a problem onto the physical qubits and couplers of the quantum chip, from which a solution is read after a process called quantum annealing. However, this process is subject to multiple sources of biases, including poor calibration, leakage between adjacent qubits, control biases, etc., which might negatively influence the quality of the annealing results. In this work, we aim at mitigating the effect of such biases for solving constrained optimization problems, by offering a two-step method, and apply it to Graph Partitioning. In the first step, we measure and reduce any biases that result from implementing the constraints of the problem. In the second, we add the objective function to the resulting bias-corrected implementation of the constraints, and send the problem to the quantum annealer. We apply this concept to Graph Partitioning, an important NP-hard problem, which asks to find a partition of the vertices of a graph that is balanced (the constraint) and minimizes the cut size (the objective). We first quantify the bias of the implementation of the constraint on the quantum annealer, that is, we require, in an unbiased implementation, that any two vertices have the same likelihood of being assigned to the same or to different parts of the partition. We then propose an iterative method to correct any such biases. We demonstrate that, after adding the objective, solving the resulting bias-corrected Ising problem on the quantum annealer results in a higher solution accuracy.
△ Less
Submitted 8 March, 2021;
originally announced March 2021.
-
Optimizing embedding-related quantum annealing parameters for reducing hardware bias
Authors:
Aaron Barbosa,
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Quantum annealers have been designed to propose near-optimal solutions to NP-hard optimization problems. However, the accuracy of current annealers such as the ones of D-Wave Systems, Inc., is limited by environmental noise and hardware biases. One way to deal with these imperfections and to improve the quality of the annealing results is to apply a variety of pre-processing techniques such as spi…
▽ More
Quantum annealers have been designed to propose near-optimal solutions to NP-hard optimization problems. However, the accuracy of current annealers such as the ones of D-Wave Systems, Inc., is limited by environmental noise and hardware biases. One way to deal with these imperfections and to improve the quality of the annealing results is to apply a variety of pre-processing techniques such as spin reversal (SR), anneal offsets (AO), or chain weights (CW). Maximizing the effectiveness of these techniques involves performing optimizations over a large number of parameters, which would be too costly if needed to be done for each new problem instance. In this work, we show that the aforementioned parameter optimization can be done for an entire class of problems, given each instance uses a previously chosen fixed embedding. Specifically, in the training phase, we fix an embedding E of a complete graph onto the hardware of the annealer, and then run an optimization algorithm to tune the following set of parameter values: the set of bits to be flipped for SR, the specific qubit offsets for AO, and the distribution of chain weights, optimized over a set of training graphs randomly chosen from that class, where the graphs are embedded onto the hardware using E. In the testing phase, we estimate how well the parameters computed during the training phase work on a random selection of other graphs from that class. We investigate graph instances of varying densities for the Maximum Clique, Maximum Cut, and Graph Partitioning problems. Our results indicate that, compared to their default behavior, substantial improvements of the annealing results can be achieved by using the optimized parameters for SR, AO, and CW.
△ Less
Submitted 1 December, 2020; v1 submitted 1 November, 2020;
originally announced November 2020.
-
Decomposition algorithms for solving NP-hard problems on a quantum annealer
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
NP-hard problems such as the maximum clique or minimum vertex cover problems, two of Karp's 21 NP-hard problems, have several applications in computational chemistry, biochemistry and computer network security. Adiabatic quantum annealers can search for the optimum value of such NP-hard optimization problems, given the problem can be embedded on their hardware. However, this is often not possible…
▽ More
NP-hard problems such as the maximum clique or minimum vertex cover problems, two of Karp's 21 NP-hard problems, have several applications in computational chemistry, biochemistry and computer network security. Adiabatic quantum annealers can search for the optimum value of such NP-hard optimization problems, given the problem can be embedded on their hardware. However, this is often not possible due to certain limitations of the hardware connectivity structure of the annealer. This paper studies a general framework for a decomposition algorithm for NP-hard graph problems aiming to identify an optimal set of vertices. Our generic algorithm allows us to recursively divide an instance until the generated subproblems can be embedded on the quantum annealer hardware and subsequently solved. The framework is applied to the maximum clique and minimum vertex cover problems, and we propose several pruning and reduction techniques to speed up the recursive decomposition. The performance of both algorithms is assessed in a detailed simulation study.
△ Less
Submitted 10 September, 2020;
originally announced September 2020.
-
Inferring the Dynamics of the State Evolution During Quantum Annealing
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
To solve an optimization problem using a commercial quantum annealer, one has to represent the problem of interest as an Ising or a quadratic unconstrained binary optimization (QUBO) problem and submit its coefficients to the annealer, which then returns a user-specified number of low-energy solutions. It would be useful to know what happens in the quantum processor during the anneal process so th…
▽ More
To solve an optimization problem using a commercial quantum annealer, one has to represent the problem of interest as an Ising or a quadratic unconstrained binary optimization (QUBO) problem and submit its coefficients to the annealer, which then returns a user-specified number of low-energy solutions. It would be useful to know what happens in the quantum processor during the anneal process so that one could design better algorithms or suggest improvements to the hardware. However, existing quantum annealers are not able to directly extract such information from the processor. Hence, in this work we propose to use advanced features of D-Wave 2000Q to indirectly infer information about the dynamics of the state evolution during the anneal process. Specifically, D-Wave 2000Q allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the anneal. Using this feature, we design a set of modified anneal schedules whose outputs can be used to generate information about the states of the system at user-defined time points during a standard anneal. With this process, called "slicing", we obtain approximate distributions of lowest-energy anneal solutions as the anneal time evolves. We use our technique to obtain a variety of insights into the annealer, such as the state evolution during annealing, when individual bits in an evolving solution flip during the anneal process and when they stabilize, and we introduce a technique to estimate the freeze-out point of both the system as well as of individual qubits.
△ Less
Submitted 19 November, 2020; v1 submitted 10 September, 2020;
originally announced September 2020.
-
Advanced unembedding techniques for quantum annealers
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
The D-Wave quantum annealers make it possible to obtain high quality solutions of NP-hard problems by mapping a problem in a QUBO (quadratic unconstrained binary optimization) or Ising form to the physical qubit connectivity structure on the D-Wave chip. However, the latter is restricted in that only a fraction of all pairwise couplers between physical qubits exists. Modeling the connectivity stru…
▽ More
The D-Wave quantum annealers make it possible to obtain high quality solutions of NP-hard problems by mapping a problem in a QUBO (quadratic unconstrained binary optimization) or Ising form to the physical qubit connectivity structure on the D-Wave chip. However, the latter is restricted in that only a fraction of all pairwise couplers between physical qubits exists. Modeling the connectivity structure of a given problem instance thus necessitates the computation of a minor embedding of the variables in the problem specification onto the logical qubits, which consist of several physical qubits "chained" together to act as a logical one. After annealing, it is however not guaranteed that all chained qubits get the same value (-1 or +1 for an Ising model, and 0 or 1 for a QUBO), and several approaches exist to assign a final value to each logical qubit (a process called "unembedding"). In this work, we present tailored unembedding techniques for four important NP-hard problems: the Maximum Clique, Maximum Cut, Minimum Vertex Cover, and Graph Partitioning problems. Our techniques are simple and yet make use of structural properties of the problem being solved. Using Erdős-Rényi random graphs as inputs, we compare our unembedding techniques to three popular ones (majority vote, random weighting, and minimize energy). We demonstrate that our proposed algorithms outperform the currently available ones in that they yield solutions of better quality, while being computationally equally efficient.
△ Less
Submitted 27 December, 2022; v1 submitted 10 September, 2020;
originally announced September 2020.
-
Advanced anneal paths for improved quantum annealing
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
Advances in quantum annealing technology make it possible to obtain high quality approximate solutions of important NP-hard problems. With the newer generations of the D-Wave annealer, more advanced features are available which allow the user to have greater control of the anneal process. In this contribution, we study how such features can help in improving the quality of the solutions returned b…
▽ More
Advances in quantum annealing technology make it possible to obtain high quality approximate solutions of important NP-hard problems. With the newer generations of the D-Wave annealer, more advanced features are available which allow the user to have greater control of the anneal process. In this contribution, we study how such features can help in improving the quality of the solutions returned by the annealer. Specifically, we focus on two of these features: reverse annealing and h-gain. Reverse annealing (RA) was designed to allow refining a known solution by backward annealing from a classical state representing the solution to a mid-anneal point where a transverse field is present, followed by an ordinary forward anneal, which is hoped to improve on the previous solution. The h-gain (HG) feature stands for time-dependent gain in Hamiltonian linear ($h$) biases and was originally developed to help study freezeout times and phase transitions in spin glasses. Here we apply HG to bias the quantum state in the beginning of the annealing process towards the known solution as in the RA case, but using a different apparatus. We also investigate a hybrid reverse annealing/h-gain schedule, which has a backward phase resembling an RA step and whose forward phase uses the HG idea. To optimize the parameters of the schedules, we employ a Bayesian optimization framework. We test all techniques on a variety of input problems including the weighted Maximum Cut problem and the weighted Maximum Clique problem. Our results show that each technique may dominate the others depending on the input instance, and that the HG technique is a viable alternative to RA for some problems.
△ Less
Submitted 23 October, 2020; v1 submitted 10 September, 2020;
originally announced September 2020.
-
Siblings of countable cographs
Authors:
Gena Hahn,
Maurice Pouzet,
Robert Woodrow
Abstract:
We show that every countable cograph has either one or infinitely many siblings. This answers, very partially, a conjecture of Thomassé. The main tools are the notion of well quasi ordering and the correspondence between cographs and some labelled ordered trees.
We show that every countable cograph has either one or infinitely many siblings. This answers, very partially, a conjecture of Thomassé. The main tools are the notion of well quasi ordering and the correspondence between cographs and some labelled ordered trees.
△ Less
Submitted 26 April, 2020;
originally announced April 2020.
-
Influence of the Yarkovsky force on Jupiter Trojan asteroids
Authors:
Stephan Hellmich,
Stefano Mottola,
Gerhard Hahn,
Ekkehard Kührt,
Detlef de Niem
Abstract:
Aims. We investigate the influence of the Yarkovsky force on the long-term orbital evolution of Jupiter Trojan asteroids. Methods. Clones of the observed population with different sizes and different thermal properties were numerically integrated for 1 Gyr with and without the Yarkovsky effect. The escape rate of these objects from the Trojan region as well as changes in the libration amplitude, e…
▽ More
Aims. We investigate the influence of the Yarkovsky force on the long-term orbital evolution of Jupiter Trojan asteroids. Methods. Clones of the observed population with different sizes and different thermal properties were numerically integrated for 1 Gyr with and without the Yarkovsky effect. The escape rate of these objects from the Trojan region as well as changes in the libration amplitude, eccentricity, and inclination were used as a metric of the strength of the Yarkovsky effect on the Trojan orbits. Results. Objects with radii $R\leq$1 km are significantly influenced by the Yarkovsky force. The effect causes a depletion of these objects over timescales of a few hundred million years. As a consequence, we expect the size-frequency distribution of small Trojans to show a shallower slope than that of the currently observable population ($R$ $\gtrsim$ 1 km), with a turning point between $R$ = 100 m and $R$ = 1 km. The effect of the Yarkovsky acceleration on the orbits of Trojans depends on the sense of rotation in a complex way. The libration amplitude of prograde rotators decreases with time while the eccentricity increases. Retrograde rotators experience the opposite effect, which results in retrograde rotators being ejected faster from the 1:1 resonance region. Furthermore, for objects affected by the Yarkovsky force, we find indications that the effect tends to smooth out the differences in the orbital distribution between the two clouds.
△ Less
Submitted 19 September, 2019; v1 submitted 18 September, 2019;
originally announced September 2019.
-
A panchromatic spatially resolved study of the inner 500pc of NGC1052 -- II: Gas excitation and kinematics
Authors:
Luis G. Dahmer-Hahn,
Rogério Riffel,
Tiago V. Ricci,
João E. Steiner,
Thaisa Storchi-Bergmann,
Rogemar A. Riffel,
Roberto B. Menezes,
Natacha Z. Dametto,
Marlon R. Diniz,
Juliana C. Motter,
Daniel Ruschel-Dutra
Abstract:
We map the optical and near-infrared (NIR) emission-line flux distributions and kinematics of the inner 320$\times$535pc$^2$ of the elliptical galaxy NGC1052. The integral field spectra were obtained with the Gemini Telescope using the GMOS-IFU and NIFS instruments, with angular resolutions of 0''88 and 0''1 in the optical and NIR, respectively. We detect five kinematic components: (1 and 2) Two s…
▽ More
We map the optical and near-infrared (NIR) emission-line flux distributions and kinematics of the inner 320$\times$535pc$^2$ of the elliptical galaxy NGC1052. The integral field spectra were obtained with the Gemini Telescope using the GMOS-IFU and NIFS instruments, with angular resolutions of 0''88 and 0''1 in the optical and NIR, respectively. We detect five kinematic components: (1 and 2) Two spatially unresolved components, being a broad line region visible in H$α$, with a FWHM of $\sim$3200km s$^{-1}$ and an intermediate-broad component seen in the [OIII]$λλ$4959,5007 doublet; (3) an extended intermediate-width component with 280<FWHM<450km s$^{-1}$ and centroid velocities up to 400km s$^{-1}$, which dominates the flux in our data, attributed either to a bipolar outflow related to the jets, rotation in an eccentric disc or a combination of a disc and large-scale gas bubbles; (4 and 5) two narrow (FWHM<150km s$^{-1}$) components, one visible in [OIII], and one visible in the other emission lines, extending beyond the field-of-view of our data, which is attributed to large-scale shocks. Our results suggest that the ionization within the observed field of view cannot be explained by a single mechanism, with photoionization being the dominant mechanism in the nucleus with a combination of shocks and photoionization responsible for the extended ionization.
△ Less
Submitted 4 September, 2019;
originally announced September 2019.
-
Peering into the Anneal Process of a Quantum Annealer
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
Commercial adiabatic quantum annealers have the potential to solve important NP-hard optimization problems efficiently. The newest generation of those machines additionally allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the annealing. In this work we use the aforementioned feature of the D-Wave 2000Q to…
▽ More
Commercial adiabatic quantum annealers have the potential to solve important NP-hard optimization problems efficiently. The newest generation of those machines additionally allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the annealing. In this work we use the aforementioned feature of the D-Wave 2000Q to attempt to monitor how the anneal solution evolves during the anneal process. This process we call slicing: at each time slice during the anneal, we are able to obtain an approximate distribution of anneal solutions. We use our technique to obtain a variety of insights into the D-Wave 2000Q. For example, we observe when individual bits flip during the anneal process and when they stabilize, which allows us to determine the freeze-out point for each qubit individually. We highlight our results using both random QUBO (quadratic unconstrained binary optimization) instances and, for better visualization, instances which we specifically optimize (using our own genetic algorithm) to exhibit a pronounced evolution of its solution during the anneal.
△ Less
Submitted 11 September, 2019; v1 submitted 7 August, 2019;
originally announced August 2019.
-
Using Graph Partitioning for Scalable Distributed Quantum Molecular Dynamics
Authors:
Hristo N. Djidjev,
Georg Hahn,
Susan M. Mniszewski,
Christian F. A. Negre,
Anders M. N. Niklasson
Abstract:
The simulation of the physical movement of multi-body systems at an atomistic level, with forces calculated from a quantum mechanical description of the electrons, motivates a graph partitioning problem studied in this article. Several advanced algorithms relying on evaluations of matrix polynomials have been published in the literature for such simulations. We aim to use a special type of graph p…
▽ More
The simulation of the physical movement of multi-body systems at an atomistic level, with forces calculated from a quantum mechanical description of the electrons, motivates a graph partitioning problem studied in this article. Several advanced algorithms relying on evaluations of matrix polynomials have been published in the literature for such simulations. We aim to use a special type of graph partitioning in order to efficiently parallelize these computations. For this, we create a graph representing the zero-nonzero structure of a thresholded density matrix, and partition that graph into several components. Each separate submatrix (corresponding to each subgraph) is then substituted into the matrix polynomial, and the result for the full matrix polynomial is reassembled at the end from the individual polynomials. This paper starts by introducing a rigorous definition as well as a mathematical justification of this partitioning problem. We assess the performance of several methods to compute graph partitions with respect to both the quality of the partitioning and their runtime.
△ Less
Submitted 7 September, 2019; v1 submitted 26 June, 2019;
originally announced June 2019.
-
Optimizing the spin reversal transform on the D-Wave 2000Q
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo Djidjev
Abstract:
Commercial quantum annealers from D-Wave Systems make it possible to obtain approximate solutions of high quality for certain NP-hard problems in nearly constant time. Before solving a problem on D-Wave, several pre-processing methods can be applied, one of them being the so-called spin reversal or gauge transform. The spin reversal transform flips the sign of selected variables and coefficients o…
▽ More
Commercial quantum annealers from D-Wave Systems make it possible to obtain approximate solutions of high quality for certain NP-hard problems in nearly constant time. Before solving a problem on D-Wave, several pre-processing methods can be applied, one of them being the so-called spin reversal or gauge transform. The spin reversal transform flips the sign of selected variables and coefficients of the Ising or QUBO (quadratic unconstrained binary optimization) representation of the problem that D-Wave minimizes. The spin reversal transform leaves the ground state of the Ising model invariant, but can average out the biases induced through analog and systematic errors on the device, thus improving the quality of the solution that D-Wave returns. This work investigates the effectiveness of the spin reversal transform for D-Wave 2000Q. We consider two important NP-hard problems, the Maximum Clique and the Minimum Vertex Cover problems, and show on a variety of input problem graphs that using the spin reversal transform can yield substantial improvements in solution quality. In contrast to the native spin reversal built into D-Wave, we consider more general ways to reverse individual spins and we investigate the dependence on the problem type, on the spin reversal probability, and possible advantages of carrying out reversals on the qubit instead of the chain level. Most importantly, for a given individual problem, we use our findings to optimize the spin reversal transform using a genetic optimization algorithm.
△ Less
Submitted 7 September, 2019; v1 submitted 26 June, 2019;
originally announced June 2019.
-
Solving large Minimum Vertex Cover problems on a quantum annealer
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
We consider the minimum vertex cover problem having applications in e.g. biochemistry and network security. Quantum annealers can find the optimum solution of such NP-hard problems, given they can be embedded on the hardware. This is often infeasible due to limitations of the hardware connectivity structure. This paper presents a decomposition algorithm for the minimum vertex cover problem: The al…
▽ More
We consider the minimum vertex cover problem having applications in e.g. biochemistry and network security. Quantum annealers can find the optimum solution of such NP-hard problems, given they can be embedded on the hardware. This is often infeasible due to limitations of the hardware connectivity structure. This paper presents a decomposition algorithm for the minimum vertex cover problem: The algorithm recursively divides an arbitrary problem until the generated subproblems can be embedded and solved on the annealer. To speed up the decomposition, we propose several pruning and reduction techniques. The performance of our algorithm is assessed in a simulation study.
△ Less
Submitted 30 April, 2019; v1 submitted 29 March, 2019;
originally announced April 2019.
-
Solving large Maximum Clique problems on a quantum annealer
Authors:
Elijah Pelofske,
Georg Hahn,
Hristo N. Djidjev
Abstract:
Commercial quantum annealers from D-Wave Systems can find high quality solutions of quadratic unconstrained binary optimization problems that can be embedded onto its hardware. However, even though such devices currently offer up to 2048 qubits, due to limitations on the connectivity of those qubits, the size of problems that can typically be solved is rather small (around 65 variables). This limi…
▽ More
Commercial quantum annealers from D-Wave Systems can find high quality solutions of quadratic unconstrained binary optimization problems that can be embedded onto its hardware. However, even though such devices currently offer up to 2048 qubits, due to limitations on the connectivity of those qubits, the size of problems that can typically be solved is rather small (around 65 variables). This limitation poses a problem for using D-Wave machines to solve application-relevant problems, which can have thousands of variables. For the important Maximum Clique problem, this article investigates methods for decomposing larger problem instances into smaller ones, which can subsequently be solved on D-Wave. During the decomposition, we aim to prune as many generated subproblems that do not contribute to the solution as possible, in order to reduce the computational complexity. The reduction methods presented in this article include upper and lower bound heuristics in conjunction with graph decomposition, vertex and edge extraction, and persistency analysis.
△ Less
Submitted 9 May, 2019; v1 submitted 22 January, 2019;
originally announced January 2019.
-
A panchromatic spatially resolved study of the inner 500 pc of NGC 1052 - I: Stellar population
Authors:
Luis Gabriel Dahmer-Hahn,
Rogério Riffel,
João Evangelista Steiner,
Rogemar André Riffel,
Roberto Bertoldo Menezes,
Tiago Vecchi Ricci,
Natacha Zanon Dametto,
Thaisa Storchi-Bergmann,
Marlon Rodrigo Diniz
Abstract:
We map optical and near-infrared (NIR) stellar population properties of the inner 320$\times$535pc$^2$ of the elliptical galaxy NGC1052. The optical and NIR spectra were obtained using the Gemini Integral Field Units of the GMOS instrument and NIFS, respectively. By performing stellar population synthesis in the optical alone, we find that this region of the galaxy is dominated by old (t$>$10Gyr)…
▽ More
We map optical and near-infrared (NIR) stellar population properties of the inner 320$\times$535pc$^2$ of the elliptical galaxy NGC1052. The optical and NIR spectra were obtained using the Gemini Integral Field Units of the GMOS instrument and NIFS, respectively. By performing stellar population synthesis in the optical alone, we find that this region of the galaxy is dominated by old (t$>$10Gyr) stellar populations. Using the NIR, we find the nucleus to be dominated by old stellar populations, and a circumnuclear ring with younger ($\sim$2.5Gyr) stars. We also combined the optical and NIR datacubes and performed a panchromatic spatially resolved stellar population synthesis, which resulted in a dominance of older stellar populations, in agreement with optical results. We argue that the technique of combining optical and NIR data might be useful to isolate the contribution of stellar population ages with strong NIR absorption bands. We also derive the stellar kinematics and find that the stellar motions are dominated by a high ($\sim$240km$\cdot$s$^{-1}$) velocity dispersion in the nucleus, with stars also rotating around the center. Lastly, we measure the absorption bands, both in the optical and in the NIR, and find a nuclear drop in their equivalent widths. The favored explanation for this drop is a featureless continuum emission from the low luminosity active galactic nucleus.
△ Less
Submitted 7 November, 2018;
originally announced November 2018.
-
On the expected runtime of multiple testing algorithms with bounded error
Authors:
Georg Hahn
Abstract:
Consider testing multiple hypotheses in the setting where the p-values of all hypotheses are unknown and thus have to be approximated using Monte Carlo simulations. One class of algorithms published in the literature for this scenario provides guarantees on the correctness of their testing result through the computation of confidence statements on all approximated p-values. This article focuses on…
▽ More
Consider testing multiple hypotheses in the setting where the p-values of all hypotheses are unknown and thus have to be approximated using Monte Carlo simulations. One class of algorithms published in the literature for this scenario provides guarantees on the correctness of their testing result through the computation of confidence statements on all approximated p-values. This article focuses on the expected runtime of such algorithms and derives a variety of finite and infinite expected runtime results.
△ Less
Submitted 12 June, 2020; v1 submitted 27 July, 2018;
originally announced July 2018.
-
Determination of lattice functions and 2nd-order transfer matrix for High Power Cyclotron
Authors:
B-H. Oh,
C. Choi,
G-R. Hahn,
M. Chung,
S. Shin
Abstract:
This paper describes the development of beam dynamic simulation code for cyclotron. Starting from a description of beam dynamics in the cyclotron, lattice functions were determined and the solutions for the 2nd-order nonlinear Hamiltonian were revised. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design.
This paper describes the development of beam dynamic simulation code for cyclotron. Starting from a description of beam dynamics in the cyclotron, lattice functions were determined and the solutions for the 2nd-order nonlinear Hamiltonian were revised. Based on the description of beam dynamics in the cyclotron, simulation code was also developed for cyclotron design.
△ Less
Submitted 3 July, 2018;
originally announced July 2018.
-
Probing Evolutionary Population Synthesis Models in the Near Infrared with Early Type Galaxies
Authors:
Luis Gabriel Dahmer-Hahn,
Rogério Riffel,
Alberto Rodríguez-Ardila,
Lucimara P. Martins,
Carolina Kehrig,
Timothy M. Heckman,
Miriani G. Pastoriza,
Natacha Z. Dametto
Abstract:
We performed a near-infrared (NIR, $\sim$1.0$μ$m-2.4$μ$m) stellar population study in a sample of early type galaxies. The synthesis was performed using five different evolutionary population synthesis libraries of models. Our main results can be summarized as follows: low spectral resolution libraries are not able to produce reliable results when applied to the NIR alone, with each library findin…
▽ More
We performed a near-infrared (NIR, $\sim$1.0$μ$m-2.4$μ$m) stellar population study in a sample of early type galaxies. The synthesis was performed using five different evolutionary population synthesis libraries of models. Our main results can be summarized as follows: low spectral resolution libraries are not able to produce reliable results when applied to the NIR alone, with each library finding a different dominant population. The two newest higher resolution models, on the other hand, perform considerably better, finding consistent results to each other and to literature values. We also found that optical results are consistent with each other even for lower resolution models. We also compared optical and NIR results, and found out that lower resolution models tend to disagree in the optical and in the NIR, with higher fraction of young populations in the NIR and dust extinction $\sim$1 magnitude higher than optical values. For higher resolution models, optical and NIR results tend do aggree much better, suggesting that a higher spectral resolution is fundamental to improve the quality of the results.
△ Less
Submitted 22 February, 2018;
originally announced February 2018.
-
Task-based Parallel Computation of the Density Matrix in Quantum-based Molecular Dynamics using Graph Partitioning
Authors:
Purnima Ghale,
Matthew P. Kroonblawd,
Susan M. Mniszewski,
Christian F. A. Negre,
Robert Pavel,
Sergio Pino,
Vivek B. Sardeshmukh,
Guangjie Shi,
Georg Hahn
Abstract:
Quantum-based molecular dynamics (QMD) is a highly accurate and transferable method for material science simulations. However, the time scales and system sizes accessible to QMD are typically limited to picoseconds and a few hundred atoms. These constraints arise due to expensive self-consistent ground-state electronic structure calculations that can often scale cubically with the number of atoms.…
▽ More
Quantum-based molecular dynamics (QMD) is a highly accurate and transferable method for material science simulations. However, the time scales and system sizes accessible to QMD are typically limited to picoseconds and a few hundred atoms. These constraints arise due to expensive self-consistent ground-state electronic structure calculations that can often scale cubically with the number of atoms. Linearly scaling methods depend on computing the density matrix P from the Hamiltonian matrix H by exploiting the sparsity in both matrices. The second-order spectral projection (SP2) algorithm is an O(N) algorithm that computes P with a sequence of 40-50 matrix-matrix multiplications. In this paper, we present task-based implementations of a recently developed data-parallel graph-based approach to the SP2 algorithm, G-SP2. We represent the density matrix P as an undirected graph and use graph partitioning techniques to divide the computation into smaller independent tasks. The partitions thus obtained are generally not of equal size and give rise to undesirable load imbalances in standard MPI-based implementations. This load-balancing challenge can be mitigated by dynamically scheduling parallel computations at runtime using task-based programming models. We develop task-based implementations of the data-parallel G-SP2 algorithm using both Intel's Concurrent Collections (CnC) as well as the Charm++ programming model and evaluate these implementations for future use. Scaling and performance results of our implementations are investigated for representative segments of QMD simulations for solvated protein systems containing more than 10,000 atoms.
△ Less
Submitted 25 January, 2018;
originally announced January 2018.
-
Efficient Combinatorial Optimization Using Quantum Annealing
Authors:
Hristo N. Djidjev,
Guillaume Chapuis,
Georg Hahn,
Guillaume Rizk
Abstract:
The recent availability of the first commercial quantum computers has provided a promising tool to tackle NP hard problems which can only be solved heuristically with present techniques. However, it is unclear if the current state of quantum computing already provides a quantum advantage over the current state of the art in classical computing. This article assesses the performance of the D-Wave 2…
▽ More
The recent availability of the first commercial quantum computers has provided a promising tool to tackle NP hard problems which can only be solved heuristically with present techniques. However, it is unclear if the current state of quantum computing already provides a quantum advantage over the current state of the art in classical computing. This article assesses the performance of the D-Wave 2X quantum annealer on two NP hard graph problems, in particular clique finding and graph partitioning. For this, we provide formulations as Qubo and Ising Hamiltonians suitable for the quantum annealer and compare a variety of quantum solvers (Sapi, QBSolv, QSage provided by D-Wave Sys, Inc.) to current classical algorithms (METIS, Simulated Annealing, third-party clique finding and graph splitting heuristics) on certain test sets of graphs. We demonstrate that for small graph instances, classical methods still outperform the quantum annealer in terms of computing time, even though the quality of the best solutions obtained is comparable. Nevertheless, due to the limited problem size which can be embedded on the D-Wave 2X chip, the aforementioned finding applies to most of problems of general nature solvable on the quantum annealer. For instances specifically designed to fit the D-Wave 2X architecture, we observe substantial speed-ups in computing time over classical approaches.
△ Less
Submitted 29 January, 2018; v1 submitted 25 January, 2018;
originally announced January 2018.
-
Reducing Binary Quadratic Forms for More Scalable Quantum Annealing
Authors:
Georg Hahn,
Hristo N. Djidjev
Abstract:
Recent advances in the development of commercial quantum annealers such as the D-Wave 2X allow solving NP-hard optimization problems that can be expressed as quadratic unconstrained binary programs. However, the relatively small number of available qubits (around 1000 for the D-Wave 2X quantum annealer) poses a severe limitation to the range of problems that can be solved. This paper explores the…
▽ More
Recent advances in the development of commercial quantum annealers such as the D-Wave 2X allow solving NP-hard optimization problems that can be expressed as quadratic unconstrained binary programs. However, the relatively small number of available qubits (around 1000 for the D-Wave 2X quantum annealer) poses a severe limitation to the range of problems that can be solved. This paper explores the suitability of preprocessing methods for reducing the sizes of the input programs and thereby the number of qubits required for their solution on quantum computers. Such methods allow us to determine the value of certain variables that hold in either any optimal solution (called strong persistencies) or in at least one optimal solution (weak persistencies). We investigate preprocessing methods for two important NP-hard graph problems, the computation of a maximum clique and a maximum cut in a graph. We show that the identification of strong and weak persistencies for those two optimization problems is very instance-specific, but can lead to substantial reductions in the number of variables.
△ Less
Submitted 29 January, 2018; v1 submitted 25 January, 2018;
originally announced January 2018.
-
Finding Maximum Cliques on the D-Wave Quantum Annealer
Authors:
Guillaume Chapuis,
Hristo N. Djidjev,
Georg Hahn,
Guillaume Rizk
Abstract:
This paper assesses the performance of the D-Wave 2X (DW) quantum annealer for finding a maximum clique in a graph, one of the most fundamental and important NP-hard problems. Because the size of the largest graphs DW can directly solve is quite small (usually around 45 vertices), we also consider decomposition algorithms intended for larger graphs and analyze their performance. For smaller graphs…
▽ More
This paper assesses the performance of the D-Wave 2X (DW) quantum annealer for finding a maximum clique in a graph, one of the most fundamental and important NP-hard problems. Because the size of the largest graphs DW can directly solve is quite small (usually around 45 vertices), we also consider decomposition algorithms intended for larger graphs and analyze their performance. For smaller graphs that fit DW, we provide formulations of the maximum clique problem as a quadratic unconstrained binary optimization (QUBO) problem, which is one of the two input types (together with the Ising model) acceptable by the machine, and compare several quantum implementations to current classical algorithms such as simulated annealing, Gurobi, and third-party clique finding heuristics. We further estimate the contributions of the quantum phase of the quantum annealer and the classical post-processing phase typically used to enhance each solution returned by DW. We demonstrate that on random graphs that fit DW, no quantum speedup can be observed compared with the classical algorithms. On the other hand, for instances specifically designed to fit well the DW qubit interconnection network, we observe substantial speed-ups in computing time over classical approaches.
△ Less
Submitted 24 April, 2018; v1 submitted 25 January, 2018;
originally announced January 2018.
-
Bypass rewiring and extreme robustness of Eulerian networks
Authors:
Junsang Park,
Seungwon Shin,
Sang Geun Hahn
Abstract:
Bypass rewiring improves connectivity and robustness of networks against removal of nodes including failures and attacks. A concept of bypass rewiring on directed networks is proposed, and random bypass rewiring on infinite directed random networks is analytically and numerically investigated with simulations. As a result, random bypass rewiring makes infinite directed (undirected) random networks…
▽ More
Bypass rewiring improves connectivity and robustness of networks against removal of nodes including failures and attacks. A concept of bypass rewiring on directed networks is proposed, and random bypass rewiring on infinite directed random networks is analytically and numerically investigated with simulations. As a result, random bypass rewiring makes infinite directed (undirected) random networks extremely robust for arbitrary occupation probabilities if and only if in-degree of every node except a fixed number of nodes is equal to the out-degree (every node except a finite number of nodes has even degree); bypass rewiring makes the percolation threshold 0. Consequently, a finite network has a strongly connected spanning subnetwork which has an Eulerian path or cycle if and only if there exists an way of bypass rewiring to make the finite network extremely robust for every combination of removed nodes; Eulerian networks are extremely robust for every combination of removed nodes.
△ Less
Submitted 14 December, 2017; v1 submitted 31 October, 2017;
originally announced October 2017.
-
A lower bound on the size of an absorbing set in an arc-coloured tournament
Authors:
Laurent Beaudou,
Luc Devroye,
Gena Hahn
Abstract:
Bousquet, Lochet and Thomassé recently gave an elegant proof that for any integer $n$, there is a least integer $f(n)$ such that any tournament whose arcs are coloured with $n$ colours contains a subset of vertices $S$ of size $f(n)$ with the property that any vertex not in $S$ admits a monochromatic path to some vertex of $S$. In this note we provide a lower bound on the value $f(n)$.
Bousquet, Lochet and Thomassé recently gave an elegant proof that for any integer $n$, there is a least integer $f(n)$ such that any tournament whose arcs are coloured with $n$ colours contains a subset of vertices $S$ of size $f(n)$ with the property that any vertex not in $S$ admits a monochromatic path to some vertex of $S$. In this note we provide a lower bound on the value $f(n)$.
△ Less
Submitted 30 August, 2017; v1 submitted 29 August, 2017;
originally announced August 2017.
-
Implementing Monte Carlo Tests with P-value Buckets
Authors:
Axel Gandy,
Georg Hahn,
Dong Ding
Abstract:
Software packages usually report the results of statistical tests using p-values. Users often interpret these by comparing them to standard thresholds, e.g. 0.1%, 1% and 5%, which is sometimes reinforced by a star rating (***, **, *). We consider an arbitrary statistical test whose p-value p is not available explicitly, but can be approximated by Monte Carlo samples, e.g. by bootstrap or permutati…
▽ More
Software packages usually report the results of statistical tests using p-values. Users often interpret these by comparing them to standard thresholds, e.g. 0.1%, 1% and 5%, which is sometimes reinforced by a star rating (***, **, *). We consider an arbitrary statistical test whose p-value p is not available explicitly, but can be approximated by Monte Carlo samples, e.g. by bootstrap or permutation tests. The standard implementation of such tests usually draws a fixed number of samples to approximate p. However, the probability that the exact and the approximated p-value lie on different sides of a threshold (the resampling risk) can be high, particularly for p-values close to a threshold. We present a method to overcome this. We consider a finite set of user-specified intervals which cover [0,1] and which can be overlapping. We call these p-value buckets. We present algorithms that, with arbitrarily high probability, return a p-value bucket containing p. We prove that for both a bounded resampling risk and a finite runtime, overlapping buckets need to be employed, and that our methods both bound the resampling risk and guarantee a finite runtime for such overlapping buckets. To interpret decisions with overlapping buckets, we propose an extension of the star rating system. We demonstrate that our methods are suitable for use in standard software, including for low p-value thresholds occurring in multiple testing settings, and that they can be computationally more efficient than standard implementations.
△ Less
Submitted 4 November, 2019; v1 submitted 27 March, 2017;
originally announced March 2017.
-
A simple method for implementing Monte Carlo tests
Authors:
Dong Ding,
Axel Gandy,
Georg Hahn
Abstract:
We consider a statistical test whose p-value can only be approximated using Monte Carlo simulations. We are interested in deciding whether the p-value for an observed data set lies above or below a given threshold such as 5%. We want to ensure that the resampling risk, the probability of the (Monte Carlo) decision being different from the true decision, is uniformly bounded. This article introduce…
▽ More
We consider a statistical test whose p-value can only be approximated using Monte Carlo simulations. We are interested in deciding whether the p-value for an observed data set lies above or below a given threshold such as 5%. We want to ensure that the resampling risk, the probability of the (Monte Carlo) decision being different from the true decision, is uniformly bounded. This article introduces a simple open-ended method with this property, the confidence sequence method (CSM). We compare our approach to another algorithm, SIMCTEST, which also guarantees an (asymptotic) uniform bound on the resampling risk, as well as to other Monte Carlo procedures without a uniform bound. CSM is free of tuning parameters and conservative. It has the same theoretical guarantee as SIMCTEST and, in many settings, similar stopping boundaries. As it is much simpler than other methods, CSM is a useful method for practical applications.
△ Less
Submitted 9 October, 2019; v1 submitted 5 November, 2016;
originally announced November 2016.
-
Closure properties of classes of multiple testing procedures
Authors:
Georg Hahn
Abstract:
Statistical discoveries are often obtained through multiple hypothesis testing. A variety of procedures exists to evaluate multiple hypotheses, for instance the ones of Benjamini-Hochberg, Bonferroni, Holm or Sidak. We are particularly interested in multiple testing procedures with two desired properties: (solely) monotonic and well-behaved procedures. This article investigates to which extent the…
▽ More
Statistical discoveries are often obtained through multiple hypothesis testing. A variety of procedures exists to evaluate multiple hypotheses, for instance the ones of Benjamini-Hochberg, Bonferroni, Holm or Sidak. We are particularly interested in multiple testing procedures with two desired properties: (solely) monotonic and well-behaved procedures. This article investigates to which extent the classes of (monotonic or well-behaved) multiple testing procedures, in particular the subclasses of so-called step-up and step-down procedures, are closed under basic set operations, specifically the union, intersection, difference and the complement of sets of rejected or non-rejected hypotheses. The present article proves two main results: First, taking the union or intersection of arbitrary (monotonic or well-behaved) multiple testing procedures results in new procedures which are monotonic but not well-behaved, whereas the complement or difference generally preserves neither property. Second, the two classes of (solely monotonic or well-behaved) step-up and step-down procedures are closed under taking the union or intersection, but not the complement or difference.
△ Less
Submitted 29 April, 2017; v1 submitted 14 July, 2016;
originally announced July 2016.
-
Bypass rewiring and robustness of complex networks
Authors:
Junsang Park,
Sang Geun Hahn
Abstract:
A concept of bypass rewiring is introduced and random bypass rewiring is analytically and numerically investigated with simulations. Our results show that bypass rewiring makes networks robust against removal of nodes including random failures and attacks. In particular, random bypass rewiring connects all nodes except the removed nodes on an even degree infinite network and makes the percolation…
▽ More
A concept of bypass rewiring is introduced and random bypass rewiring is analytically and numerically investigated with simulations. Our results show that bypass rewiring makes networks robust against removal of nodes including random failures and attacks. In particular, random bypass rewiring connects all nodes except the removed nodes on an even degree infinite network and makes the percolation threshold 0 for arbitrary occupation probabilities. In our example, the even degree network is more robust than the original network with random bypass rewiring while the original network is more robust than the even degree networks without random bypass. We propose a greedy bypass rewiring algorithm which guarantees the maximum size of the largest component at each step, assuming which node will be removed next is unknown. The simulation result shows that the greedy bypass rewiring algorithm improves the robustness of the autonomous system of the Internet under attacks more than random bypass rewiring.
△ Less
Submitted 26 September, 2017; v1 submitted 30 June, 2016;
originally announced June 2016.
-
Orbit Error Correction on the High Energy Beam Transport Line at the KHIMA Accelerator System
Authors:
Chawon Park,
Heejoong Yim,
Garam Hahn,
Dong Hyun An
Abstract:
For the purpose of treatment of various cancer and medical research, the synchrotron based medical machine under the Korea Heavy Ion Medical Accelerator (KHIMA) project have been conducted and is going to treat the patient at the beginning of 2018. The KHIMA synchrotron is designed to accelerate and extract the carbon ion (proton) beam with various energy range, 110 up to 430 MeV/u (60 up to 230 M…
▽ More
For the purpose of treatment of various cancer and medical research, the synchrotron based medical machine under the Korea Heavy Ion Medical Accelerator (KHIMA) project have been conducted and is going to treat the patient at the beginning of 2018. The KHIMA synchrotron is designed to accelerate and extract the carbon ion (proton) beam with various energy range, 110 up to 430 MeV/u (60 up to 230 MeV). A lattice design and beam optics studies for the High Energy Beam Transport (HEBT) line at the KHIMA accelerator system have been carried out with WinAgile and the MAD-X codes. Because the magnetic eld errors and the mis-alignments introduce to the deviations from the design parameters, these error sources should be treated explicitly and the sensitivity of the machine's lattice to di erent individual error sources is considered. Various types of errors which are static and dynamic one have been taken into account and have been consequentially corrected with a dedicated correction algorithm by using the MAD-X program. As a result, the tolerances for the diverse error contributions have been speci ed for the dedicated lattice components in the whole HEBT lines.
△ Less
Submitted 8 June, 2016;
originally announced June 2016.
-
Graph Partitioning Methods for Fast Parallel Quantum Molecular Dynamics
Authors:
Hristo N. Djidjev,
Georg Hahn,
Susan M. Mniszewski,
Christian F. A. Negre,
Anders M. N. Niklasson,
Vivek B. Sardeshmukh
Abstract:
We study a graph partitioning problem motivated by the simulation of the physical movement of multi-body systems on an atomistic level, where the forces are calculated from a quantum mechanical description of the electrons. Several advanced algorithms have been published in the literature for such simulations that are based on evaluations of matrix polynomials. We aim at efficiently parallelizing…
▽ More
We study a graph partitioning problem motivated by the simulation of the physical movement of multi-body systems on an atomistic level, where the forces are calculated from a quantum mechanical description of the electrons. Several advanced algorithms have been published in the literature for such simulations that are based on evaluations of matrix polynomials. We aim at efficiently parallelizing these computations by using a special type of graph partitioning. For this, we represent the zero-nonzero structure of a thresholded matrix as a graph and partition that graph into several components. The matrix polynomial is then evaluated for each separate submatrix corresponding to the subgraphs and the evaluated submatrix polynomials are used to assemble the final result for the full matrix polynomial. The paper provides a rigorous definition as well as a mathematical justification of this partitioning problem. We use several algorithms to compute graph partitions and experimentally evaluate their performance with respect to the quality of the partition obtained with each method and the time needed to produce it.
△ Less
Submitted 25 January, 2018; v1 submitted 3 May, 2016;
originally announced May 2016.
-
Cops and Robbers ordinals of cop-win trees
Authors:
Anthony Bonato,
Przemysław Gordinowicz,
Gena Hahn
Abstract:
A relational characterization of cop-win graphs was provided by Nowakowski and Winkler in their seminal paper on the game of Cops and Robbers. As a by-product of that characterization, each cop-win graph is assigned a unique ordinal, which we refer to as a CR-ordinal. For finite graphs, CR-ordinals correspond to the length of the game assuming optimal play, with the cop beginning the game in a lea…
▽ More
A relational characterization of cop-win graphs was provided by Nowakowski and Winkler in their seminal paper on the game of Cops and Robbers. As a by-product of that characterization, each cop-win graph is assigned a unique ordinal, which we refer to as a CR-ordinal. For finite graphs, CR-ordinals correspond to the length of the game assuming optimal play, with the cop beginning the game in a least favourable initial position. For infinite graphs, however, the possible values of CR-ordinals have not been considered in the literature until the present work. We classify the CR-ordinals of cop-win trees as either a finite ordinal, or those of the form $α+ ω$, where $α$ is a limit ordinal. For general infinite cop-win graphs, we provide an example whose CR-ordinal is not of this form. We finish with some problems on characterizing the CR-ordinals in the general case of cop-win graphs.
△ Less
Submitted 11 January, 2017; v1 submitted 14 March, 2016;
originally announced March 2016.
-
Optimal allocation of Monte Carlo simulations to multiple hypothesis tests
Authors:
Georg Hahn
Abstract:
Multiple hypothesis tests are often carried out in practice using p-value estimates obtained with bootstrap or permutation tests since the analytical p-values underlying all hypotheses are usually unknown. This article considers the allocation of a pre-specified total number of Monte Carlo simulations $K \in \mathbb{N}$ (i.e., permutations or draws from a bootstrap distribution) to a given number…
▽ More
Multiple hypothesis tests are often carried out in practice using p-value estimates obtained with bootstrap or permutation tests since the analytical p-values underlying all hypotheses are usually unknown. This article considers the allocation of a pre-specified total number of Monte Carlo simulations $K \in \mathbb{N}$ (i.e., permutations or draws from a bootstrap distribution) to a given number of $m \in \mathbb{N}$ hypotheses in order to approximate their p-values $p \in [0,1]^m$ in an optimal way, in the sense that the allocation minimises the total expected number of misclassified hypotheses. A misclassification occurs if a decision on a single hypothesis, obtained with an approximated p-value, differs from the one obtained if its p-value was known analytically. The contribution of this article is threefold: Under the assumption that $p$ is known and $K \in \mathbb{R}$, and using a normal approximation of the Binomial distribution, the optimal real-valued allocation of $K$ simulations to $m$ hypotheses is derived when correcting for multiplicity with the Bonferroni correction, both when computing the p-value estimates with or without a pseudo-count. Computational subtleties arising in the former case will be discussed. Second, with the help of an algorithm based on simulated annealing, empirical evidence is given that the optimal integer allocation is likely of the same form as the optimal real-valued allocation, and that both seem to coincide asympotically. Third, an empirical study on simulated and real data demonstrates that a recently proposed sampling algorithm based on Thompson sampling asympotically mimics the optimal (real-valued) allocation when the p-values are unknown and thus estimated at runtime.
△ Less
Submitted 5 October, 2019; v1 submitted 27 February, 2015;
originally announced February 2015.