ACCELERATED PATTERN MATCHING METHOD ON A QUANTUM COMPUTING SYSTEM

Information

  • Patent Application
  • 20240013078
  • Publication Number
    20240013078
  • Date Filed
    September 18, 2023
    8 months ago
  • Date Published
    January 11, 2024
    4 months ago
Abstract
A method of determining a pattern in a sequence of bits using a quantum computing system includes setting a first register of a quantum processor in a superposition of a plurality of string index states, encoding a bit string in a second register of the quantum processor, encoding a bit pattern in a third register of the quantum processor, circularly shifting qubits of the second register conditioned on the first register, amplifying an amplitude of a state combined with the first register in which the circularly shifted qubits of the second register matches qubits of the third register, measuring an amplitude of the first register and determining a string index state of the plurality of string index states associated with the amplified state, and outputting, by use of a classical computer, a string index associated with the first register in the measured state.
Description
BACKGROUND
Field

The present disclosure generally relates to a method of performing computation in a quantum computing system, and more specifically, to a method of searching for a bit pattern in a bit string using a computing system that are a group of trapped ions.


Description of the Related Art

In quantum computing, quantum bits or qubits, which are analogous to bits representing a “0” and a “1” in a classical (digital) computer, are required to be prepared, manipulated, and measured (read-out) with near perfect control during a computation process.


Among physical systems upon which it is proposed to build large-scale quantum computers, is a group of ions (e.g., charged atoms), which are trapped and suspended in vacuum by electromagnetic fields. The ions have internal hyperfine states which are separated by frequencies in the several GHz range and can be used as the computational states of a qubit (referred to as “qubit states”). These hyperfine states can be controlled using radiation provided from a laser, or sometimes referred to herein as the interaction with laser beams. The ions can be cooled to near their motional ground states using such laser interactions. The ions can also be optically pumped to one of the two hyperfine states with high accuracy (preparation of qubits), manipulated between the two hyperfine states (single-qubit gate operations) by laser beams, and their internal hyperfine states detected by fluorescence upon application of a resonant laser beam (read-out of qubits). A pair of ions can be controllably entangled (two-qubit gate operations) by qubit-state dependent force using laser pulses that couple the ions to the collective motional modes of a group of trapped ions, which arise from their Coulombic interaction between the ions. In general, entanglement occurs when pairs or groups of ions (or particles) are generated, interact, or share spatial proximity in ways such that the quantum state of each ion cannot be described independently of the quantum state of the others, even when the ions are separated by a large distance.


Pattern matching is ubiquitously used in various areas, such as image processing, study of DNA sequences, and data compression and statistics. Computational complexity of a pattern matching process (i.e., the amount of resources required to run the pattern matching process, in particular, time and memory requirements) varies as the size of input varies. For example, in a string matching process that is the simplest pattern matching process, in which a bit pattern custom-character of length M is searched within a bit string custom-character of length N (M≤N), the amount of time resources scales as Θ(N+M) in units of elementary logic operations (which takes a constant amount of time on a given computer) and the amount of memory resources scales as custom-character(N+M) in units of bits, by the best known classical methods.


Therefore, there is a need for method for using a quantum processor to speed up pattern matching in a quantum computing system.


SUMMARY

Embodiments described herein provide a method of determining a pattern in a sequence of bits using a quantum computing system. The method includes setting a first register of a quantum processor in a superposition of a plurality of string index states, encoding a bit string in a second register of the quantum processor, encoding a bit pattern in a third register of the quantum processor, circularly shifting qubits of the second register conditioned on the first register, amplifying an amplitude of a state combined with the first register in which the circularly shifted qubits of the second register matches qubits of the third register, measuring an amplitude of the first register and determining a string index state of the plurality of string index states associated with the amplified state, and outputting, by use of a classical computer, a string index associated with the first register in the measured state.


Embodiments described herein also provides a quantum computing system for determining a pattern in a sequence of bits. The quantum computing system includes a quantum processor including a group of trapped ions, each of the trapped ions having two frequency-separated states, one or more lasers configured to emit a laser beam, which is provided to trapped ions in the quantum processor, a classical computer, and a system controller. The classical computer is configured to select a bit string and a bit pattern to be searched within the bit string, each bit of the bit string having a string index. The system controller is configured to set a first register of the quantum processor in a superposition of a plurality of string index states, each of which is associated with a string index, encode the bit string in a second register of the quantum processor, encode the bit pattern in a third register of the quantum processor, circularly shift qubits of the second register conditioned on the first register in each string index state in the superposition of the plurality of string index states, amplify an amplitude of a state combined with the first register in a string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of the third register, and measure an amplitude of the first register and determining the string index state of the plurality of string index states associated with the amplified state. The classical computer is further configured to output the string index associated with the first register in the measured state.


Embodiments described herein further provide a quantum computing system including non-volatile memory having a number of instructions stored therein. The instruction, when executed by one or more processors, cause the quantum computing system to perform operations including setting a first register of the quantum processor in a superposition of a plurality of string index states, each of which is associated with a string index, where the quantum processor includes a plurality of qubits, encoding a bit string in a second register of the quantum processor, encoding a bit pattern in a third register of the quantum processor, circularly shifting qubits of the second register conditioned on the first register in each string index state in the superposition of the plurality of string index states, amplifying an amplitude of a state combined with the first register in a string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of the third register, measuring an amplitude of the first register and determining the string index state of the plurality of string index states associated with the amplified state, and outputting, by use of a classical computer, the string index associated with the first register in the measured state. The quantum computing system further includes selecting, by use of the classical computer, the bit string and the bit pattern to be searched within the bit string, each bit of the bit string having a string index, where each qubit of the plurality of qubits includes a trapped ion having two frequency-separated states. The quantum computing system further includes preparing the quantum processor in an initial state by setting, by a system controller, each trapped ion in the quantum processor in the lower energy state of the two frequency-separated states. Setting the first register of the quantum processor in the superposition of the plurality of string index states includes transferring, by the system controller, each trapped ion in the first register in a superposition of the two frequency-separated states. Encoding the bit string in the second register of the quantum processor includes applying, by the system controller, a combination of single-qubit gate operations to the second register of the quantum processor. Encoding the bit pattern in the third register of the quantum processor includes applying, by the system controller, a combination of single-qubit gate operations to the third register of the quantum processor. Circularly shifting qubits of the second register conditioned on the first register includes applying, by the system controller, a combination of single-qubit gate operations and two-qubit gate operations to the first and second registers of the quantum processor. The two-qubit gate operations includes controlled-SWAP operations applied to the second register conditioned on the first register. The two-qubit gate operations includes controlled-SWAP operations applied to the second register conditioned on a plurality of ancillary qubits in which the first register is copied by fan-out CNOT operations applied on the plurality of ancillary qubits conditioned on the first register. Amplifying the state combined with the first register in the string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of the third register includes applying, by the system controller, a combination of single-qubit gate operations and two-qubit gate operations to the second and third registers of the quantum processor.





BRIEF DESCRIPTION OF THE DRAWINGS

So that the manner in which the above-recited features of the present disclosure can be understood in detail, a more particular description of the disclosure, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this disclosure and are therefore not to be considered limiting of its scope, for the disclosure may admit to other equally effective embodiments.



FIG. 1 is a schematic partial view of an ion trap quantum computing system according to one embodiment.



