ORIGINAL RESEARCH article
Sec. Quantum Engineering and Technology
Efficient Criteria of Quantumness for a Large System of Qubits
- 1Department of Physics, Faculty of Science Division I, Tokyo University of Science, Shinjuku, Japan
- 2Division of Nano-Quantum Information Science and Technology, Research Institute for Science and Technology, Tokyo University of Science, Shinjuku, Japan
- 3Department of Physics, University of Notre Dame, Notre Dame, IN, United States
- 4Notre Dame Institute for Advanced Study, University of Notre Dame, Notre Dame, IN, United States
- 5Research Center for Emerging Computing Technologies (RCECT), National Institute of Advanced Industrial Science and Technology (AIST), Ibaraki, Japan
- 6Department of Physics, Loughborough University, Loughborough, United Kingdom
In order to model and evaluate large-scale quantum systems, e.g., quantum computer and quantum annealer, it is necessary to quantify the “quantumness” of such systems. In this paper, we discuss the dimensionless combinations of basic parameters of large, partially quantum coherent systems, which could be used to characterize their degree of quantumness. Based on analytical and numerical calculations, we suggest one such number for a system of qubits undergoing adiabatic evolution, i.e., the accessibility index. Applying it to the case of D-Wave One superconducting quantum annealing device, we find that its operation as described falls well within the quantum domain.
One of the key obstacles in the way to the full development of quantum technologies 2.0  is the same circumstance which stimulated their development in the first place: the fundamental impossibility of an efficient simulation of a large enough, quantum coherent structure with classical means. In practice “large enough” turned out to be a system comprising a hundred or so quantum bits, which is still too small to form a quantum computer capable of simulating other “large enough” quantum systems. On the other hand, artificial quantum coherent systems comprising thousands of qubits are being fabricated  and even successfully used, like commercial quantum annealers [3, 4]. Arrays of superconducting qubits are also being considered as microwave range detectors capable of exceeding the standard quantum limit (in such application as, e.g., search for galactic axions ). Quantum coherence of the array is the key element of the detection mechanism. This “quantum capacity gap”  needs to be bridged, in order to allow a systematic progress towards the development of the full potential of quantum technologies 2.0, such as noisy intermediate-scale quantum (NISQ) devices  and universal fault tolerant quantum computers.
The impossibility of an efficient classical simulation of a large quantum system is not absolute, in the sense that it concerns the simulation of an arbitrary evolution of such a system, whereby its state vector can reach all of its (exponentially high-dimensional) Hilbert space and has potentially infinite time to do so. The Margolus-Levitin theorem and its generalizations [8–13] put a limit on the speed of such evolution, thus restricting the accessible part of the Hilbert space for any finite time interval. This agrees with a proof  that the manifold of all quantum many-body states that can be generated by arbitrary time-dependent local Hamiltonians in a time that scales polynomially in the system size occupies an exponentially small volume in its Hilbert space. (This is a literally correct statement, since the Hilbert space of a system of qubits is a finite-dimensional complex projective space; that is, it is compact and, moreover, it has a unitary invariant Fubini-Study metric ). Numerical and analytical studies also indicate that the number of independent constraints describing quantum evolution may be much less than the dimensionality of the Hilbert space . It is therefore reasonable to suggest that a “general case” evolution of a large quantum coherent quantum structure can be characterized by a non-exponentially large set of dimensionless parameters, which correspond to qualitatively different regimes of evolution of this structure. Our recent numerical simulations indicate the existence of such regimes in a set of qubits with pumping and dissipation .
Such dimensionless parameters, if exist, will be combinations of fundamental physical constants and parameters, which characterize the system. (We will only need the Planck constant, since the speed of light and gravity constant are not relevant for the currently feasible devices). For example, in the standard approximation, a system of qubits is described by a quantum Ising Hamiltonian,
and a set of Lindblad operators responsible for dephasing and relaxation of separate qubits, with characteristic times, respectively, tϕ and tr, and their combination, the decoherence time tD. In the case of an adiabatic quantum processor [3, 4], the time dependence of the Hamiltonian parameters Jij, hj and Δj (except that induced by the ambient noise) is determined by that of the adiabatic parameter λ (t) (in the simplest case, by the time of adiabatic evolution, tf). Then, the dimensionless characteristics of the system should be the combinations of ℏ with the following quantities:
(1) Dimensionless: N (number of qubits); ⟨z⟩ (connectivity of the network: average number of couplings per qubit); ⟨δz2⟩ (its dispersion); ⟨zizj⟩ (its correlation function); …;
(2) Powers of energy: ⟨E⟩ (average qubit excitation energy); ⟨δE2⟩; ⟨EiEj⟩; ⟨J⟩ (average coupling strength); ⟨δJ2⟩; ⟨JiJk⟩; kBT; …;
(3) Powers of time: tϕ; tr; tD; tf;
and such additional parameters as, e.g., the spectral density of ambient noise SA (f). Note that all these parameters can be efficiently obtained by either direct measurements or straightforward calculations.
The field of several dozens of independent dimensionless combinations of the above parameters is narrowed for a particular quantum system and the mode of its operation. Here we concentrate on adiabatic quantum computing (quantum annealing). The Hamiltonian of a quantum system is here manipulated in such a way that its ground state changes from the easily accessible one to the one encoding a solution to the desired problem (and presumably having a very complex structure, so that reaching it by annealing would be improbable). In the case of a slow enough evolution of the Hamiltonian H (t), the system initialized in the ground state of H (0) will evolve into the ground state of H (tf) by virtue of the adiabatic theorem . If the system is totally insulated, the quantum speed and accessibility theorems [8–14] do not put fundamental constraints on an adiabatic quantum computer. Nevertheless, in a realistic case, the computation time is limited by interactions with the outside world leading to nonunitarity , and the question arises whether the evolution from the initial to the desired final state of the system is possible. As we see, the question of accessibility of different regions of the Hilbert space is especially relevant in this case.
An intriguing twist is added by the fact that the operation of D-Wave processors demonstrated what looked convincingly like quantum annealing [20, 21] despite the large discrepancy between the adiabatic evolution time, tf (microseconds ), and the qubit decoherence time, tD ≪ tf (tens of nanoseconds ), in the absence of any quantum error correction. On the second thought, it is not so surprising. The quantum state of an adiabatic quantum computer evolves starting at t = 0 from a factorized state, |in⟩. The computation is successful, if at the time t = tf there is a sufficient ratio of quantum trajectories ending in the state |out⟩, which is also factorized by design. Decoherence tends to disrupt quantum correlations between different qubits, and thus constrain the trajectories to partially factorized submanifolds of the Hilbert space. Nevertheless the success does not necessarily require that these trajectories pass through globally entangled states, and thus certain degree of decoherence may not necessarily make the proper operation of the device impossible.
While it cannot be predicted whether the evolution of a given quantum system can take it from the given initial to the desired final state, the average of the maximal distance between some initial and some final state of the system, for given values of tf and other system parameters, may serve as a heuristic indicator of success. This distance can be naturally determined via the Fubini-Study metric , in which the distance s (ϕ, ψ) between states |ϕ⟩ and |ψ⟩ is given by
The maximal Fubini-Study distance in the Hilbert space is π/2, the distance between mutually orthogonal states. We will therefore choose the quantity
as an ad hoc parameter, which characterizes the ability of an adiabatically evolving quantum device to reach its desired quantum state. The bar denotes the averaging over all initial states and over all quantum trajectories accessible to the system, which connect them to the final states maximally removed from them (in terms of the Fubini-Study distance).
2 Random Walk Model: Heuristic Treatment
The evolution of the state vector of a quantum system can be modeled by a series of random collapses to one of its instantaneous eigenstates at the moments t1, t2, …, and unitary evolutions under the Hamiltonian
For our purpose it is essential to keep the decoherence time as an explicit parameter. Therefore, the model we use is a random walk in the Hilbert space of the system, with the step (Fubini-Study) lengths Δsj dependent on Δt. For small time intervals of unitary evolution
Substituting here the maximal possible value of s = π/2, we find
This coincides with the rigorous Mandelstam-Tamm expression for the minimal time necessary to evolve from an initial state to a state orthogonal to it .
Now consider the random walk of M ≫ 1 steps of identical duration Δt, controlled by independently distributed random Hamiltonians
Here the summation is taken over the states |χ⟩ from the orthonormal basis of the Hilbert space of our system, which includes the state |ψ⟩.
Making a further simplification, assume that all random steps have the same length
where the bar average is over random choice of Hj-Hamiltonians and initial states |ψ⟩.
After averaging over random Hamiltonians
Here, f(N) is some function of the dimension of the Hilbert space D = 2N, which is given by
Then we obtain from Eq. 3
In particular, the condition
3 Random Walk Model: Numerical Approach
The goal of our numerical simulations is to determine the function f(N) in Eq. 13 through the relation between the Fubini-Study length of a single step, ΔS, and the number of steps, M, of a maximal random walk in the D = 2N-dimensional Hilbert space (see Figures 1A,B).
FIGURE 1. Random walk in the Hilbert space. (A) Schematics of states walking randomly in the Hilbert space. The points represent states. The state walks randomly from |Ψt=0⟩ to |Ψt=M⟩, where the stride of each step is given by the Fubini-Study distance ΔS. (B) Schematic of the random walk in the Hilbert space in the Bloch sphere. The red points indicate the states, which walks with the random direction with a fixed stride ΔS. The random walker on the Bloch sphere starts from the initial state |0⟩ at the north pole, and finally reaches the state |1⟩ at the south pole. (C) For considering the uniform random walk in the Hilbert space, we generate a random vector in the local orthogonal coordinate (spanned by the bases
A D-dimensional quantum state can be parameterized by 2 (D—1) real parameters
We consider the uniform random walk in the Hilbert space using the form Eq. 14, where each step length is fixed with the Fubini-Study distance ΔS. First, we randomly prepare an initial state |Ψt=0⟩ in the D-dimensional Hilbert space. The state |Ψt(θ)⟩ at a step t is parameterized by θ. The state at t + 1 is updated as |Ψt+1 (θ + dθ)⟩.
The Fubini-Study distance between these states |Ψ (θ + dθ)⟩ and |Ψ (θ)⟩ can be written as
where the Fubini-Study metric is given by
The span of the random walk from t = 0 to M is then given by
In our simulations we fix the number of steps M and search for such a value of ΔS (“critial value”) that the final state |Ψt=M⟩ satisfies the condition |π/2—s (Ψ0, ΨM)| ≤ ϵ for some small ϵ.
Figure 1D shows the critical value of ΔS as a function of M. Fitting yields
where A(N) = 0.309 (8)N0.76(2) and B(N) = 0.472 7 (8)N0.020(1) (see Figures 1E,F). In these numerical simulations, we used
4 Accessibility Index, “Quantumness” Criterion and Comparison to Experiment
The quantities ΔS and M do not have a direct experimental significance. In our approximate treatment we can relate them instead to the decoherence time (during which the unitary evolution takes place), tD, and the total time of adiabatic evolution, tf, via
where J is the typical coupling between qubits, so that J is a reasonable measure of the uncertainty of the N-qubit system’s energy during the adiabatic evolution (here we use the Mandelstam-Tamm expression). Then from Eq. 18 we find, that the necessary “quantumness” criterion for the adiabatic evolution is
where C ≈ 0.31. Given all the approximations we have made, we can as well take C = 1. Then the “quantumness” condition can be written as
where the accessibility index for a system of N qubits with average coupling strength J is
and the necessary “quantumness” condition was satisfied. This indicates that the results of Refs. [20, 22] are consistent with quantum annealing. From Eq. 21, we can evaluate the maximal size of a quantum processor for which the “quantumness” condition holds, other things being the same as in [20, 22]:
Note that the condition
We have investigated the generic behavior of a partially coherent system of qubits undergoing adiabatic evolution. Our aim was finding a convenient dimensionless parameter, which could characterize the degree of “quantumness” of our system. Basing our analysis on a random walk model of quantum evolution in the Hilbert space, we found a parameter allowing to evaluate the likelihood of a successful quantum transition between the initial and desired final states of the system. It is the accessibility index, expressed through the qubit decoherence time, time of evolution and the qubit coupling strength. Applying it to the case of 128-qubit D-Wave processors, we found that their evolution was consistent with quantum adiabatic transitions despite the qubit decoherence time being much smaller than the evolution time. An important future study is to assess the effectiveness of the accessibility index in more detail by considering a specific Hamiltonian system coupled to an environment.
Data Availability Statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
AZ: posing the problem; SW and AZ: developing the methodology and performing analytical calculations; SW and MS: numerical calculations; SW, SK, and AZ: interpretation and analysis on the results. All authors made intellectual contribution to the work, contributed to writing and editing of the manuscript and approved the submitted version.
This paper is partly based on results obtained from a project, JPNP16007, commissioned by the New Energy and Industrial Technology Development Organization (NEDO), Japan. SW was supported by Nanotech CUPAL, National Institute of Advanced Industrial Science and Technology (AIST). AZ was partially supported by NDIAS Residential Fellowship and by the SUPERGALAX project which has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 863313.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2021.773128/full#supplementary-material
2. Kakuyanagi K, Matsuzaki Y, Déprez C, Toida H, Semba K, Yamaguchi H, et al. Observation of Collective Coupling between an Engineered Ensemble of Macroscopic Artificial Atoms and a Superconducting Resonator. Phys Rev Lett (2016) 117:210503. doi:10.1103/physrevlett.117.210503
3.D‐Wave Systems Introduction to the D-Wave Quantum Hardware (2019) company manual. https://www.dwavesys.com/tutorials/background-reading-series/introduction-d-wave-quantum-hardware.
4. Bunyk PI, Hoskinson EM, Johnson MW, Tolkacheva E, Altomare F, Berkley AJ, et al. Architectural Considerations in the Design of a Superconducting Quantum Annealing Processor. IEEE Trans Appl Supercond (2014) 24:1–10. doi:10.1109/tasc.2014.2318294
14. Poulin D, Qarry A, Somma R, Verstraete F. Quantum Simulation of Time-dependent Hamiltonians and the Convenient Illusion of Hilbert Space. Phys Rev Lett (2011) 106:170501. doi:10.1103/physrevlett.106.170501
17. Andreev AV, Balanov AG, Fromhold TM, Greenaway MT, Hramov AE, Li W, et al. Emergence and Control of Complex Behaviors in Driven Systems of Interacting Qubits with Dissipation. Npj Quan Inf (2021) 7:1. doi:10.1038/s41534-020-00339-1
Keywords: quantumness, large system of qubits, quantum annealing (QA), quantumness criterion, random walk
Citation: Watabe S, Serikow MZ, Kawabata S and Zagoskin A (2022) Efficient Criteria of Quantumness for a Large System of Qubits. Front. Phys. 9:773128. doi: 10.3389/fphy.2021.773128
Received: 09 September 2021; Accepted: 22 November 2021;
Published: 05 January 2022.
Edited by:Alexandre M. Souza, Centro Brasileiro de Pesquisas Físicas, Brazil
Reviewed by:Rogério De Sousa, University of Victoria, Canada
Gustavo Rigolin, Federal University of São Carlos, Brazil
Copyright © 2022 Watabe, Serikow, Kawabata and Zagoskin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shohei Watabe, firstname.lastname@example.org