Arxiv:2104.14310V3 [Quant-Ph] 2 Jun 2021

Total Page:16

File Type:pdf, Size:1020Kb

Arxiv:2104.14310V3 [Quant-Ph] 2 Jun 2021 Preparing Dicke states in a spin ensemble using phase estimation Yang Wang QuTech, Delft University of Technology, P.O. Box 5046, 2600 GA Delft, The Netherlands Barbara M. Terhal QuTech, Delft University of Technology, P.O. Box 5046, 2600 GA Delft, The Netherlands and JARA Institute for Quantum Information, Forschungszentrum Juelich, D-52425 Juelich, Germany (Dated: September 13, 2021) We present a Dicke state preparation scheme which uses global control of N spin qubits: our scheme is based on the standard phase estimation algorithm, which estimates the eigenvalue of a unitary operator. The scheme prepares a Dicke state non-deterministically by collectively coupling the spins to an ancilla qubit via a ZZ-interaction, using dlog2 Ne + 1 ancilla qubit measurements. The preparation of such Dicke states can be useful if the spins in the ensemble are used for magnetic sensing: we discuss a possible realization using an ensemble of electronic spins located at diamond nitrogen-vacancy centers coupled to a single superconducting flux qubit. We also analyze the effect of noise and limitations in our scheme. I. INTRODUCTION preparing Dicke states which are based on the full ad- dressability of the qubits [25, 26]. To address this issue, A promising application of the emerging quantum tech- some work has been dedicated to schemes which require nology is quantum-enhanced sensing, sometimes referred only a global control of the spin ensemble, such as using to as quantum metrology [1,2]. Using entangled states, steady state evolution [27], repeated energy transfer [7], continuous weak measurements [28], and the use of ge- one can, in principle, improve the measurementp sensi- tivity from the standard quantum limit (1= N) to the ometric phase gates [29]. Unfortunately, these methods Heisenberg limit (1=N)[1,3,4], where N is the number are still demanding currently when N is large, as they of- of probes or repetitions. However, preserving this quan- ten need complicated measurement-based feedback, high tum advantage is difficult in the presence of decoherence fidelity control and long preparation times. For example, [5]. For instance, a single-qubit Pauli Z-error can to- the optimized scheme in Ref. [7] uses O(N) rounds of ini- tally dephase a N-qubit GHZ-state, which would obtain tialization and evolution of an ancilla qubit. Our goal is Heisenberg-limited sensitivity in the noiseless case [6]. to improve the scaling with N so that one could possibly N-qubit Dicke states form a class of entangled states handle a larger error rate on the ancilla qubit. which are interesting for metrology [6{10]. Compared to In this paper, we present a Dicke state preparation other states used in quantum sensing, Dicke states have scheme that uses standard phase estimation [30], which been argued to be more robust to various noise sources prepares an eigenstate of a unitary operator by estimat- such as spin dephasing, spin damping, and spin number ing its eigenvalue. This algorithm is based on executing fluctuations [9]. Recent work has demonstrated a scheme projective measurements on the spin ensemble using an to use Dicke states for detecting the magnetic field in- ancilla qubit, and it will prepare a random Dicke state. duced by a single spin [7]. Another distinctive feature The scheme requires a ZZ-coupling between each spin of the use of Dicke states is that the optimal sensitivity in the ensemble to a single ancilla qubit. In Section III can be obtained through only global control on the set we detail how this coupling could be realized between an of spins [9, 10]. This is relevant for realizing practical ensemble of NV electronic spins and a superconducting quantum sensing using entangled states, as precise indi- flux qubit as ancilla. vidual spin qubit control can be difficult. Furthermore, Our scheme is efficient with respect to the number of arXiv:2104.14310v4 [quant-ph] 10 Sep 2021 superpositions of Dicke states can be used for quantum operations. It uses only dlog2 Ne + 1 rounds of phase error correction [11, 12]. estimation for preparing a random N-spin Dicke state. Dicke state preparation has been experimentally re- Each round of phase estimation measures a global opera- alized using photons [13, 14] and trapped-ion qubits tor of the spins, i.e., it applies an ancilla-qubit controlled 1 PN [15, 16], and there also exist many theoretical prepara- global Jz = 2 i=1 Zi rotation followed by ancilla qubit tion proposals suitable for a few qubits, see Ref. [17{ readout. The total time for performing the controlled 20] for example. However, it remains a challenge for rotations is upper-bounded by a constant and the prepa- large spin ensembles like N > O(100) diamond Nitrogen- ration timep thus scales as O(log2 N). With a probabil- Vacancy centers (negatively-charged NV [21]), each host- ity ∼ O(1= N), the prepared Dicke state would obtain ing an electronic spin S = 1. Since these NV center spins Heisenberg-limited sensitivity using only global control. are rather isolated from each other, it is costly to per- Besides the efficiency, our scheme also has some noise- form entangling gates between the electronic spins [22{ resilience: phase estimation can be realized with inte- 24]. This limitation excludes quantum algorithms for grated dynamical decoupling, which provides robustness 2 to the dephasing of the ancilla qubit as well as the de- The variance of θ can be calculated by the error propa- phasing of the spins in the ensemble. Furthermore, by gation formula repeating the projective measurements and performing a 2 simple majority vote, the effects of ancilla qubit decay 2 (∆M(θ)) (∆θ) = 2 ; (2) and flipped measurements (due to ancilla qubit dephas- j@θ hM(θ)ij ing or imperfect measurement) can be mitigated. This paper is organized as follows. In Sec.IA, we where the expectation value is with respect to the ini- briefly review Dicke states and Heisenberg-limited sens- tial state jN; mzi and M(θ) is the Heisenberg-evolved ing. In Sec. II we present the idea of using phase esti- operator. If we were to measure M = αJx + βJz, then mation to prepare Dicke states. In Sec. III we discuss hM(θ)i = hJzi (β cos(θ) + α sin(θ)) the Hamiltonian and a possible experimental setup with (3) multiple NV centers coupled to a flux qubit. In Sec.IV ≈ mz(β + θα); we numerically consider the performance of the scheme given the dominant noise sources. Finally, we discuss the for small θ (note that hN; mzj Jx jN; mzi = 0). We mea- results in Sec.V. sure Jx by choosing β = 0; α = 1, its expectation value hJx(θ)i has an optimal dependence on θ when mz is large. 2 However, the variance (∆Jx(θ)) will be large in a rotated A. Dicke states Dicke state, precluding any Heisenberg gains. 2 The proposal is instead to measure M = Jz , so that For simplicity, we assume even spin number N the variance is given by (see details in Ref. [10]): throughout this paper (odd spin-number N can be 2 2 2 2 2 treated similarly). The N-spin (or qubit) Dicke (∆θ) = [(∆Jx ) f(θ) + 4hJx i − 3hJy i state jN; m i with m 2 {− N ;:::; N g is a uniform, 2 2 2 2 2 2 z z 2 2 − 2hJz i × (1 + hJx i) + 6hJzJx Jzi] = [4(hJx i − hJz i) ] permutation-symmetric, superposition of N-bit strings (4) jxi where all bit-strings have N=2 + mz spins in j0i, 2 2 i.e., their Hamming weight is N=2 − mz. For example, (∆Jz ) 2 1 with f(θ) = 2 2 2 + tan (θ). The minimal vari- jN = 4; m = 0i = p (j0011i+j0101i+j0110i+j1001i+ (∆Jx) tan (θ) z 6 ance is obtained when tan2(θ) = p(∆J 2)2=(∆J 2)2. For j1010i + j1100i). A Dicke state jN; mzi is an eigenstate z x of the collective spin operator Dicke state jN; mzi the minimal variance (obtained at θ ≈ 0) is N 1 X Jz = Zi; (1) 2 4 2 2 2 2mz + 2 64mz − 16mz i=1 (∆θmin) = 2 2 + 2 2 2 : N + 2N − 12mz (N + 2N − 12mz) with eigenvalue mz. Here Zi is the Pauli Z operator on (5) 1 PN the spin labeled i. In addition, we have Jx = 2 i=1 Xi 1 PN Note that the sensitivity can surpass the standard quan- and Jy = 2 i=1 Yi. p To use such states for metrology, one imagines that tum limit when mz ∼ O( N) and is Heisenberg-limited 2 the prepared quantum state is transformed by e−iθJy when mz ∼ O(1). In addition, when mz = 0, (∆θmin) = 2 and the goal is to estimate the rotation angle θ which N(N+2) saturates the quantum Cram´er-Raobound [6]. 2 is assumed to be small. A standard metrological method The expectation value hJz i can in principle be obtained (for NV centers, limited by T2 and optical measurement by measuring Jz, squaring its outcome and gathering suf- accuracy) is Ramsey spectroscopy [31] using a single ficient statistics by repeating the measurements. We are qubit state repeatedly (or, equivalently, using a prod- thus especially interested in Dicke states close to jN; 0i, uct state of multiple qubits). In this context, the Ram- i.e., jN; mzi with mz ∼ O(1). Other than this motiva- sey method corresponds to preparing a simple prod- tion, we do not focus on aspects of using a (noisy) Dicke i π J uct state e 2 y j00 ::: 0i and letting it thus evolve to state for metrology in this paper.
Recommended publications
  • Towards the First Practical Applications of Quantum Computers
    Towards the First Practical Applications of Quantum Computers by Kevin J. Sung A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Computer Science and Engineering) in the University of Michigan 2020 Doctoral Committee: Professor Christopher Peikert, Co-Chair Adjunct Professor Yaoyun Shi, Co-Chair Dr. Ryan Babbush, Google Professor John P. Hayes Professor Emeritus Kim Winick Kevin J. Sung kevjsung@umich.edu ORCID iD: 0000-0001-6459-6374 ©2020 Kevin J. Sung To my family ii Acknowledgments This thesis was made possible by my supportive advisors Yaoyun Shi and Christopher Peikert, my friends, and my family. I would like to thank the Google AI Quantum team for hosting me while I performed much of the research in this thesis. I cannot express how lucky I feel to have had ac- cess to their state-of-the-art quantum computing hardware. The experiments on this hardware that I present in this thesis were the result of large team-wide collaborations. iii TABLE OF CONTENTS Dedication ....................................... ii Acknowledgments ................................... iii List of Figures ..................................... vii List of Tables ...................................... xiii List of Algorithms ................................... xv List of Appendices ................................... xvi Abstract ......................................... xvii Chapter 1 Introduction ..................................... 1 1.1 Quantum computing in the NISQ era....................1 1.1.1 The
    [Show full text]
  • Arxiv:1812.09167V1 [Quant-Ph] 21 Dec 2018 It with the Tex Typesetting System Being a Prime Example
    Open source software in quantum computing Mark Fingerhutha,1, 2 Tomáš Babej,1 and Peter Wittek3, 4, 5, 6 1ProteinQure Inc., Toronto, Canada 2University of KwaZulu-Natal, Durban, South Africa 3Rotman School of Management, University of Toronto, Toronto, Canada 4Creative Destruction Lab, Toronto, Canada 5Vector Institute for Artificial Intelligence, Toronto, Canada 6Perimeter Institute for Theoretical Physics, Waterloo, Canada Open source software is becoming crucial in the design and testing of quantum algorithms. Many of the tools are backed by major commercial vendors with the goal to make it easier to develop quantum software: this mirrors how well-funded open machine learning frameworks enabled the development of complex models and their execution on equally complex hardware. We review a wide range of open source software for quantum computing, covering all stages of the quantum toolchain from quantum hardware interfaces through quantum compilers to implementations of quantum algorithms, as well as all quantum computing paradigms, including quantum annealing, and discrete and continuous-variable gate-model quantum computing. The evaluation of each project covers characteristics such as documentation, licence, the choice of programming language, compliance with norms of software engineering, and the culture of the project. We find that while the diversity of projects is mesmerizing, only a few attract external developers and even many commercially backed frameworks have shortcomings in software engineering. Based on these observations, we highlight the best practices that could foster a more active community around quantum computing software that welcomes newcomers to the field, but also ensures high-quality, well-documented code. INTRODUCTION Source code has been developed and shared among enthusiasts since the early 1950s.
    [Show full text]
  • Arxiv:1302.3247V5 [Quant-Ph] 10 Dec 2014 Pairwise Interaction
    Permutation-invariant quantum codes Yingkai Ouyang University of Waterloo, Waterloo, Ontario, Canada and Singapore University of Technology and Design, Singapore∗ A quantum code is a subspace of a Hilbert space of a physical system chosen to be correctable against a given class of errors, where information can be encoded. Ideally, the quantum code lies within the ground space of the physical system. When the physical model is the Heisenberg ferromagnet in the absence of an external magnetic field, the corresponding ground-space contains all permutation-invariant states. We use techniques from combinatorics and operator theory to construct families of permutation-invariant quantum codes. These codes have length proportional to t2; one family of codes perfectly corrects arbitrary weight t errors, while the other family of codes approximately correct t spontaneous decay errors. The analysis of our codes' performance with respect to spontaneous decay errors utilizes elementary matrix analysis, where we revisit and extend the quantum error correction criterion of Knill and Laflamme, and Leung, Chuang, Nielsen and Yamamoto. PACS numbers: 03.65.Aa,03.67.Pp,05.30.-d,75.10.Pq I. INTRODUCTION label the particles in the system, and e = fi; jg labels the exchange interactions in the system. Here Je and A quantum bit (qubit) is a fundamental resource Si denote the exchange constants and the vector of required in many quantum information theoretic tasks, spin operators respectively. Since exchange operators such as in quantum cryptographic protocols [1] and in essentially swap particles (see ex. 1.9 of Ref. [12] quantum computers [2]. To combat decoherence, a two- or Ref.
    [Show full text]
  • Low-Density Parity-Check Code-Based Cryptographic Systems Specification Revision 3.0 – April, 2020 Name of the Proposed Cryptosystem
    LEDAcrypt: Low-dEnsity parity-check coDe-bAsed cryptographic systems Specification revision 3:0 { April, 2020 Name of the proposed cryptosystem LEDAcrypt (Low-dEnsity parity-check coDe-bAsed cryptographic systems) Submitters This submission is from the following team, listed in alphabetical order: • Marco Baldi, Universit`aPolitecnica delle Marche, Ancona, Italy • Alessandro Barenghi, Politecnico di Milano, Milano, Italy • Franco Chiaraluce, Universit`aPolitecnica delle Marche, Ancona, Italy • Gerardo Pelosi, Politecnico di Milano, Milano, Italy • Paolo Santini, Universit`aPolitecnica delle Marche, Ancona, Italy E-mail addresses: m.baldi@univpm.it, alessandro.barenghi@polimi.it, f.chiaraluce@univpm.it, gerardo.pelosi@polimi.it, p.santini@pm.univpm.it. Contact telephone and address Marco Baldi (phone: +39 071 220 4894), Universit`aPolitecnica delle Marche, Dipartimento di Ingegneria dell'Informazione (DII), Via Brecce Bianche 12, I-60131, Ancona, Italy. Names of auxiliary submitters There are no auxiliary submitters. The principal submitter is the team listed above. Name of the inventors/developers of the cryptosystem Same as submitters. Name of the owner, if any, of the cryptosystem Same as submitters. Backup contact telephone and address Gerardo Pelosi (phone: +39 02 2399 3476), Politecnico di Milano, Dipartimento di Elettronica, Informazione e Bioingegneria (DEIB), Via G. Ponzio 34/5, I-20133, Milano, Italy. Signature of the submitter1 × 1See also printed version of \Statement by Each Submitter". LEDAcrypt Page 3 Contents Foreword 9 1 Complete written specification 11 1.1 Preliminaries . 11 1.1.1 Finite fields and circulant matrix algebra . 11 1.1.2 Quasi-cyclic low-density parity-check codes and their efficient decoding . 15 1.1.3 Classic code-based cryptosystems .
    [Show full text]
  • Arxiv:2012.11382V2 [Quant-Ph] 11 Jan 2021
    Quantum Integer Programming (QuIP) 47-779: Lecture Notes David E. Bernal, Sridhar Tayur, Davide Venturelli Fall 2020 Abstract This lecture series on Quantum Integer Programming (QuIP) { created by Professor Sridhar Tayur, David E. Bernal and Dr. Davide Venturelli, a collaboration between CMU and USRA, with the support from Amazon Braket during Fall 2020 { is intended for students and researchers interested in Integer Programming and the potential of near term quantum and quantum inspired computing in solving optimization problems. Originally created for Tepper School of Business course 47-779 (at CMU), these were also used for the course ID5840 (at IIT-Madras, by Professors Anil Prabhakar and Prabha Mandayam) whose students (listed at the beginning of each lecture) were scribes. Dr. Vikesh Siddhu, post-doc in CMU Quantum Computing Group, assisted during the lectures, student projects and with proof-reading this scribe. Through these lectures one will learn to formulate a problem and map it to a Quadratic Unconstrained Binary Optimization (QUBO) problem, understand various mapping and techniques like the Ising model, Graver Augmented Multiseed Algorithm (GAMA), Simulated or Quantum Annealing and QAOA, and ideas on how to solve these Integer problems using these quantum and classical methods. The course website (with lecture videos and colab notebooks): https://bernalde.github.io/QuIP/ Keywords: Ising model, Integer Programming, Computational Algebraic Geometry, Graver Basis, Quan- tum Annealing, Simulated Annealing, Combinatorial Optimization, Graph coloring, discrete nonlinear opti- mization. arXiv:2012.11382v2 [quant-ph] 11 Jan 2021 0-1 47-779. Quantum Integer Programming Mini-1, Fall 2020 Room: Zoom Online Time: Tuesday and Thursday 5:20pm-7:10pm Instructors: Sridhar Tayur Email: stayur@cmu.edu Office: 4216 Tepper Quad Davide Venturelli Email: DVenturelli@usra.edu Office: Online David E.
    [Show full text]
  • Solving Constrained Quadratic Binary Problems Via Quantum Adiabatic Evolution
    Solving Constrained Quadratic Binary Problems via Quantum Adiabatic Evolution Pooya Ronagh, Brad Woods, Ehsan Iranmanesh Abstract Quantum adiabatic evolution is perceived as useful for binary quadratic programming problems that are a priori unconstrained. For constrained problems, it is a common practice to relax linear equality constraints as penalty terms in the objective function. However, there has not yet been proposed a method for efficiently dealing with inequality constraints using the quantum adiabatic approach. In this paper, we give a method for solving the Lagrangian dual of a binary quadratic programming (BQP) problem in the presence of inequality constraints and employ this procedure within a branch-and-bound framework for constrained BQP (CBQP) problems. Keywords: Adiabatic quantum computation, Constrained integer programming, Branch and bound, Lagrangian duality 1 Introduction The unconstrained binary quadratic programming (UBQP) problem is defined by Minimize xT Qx subject to x 2 f0; 1gn; where, without loss of generality, Q 2 Zn×n. Recent advancements in quantum computing technology [6, 20, 22] have raised hopes of the production of computing systems that are capable of solving UBQP problems, showing quantum speedup. The stochastic nature of such systems, together with all sources of noise and error, are challenges yet to be overcome in achieving scalable quantum computing systems of this type. This paper is regardless motivated by the assumption of the existence of systems that can solve UBQP problems efficiently and to optimality, or at leastin conjunction with a framework of noise analysis of the suboptimal results. We call such a computing system a UBQP oracle. Many NP-hard combinatorial optimization problems arise naturally or can easily be reformulated as UBQP problems, such as the quadratic assignment problem, the maximum cut problem, the maximum clique problem, the set packing problem, and the graph colouring problem (see, for instance, Boros and Prékopa [9], Boros and Hammer [8], Bourjolly et al.
    [Show full text]
  • Master Thesis & Internship Projects @ Imec
    Master Thesis & Internship Projects @ imec 2018 Topic Guide 1 Information.............................................................................................................................................................................................1 I. CMOS & beyond CMOS ................................................................................................... 7 Machine learning based computational lithography........................................................................................................... 7 Ferroelectric device electrical characterization ................................................................................................................. 7 Selective low-k dielectric growth using metal nitride conversion approach .............................................................. 8 Impact of thermal budget on spin-orbit torque and SOT-MRAM performances ..................................................... 9 Low temperature epitaxial growth of Ga doped SiGe for future device applications ........................................... 10 Study of an optical metrology technique for advanced semiconductor process control...................................... 10 Spin wave-like excitations in low-dimensional ferromagnets ....................................................................................... 11 Time-dependent Ginzburg-Landau simulations for superconducting rings ............................................................... 12 Modeling of Transmon Qubits.............................................................................................................................................
    [Show full text]
  • Enfield: an Openqasm Compiler
    Enfield: An OpenQASM Compiler Marcos Yukio Siraichi1, Caio Henrique Segawa Tonetti1 1Departamento de Cienciaˆ da Computac¸ao˜ – UFMG Av. Antonioˆ Carlos, 6627 – 31.270-010 – Belo Horizonte – MG – Brazil fyukio.siraichi,caiotonettig@dcc.ufmg.br Abstract. In 2016, IBM made quantum computing accessible to everyone by deploying a cloud-based quantum computer. It elicited much enthusiasm for quantum computing research and testing. When writing programs for execut- ing in this environment, the programmer had to manually tune the program so that it complied to some resource restrictions specific to IBM’s machine. Hence, we generalized and automated this process, building Enfield: a source-to-source compiler for transforming quantum programs, allocating resources and optimiz- ing. It features state-of-the-art implementations and the extensibility necessary for implementing and comparing new solutions. We illustrate its effectiveness with benchmarks from IBM that were used in previous research and show that our tool performs at least 22% better. Video Link: https://youtu.be/u5OgQAlB3L0 1. Introduction With the recent popularization of quantum computing and the introduction of experi- mental prototypes available for the general public, the possibility of building and pro- gramming quantum experiments has elicited much enthusiasm [Devitt 2016]. In this plat- form, tools that aid the creation of programs are available, such as visual circuit repre- sentation and a gate-level syntax based on the Open Quantum Assembly (OpenQASM) [Cross et al. 2017]. Visual circuits are good for learning the basics of quantum program- ming, but it quickly becomes a tiresome activity for larger programs. OpenQASM, on the other hand, offers a low level of abstraction, demanding a deeper understanding of each gate and the specific constraints for the machine to be executed on [Haner¨ et al.
    [Show full text]
  • Verifying Results of the IBM Qiskit Quantum Circuit Compilation Flow
    Verifying Results of the IBM Qiskit Quantum Circuit Compilation Flow Lukas Burgholzer∗ Rudy Raymondy Robert Wille∗z ∗Institute for Integrated Circuits, Johannes Kepler University Linz, Austria yIBM Research – Tokyo, Japan zSoftware Competence Center Hagenberg GmbH (SCCH), Hagenberg, Austria lukas.burgholzer@jku.at rudyhar@jp.ibm.com robert.wille@jku.at https://iic.jku.at/eda/research/quantum/ Abstract—Realizing a conceptual quantum algorithm on an two quantum circuits G and G0 boils down to comparing actual physical device necessitates the algorithm’s quantum the corresponding unitary matrices U and U 0, respectively. circuit description to undergo certain transformations in order to However, since these matrices U and U 0 are exponentially adhere to all constraints imposed by the hardware. In this regard, the individual high-level circuit components are first synthesized large, many existing approaches to this problem (e.g., [5]–[9]) to the supported low-level gate-set of the quantum computer, remain unsatisfactory and can hardly be employed on a large before being mapped to the target’s architecture—utilizing sev- scale. eral optimizations in order to improve the compilation result. Recently, a promising solution to address this problem Specialized tools for this complex task exist, e.g., IBM’s Qiskit, Google’s Cirq, Microsoft’s QDK, or Rigetti’s Forest. However, to has been proposed in [10]. There, the use of dedicated date, the circuits resulting from these tools are hardly verified, data-structures, in particular decision diagrams [9], [11], [12], which is mainly due to the immense complexity of checking if has been proposed which allow for a non-exponential rep- two quantum circuits indeed realize the same functionality.
    [Show full text]
  • Computer Physics Communications Massively Parallel Quantum Computer Simulator, Eleven Years Later
    Computer Physics Communications 237 (2019) 47–61 Contents lists available at ScienceDirect Computer Physics Communications journal homepage: www.elsevier.com/locate/cpc Massively parallel quantum computer simulator, eleven years later Hans De Raedt a, Fengping Jin b, Dennis Willsch b,c, Madita Willsch b,c, Naoki Yoshioka d, ∗ ∗∗ Nobuyasu Ito d,e, Shengjun Yuan f, , Kristel Michielsen b,c, a Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, NL-9747 AG Groningen, The Netherlands b Institute for Advanced Simulation, Jülich Supercomputing Center, Forschungzentrum Jülich, D-52425 Jülich, Germany c RWTH Aachen University, D-52056 Aachen, Germany d RIKEN Center for Computational Science, 7-1-26 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan e Department of Applied Physics, School of Engineering, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan f School of Physics and Technology, Wuhan University, Wuhan 430072, China article info a b s t r a c t Article history: A revised version of the massively parallel simulator of a universal quantum computer, described in this Received 12 May 2018 journal eleven years ago, is used to benchmark various gate-based quantum algorithms on some of the Received in revised form 30 October 2018 most powerful supercomputers that exist today. Adaptive encoding of the wave function reduces the Accepted 8 November 2018 memory requirement by a factor of eight, making it possible to simulate universal quantum computers Available online 29 November 2018 with up to 48 qubits on the Sunway TaihuLight and on the K computer. The simulator exhibits close-to- Keywords: ideal weak-scaling behavior on the Sunway TaihuLight, on the K computer, on an IBM Blue Gene/Q, and Quantum computation on Intel Xeon based clusters, implying that the combination of parallelization and hardware can track the Computer simulation exponential scaling due to the increasing number of qubits.
    [Show full text]
  • Massively Parallel Quantum Computer Simulator, Eleven Years Later
    Computer Physics Communications 237 (2019) 47–61 Contents lists available at ScienceDirect Computer Physics Communications journal homepage: www.elsevier.com/locate/cpc Massively parallel quantum computer simulator, eleven years later Hans De Raedt a, Fengping Jin b, Dennis Willsch b,c, Madita Willsch b,c, Naoki Yoshioka d, ∗ ∗∗ Nobuyasu Ito d,e, Shengjun Yuan f, , Kristel Michielsen b,c, a Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, NL-9747 AG Groningen, The Netherlands b Institute for Advanced Simulation, Jülich Supercomputing Center, Forschungzentrum Jülich, D-52425 Jülich, Germany c RWTH Aachen University, D-52056 Aachen, Germany d RIKEN Center for Computational Science, 7-1-26 Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan e Department of Applied Physics, School of Engineering, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan f School of Physics and Technology, Wuhan University, Wuhan 430072, China article info a b s t r a c t Article history: A revised version of the massively parallel simulator of a universal quantum computer, described in this Received 12 May 2018 journal eleven years ago, is used to benchmark various gate-based quantum algorithms on some of the Received in revised form 30 October 2018 most powerful supercomputers that exist today. Adaptive encoding of the wave function reduces the Accepted 8 November 2018 memory requirement by a factor of eight, making it possible to simulate universal quantum computers Available online 29 November 2018 with up to 48 qubits on the Sunway TaihuLight and on the K computer. The simulator exhibits close-to- Keywords: ideal weak-scaling behavior on the Sunway TaihuLight, on the K computer, on an IBM Blue Gene/Q, and Quantum computation on Intel Xeon based clusters, implying that the combination of parallelization and hardware can track the Computer simulation exponential scaling due to the increasing number of qubits.
    [Show full text]
  • Free Software for Quantum Computation
    Free software for quantum computation Piotr Gawron Institute of Theoretical and Applied Informatics Polish Academy of Sciences KQIS AGH 13 March 2019 Kraków Outline Introduction—Free software Free software actions for quantum computation Quantum Open Source Foundation Fosdem 19 Quantum computing track Xanadu.ai Julia and quantum computing Introduction to Julia QuantumInformation.jl Programming D-Wave Annealer Quantum annealing D-Wave annealer D-Wave software stack Future work—interesting challenges Interesting goals to pursue Outline Introduction—Free software Free software actions for quantum computation Quantum Open Source Foundation Fosdem 19 Quantum computing track Xanadu.ai Julia and quantum computing Introduction to Julia QuantumInformation.jl Programming D-Wave Annealer Quantum annealing D-Wave annealer D-Wave software stack Future work—interesting challenges Interesting goals to pursue Free software IT Giants conference AGH 2009 Free software A program is free software if the program’s users have the four essential freedoms: I The freedom to run the program as you wish, for any purpose (freedom 0). I The freedom to study how the program works, and change it so it does your computing as you wish (freedom 1). Access to the source code is a precondition for this. I The freedom to redistribute copies so you can help others (freedom 2). I The freedom to distribute copies of your modified versions to others (freedom 3). By doing this you can give the whole community a chance to benefit from your changes. Access to the source code is a
    [Show full text]