FIG. 2 depicts a schematic view of an ion trap for confining ions in a group according to one embodiment.



FIG. 3 depicts a schematic energy diagram of each ion in a group of trapped ions according to one embodiment.



FIG. 4 depicts a qubit state of an ion represented as a point on a surface of the Bloch sphere.



FIGS. 5A, 5B, and 5C depict a few schematic collective transverse motional mode structures of a group of five trapped ions.



FIGS. 6A and 6B depict schematic views of motional sideband spectrum of each ion and a motional mode according to one embodiment.



FIG. 7 depicts a flowchart illustrating a method of searching for a bit pattern within a bit string according to one embodiment.





To facilitate understanding, identical reference numerals have been used, where possible, to designate identical elements that are common to the figures. In the figures and the following description, an orthogonal coordinate system including an X-axis, a Y-axis, and a Z-axis is used. The directions represented by the arrows in the drawing are assumed to be positive directions for convenience. It is contemplated that elements disclosed in some embodiments may be beneficially utilized on other implementations without specific recitation.


DETAILED DESCRIPTION

Embodiments described herein are generally related to a method of performing computation in a quantum computing system, and more specifically, to a method of searching for a bit pattern in a bit string using a computing system that are a group of trapped ions.


A quantum computing system that is able to solve a pattern matching problem may include a classical computer, a system controller, and a quantum processor. The classical computer performs supporting and system control tasks including selecting a bit string and a bit pattern to be searched within the bit string by use of a user interface. The classical computer further converts a series of logic gates into laser pulses to be applied on the quantum processor by the system processor. A software program for performing the tasks is stored in a non-volatile memory within the classical computer. The quantum processor includes trapped ions that are coupled by use of various hardware, including lasers to manipulate internal hyperfine states (qubit states) of the trapped ions and photomultiplier tubes (PMTs) to read-out the internal hyperfine states (qubit states) of the trapped ions. The system controller receives from the classical computer instructions for controlling the quantum processor, and controls various hardware associated with controlling any and all aspects of the process of controlling the quantum processor, and returns a read-out of the quantum processor and thus output of results of the computation(s) to the classical computer.


The methods and systems described herein include an efficient pattern matching routine that provides improvements over conventional pattern matching routines.


General Hardware Configurations


FIG. 1 is a schematic partial view of an ion trap quantum computing system, or system 100, according to one embodiment. The system 100 includes a classical (digital) computer 102, a system controller 104 and a quantum processor that is a group 106 of trapped ions (i.e., five shown) that extend along the Z-axis. Each ion in the group 106 of trapped ions is an ion having a nuclear spin I and an electron spin S such that a difference between the nuclear spin I and the electron spin S is zero, such as a positive ytterbium ion, 171Yb+, a positive barium ion 133Ba+, a positive cadmium ion 111cd+ or 113Cd+, which all have a nuclear spin I=½ and the 2S1/2 hyperfine states. In some embodiments, all ions in the group 106 of trapped ions are the same species and isotope (e.g., 171Yb+.) In some other embodiments, the group 106 of trapped ions includes one or more species or isotopes (e.g., some ions are 171Yb+ and some other ions are 133Ba+). In yet additional embodiments, the group 106 of trapped ions may include various isotopes of the same species (e.g., different isotopes of Yb, different isotopes of Ba). The ions in the group 106 of trapped ions are individually addressed with separate laser beams. The classical computer 102 includes a central processing unit (CPU), memory, and support circuits (or I/O). The memory is connected to the CPU, and may be one or more of a readily available memory, such as a read-only memory (ROM), a random access memory (RAM), floppy disk, hard disk, or any other form of digital storage, local or remote. Software instructions, algorithms and data can be coded and stored within the memory for instructing the CPU. The support circuits (not shown) are also connected to the CPU for supporting the processor in a conventional manner. The support circuits may include conventional cache, power supplies, clock circuits, input/output circuitry, subsystems, and the like.


An imaging objective 108, such as an objective lens with a numerical aperture (NA), for example, of 0.37, collects fluorescence along the Y-axis from the ions and maps each ion onto a multi-channel photo-multiplier tube (PMT) 110 for measurement of individual ions. Non-copropagating Raman laser beams from a laser 112, which are provided along the X-axis, perform operations on the ions. A diffractive beam splitter 114 creates an array of static Raman beams 116 that are individually switched using a multi-channel acousto-optic modulator (AOM) 118 and is configured to selectively act on individual ions. A global Raman laser beam 120 illuminates all ions at once. In some embodiments, individual Raman laser beams (not shown) each illuminate individual ions. The system controller (also referred to as a “RF controller”) 104 controls the AOM 118 and thus controls laser pulses to be applied to trapped ions in the group 106 of trapped ions. The system controller 104 includes a central processing unit (CPU) 122, a read-only memory (ROM) 124, a random access memory (RAM) 126, a storage unit 128, and the like. The CPU 122 is a processor of the system controller 104. The ROM 124 stores various programs and the RAM 126 is the working memory for various programs and data. The storage unit 128 includes a nonvolatile memory, such as a hard disk drive (HDD) or a flash memory, and stores various programs even if power is turned off. The CPU 122, the ROM 124, the RAM 126, and the storage unit 128 are interconnected via a bus 130. The system controller 104 executes a control program which is stored in the ROM 124 or the storage unit 128 and uses the RAM 126 as a working area. The control program will include software applications that include program code that may be executed by processor in order to perform various functionalities associated with receiving and analyzing data and controlling any and all aspects of the methods and hardware used to create the ion trap quantum computer system 100 discussed herein.



FIG. 2 depicts a schematic view of an ion trap 200 (also referred to as a Paul trap) for confining ions in the group 106 according to one embodiment. The confining potential is exerted by both static (DC) voltage and radio frequency (RF) voltages. A static (DC) voltage V s is applied to end-cap electrodes 210 and 212 to confine the ions along the Z-axis (also referred to as an “axial direction” or a “longitudinal direction”). The ions in the group 106 are nearly evenly distributed in the axial direction due to the Coulomb interaction between the ions. In some embodiments, the ion trap 200 includes four hyperbolically-shaped electrodes 202, 204, 206, and 208 extending along the Z-axis.


During operation, a sinusoidal voltage V1 (with an amplitude VRF/2) is applied to an opposing pair of the electrodes 202, 204 and a sinusoidal voltage V2 with a phase shift of 180° from the sinusoidal voltage V1 (and the amplitude VRF/2) is applied to the other opposing pair of the electrodes 206, 208 at a driving frequency ωRF, generating a quadrupole potential. In some embodiments, a sinusoidal voltage is only applied to one opposing pair of the electrodes 202, 204, and the other opposing pair 206, 208 is grounded. The quadrupole potential creates an effective confining force in the X-Y plane perpendicular to the Z-axis (also referred to as a “radial direction” or “transverse direction”) for each of the trapped ions, which is proportional to a distance from a saddle point (i.e., a position in the axial direction (Z-direction)) at which the RF electric field vanishes. The motion in the radial direction (i.e., direction in the X-Y plane) of each ion is approximated as a harmonic oscillation (referred to as secular motion) with a restoring force towards the saddle point in the radial direction and can be modeled by spring constants kx and ky, respectively, as is discussed in greater detail below. In some embodiments, the spring constants in the radial direction are modeled as equal when the quadrupole potential is symmetric in the radial direction. However, undesirably in some cases, the motion of the ions in the radial direction may be distorted due to some asymmetry in the physical trap configuration, a small DC patch potential due to inhomogeneity of a surface of the electrodes, or the like and due to these and other external sources of distortion the ions may lie off-center from the saddle points.



FIG. 3 depicts a schematic energy diagram 300 of each ion in the group 106 of trapped ions according to one embodiment. Each ion in the group 106 of trapped ions is an ion having a nuclear spin I and an electron spin S such that a difference between the nuclear spin I and the electron spin S is zero. In one example, each ion may be a positive Ytterbium ion, 171Yb+, which has a nuclear spin I=½ and the 2S1/2 hyperfine states (i.e., two electronic states) with an energy split corresponding to a frequency difference (referred to as a “carrier frequency”) of ω01/2π=12.642821 GHz. In other examples, each ion may be a positive barium ion 133Ba+, a positive cadmium ion ca or 113Cd+, which all have a nuclear spin I=½ and the 2S1/2 hyperfine states. A qubit is formed with the two hyperfine states, denoted as |0custom-character and |1custom-character, where the hyperfine ground state (i.e., the lower energy state of the 2S1/2 hyperfine states) is chosen to represent |0custom-character. Hereinafter, the terms “hyperfine states,” “internal hyperfine states,” and “qubits” may be interchangeably used to represent |0custom-character and |1custom-character. Each ion may be cooled (i.e., kinetic energy of the ion may be reduced) to near the motional ground state |0custom-characterm for any motional mode m with no phonon excitation (i.e., nph=0) by known laser cooling methods, such as Doppler cooling or resolved sideband cooling, and then the qubit state prepared in the hyperfine ground state |0custom-character by optical pumping. Here, |0custom-character represents the individual qubit state of a trapped ion whereas |0custom-characterm with the subscript m denotes the motional ground state for a motional mode m of a group 106 of trapped ions.


An individual qubit state of each trapped ion may be manipulated by, for example, a mode-locked laser at 355 nanometers (nm) via the excited 2P1/2 level (denoted as |ecustom-character). As shown in FIG. 3, a laser beam from the laser may be split into a pair of non-copropagating laser beams (a first laser beam with frequency ω1 and a second laser beam with frequency ω2) in the Raman configuration, and detuned by a one-photon transition detuning frequency Δ=ω1−ω0e with respect to the transition frequency ω0e between |0custom-character and |ecustom-character, as illustrated in FIG. 3. A two-photon transition detuning frequency δ includes adjusting the amount of energy that is provided to the trapped ion by the first and second laser beams, which when combined is used to cause the trapped ion to transfer between the hyperfine states |0custom-character and |1custom-character. When the one-photon transition detuning frequency Δ is much larger than a two-photon transition detuning frequency (also referred to simply as “detuning frequency”) δ=ω1−φ2−ω01 (hereinafter denoted as ±μ, μ being a positive value), single-photon Rabi frequencies Ω0e(t) and Ω1e(t) (which are time-dependent, and are determined by amplitudes and phases of the first and second laser beams), at which Rabi flopping between states |0custom-character and |ecustom-character and between states |1custom-character and |ecustom-character respectively occur, and a spontaneous emission rate from the excited state |ecustom-character, Rabi flopping between the two hyperfine states |0custom-character and |1custom-character (referred to as a “carrier transition”) is induced at the two-photon Rabi frequency Ω(t). The two-photon Rabi frequency Ω(t) has an intensity (i.e., absolute value of amplitude) that is proportional to Ω0eΩ1e/2Δ, where Ω0e and Ω1e are the single-photon Rabi frequencies due to the first and second laser beams, respectively. Hereinafter, this set of non-propagating laser beams in the Raman configuration to manipulate internal hyperfine states of qubits (qubit states) may be referred to as a “composite pulse” or simply as a “pulse,” and the resulting time-dependent pattern of the two-photon Rabi frequency Ω(t) may be referred to as an “amplitude” of a pulse or simply as a “pulse,” which are illustrated and further described below. The detuning frequency δ=ω1−ω2−ω01 may be referred to as detuning frequency of the composite pulse or detuning frequency of the pulse. The amplitude of the two-photon Rabi frequency Ω(t), which is determined by amplitudes of the first and second laser beams, may be referred to as an “amplitude” of the composite pulse.


It should be noted that the particular atomic species used in the discussion provided herein is just one example of atomic species which have stable and well-defined two-level energy structures when ionized and an excited state that is optically accessible, and thus is not intended to limit the possible configurations, specifications, or the like of an ion trap quantum computer according to the present disclosure. For example, other ion species include alkaline earth metal ions (Be+, Ca+, Sr+, Mg+, and Ba+) or transition metal ions (Zn+, Hg+, Cd+).



FIG. 4 is provided to help visualize a qubit state of an ion is represented as a point on a surface of the Bloch sphere 400 with an azimuthal angle ϕ and a polar angle θ. Application of the composite pulse as described above, causes Rabi flopping between the qubit state |0custom-character (represented as the north pole of the Bloch sphere) and |1custom-character (the south pole of the Bloch sphere) to occur. Adjusting time duration and amplitudes of the composite pulse flips the qubit state from |0custom-character to |1custom-character (i.e., from the north pole to the south pole of the Bloch sphere), or the qubit state from |1custom-character to |0custom-character (i.e., from the south pole to the north pole of the Bloch sphere). This application of the composite pulse is referred to as a “π-pulse”. Further, by adjusting time duration and amplitudes of the composite pulse, the qubit state |0custom-character may be transformed to a superposition state |0custom-character+|1custom-character, where the two qubit states |0custom-character and |1custom-character are added and equally-weighted in-phase (a normalization factor of the superposition state is omitted hereinafter for convenience) and the qubit state |1custom-character to a superposition state |0custom-character−|1custom-character, where the two qubit states |0custom-character and |1custom-character are added equally-weighted but out of phase. This application of the composite pulse is referred to as a “π/2-pulse”. More generally, a superposition of the two qubits states |0custom-character and |1custom-character that are added and equally-weighted is represented by a point that lies on the equator of the Bloch sphere. For example, the superposition states |0custom-character±|1custom-character correspond to points on the equator with the azimuthal angle ϕ being zero and π, respectively. The superposition states that correspond to points on the equator with the azimuthal angle ϕ are denoted as |0custom-character+e1custom-character (e.g., |0custom-character±i|1custom-character for ϕ=±π/2). Transformation between two points on the equator (i.e., a rotation about the Z-axis on the Bloch sphere) can be implemented by shifting phases of the composite pulse.


Entanglement Formation



FIGS. 5A, 5B, and 5C depict a few schematic structures of collective transverse motional modes (also referred to simply as “motional mode structures”) of a group 106 of five trapped ions, for example. Here, the confining potential due to a static voltage V s applied to the end-cap electrodes 210 and 212 is weaker compared to the confining potential in the radial direction. The collective motional modes of the group 106 of trapped ions in the transverse direction are determined by the Coulomb interaction between the trapped ions combined with the confining potentials generated by the ion trap 200. The trapped ions undergo collective transversal motions (referred to as “collective transverse motional modes,” “collective motional modes,” or simply “motional modes”), where each mode has a distinct energy (or equivalently, a frequency) associated with it. A motional mode having the m-th lowest energy is hereinafter referred to as |nphcustom-characterm, where nph denotes the number of motional quanta (in units of energy excitation, referred to as phonons) in the motional mode, and the number of motional modes M in a given transverse direction is equal to the number of trapped ions in the group 106. FIGS. 5A-5C schematically illustrates examples of different types of collective transverse motional modes that may be experienced by five trapped ions that are positioned in a group 106. FIG. 5A is a schematic view of a common motional mode |nphcustom-characterm having the highest energy, where M is the number of motional modes. In the common motional mode |ncustom-characterm, all ions oscillate in phase in the transverse direction. FIG. 5B is a schematic view of a tilt motional mode |nphcustom-characterm-1 which has the second highest energy. In the tilt motional mode, ions on opposite ends move out of phase in the transverse direction (i.e., in opposite directions). FIG. 5C is a schematic view of a higher-order motional mode |nphcustom-characterm-3 which has a lower energy than that of the tilt motional mode |nphcustom-characterm-1, and in which the ions move in a more complicated mode pattern.


It should be noted that the particular configuration described above is just one among several possible examples of a trap for confining ions according to the present disclosure and does not limit the possible configurations, specifications, or the like of traps according to the present disclosure. For example, the geometry of the electrodes is not limited to the hyperbolic electrodes described above. In other examples, a trap that generates an effective electric field causing the motion of the ions in the radial direction as harmonic oscillations may be a multi-layer trap in which several electrode layers are stacked and an RF voltage is applied to two diagonally opposite electrodes, or a surface trap in which all electrodes are located in a single plane on a chip. Furthermore, a trap may be divided into multiple segments, adjacent pairs of which may be linked by shuttling one or more ions, or coupled by photon interconnects. A trap may also be an array of individual trapping regions arranged closely to each other on a micro-fabricated ion trap chip. In some embodiments, the quadrupole potential has a spatially varying DC component in addition to the RF component described above.


In an ion trap quantum computer, the motional modes may act as a data bus to mediate entanglement between two qubits and this entanglement is used to perform an XX gate operation. That is, each of the two qubits is entangled with the motional modes, and then the entanglement is transferred to an entanglement between the two qubits by using motional sideband excitations, as described below. FIGS. 6A and 6B schematically depict views of a motional sideband spectrum for an ion in the group 106 in a motional mode |nphcustom-characterM having frequency ωm according to one embodiment. As illustrated in FIG. 6B, when the detuning frequency of the composite pulse is zero (i.e., a frequency difference between the first and second laser beams is tuned to the carrier frequency, δ=ω1−ω2−ω01=0), simple Rabi flopping between the qubit states |0custom-character and |1custom-character (carrier transition) occurs. When the detuning frequency of the composite pulse is positive (i.e., the frequency difference between the first and second laser beams is tuned higher than the carrier frequency, δ=ω1−ω2−ω01=μ>0, referred to as a blue sideband), Rabi flopping between combined qubit-motional states |0custom-character|nphcustom-characterm and |1custom-character|nph+1custom-characterm occurs (i.e., a transition from the m-th motional mode with n-phonon excitations denoted by |nphcustom-characterm to the m-th motional mode with (nph+1)-phonon excitations denoted by |nph+1custom-characterm occurs when the qubit state |0custom-character flips to |1custom-character). When the detuning frequency of the composite pulse is negative (i.e., the frequency difference between the first and second laser beams is tuned lower than the carrier frequency by the frequency ωm of the motional mode |nphcustom-characterm, δ=ω1−ω2−ω01=−μ<0, referred to as a red sideband), Rabi flopping between combined qubit-motional states |0custom-character|nphcustom-characterm and |1custom-character|nph−1custom-characterm occurs (i.e., a transition from the motional mode |nphcustom-characterm to the motional mode |nph−1custom-characterm with one less phonon excitations occurs when the qubit state |0custom-character flips to |1custom-character). A π/2-pulse on the blue sideband applied to a qubit transforms the combined qubit-motional state |0custom-character|nphcustom-characterm into a superposition of |0custom-character|nphcustom-characterm and |1custom-character|nph+1custom-characterm. A π/2-pulse on the red sideband applied to a qubit transforms the combined qubit-motional |0custom-character|nphcustom-characterm into a superposition of |0custom-character|nphcustom-characterm and |1custom-character|nph−1custom-characterm. When the two-photon Rabi frequency Ω(t) is smaller as compared to the detuning frequency δ=ω1−ω2−ω01=±μ, the blue sideband transition or the red sideband transition may be selectively driven. Thus, a qubit can be entangled with a desired motional mode by applying the right type of pulse, such as a π/2-pulse, which can be subsequently entangled with another qubit, leading to an entanglement between the two qubits that is needed to perform an XX-gate operation in an ion trap quantum computer.


By controlling and/or directing transformations of the combined qubit-motional states as described above, an XX-gate operation may be performed on two qubits (i-th and j-th qubits). In general, the XX-gate operation (with maximal entanglement) respectively transforms two-qubit states |0custom-characteri|0custom-characterj, |0custom-characteri|1custom-characterj, |1custom-characteri|0custom-characterj, and |1custom-characteri|1custom-characterjas follows:

    • |0 custom-characteri|0custom-characterj→|0custom-characteri|0custom-characterj−i|1custom-characteri|1custom-characterj
    • |0 custom-characteri|1custom-characterj→|0custom-characteri|1custom-characterj−i|1custom-characteri|0custom-characterj
    • |1custom-characteri|0custom-characterj→−i|0custom-characteri|1custom-characterj+|1custom-characteri|0custom-characterj
    • |1custom-characteri|1custom-characterj→−i|0custom-characteri|0custom-characterj+|1custom-characteri|1custom-characterj

      For example, when the two qubits (i-th and j-th qubits) are both initially in the hyperfine ground state |0custom-character (denoted as |0custom-characteri|0custom-characterj) and subsequently a π/2-pulse on the blue sideband is applied to the i-th qubit, the combined state of the i-th qubit and the motional mode |0custom-characteri|nphcustom-characterm is transformed into a superposition of |0custom-characteri|nphcustom-characterm and |1custom-characteri|nph+1custom-characterm, and thus the combined state of the two qubits and the motional mode is transformed into a superposition of |0custom-characteri|0custom-characterj|nphcustom-characterm and |1custom-characteri|0custom-characterj|nph+1custom-characterm. When a π/2-pulse on the red sideband is applied to the j-th qubit, the combined state of the j-th qubit and the motional mode |0custom-characterj|nphcustom-characterm is transformed to a superposition of |0custom-characterj|nphcustom-characterm and |1custom-characterj|nph−1custom-characterm and the combined state |0custom-characterj|nph+1custom-characterm is transformed into a superposition of |0custom-characterj|nph+1custom-characterm and |1custom-characterj|nphcustom-characterm.


Thus, applications of a π/2-pulse on the blue sideband on the i-th qubit and a π/2-pulse on the red sideband on the j-th qubit may transform the combined state of the two qubits and the motional mode |0custom-characteri|0custom-characterj|nphcustom-characterm into a superposition of |0custom-characteri|0custom-characterj|nphcustom-characterm and |1custom-characteri|1custom-characterj|nphcustom-characterm, the two qubits now being in an entangled state. For those of ordinary skill in the art, it should be clear that two-qubit states that are entangled with motional mode having a different number of phonon excitations from the initial number of phonon excitations nph (i.e., |1custom-characteri|0custom-characterj|nph+1custom-characterm and |0custom-characteri|1custom-characterj|nph−1custom-characterm) can be removed by a sufficiently complex pulse sequence, and thus the combined state of the two qubits and the motional mode after the XX-gate operation may be considered disentangled as the initial number of phonon excitations nph in the m-th motional mode stays unchanged at the end of the XX-gate operation. Thus, qubit states before and after the XX-gate operation will be described below generally without including the motional modes.


More generally, the combined state of i-th and j-th qubits transformed by the application of pulses on the sidebands for duration T (referred to as a “gate duration”), having amplitudes Ω(i) and Ω(i) and detuning frequency μ, can be described in terms of an entangling interaction x(i,j)(t) as follows:

    • |0custom-characteri|0custom-characterj→cos (2x(i,j)(t))|0custom-characteri|0custom-characterj−i sin (2x(i,j)(t))|1custom-characteri|1custom-characterj
    • |0custom-characteri|1custom-characterj→cos (2x(i,j)(t))|0custom-characteri|1custom-characterj−i sin (2x(i,j)(t))|1custom-characteri|0custom-characterj
    • |1custom-characteri|0custom-characterj→−i sin (2x(i,j)(t))|0custom-characteri|1custom-charactern,j+cos (2X(i,j)(t))|1custom-characteri|0custom-characterj
    • |1custom-characteri|1custom-characterj→−i sin (2X(i,j)(t))|0custom-characteri|0custom-characterj+cos (2X(i,j)(t))|1custom-characteri|1custom-characterj

      where,








χ

(

i
,
j

)


(
τ
)

=


-
4






m
=
1

M



η
m

(
i
)




η
m

(
j
)






0
τ



dt
2





0

t
2




dt
1




Ω

(
i
)


(

t
2

)




Ω

(
j
)


(

t
1

)



cos

(

μ


t
2


)



cos

(

μ


t
1


)



sin
[


ω
m

(


t
2

-

t
1


)

]












and nm(i) is the Lamb-Dicke parameter that quantifies the coupling strength between the i-th ion and the m-th motional mode having the frequency Wm, and M is the number of the motional modes (equal to the number N of ions in the group 106).


The entanglement interaction between two qubits described above can be used to perform an XX-gate operation. The XX-gate operation (X(gate) along with single-qubit operations (R gates) forms a set of elementary logic operations {R,XX} that can be used to build a quantum computer that is configured to perform desired computational processes. Among several known sets of elementary logic operations by which any quantum algorithm can be decomposed, a set of elementary logic operations, commonly denoted as {R, XX}, is native to a quantum computing system of trapped ions described herein. Here, the R gate corresponds to manipulation of individual qubit states of trapped ions, and the XX gate (also referred to as an “entangling gate”) corresponds to manipulation of the entanglement of two trapped ions.


To perform an XX-gate operation between i-th and j-th qubits, pulses that satisfy the condition x(i,j)(t)=θ(i,j)(0<θ(i,j)≤π/8) (i.e., the entangling interaction x(i,j)(t) has a desired value OM, referred to as condition for a non-zero entanglement interaction) are constructed and applied to the i-th and the j-th qubits. The transformations of the combined state of the i-th and the j-th qubits described above corresponds to the XX-gate operation with maximal entanglement when θ(i,j)=π/8. Amplitudes Ω(i)(t) and Ω(i)(t) of the pulses to be applied to the i-th and the j-th qubits are control parameters that can be adjusted to ensure a non-zero tunable entanglement of the i-th and the j-th qubits to perform a desired XX gate operation on i-th and j-th qubits.


Other two-qubit gate operations, such as a controlled-NOT (referred to as “CNOT”) gate and a SWAP gate, may be used in the method of pattern matching described below and implemented by use of a combination of the XX gate operations and single-qubit gate operations.


In general, a CNOT gate operation conditioned on i-th qubit (control bit) and targeted on j-th qubit (target bit) inverts the j-th qubit (target bit) if the i-th qubit (control bit) in state |1custom-character, and leaves both the i-th and j-th qubits unchanged otherwise, thus transforming a two-qubit state |xcustom-characteri|ycustom-characterj(x,y={0,1}) to a two-qubit state |xcustom-characteri|x⊕ycustom-character:




















"\[LeftBracketingBar]"

x



i





"\[LeftBracketingBar]"

y




j



CNOT




"\[LeftBracketingBar]"

x




i





"\[LeftBracketingBar]"


x


y






j

.




The output of the j-th qubit (target bit) is the exclusive OR (XOR, denoted by “⊕”) of the i-th and j-th qubits. A SWAP gate operation between i-th and j-th qubits swaps the i-th and j-th qubits, thus transforming a two-qubit state |xcustom-characteri|ycustom-characterj(x,y={1,1}) to a two-qubit state


























"\[LeftBracketingBar]"

y



i





"\[LeftBracketingBar]"

x




j

:




"\[LeftBracketingBar]"

x




i





"\[LeftBracketingBar]"

y




j



SWAP




"\[LeftBracketingBar]"

y




i





"\[LeftBracketingBar]"

x




j

.




Three-qubit gate operations, such as a controlled-controlled-NOT gate (also referred to as “Toffoli gate”) and a controlled-SWAP gate (also referred to as “Fredkin gate”), may also be used in the method of pattern matching described below and implemented by use of a combination of the XX gate operations and single-qubit gate operations. A Toffoli gate operation applied to two control bits and one target bit inverts the target bit only if both of the control bits are in state II), and leaves all three qubits unchanged otherwise, thus transforming a three-qubit state |xcustom-character|ycustom-character|zcustom-character(x,y,z={0,1}) to a three-qubit state |xcustom-character|ycustom-character|(x y)⊕z):




















"\[LeftBracketingBar]"

x







"\[LeftBracketingBar]"

y








"\[LeftBracketingBar]"

z






Toffoli




"\[LeftBracketingBar]"

x








"\[LeftBracketingBar]"

y








"\[LeftBracketingBar]"



(

x
·
y

)


z





.




The output of the target bit includes XOR (⊕) of the target bit and AND (denoted by “·”) of the control qubits. A Fredkin gate operation applied to one control bit and two target bits swaps the target bits if the control bit is state |1custom-character, and laves all three qubits unchanged otherwise, thus transforming a three-qubit state |xcustom-character|ycustom-character|zcustom-character(x,y,z={0,1}) to a three-qubit state |xcustom-character|(x·Z)⊕(x·Y))|(X·Y)⊕(x·Z)):




















"\[LeftBracketingBar]"

x







"\[LeftBracketingBar]"

y








"\[LeftBracketingBar]"

z






Fredkin




"\[LeftBracketingBar]"

x








"\[LeftBracketingBar]"



(

x
·
z

)



(


x
_

·
y

)










"\[LeftBracketingBar]"



(

x
·
y

)



(


x
_

·
z

)






,




where the output of the target bits include NOT of the target bit (denoted by “−”).


Accelerated Pattern Matching Method on a Quantum Computing System

In a quantum computing system, a quantum computer is a domain-specific accelerator (also referred to as a “quantum processor” hereinafter) that is able to accelerate certain computational tasks beyond the reach of classical computers. Examples of such computational tasks include searching a given set of data that can be used for pattern matching. Pattern matching is used in a wide range of applications, such image processing, study of DNA sequences, and data compression and statistics. However, complexity of pattern matching increases drastically as a problem size (e.g., the lengths of a string and a pattern to be searched within the string) increases and may be unsolvable or be too complex to complete in a reasonable amount of time by a classical computer alone. Thus, methods for accelerated pattern matching on a quantum computing system are described herein.


It should be noted that the particular example embodiments described herein are just some possible examples of a quantum computing system according to the present disclosure and do not limit the possible configurations, specifications, or the like of quantum computer systems according to the present disclosure. For example, a quantum computing system according to the present disclosure can be applied to other types of computational task in which searching of data contributes to the computational complexity and can be accelerated by use of a quantum processor.



FIG. 7 depicts a flowchart illustrating a method 700 of searching for a bit pattern within a bit string. A bit string and a bit pattern are sequences of bits, each of which has single binary value, either 0 or 1. In the example described herein, the quantum processor is the group 106 of trapped ions, in which the two hyperfine states of each of the trapped ions form a qubit that has binary qubit states |0custom-character and |1custom-character.


In block 702, by the classical computer 102, a pattern matching problem is selected, for example, by use of a user interface of the classical computer 102, or retrieved from the memory of the classical computer 102. Specifically, a bit string custom-character (t0t1, . . . tN−1) of length N and a bit pattern custom-character(p0p1 . . . PM−1) of length M with M≤N are selected. The pattern matching problem is to search and locate the bit pattern custom-character within the bit string custom-character. Pattern matching may be exact match (i.e., a M-bit long sequence contained in the bit string custom-character exactly matches the bit pattern custom-character) or fuzzy match (i.e., (M−d) bits of a M-bit long sequence contained in the bit string custom-character matches the bit pattern custom-character, where d is a predetermined number indicating fuzziness in the matching).


In block 704, by the system controller 104, the quantum processor 106 is set in an initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterP. The first register (denoted by the subscript “I” and referred to also as an “index register” hereinafter) is formed of n=log2/N qubits |k0custom-character|k1custom-character. . . |kn−1custom-character=⊗i=0n−1|kicustom-character to be used to encode a string index k(=20k0+21k1+ . . . +2n−1kn−1=E1=0n−12jKj=0,1, . . . ,N−1). In the initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterp, all qubits of the index register are prepared in state |0custom-character(i.e., |0custom-characterI=|0custom-character⊗n), for example, the hyperfine ground state |0custom-character, by optical pumping in an exemplary quantum computer with trapped ions. The second register (denoted by the subscript “T” and referred to also as a “string register” hereinafter) is formed of N qubits to be used to encode the bit string custom-character(t0t1 . . . tN−1). In the initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterp, all qubits of the string register are prepared in state |0custom-character(i.e., |0custom-characterT=|0custom-character⊗N). The third register (denoted by the subscript “P” and referred to as a “pattern register” hereinafter) is formed of M qubits to be used to encode the bit pattern custom-character. In the initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterp, all qubits of the pattern register are prepared in state |0custom-character(i.e., |0custom-characterp=|0custom-characterεM).


In block 706, by the system controller 104, each qubit of the index register is transformed from state |0custom-character to a linear superposition state of |0custom-character and |1custom-character, (|0custom-character+|1custom-character)/√{square root over (2)}, thus transforming the index register from the initial state |0custom-characterI=|0custom-character on to a superposition state








(







"\[LeftBracketingBar]"

0



+



"\[LeftBracketingBar]"

1





2


)



n


=




j
=
0


n
-
1



(








k
j

=

{

0
,
1

}









"\[LeftBracketingBar]"


k
j





2



)






in the quantum processor 106. The index register is now in a superposition of 2n(=N) states (|0custom-character|0custom-character. . . |0custom-character,|0custom-character|0custom-character. . . |1custom-character, . . . , |1custom-character|1custom-character. . . |1custom-character) each associated with a string index k(=0, 1, . . . , N−1) (thus, referred to as string index states), and may be denoted as










1

N









k
=
0


N
-
1






"\[LeftBracketingBar]"

k




I

,




where |kcustom-characterI is a stale formed of n (=log2N) qubits defined as |kcustom-characterI=|k0custom-character|k1custom-character. . . |kn−1custom-character(k=Σj=0n−12jkj=0,1, . . . , N−1). By this transformation, the quantum processor 106 is transformed from the initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterp to an initial superposition state


















"\[LeftBracketingBar]"


ψ
S




=


1

N









k
=
0


N
-
1






"\[LeftBracketingBar]"

k





I





"\[LeftBracketingBar]"

0




T





"\[LeftBracketingBar]"

0




P




(i.e., a superposition of combined states of the index register, the string register, and the pattern register). This transformation may be implemented by application of a proper combination of single-qubit operations to the n=(log2N) qubits of the index register in the initial state |Ψ0custom-character=|0custom-characterI|0custom-characterT|0custom-characterp in a single time step in units of elementary logic operations.


In block 708, by the system controller 104, the bit string custom-character(t0t1 . . . tN−1) is encoded in the string register, and the bit pattern custom-character(p0p1 . . . pM−1) is encoded in the pattern register in the quantum processor 106. Specifically, the i-th bit information ti(i=0,1, . . . N−1) of the bit string custom-character(t0t1 . . . tN−1) is encoded in the i-th qubit of the string register as |custom-charactercustom-characterT=|t0custom-character|t1custom-character. . . |tN−1custom-character=⊗i=0N−1|ticustom-character(i=0,1, . . . N−1) (referred to as a string encoded state). That is, the string register is transformed form the initial state of the string register |0custom-characterT to the string encoded state |custom-charactercustom-characterT. Similarly, the i-th bit information pi(i=0,1, . . . M−1) of the bit pattern custom-character(p0p1 . . . PM) is encoded in the i-th qubit of the pattern register as |custom-charactercustom-characterp=|p0custom-character|p1custom-character. . . |pM−1custom-character=⊗i=0M−1|picustom-character(i=0, 1, . . . M−1) (referred to as a pattern encoded state). That is, the pattern register is transformed from the initial state of the pattern register |0custom-characterp to the pattern encoded state βcustom-charactercustom-characterp. These operations may be implemented by application of a proper combination of single-qubit operations to the N qubits of the string register in the initial state |0custom-characterT and to the M qubits of the pattern register in the initial state |0custom-characterp in a single time step in units of elementary logic operations.


In block 710, by the system controller 104, a controlled-circular-shift operation S is applied on the string register in the string encoded state |custom-charactercustom-characterT=⊗i=0N−1|ticustom-character conditioned on the index register in the initial superposition state









1

N









k
=
0


N
-
1






"\[LeftBracketingBar]"

k




I




in the quantum processor 106. That is, all qubits of the string register are circularly shifted (e.g., left-circularly shifted) by k bit positions in a combined state having the index register in state |kcustom-characterI (referred to as an operation Sk):






custom-character[|kcustom-characterIi=0M−1|ticustom-character]=|kcustom-characterIi=0N−1custom-characterk|ticustom-character=|kcustom-characterIi=0N−1|ti+k mod Ncustom-character.


Thus, the operation Sk transforms the string register in the state having the index register in |kcustom-characterI from the string encoded state |custom-charactercustom-characterT=⊗i=0N−1|ticustom-character to a circularly shifted state |custom-charactercustom-characterT(k)=⊗i=0N−1|ti+k mod Ncustom-character.


The controlled-circular-shift operation Sk conditioned on the n-qubit long index register |kcustom-characterI may be implemented as a series of bit-wise circular shift operation custom-character2j(k) conditioned on the j-th qubit |kicustom-character of the index register that left-circularly shifts all qubits of the string register by 2j bit positions: |kcustom-characterIi=0N−1custom-characterk|ticustom-character=(⊗j=0n−1|kjcustom-character)⊗i=0N−1j=0n−1custom-character2j(kj))|ticustom-character.


The bit-wise controlled-circular-shift operator S2j(kj) may be implemented as a combination of controlled-SWAP operations targeted on the string register (each of which interchanges a pair among the N qubits of the string register) conditioned on the j-th qubit |kjcustom-character of the index register. Among the required controlled-SWAP operations, N/2 controlled-SWAP operations to circularly shift N/2 qubits among the total N qubits by 2j bit positions may be applied in a single time step in units of elementary logic operations as described in detail below. In a subsequent time step, N/4 controlled-SWAP operations are applied to circularly shift a half of the N/2 remaining qubits by 2j bit positions. At each subsequent time step, the number of qubits that need to be swapped decreased by half. Therefore, the bit-wise controlled-circular-shift operator custom-character2j(kj) targeted on N qubits may be implemented in custom-character(log(N)) time steps in units of elementary logic operations using parallel controlled-SWAP operations.


To implement, in one time step in units of elementary logic operations, at most N/2 controlled-SWAP operations controlled by the same control qubit, N/2 ancillary qubits, each initialized to state |0custom-character, may be used. Using a fan-out CNOT operation on the ancillary qubits with the control qubit, the control qubit state may be copied into the ancillary qubits. With the N/2 available copies of the control qubit state, each controlled-SWAP operations may then be implemented in parallel since each one of at most N/2 controlled-SWAP operations may use each one of the available control state copies encoded in the ancillary qubits. Once all of controlled-SWAP operations in a given bit-wise i) controlled-circular-shift operator custom-character2j(kj) are implemented, the ancillary qubits may be restored to state |0custom-character by subsequently applying another fan-out CNOT operation as described above. The fan-out operation may be implemented in custom-character(log(N)) time steps in units of elementary logic operations.


The CNOT gate operations and the controlled-SWAP operations used to perform the controlled-circular-shift operation custom-character may be each implemented by application of a proper combination of single-qubit operations and two-qubit gate operations to the n=(log2N) qubits of the index register in the initial superposition state










1

N









k
=
0


N
-
1






"\[LeftBracketingBar]"

k




I

,




the N qubits of the string register in the string encoded state |custom-charactercustom-characterT=⊗i=0N−1|ticustom-character, and N/2 ancillary qubits, each initialized to state |0custom-character.


In block 712, by the system controller 104, a bit-wise CNOT gate operation is applied on the pattern register in the pattern encoded state |custom-charactercustom-characterp=⊗i=0M−1|picustom-character conditioned on the first M bits of the string register in the string encoded state |custom-charactercustom-characterT(k)=⊗i=0N−1|ti+k mod Ncustom-character. That is, the bit-wise CNOT gate operation brings the pattern register from the pattern encoded state |custom-charactercustom-characterp=⊗i=0M−1|picustom-character to a string-pattern matching evaluation state (⊗i=0M−1|pi⊗ti+k mod Ncustom-character) After this transformation, if the first M bits of the string register in the circularly shifted state |custom-charactercustom-characterT(k)=⊗i=0N−1|ti+k mod N)exactly match all qubits of the pattern register in the pattern encoded state |custom-charactercustom-characterp=⊗i=0M−1|picustom-character (i.e., exact matching), all qubits |pi⊗ti+k mod Ncustom-character of the pattern register in the string-pattern matching evaluation state (⊗i=0M−1|pi⊗ti+k mod Ncustom-character) are in state |0custom-character. If the first M qubits except for d qubits of the string register in the circularly shifted state |custom-charactercustom-characterT(k)=⊗i=0N−1|ti+k mod Ncustom-character match the qubits of the pattern register in the pattern encoded state |custom-charactercustom-characterp=⊗i=0M−1|picustom-character (i.e., fuzzy matching), (M−d) qubits |pi⊗ti+k mod N) of the pattern register in the string-pattern matching evaluation state (⊗i=0M−1|pi⊗ti+k mod Ncustom-character) are in state |0custom-character and d bits are in state |1custom-character.


The bit-wise CNOT gate operations can be implemented by the application of a combination of single-qubit gate operations and two-qubit gate operations by the system controller 104 in a single time step in units of elementary logic operations.


In block 714, by the system controller 104, the string-pattern matching evaluation state (⊗i=0M−1|pi⊗ti+k mod Ncustom-character) having (M−d) or more qubits in state |0custom-character (referred to as a matched state) is amplified and measured. The exact matching corresponds to d=0. The fuzzy matching corresponds to d≥0. The string index k specifies the location of the matched string. The amplification of the matched state may be performed by applying an oracle operation U W on the pattern register, which reverses a phase of the matched state:











i
=
0


M
-
1





"\[LeftBracketingBar]"


x
i








U
w



{







(

-
1

)




i
=
0


M
-
1





"\[LeftBracketingBar]"


x
i









if






i
=
0


M
-
1



x
i




d








(

+
1

)




i
=
0


M
-
1





"\[LeftBracketingBar]"


x
i









if






i
=
0


M
-
1



x
i




d




.






Oracle operations, such as the oracle operation UW described herein, may be implemented in custom-character(log(M)) time steps in units of elementary logic operations using custom-character(M) ancillary qubits in a quantum processor that has a long-range interactions, such as the exemplary quantum computer with trapped ions and needs to be repeated custom-character(√{square root over (N)}) times for successful amplification. Once an amplitude of the index register (in the superposition of string index states) is measured, the matched state in which the string-pattern matching evaluation state (⊗i=0M−1|pi⊗ti+k mod Ncustom-character) having (M−d) or more qubits in state |0custom-character has the largest amplitude. The string index k associated with the measured matched state corresponds to the location of the pattern custom-character within the bit string custom-character. Subsequently, the classical computer 102 generally outputs the measured result including whether or not the bit string custom-character contains the bit pattern custom-character and if it does, the location of the pattern custom-character within the bit string custom-character.


In the embodiments described herein, the methods for a more efficient and accelerated pattern matching process is provided. The amount of time resources required to run a matching process for searching a bit pattern custom-character of length M within a bit string custom-character of length N by the method described herein scales as custom-character(√{square root over (N)}(log(N)2+log(M))) time steps in units of elementary logic operations. As described above, transformation of the index register to the superposition state in block 706 and encoding of the bit string T and the pattern custom-character in block 708 may be each implemented in a single time step in units of elementary logic operations. In block 710, each bit-wise controlled-circular-shift operator custom-character2j(kj) may be implemented in custom-character(log(N)) time steps in units of elementary logic operations, and thus the total controlled-circular-shift operation custom-character(that is n=(log2N) of the bit-wise controlled-circular shift operations) may be implemented in custom-character(log(N)2) time steps in units of elementary logic operations. The application of bit-wise CNOT gate operation in block 712 may be implemented in a single time step in units of elementary logic operations. The evaluation of matching of the bit string custom-character and the pattern custom-character using the oracle operation UW in block 708 may be implemented in custom-character(log(M)) time steps in units of elementary logic operations. For amplitude amplification, this process is effectively repeated custom-character(√{square root over (N)}) times. Thus, in total, the pattern matching by the method described herein may be performed in custom-character(√{square root over (N)}(log(N)2+log(M))) time steps in units of elementary logic operations.


The amount of memory resources required to run a matching process for searching a bit pattern custom-character of length M within a bit string custom-character of length N by the method described herein scales as custom-character(N+M) qubits. As described above, the amount of memory resources includes N to encode the bit string custom-character(t0t1 . . . tN−1) of length N, M to encode the bit pattern custom-character(p0p1 . . . PM−1) of length M, n=log2N to encode the string index, N/2 of ancillary qubits to implement parallel controlled-SWAP operations, and custom-character(M) of ancillary qubits to implement the oracle operation Uw. Thus, the method requires custom-character(N+M) qubits in total. Therefore, the method described herein provides improvement in the amount of time required to run a matching process by use of a quantum processor over a purely classical computational method, by allowing simultaneous search of a bit pattern in multiple locations in a bit string due to the use of the multiple states created by a quantum superposition.


While the foregoing is directed to specific embodiments, other and further embodiments may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.

Claims
  • 1. A method of determining a pattern in a sequence of bits using a quantum computing system comprising a classical computer and a quantum processor comprising a plurality of qubits, the method comprising: circularly shifting qubits of a second register conditioned on a first register in each string index state in a superposition of a plurality of string index states, each of which is associated with a string index;amplifying an amplitude of a state combined with the first register in a string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of a third register; andmeasuring an amplitude of the first register and determining the string index state of the plurality of string index states associated with the amplified state.
  • 2. The method according to claim 1, further comprising: prior to the circularly shifting of the qubits of the second register, setting the first register of the quantum processor in the superposition of the plurality of string index states;encoding a bit string in the second register of the quantum processor; andencoding a bit pattern in the third register of the quantum processor, andsubsequent to the measuring of the amplitude of the first register, outputting, by use of the classical computer, the string index associated with the first register in the measured state.
  • 3. The method according to claim 2, wherein encoding the bit string in the second register of the quantum processor comprises applying a combination of single-qubit gate operations to the second register of the quantum processor.
  • 4. The method according to claim 2, wherein encoding the bit pattern in the third register of the quantum processor comprises applying a combination of single-qubit gate operations to the third register of the quantum processor.
  • 5. The method according to claim 1, wherein each qubit of the plurality of qubits comprises a trapped ion having two frequency-separated states.
  • 6. The method according to claim 5, further comprising: preparing the quantum processor in an initial state by setting each trapped ion in the quantum processor in the lower energy state of the two frequency-separated states.
  • 7. The method according to claim 6, wherein setting the first register of the quantum processor in the superposition of the plurality of string index states comprises transferring each trapped ion in the first register in a superposition of the two frequency-separated states.
  • 8. The method according to claim 6, wherein circularly shifting qubits of the second register conditioned on the first register comprises applying a combination of single-qubit gate operations and two-qubit gate operations to the first and second registers of the quantum processor.
  • 9. The method according to claim 8, wherein the two-qubit gate operations comprise controlled-SWAP operations applied to the second register conditioned on the first register.
  • 10. The method according to claim 8, wherein the two-qubit gate operations comprise controlled-SWAP operations applied to the second register conditioned on a plurality of ancillary qubits in which the first register is copied by fan-out CNOT operations applied on the plurality of ancillary qubits conditioned on the first register.
  • 11. The method according to claim 5, wherein amplifying the state combined with the first register in the string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of the third register comprises applying a combination of single-qubit gate operations and two-qubit gate operations to the second and third registers of the quantum processor.
  • 12. A quantum computing system, comprising: a quantum processor comprising a plurality of qubits, each qubit of the plurality of qubits comprising a trapped ion having two frequency-separated states;one or more lasers configured to emit a laser beam, which is provided to trapped ions in the quantum processor;a classical computer; anda system controller configured to: circularly shift qubits of a second register conditioned on a first register in each string index state in a superposition of a plurality of string index states, each of which is associated with a string index;amplify an amplitude of a state combined with the first register in a string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of a third register; andmeasure an amplitude of the first register and determining the string index state of the plurality of string index states associated with the amplified state.
  • 13. The quantum computing system according to claim 12, wherein the system controller is further configured to: prior to the circularly shifting of the qubits of the second register, set a first register of the quantum processor in a superposition of a plurality of string index states, each of which is associated with a string index, wherein the quantum processor comprises a plurality of qubits;encode a bit string in a second register of the quantum processor; andencode a bit pattern in the third register of the quantum processor, andthe classical computer is further configured to: subsequent to the measuring of the amplitude of the first register, output the string index associated with the first register in the measured state.
  • 14. The quantum computing system according to claim 13, wherein encoding the bit string in the second register of the quantum processor comprises applying a combination of single-qubit gate operations to the second register of the quantum processor, andencoding the bit pattern in the third register of the quantum processor comprises applying a combination of single-qubit gate operations to the third register of the quantum processor.
  • 15. The quantum computing system according to claim 12, wherein the system controller is further configured to: prepare the quantum processor in an initial state by setting each trapped ion in the quantum processor in the lower energy state of the two frequency-separated states.
  • 16. The quantum computing system according to claim 12, wherein setting the first register of the quantum processor in the superposition of the plurality of string index states comprises transferring each trapped ion in the first register in a superposition of the two frequency-separated states.
  • 17. The quantum computing system according to claim 12, wherein circularly shifting qubits of the second register conditioned on the first register comprises applying, by the system controller, a combination of single-qubit gate operations and two-qubit gate operations to the first and second registers of the quantum processor.
  • 18. The quantum computing system according to claim 17, wherein the two-qubit gate operations comprise controlled-SWAP operations applied to the second register conditioned on the first register.
  • 19. The quantum computing system according to claim 17, wherein the two-qubit gate operations comprise controlled-SWAP operations applied to the second register conditioned on a plurality of ancillary qubits in which the first register is copied by fan-out CNOT operations applied on the plurality of ancillary qubits conditioned on the first register.
  • 20. The quantum computing system according to claim 12, wherein amplifying the state combined with the first register in the string index state of the plurality of string index states in which the circularly shifted qubits of the second register matches qubits of the third register comprises applying, by the system controller, a combination of single-qubit gate operations and two-qubit gate operations to the second and third registers of the quantum processor.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation application of U.S. patent application Ser. No. 17/313,671, filed May 6, 2021, which claims the benefit to U.S. Provisional Application No. 63/031,385, filed May 28, 2020, both of which are incorporated by reference herein.

GOVERNMENT LICENSE RIGHTS

This invention was made with Government support under DESC0019040 and DE-SC0020312 awarded by the Department of Energy. The Government has certain rights in the invention.

Provisional Applications (1)
Number Date Country
63031385 May 2020 US
Continuations (1)
Number Date Country
Parent 17313671 May 2021 US
Child 18369327 US