Next Article in Journal
A New Hyperchaotic 4D-FDHNN System with Four Positive Lyapunov Exponents and Its Application in Image Encryption
Next Article in Special Issue
Gaussian Amplitude Amplification for Quantum Pathfinding
Previous Article in Journal
State-of-the-Art Statistical Approaches for Estimating Flood Events
Previous Article in Special Issue
Quantum Simulation of Pseudo-Hermitian-φ-Symmetric Two-Level Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A NISQ Method to Simulate Hermitian Matrix Evolution

1
Peng Cheng Laboratory, Shenzhen 518055, China
2
Beijing Academy of Quantum Information Sciences, Beijing 100193, China
*
Authors to whom correspondence should be addressed.
Entropy 2022, 24(7), 899; https://doi.org/10.3390/e24070899
Submission received: 23 May 2022 / Revised: 21 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022
(This article belongs to the Special Issue Quantum Computing for Complex Dynamics)

Abstract

:
As a universal quantum computer requires millions of error-corrected qubits, one of the current goals is to exploit the power of noisy intermediate-scale quantum (NISQ) devices. Based on a NISQ module–layered circuit, we propose a heuristic protocol to simulate Hermitian matrix evolution, which is widely applied as the core for many quantum algorithms. The two embedded methods, with their own advantages, only require shallow circuits and basic quantum gates. Capable to being deployed in near future quantum devices, we hope it provides an experiment-friendly way, contributing to the exploitation of power of current devices.

1. Introduction

Building up a large-scale error-corrected quantum computer is to be one of the greatest scientific and engineering achievements [1,2,3,4]. However, stringent requirements such as millions of qubits with high accuracy are far to meet. Preskill coined “Noisy Intermediate-Scale Quantum” (NISQ) to describe this era, where NISQ devices represent the current state of the art in the fabrication of quantum devices [5]. The leading quantum computers contain up to a few hundred physical qubits, but provide rare practical applications as error correction is missing [6,7]. Therefore, while polishing the hardware-related techniques, one present-day goal is to exploit the power of current machines.
Matrix evolution is computationally hard in numerical mathematics, as O ( d 3 ) operations are required for an unstructured d × d matrix [8]. With advent of quantum algorithms, this can be solved to some extent by instinct priority of fault-tolerant quantum computation on matrix multiplications. For example, complexity of t-time analog Hamiltonian evolution is O ( t ) . Furthermore, certain digital algorithms were proposed for hermitian matrices’ evolution, which produce quantum speedups in many scenarios, such as simulation algorithms, quantum principal component analysis, quantum matrix inversion, and their generalizations [9,10,11,12,13,14,15]. However, deep quantum circuits and inaccessible oracles are required, which hinder their applications on near term devices. Accordingly, one question need be addressed: How to realize matrix evolution on current NISQ devices?
NISQ algorithms are a class of algorithms with no explicit requirements for error correction, promising to be deployed on NISQ hardware [5,16]. In regard to matrix evolution, we explore its near future application by introducing one typical NISQ module–layered circuit. Consequently, in Section 2, a heuristic protocol is proposed to employ a layered circuit to simulate hermitian matrix evolution, which is experimental friendly and can be employed in near-term applications of algorithms. To generate proper layered circuits, two methods are embedded in this protocol. The first is inspired by optimal control theory, which finds the simulating circuit directly, but with no scalability. The second, averting the scalability problem, generates the simulating circuit by a hybrid quantum-classical paradigm [17]. Both simulating circuits are with basic quantum gates and a pre-set depth, with the consumption of generating those circuits analyzed in Section 3. To support feasibility, simulating circuits are generated and validated numerically in Section 4, where hermitian matrices are set as density matrix of Bell state, GHZ state, and Hamiltonian of Crotonic acid molecular [18]. Compared with the method such as product formula and density matrix evolution, simulating circuits by our methods are with shallower depth, and more friendly to the experimental realizations. Furthermore, as a generalization of linear combination unitaries and layered circuit [19,20,21], ancillary layered circuit is proposed in our protocol, which not only serves as a quantum compiler here, but also an essential subroutine for NISQ algorithms.

2. Result

For hermitian matrices evolution, Lie–Trotter products and density matrix evolution provide solutions. We briefly review them here.
In the first case, for general Hamiltonian ρ = i ρ i with ρ i be local interactions, e i ρ t can be simulated by the lowest order Lie–Trotter product formula [22],
( i e i ρ i Δ t ) n e i ρ t + O ( | | ρ | | 2 t 2 / n ) ,
with Δ t = t / n , as shown in the schematic process in Figure 1a. n = O ( t 2 ϵ 1 | | ρ | | 2 ) denotes the repetition times, and therefore the circuits size and depth for full e i ρ t with desired accuracy ϵ . Though repeated applications of simulated circuit is feasible in theory, a low-depth quantum circuit is preferred by current NISQ devices as the limited coherence time. Therefore, under the circumstance of a tolerant accuracy, the implementation appears unfriendly to current or near term devices.
In the second case, ρ is not only Hermitian but also positive semi-definite and unit trace, that is, ρ is a physical quantum state, e i ρ t can be realized by multiple copies of ρ and infinitesimal swap operations [10]. Assuming that σ is another quantum state that e i ρ t act on, the infinitesimal swap operation has such effect,
tr p [ e i S Δ t ρ σ e i S Δ t ] = σ i Δ t ρ , σ + O ( Δ t 2 ) ,
where tr p is the partial trace over ρ and S is the swap operator. Shown in Figure 1b, density matrix evolution with respect to σ can be constructed by repeated applications of (2) with n = t / Δ t copies of ρ . Therefore, if the swap S and its infinitesimal exponential operation e i S Δ t can be implemented in a single layer circuit, both size and depth of using Equation (2) are O ( t 2 ϵ 1 ) . For accurate simulation, this depth and size of the circuit still dissatisfy the characteristics of current devices, and thus hinder its near term application.
One typical NISQ module is the layered circuit, which concretely implement near term applications [23,24]. Specifically, a m layered circuit U ( θ ) = i = 1 m U i ( θ i ) is presented in Figure 2a, which is parameterized by θ . As NISQ devices are with characteristics such as, limited size, short coherence time, and basic quantum operations, a shallow layered circuit seems perfect to undertake a NISQ applications. Remarkably, although shallow circuits and basic operations are utilized, the expressivity of layered circuits can be nontrivial and has been investigated in some recent papers [25,26]. Therefore, in our work, a layered circuit is employed to approach the simulation of matrix evolution, mathematically,
U m ( θ m ) U i ( θ 1 ) e i t ρ ,
where ρ is target hermitian matrix.
A layered circuit can be fully determined by its structure and parameters. For most applications, structure is configured previously, which depends on the tasks at hand. The widely used structures include quantum alternating operator ansatz, variational Hamiltonian ansatz and unitary coupled clustered ansatz [27,28,29,30]. In this study, according to hardware efficiency, we employed an m = O ( n ) layered circuit for following n-qubit tasks. All two body interactions are involved and gates in one layer are commuted with each others, which aims at employing the two-body interactions of the devices with a limited depth circuit. This circuit is problem-agnostic, which is given in Appendix B For specific problem faced, the circuit should be re-designed, which exploits both expressibility and trainability.
Furthermore, before stepping into parameter determination, ancillary layered circuit is introduced as preliminary. It is schematically presented in Figure 2b. Besides the principal system which is denoted as p, the ancillary register a is added with the same size as Chio-Jamiolkowski isomorphism is employed. The entire circuit includes three parts: encoding circuit ( U e , colored blue), layered circuit part ( U l I , colored pink), and decoding circuit ( U d , colored green). The details of dynamics in Figure 2b are shown as follows.
For initialization, two registers are jointly prepared on Ω 0 = | ψ 0 ψ 0 | , where
| ψ 0 = | 0 p | 0 a .
Encoding circuit U e is supposed to evolve the system into Ω 1 = | ψ 1 ψ 1 | , where | ψ 1 = i | i p | i a are pairs of Bell states. This step can be realized by a bunch of control-z gates, C z and Hadamard gates, H a ,
Ω 1 = U e Ω 0 U e , U e = I p H a · C z · H p H a .
Then, U l is applied on p, driving the system into Ω 2 , i.e., the Choi matrix of U l ,
Ω 2 = | ψ 2 ψ 2 | = i , j U l ( | i j | ) p U l | i j | a .
Noted that, arbitrary state ρ with the application of U l has a state-channel duality,
U l ( ρ ) = tr a [ Ω 2 · I ρ * ] .
Thus, knowing Ω 2 is sufficient to completely determine U l , i.e., converting quantum channel characterization to state characterization. Choi-Jamiolkowski isomorphism is exactly the correspondences in Equations (6) and (7).
Finally, Ω 3 is the result generated by steering Ω 2 through the decoding circuit U d = U e , which disentangles the system and is essential in our methods.
For parameter determination, it can be solved by two optimization methods, which minimize the distance,
U m ( θ m ) U i ( θ 1 ) e i t ρ ,
where ρ is the target hermitian matrix. Specifically, we substitute U i ( θ i ) with e i t θ i ρ i , where ρ i are tensor products of Pauli operators. In this configuration, the problem states as approaching e i t ρ by
e i t θ m ρ m e i t θ 1 ρ 1 .
The formalism of ρ i determines the structure of layered circuit. It is an empirical task and is set in advance, which we have depicted before.
From the point of view of state, for a specified time t, ρ drives an ideal arbitrary system, which is labelled as σ , to
σ 0 ( t ) = e i t ρ σ e i t ρ .
Simultaneously, with the application of Equation (9), the real system is steered to σ ( t ) , where
σ ( t ) = i = 1 m e i t θ i ρ i σ i = 1 m e i t θ i ρ i .
The objective function f is thus defined as the overlap which is measured by the standard inner product
f ( θ ) = tr ( σ 0 ( t ) σ ( t ) ) = tr ( μ k · ν k ) ,
where
μ k = i = k + 1 m e i t θ i ρ i e i t ρ σ e i t ρ i = m k + 1 e i t θ i ρ i , ν k = i = k 1 e i t θ i ρ i σ i = 1 k e i t θ i ρ i .
Obviously, f ( θ ) 1 , the maximum situation is satisfied when the outputs of i = 1 m e i t θ i ρ i and e i t ρ are the same. This is a state to state situation. To approach the dynamics, Choi matrix is employed, which can be produced by ancillary layered circuit. In this situation, the notations are redefined,
σ = U e Ω 0 U e , μ k = i = k + 1 m ( e i t θ i ρ i · e i t ρ ) I · σ · i = m k + 1 ( e i t ρ · e i t θ i ρ i ) I , ν k = i = k 1 e i t θ i ρ i I · σ · i = 1 k e i t θ i ρ i I .
Method 1 is a traditional method of gradient, where the partial derivative with respect to the parameter θ k are
f θ k = i t μ k | [ ρ k , ν k ] .
· means trace here. Therefore, by updating θ k by θ k with a learning rate η ,
θ k = θ k + η f θ k ,
the objective function would increase along the direction of gradient until converging into the local maximum. If f 1 is iteratively achieved, this layered circuit is said to approach the application of e i t ρ .
However, one vital problem is at calculating the gradient. As obtaining ν k and μ k is in general inefficient for tremendous time-consuming. Therefore, Method 2 is reported.
Method 2 is outlined as follows, (i) implementing e i ρ Δ t with a layered circuit, where Δ t is a sufficiently short period of time; (ii) implementing e i ρ t with (i) as the starting point step by step, where ancillary layered circuit is used as a compiler and O ( log ( t / Δ t ) ) steps are required. This idea is extended from trotter decomposition, which induces a trotter error with ϵ t O ( t Δ t ) . We illustrate the details as follows.
The first step is to realize an approximation,
i e i Δ t θ i ρ i e i Δ t ρ ,
which requires us realizing following equation if only lowest order is considered,
θ i ρ i = ρ .
This is coincident with Lie–Trotter decomposition, where sufficiently small Δ t is adopted. Related details for derivation can be found in Appendix A. Accordingly, assuming that Δ t is small enough, e i ρ Δ t can be simulated efficiently by a layered circuit, where structure and parameters are given by Equation (18).
The second step is to simulate e i ρ t . On basis of (i), this can be achieved by repeatedly applying ancillary layered circuit as compiler for O ( log ( t / Δ t ) ) times. We depict the i-th iteration in Figure 3, where U l consists of U i + 1 ( θ ) and U i 1 followed by n c times. The ancillary layered circuit serves as a compiler, learning repeated U i by U i + 1 ( θ ) . n c represents the compression efficiency, U i is known, and the structure of U i + 1 is set previously. Specifically, the task for i-th iteration is to generate U i + 1 ( θ ) with an explicitly prior known U i , where U i + 1 U i n c and θ are the parameters to be optimized.
The objective function is defined as,
f = t r ( Ω 1 Ω 2 ) = t r ( Ω 0 Ω 3 ) ,
which reaches the maximum at U i + 1 ( θ ) = U i n c . As Ω 1 are pairs of Bell states, this objective function cannot be measured with the local bases directly and efficiently. We investigate the second term in Equation (19): A decoding circuit U d disentangles the target component into Ω 3 . In this situation, only with the local measurement on Ω 0 , the objective function can be estimated and optimized as conventional hybrid quantum algorithms [31]. Iteratively, parameters in U i + 1 , which are obtained in the i-th iteration, could be applied in the ( i + 1 ) -th step, where U i + 2 = U i + 1 n c is realized.
Accordingly, by log ( t / δ ) steps, e i t ρ can be implemented based on e i Δ t ρ . Compared to repeatedly applications of U i which will cause the circuit depth affordable, this strategy trades the circuit depth with repeated applications of ancillary layered circuit(shown in Figure 3), which increases the training time while compresses the repeated circuit depth.

3. Analysis

To simulate hermitian matrix evolution with layered circuit, time complexity is bounded by the depth of used circuit, which is easily analyzed if circuits are determined. For errors introduced by simulation, it is the distance between target evolution and layered circuit, which is also generated by two embedded methods. Therefore, analysis for embedded methods are important: Learning complexities and learning errors are analyzed in this section.
We list procedures of two methods in the Method 1 and Method 2 and analyze the learning complexity first.
Method 1 relies on a classical optimization to generate parameter configurations. Therefore, complexities for time and memory is based on classical resources, such as classical logic gates and registers. During the procedure, approximating the gradient by Equation (15) costs the most. For each iteration, calculating μ k and ν k lies at the heart, which requires the implementation of matrix multiplications on a classical computer. Though e i t ρ i can be efficiently simulated individually for locality of ρ i , simulating e i t ρ is in general hard for most cases and there is no even universal efficient algorithms. To simulate an n-qubit problem, 2 n × 2 n matrix should be generated and stored, with O ( 4 n ) operations for evaluating the objective function and obtaining the gradient. Thus, the time complexity is O ( 4 n ) × r for a total r iterations while memory complexity is O ( 4 n ) . As Choi states are employed, our expedition will be enlarged, qubits are doubled.
Accordingly, this method is in general inefficient. However, it would be applicable to specific circumstances where ν k and μ k can be efficiently simulated. In fact, some investigations have employed tensor network to work out certain structure simulation problem and shed the light on the middle-scale quantum systems [32].
Method 2 is a hybrid quantum classical paradigm, where ancillary layered circuits are executed on quantum computers, parameters are updated on classical computers. For tasks on quantum computers, ancillary layered circuit brings a 2 n -qubit consumption on quantum register, where n is the size of simulation quantum system. Time complexity to execute ancillary layered circuit is additional O ( n ) Hadamard and control-z gates, with a depth of O ( m ) layered circuit consisting of at most O ( m × n ) basic quantum gates. Luckily, our measurement is on Ω 0 , which is local measurement. Therefore, for single iteration, updating all parameters requires O ( m × n ) times repeatedly applying and measure ancillary layered circuit. For r repetitions, O ( m × n ) × r is required. If e i t ρ is to be realized by e i Δ t ρ , an additional multiplied factor log ( t / Δ t ) should be added. For tasks on classical computers, a storage of O ( m × n ) parameters and their numerical gradient is required. Time complexity depends only on operations of fundamental arithmetic, which is within the reach of current machines.
Accordingly, time complexity of Method 2 is O ( m × n × r × log ( t / Δ t ) ) . As in our configuration, m = O ( n ) , which leads our method acceptable with respect to efficiency.
Method 1
  Input: layered circuit i = 1 e i t θ i ρ i . θ is to be optimized, ϵ o is an optimized threshold, and δ 1 is the tolerance for improving.
  Output: θ , parameter configuration, which is optimized to approach e i t ρ via i = 1 e i t θ i ρ i .
 
1: Evaluate the objective function Equation (12) with existed θ , denoted as f. If f 1 ϵ o , go to 2, otherwise, the algorithm terminates and θ return.
2: for  k = 1 , , m  do
3:   Based on Equation (14), ν k and μ k are calculated.
4:   Evaluate f / θ k by Equation (15).
5: end for
6: Update m-element θ according to Equation (16)
7: Evaluate the objective function with the new θ
8: if f 1 ϵ o then
9:   the algorithm terminates and θ return;
10: else if Δ f δ 1 then
11:    θ are re-initialized and go to 1.
12: else
13:   go to 2.
14: end if
15:
16:
17:
18:
Method 2
  Input: layered circuit i = 1 e i θ i ρ i with known U 1 . θ is to be optimized, ϵ o is the optimized threshold, δ 2 is the tolerance for improving, and n c is the compressing factor.
  Output: θ , parameter configuration, which is optimized to approach e i t ρ via i = 1 e i t θ i ρ i
 
1: for  i = 1 , , log ( t / Δ t )  do
2:   if i = 1 then
3:     U 1 and its inversion are generated by the input.
4:   else
5:     U i are generated as the output of the last iteration.
6:   end if
7:   Evaluate the objective function in Equation (19) with existed θ and n c , which is denoted as f. If f 1 ϵ o , go to 8, otherwise, θ for U i + 1 return and go to 1.
8:    θ for U i + 1 are optimized as variant quantum algorithms, θ are updated and f is evaluated.
9:   if f 1 ϵ o then
10:     θ return and go to 1;
11:   else if Δ f δ 2 then
12:     θ are re-initialized and go to 7.
13:   else
14:    go to 8.
15:   end if
16: end for
To analyze errors of layered circuit by both methods, first, we define some notations: ϵ o is the optimized threshold which is supposed to terminate the training process; ϵ t is the deviations coming from lie-product decomposition.
If Method 1 is completed, the optimized threshold is targeted. The obtained layered circuit has an accuracy of 1 ϵ o as training process permits the error no more than ϵ o . For Method 2, the error accumulates according to a chain rule when implementing U i + 1 by U i . It ends up with an error of O ( n c s ϵ o ) , where the higher order terms are ignored and s is steps, which is in total of O ( log ( t / Δ t ) ) to realizing e i ρ t from e i ρ Δ t . Thus, the discrepancy is O ( ( t / Δ t ) ϵ o ) . Additionally, an error coming from Equation (1) is also considered as ϵ t , which is of O ( t Δ t ) . Accordingly, the total error for Method 2 is O ( ( t 2 / ϵ t ) ϵ o ) + O ( ϵ t ) .
In the end of this section, we analyze the expressive power for typical layered circuit, which is given by our appendix method. In fact, it is problem-dependent and can be replaced with smarter ansatz with tasks at hand.
Expressibility is proposed recently as a distance, which measures the output states distribution of layered circuits and the Haar [26]. To explicitly present this value, Kullback-Leibler (KL) divergence(or relative entropy) is employed to estimate this distance, which is denoted as Expr. A highly expressible circuit would produce a small Kullback-Leibler value. In this part, besides numerically calculating E x p r of circuit given in Appendix C, three other types of parametrized circuits are also studied as comparisons, which are also specified in Appendix C. Three compared circuits are repeated with up to 5 times to calculate E x p r , while our typical circuit stays unchanged. Figure 4 shows the results of Expressibility values (or KL divergences), where circle, diamond and square represents three comparison circuits and cyan cross labels the circuit given by us. The layered circuit in this work has a similar performance as multi-applications of circuit 2 and 3 in Figure A1. Remarkably, in general, repeating a circuit layer would increase the expressive power. However, as no entanglement gate exists in circuit 1 in Figure A1, this argument does not hold for that circuit. More information on expressibility and its simulation can be found in Appendix C or related work [26].

4. Numerical Experiments

In this section, numerical experiments are investigated, including applications on Hamiltonian simulation and density matrix evolution. In our protocol, ρ is taken as (a) Bell state; (b) GHZ state and (c) Hamiltonian of a liquid NMR sample, Crotonic acid, which is specified in the Appendix E.
To generate the simulating layered circuit, two learning methods are employed, with
f ( θ ) = t r ( ρ t Ω 3 )
being the objective function. θ = α , β are for Method 1 and Method 2. The target state can be expressed as,
ρ t = U d · e i ρ t I a · U e Ω 0   ( M e t h o d   1 ) , ρ t = Ω 0     ( M e t h o d   2 ) .
And additionally, Ω 3 = U d · U I a · U e Ω 0 , where
U = U 1 ( α ) = i = 1 m e i t α i ρ i   ( M e t h o d   1 ) , U = U 2 ( β ) = U i n c U i + 1 ( β )   ( M e t h o d   2 ) .
This is specified in Appendix D.
To find optimal parameter configuration, gradient-based method is employed. To begin with, random numbers by a single uniformly distribution are generated and assigned as the value of initial parameters. During the training process, the approximation of the partial derivatives can be estimated by following symmetric difference quotient
j f = f ( θ + Δ θ j ) f ( θ Δ θ j ) 2 Δ θ j ,
which can be realized by repeatedly running the training circuit with a small perturbation on the θ , where Δ θ j is naively chosen as 0.01 in our simulation. Therefore, the optimal configuration of the parameters can be obtained by repeating
θ = θ + η f ( θ ) ,
where f ( x ) = ( 1 f , 2 f , ) and the learning rate, η , is fixed at 0.02 . For simulating Method 1 and Method 2, the maximum numbers for repeating Equation (22), i.e., iterations are set as 300 and 350, respectively.
In this prototypical simulation, the program is conducted with assumptions that all one or two quantum gates are accurate. Three circuits with the depths of 3, 4 and 5 are employed for the simulation the evolution. The numbers of parameters to be optimized are thus 12, 24 and 40, which converge to O ( m × n ) , where m is the depth and n is the size of circuit. As the comparison, if ϵ is tolerant error, O ( 3 / ϵ ) , O ( 4 / ϵ ) and O ( 5 / ϵ ) depth circuits are required for simulation by lie-product decomposition. For density matrix evolution, O ( 1 / ϵ ) copies of density matrix are required as with O ( 1 / ϵ ) -depth circuit.
For data collection, strategies are different for Method 1 and Method 2. For Method 1, once the optimization completed, parameters is supposed to be optimal for simulating a t-time interval evolution. Otherwise, we need to re-run the simulation and repeat Equation (22) until the objective function satisfying the optimized threshold. For Method 2, optimizing ancillary layered circuit is not once for all. After the i-th repetition of training ancillary layered circuit with Equation (22) for 350 iterations, only U i n c U i + 1 ( β ) is learned. Therefore, O ( log ( t / Δ t ) ) repetitions are required. In our simulation, Δ t / t = 1 / 2 10 , i.e., 10 steps are required and the base number, 2, determines the compression efficiency of the circuit.
Figure 5 shows the results of our simulation, which provide supports on our protocol by converging to the target matrix evolutions. (a) (b) and (c) simulate the evolutions of Bell state, GHZ state and Hamiltonian of Crotonic acid, where the evolution timescales are set as 0.05 0.1 and 0.2 , respectively. Fidelities are shown by vertical axis and calculated by tracing the inner product of two operators. The results show that, with the accurate local quantum gate assumption, the simulating circuit gets more and more like the evolution of certain matrix evolution as with the number of iterations.

5. Limitation

In the end, defects in our protocol should be listed. First, only Hermitian operations can be dealt with. Evolution of an open quantum system is not considered in this article. With formalism for the description of open quantum systems, this target is potentially solved in future work [33]. Second, similar to the widely-used parameterized circuits architecture such as QAOA and VQE, expressive ability of ansatz in our methods determinate the upper-bound of accuracy [25,28]. Therefore, a proper ansatz will have a good performance. If the termination threshold cannot be reached, an alternative ansatz or iterative training algorithms should be resorted to. Third, optimization should be clarified, which is also the most important one. Many optimization problems are in fact NP-hard problems. The methods embedded are essentially optimization-based, and cannot get over them, too. Taking an example of gradient-based methods, they cannot avoid the local optimum problems, especially when the feasible region is complicated. Therefore, the proposed methods would fail in finding a global optimum with a bad initial guess, just as the classical cases. Although numerical results are good without considering the initialization, we have to admit a good initial guess is crucial to decrease the possibility of failure, especially for dealing with a larger problem. Investigations on the optimization method itself should attract more attention, which benefits not only the exploitation of near-future quantum devices but also for most modern technologies [34,35].

6. Conclusions

NISQ era has come in and would last for decades. Thus, finding a near-future algorithm which exploits the power of NISQ devices becomes more and more important. Simulating matrix evolution, which is an essential module for many quantum algorithms, such as density matrix evolution for quantum principal component analysis, and static or dynamic simulation for quantum systems. The existed implementation still relies on a deep quantum circuit, which is intractable with current NISQ devices.
In this paper, a heuristic layered circuit protocol is proposed for simulating e i ρ t with Hermitian ρ . To construct these circuits, two methods are given. Method 1 is with the classical optimal control theory and Method 2 is with the hybrid quantum-classical paradigm. For Method 1, learning a layered circuit requires O ( 4 n × r ) operations, where n is the size of quantum circuit and r is the number of total iterations. As computational resource may be put into its limits, we provide Method 2, which is hybridized with a quantum agent. As a comparison, time complexity is O ( m × n × r × log ( t / Δ t ) ) , with 2 n being the size of circuit and m being circuit depth. From the point of view of error, although Method 1 is in general inefficient, the error can be bounded to O ( ϵ o ) . Method 2 introduces a larger one, O ( ( t 2 / ϵ t ) ϵ o ) + O ( ϵ t ) . Only when ϵ o is sufficiently small, Method 2 can approach an accuracy as sTrotter method. Accordingly, circuits by both methods would be shallow, which is easily deployed on current devices.
Simulating matrix evolution is important for many quantum information tasks. Our protocol, which realizes hermitian matrix evolutions, is supposed to contribute to the field of NISQ algorithm. For layered circuit in our both methods, only basic quantum gates are required, which means at most 2-qubit interactions are needed. With respect to experimental technology, it is already mature for most physical platforms, such as optical lattice, spin-based, and superconducting qubit [7,36,37]. In despite of some limitations, the protocol is with an affordable computational consumption, which paves a way for possibly applications on NISQ devices. In addition, ancillary layered circuit, which serves as a compiler in this article, is promising to be an important subroutine in near future and provides a new way to exploit the power of NISQ devices.

Author Contributions

K.L. and P.G. contribute equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major Key Project OF PCL and National Natural Science Foundation of China under grant No. 11905111.

Data Availability Statement

All data for the figure and table are available on request. All other data about experiments are available upon reasonable request.

Conflicts of Interest

The authors declare no competing interests.

Abbreviations

The following abbreviations are used in this manuscript:
NISQNoisy Intermediate-Scale Quantum

Appendix A. Derivations on e i Δ t ρ in Method 2

First, if one wants to simulate e i Δ t ρ with i e i Δ t β i ρ i , the idea of optimizations can always be resorted to. An objective function can thus be defined as,
f ( β ) = t r ( e i Δ t ρ σ e i Δ t ρ i e i Δ t β i ρ i σ i e i Δ t β i ρ i ) ,
where σ is an arbitrary initial density matrix, and operators in · satisfy the permutation equality. We expand matrix exponentiation in Equation (A1) with Taylor series on Δ t ,
e i Δ t ρ = k = 0 ( i Δ t ρ ) k / k ! , i = 1 N e i Δ t β i ρ i = i = 1 N k = 0 ( i Δ t β i ρ i ) k / k ! .
With substitutions such as, H i Δ t ρ and M i i Δ t β i ρ i , ingredients in Equation (A1) can be rewritten as,
e i Δ t ρ × i = N 1 e i Δ t β i ρ i = ( H ) 0 0 ! + ( H ) 2 1 ! + ( H ) 3 3 ! + × M N 0 0 ! + M N 1 1 ! + M N 2 2 ! + M 1 0 0 ! + M 1 1 1 ! + M 1 2 2 ! + = k = 0 j i , j ρ 0 , i j i + j ρ = k ( H ) j ρ j ρ ! i = N 1 M i j i j i ! , i = 1 N e i Δ t β i ρ i × e i Δ t ρ = ( M 1 ) 0 0 ! + ( M 1 ) 1 1 ! + ( M 1 ) 2 2 ! + ( M N ) 0 0 ! + ( M N ) 1 1 ! + ( M N ) 2 2 ! + × H 0 0 ! + H 2 1 ! + H 3 3 ! + = k = 0 j i , j ρ 0 , i j i + j ρ = k i = 1 N ( M i ) j i j i ! H j ρ j ρ ! .
Under the circumstances that no approximation is adopted, the objective function is
f = ψ | k = 0 j i , j ρ 0 , i j i + j ρ = k i = 1 N ( M i ) j i j i ! H j ρ j ρ ! | ψ ψ | k = 0 j i , j ρ 0 , i j i + j ρ = k ( H ) j ρ j ρ ! i = N 1 M i j i j i ! | ψ ,
where σ = | ψ ψ | and · represents the expectation on | ψ . Obviously, it is a taylor-like polynomial with respect to Δ t . As Δ t can be a sufficiently small quantity when we simulate e i t ρ , only lowest order terms on Δ t are considered. For the objective function f, all terms which is linear with or higher than O ( Δ t 2 ) are ignored. Therefore, f can be simplified as,
( 1 + < H 2 2 ! > < H i M i > + < i M i 2 2 ! > + < i > j M i M j > ) 2 ( < H > < M i > ) 2 ,
where M i = H is the requirement for eliminating above lowest order error, which is also independent of | ψ . In the original representation, β i satisfy following correspondence
β i ρ i = ρ .
which is also the correspondence by Lie–Trotter products. Therefore, a parameterized quantum circuit for e i Δ t ρ can be obtained by Equation (A5), where ρ i is comparably easy-access to realize and benchmark in near-future quantum devices.

Appendix B. Method to Generate an N -Layer Circuit

In this section, a method to generate an N-layer circuit structure with all two-body interactions is introduced, where N is the size of the system. P 2 ( i , j ) are assumed as the products of Pauli matrices, where P 2 ( i , j ) = σ k σ l are on both i-th and j-th qubits and i , j = 1 , , N and k , l = x , y , z . All of them can be presented in an N × N triangle parameter matrix P,   
P = 1 P 2 ( 1 , 2 ) P 2 ( 1 , 3 ) P 2 ( 1 , 4 ) P 2 ( 1 , 5 ) , , P 2 ( 1 , N ) 2 P 2 ( 2 , 3 ) P 2 ( 2 , 4 ) P 2 ( 2 , 5 ) , , P 2 ( 2 , N ) 3 P 2 ( 3 , 4 ) P 2 ( 3 , 5 ) , , P 2 ( 3 , N ) 4 P 2 ( 4 , 5 ) , , P 2 ( 3 , N ) , , P 2 ( N 1 , N ) N ,
where P i i = i , P i j = P 2 ( i , j ) for i < j and P i j = 0 for i > j . Thus, the non-zeros elements in P can be divided into N group with the following Method A1, where N is supposed to be odd.
Method A1 Algorithm on grouping the Pauli words
  • Input: Parameter matrix P, N null sets, and a N × N label matrix L, where L m n = 1 for arbitrary m , n .
  • Output:N sets of Pauli words: S 1 , S 2 , S N
  • for  i = 1 , N  do
  •     for  j = i , N  do
  •         if  i = j L i j = 1  then
  •             P i j S i , L i j = 0 .
  •         else
  •            for  m = 1 , N  do
  •                if  L m i = 1 L m j = 1  then
  •                     P i j S m , L m i = 0 , L m j = 0
  •                end if
  •            end for
  •         end if
  •     end for
  • end for
With the complexity of O ( N 3 ) operations, Those Pauli words are divided into N group, where the Pauli words in each group commute with each other. Via matrix exponentiation formation multiplied by tunable parameters, they can be arranged into a parameterized circuit, with a depth of N, being linear with the size of the system. If different type Pauli combinations are considered, O ( N ) is the depth. Therefore, this algorithm provided a depth of O ( N ) parameterized circuit with all two-body interactions considered. Remarkably, this is one type of problem-agnostic circuit, which exploit the power of current quantum devices with two-body interactions available. For different problems, problem-inspired circuit can be designed, employing information about the problem, which bring better expressibility and trainability. Surely, we will study the expressibility of our proposed circuit in next section.

Appendix C. Expressibility of the Layered Circuit

A layered circuit with excellent expression is more likely to represent a target state. To measure the expressive power of a circuit, expressibility is defined in recent work as a distance of two state distributions, characterized with | | · | | , the Hilbert-Schmidt norm [26].
A = | | H a a r ( | ψ ψ | ) t d ψ θ ( | ψ θ ψ θ | ) t d ψ θ | | ,
where | ψ θ are outputs of a layered circuit with randomized configuration θ and | ψ are from a distribution according to the Haar measure. It is a generalization from the investigation of pseudo-random circuit. Thus, a highly expressible circuit would produce a small A, with A = 0 corresponding to being maximally expressive, i.e., generating a state distribution to the Haar measure.
In order to explicitly estimate the expressibility with a discrete simulated result, the Kullback-Leibler (KL) divergence, i.e., relative entropy is employed, which measures the difference between one probability distribution and a reference probability distribution. It is denoted as E x p r ,
E x p r = D K L ( P l ( f | θ , θ ) | P H a a r ( f ) ) ,
where P l ( f | θ , θ ) is the probability distribution of f= | ψ θ | ψ θ | 2 , outputs of a layered circuit with random θ and θ . For the reference probability distribution, i.e., to Haar measure, P H a a r ( f ) = ( N 1 ) ( 1 f ) N 2 is analytical, with N being size of dimension. E x p r scores the similarity to the distribution of Haar as a layered circuit with a lower value of KL divergence is a more expressible circuit.
Figure A1. A set of circuit templates considered in the study, each labeled with a circuit ID. The dashed box indicates a single circuit layer, denoted by L in the text, that can be repeated.
Figure A1. A set of circuit templates considered in the study, each labeled with a circuit ID. The dashed box indicates a single circuit layer, denoted by L in the text, that can be repeated.
Entropy 24 00899 g0a1
Accordingly, numerical experiments to compute and compare the expressibility in the strategy with three other typical layered circuits are studied, which are shown in Figure A1. They are all 4 qubit circuit. As with the original work, to construct a distribution of P l with histogram, a bin size is set as 75, and 5000 f are measured. For each type of layered circuits except for the one in our strategy, the instances where the number of layers are investigated with up to 5. Compared with 10,000 times repeatedly running the circuits and state tomography, f can be measured directly with the setup proposed in Figure A2 for 5000 time and without tomography. Only 1-qubit is added as the ancillary system A, at the end of circuit, the system is evolved as
ρ = | 0 0 | | ψ θ ψ θ | + | 1 0 | | ψ θ ψ θ | + | 0 1 | | ψ θ ψ θ | + | 1 1 | | ψ θ ψ θ | .
Figure A2. Circuit for measurement of the expressivity of the layered circuit. It costs one extra qubit as ancillary. θ and θ are two random parameter configurations for layered circuits U θ and U θ , H is a one-qubit Hadamard gate.
Figure A2. Circuit for measurement of the expressivity of the layered circuit. It costs one extra qubit as ancillary. θ and θ are two random parameter configurations for layered circuits U θ and U θ , H is a one-qubit Hadamard gate.
Entropy 24 00899 g0a2
Via the measurements on ancillary system, f can be obtained with the real and imaginary part of ψ θ | ψ θ , which are extracted as
Re = t r ( σ x I · ρ ) = 1 2 ( ψ θ | ψ θ + ψ θ | ψ θ ) Im = t r ( σ y I · ρ ) = i 2 ( ψ θ | ψ θ ψ θ | ψ θ ) .
Results are shown in the Figure 4 in main text. For three types of layered circuits which are denoted in circle, diamond and square, the instances are considered where the circuit, which is as a unit layer, are repeated with up to 5. For the circuit proposed by Method A1, the repetition is always 1.
In general, repeating a circuit layer in Figure A1 would increase the expressibility. As there is no entanglement gate in circuit 1, this argument does not hold. For circuit 2 and 3, even if a single layer for circuit 3 have pretty good expression, this argument holds. From the result, there are convergences for the metric of expressibility. For the argument whether more numbers of layers can help the expression, at least, we can not obtain more information from this metric. Additionally, the layered circuit generated in our strategy, which is labeled as cyan cross, keeps a same structure. It has a similar performance with respect to the repeated circuits 2 and 3. As a small value of expressibility imply an excellent expression, the generated circuits in the strategy performs no worse than other typical circuit.

Appendix D. Objective Functions

In this section, the explicit expression of the objective functions is derived, as well as the step-by-step illustrations of the circuit to implement the Numerical simulation. The objective function can be concluded as
f ( θ ) = t r ( ρ t Ω 3 ) ,
where θ = α , β are for Method 1 and Method 2, respectively. As for intermediate states and quantum operations of the circuit shown in original manuscript, they are depicted as,
Ω 0 = | ψ 0 ψ 0 | ,   | ψ 0 = | 0 p | 0 a , Ω 1 = | ψ 1 ψ 1 | ,   | ψ 1 = U e | ψ 0 , Ω 2 = | ψ 2 ψ 2 | ,   | ψ 2 = U l I a | ψ 1 , Ω 3 = | ψ 3 ψ 3 | ,   | ψ 3 = U d | ψ 2 .
where U e = I p H a · C z · H p H a , U d = H p H a · C z · I p H a . In addition, | ψ 1 = i = 1 d | i p | i a is the bell state pairs, where the idea of Choi matrix for assisting tomography is employed.
With respect to the U l in Method 1, it aims at finding optimal parameters α = ( α 1 , α 2 , ) to generate a layered circuit. While, for the i-th repetition of Method 2, it aims at finding appropriate parameters β = ( β 1 , β 2 , ) , generating a layered circuit which realizes repeatedly applied U i . Thus, the different U l are to be optimized
U 1 ( α ) = i = 1 m e i t α i ρ i   ( M e t h o d   1 ) , U 2 ( β ) = U i 2 U i + 1 ( β )   ( M e t h o d   2 ) .
As for target states, they are also different. Method 1 realizes an evolution by a target Hermitian matrix on system p while Method 2 realizes I p . Therefore, the density matrix ρ t for both strategies are
ρ t = U d · e i ρ t I a · U e Ω 0   ( M e t h o d   1 ) , ρ t = Ω 0     ( M e t h o d   2 ) ,
where ρ can be density matrix or Hamiltonian, for simulating hermitian matrix exponentiation.

Appendix E. Matrices for Numerical Simulation

In this section, the detailed information for ρ in numerical simulation are given. Three cases are studied. Two of them are density matrices and one is a Hamiltonian from previous experiments. The specified form are as follows.
a. Bell State
| ψ = 1 2 ( | 00 + | 11 ) ,   ρ = 1 2 1 0 0 1 0 0 0 0 0 0 0 0 1 0 0 1 .
b. GHZ State
| ψ = 1 2 ( | 000 + | 111 ) ,   ρ = 1 2 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 .
c. Hamiltonian of Crotonic Acid
The Hamiltonian simulated in our numerical experiments is from 13 C-labeled Crotonic acid dissolved in d6-acetone, which is usually employed as a four-qubit quantum system for NMR-based quantum information processing [18].
The structure of this molecule is shown in Figure A3, where C 1 to C 4 are four carbon atoms, M is one of three hydrogen atoms in a methyl group, H 1 and H 2 are two hydrogen atoms. When it is used as a four qubit quantum system, the hydrogen atoms are decoupled from the carbon atoms with a shaped radio frequency pulse. Under the circumstance of weak coupling, the internal Hamiltonian of this molecule can be expressed as,
H i n t = j = 1 4 1 2 ν j σ z j + j < k 4 π 2 J j k σ z j σ z k .
ν j and J j k are inside Hamiltonian parameters whose values are listed in the table besides the molecule. ν j , the chemical shift, are the diagonal elements and J j k , the J-coupling, are the off-diagonal elements. As for T 2 , the transverse relaxation time, is not involved in our simulation as the unitary process is assumed. All parameters can be found in related paper which is under a magnetic field of 9.4 T at room temperature (296.5 K).
Figure A3. Molecular structure of 13 C-labeled Crotonic acid. ν j and J j k are the chemical shifts and J-couplings, respectively, which are listed by the diagonal and off-diagonal elements. T 2 (in Seconds) are the relaxation time which are shown at bottom.
Figure A3. Molecular structure of 13 C-labeled Crotonic acid. ν j and J j k are the chemical shifts and J-couplings, respectively, which are listed by the diagonal and off-diagonal elements. T 2 (in Seconds) are the relaxation time which are shown at bottom.
Entropy 24 00899 g0a3

References

  1. Shor, P.W. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. SIAM Rev. 1999, 41, 303–332. [Google Scholar] [CrossRef]
  2. Grover, L.K. Quantum mechanics helps in searching for a needle in a haystack. Phys. Rev. Lett. 1997, 79, 325. [Google Scholar] [CrossRef] [Green Version]
  3. Boixo, S.; Rønnow, T.F.; Isakov, S.V.; Wang, Z.; Wecker, D.; Lidar, D.A.; Martinis, J.M.; Troyer, M. Evidence for quantum annealing with more than one hundred qubits. Nat. Phys. 2014, 10, 218–224. [Google Scholar] [CrossRef] [Green Version]
  4. Biamonte, J.; Wittek, P.; Pancotti, N.; Rebentrost, P.; Wiebe, N.; Lloyd, S. Quantum machine learning. Nature 2017, 549, 195–202. [Google Scholar] [CrossRef] [PubMed]
  5. Preskill, J. Quantum Computing in the NISQ era and beyond. Quantum 2018, 2, 79. [Google Scholar] [CrossRef]
  6. Arute, F.; Arya, K.; Babbush, R.; Bacon, D.; Bardin, J.C.; Barends, R.; Biswas, R.; Boixo, S.; Brandao, F.G.; Buell, D.A. Quantum supremacy using a programmable superconducting processor. Nature 2019, 574, 505–510. [Google Scholar] [CrossRef] [Green Version]
  7. Zhong, H.S.; Wang, H.; Deng, Y.H.; Chen, M.C.; Peng, L.C.; Luo, Y.H.; Qin, J.; Wu, D.; Ding, X.; Hu, Y.; et al. Quantum computational advantage using photons. Science 2020, 370, 1460–1463. [Google Scholar] [CrossRef]
  8. Kressner, D.; Luce, R. Fast computation of the matrix exponential for a Toeplitz matrix. SIAM J. Matrix Anal. Appl. 2018, 39, 23–47. [Google Scholar] [CrossRef] [Green Version]
  9. Low, G.H.; Chuang, I.L. Hamiltonian simulation by qubitization. Quantum 2019, 3, 163. [Google Scholar] [CrossRef] [Green Version]
  10. Lloyd, S.; Mohseni, M.; Rebentrost, P. Quantum principal component analysis. Nat. Phys. 2014, 10, 631–633. [Google Scholar] [CrossRef] [Green Version]
  11. Rebentrost, P.; Schuld, M.; Wossnig, L.; Petruccione, F.; Lloyd, S. Quantum gradient descent and Newton’s method for constrained polynomial optimization. New J. Phys. 2019, 21, 073023. [Google Scholar] [CrossRef]
  12. Gao, P.; Li, K.; Wei, S.; Gao, J.; Long, G. Quantum gradient algorithm for general polynomials. Phys. Rev. A 2021, 103, 042403. [Google Scholar] [CrossRef]
  13. Gao, P.; Li, K.; Wei, S.; Long, G.L. Quantum second-order optimization algorithm for general polynomials. Sci. China Phys. Mech. Astron. 2021, 64, 100311. [Google Scholar] [CrossRef]
  14. Harrow, A.W.; Hassidim, A.; Lloyd, S. Quantum algorithm for linear systems of equations. Phys. Rev. Lett. 2009, 103, 150502. [Google Scholar] [CrossRef] [PubMed]
  15. Rebentrost, P.; Mohseni, M.; Lloyd, S. Quantum support vector machine for big data classification. Phys. Rev. Lett. 2014, 113, 130503. [Google Scholar] [CrossRef]
  16. Bharti, K.; Cervera-Lierta, A.; Kyaw, T.H.; Haug, T.; Alperin-Lea, S.; Anand, A.; Degroote, M.; Heimonen, H.; Kottmann, J.S.; Menke, T.; et al. Noisy intermediate-scale quantum (NISQ) algorithms. arXiv 2021, arXiv:2101.08448. [Google Scholar]
  17. McClean, J.R.; Romero, J.; Babbush, R.; Aspuru-Guzik, A. The theory of variational hybrid quantum-classical algorithms. New J. Phys. 2016, 18, 023023. [Google Scholar] [CrossRef]
  18. Li, K.; Wei, S.; Gao, P.; Zhang, F.; Zhou, Z.; Xin, T.; Wang, X.; Rebentrost, P.; Long, G. Optimizing a polynomial function on a quantum processor. NPJ Quant. Inf. 2021, 7, 16. [Google Scholar] [CrossRef]
  19. Gui-Lu, L. General quantum interference principle and duality computer. Commun. Theor. Phys. 2006, 45, 825. [Google Scholar] [CrossRef] [Green Version]
  20. Long, G.; Liu, Y. Duality quantum computing. Front. Comput. Sci. China 2008, 2, 167. [Google Scholar] [CrossRef]
  21. Childs, A.M.; Wiebe, N. Hamiltonian simulation using linear combinations of unitary operations. arXiv 2012, arXiv:1202.5822. [Google Scholar] [CrossRef]
  22. Childs, A.M.; Su, Y.; Tran, M.C.; Wiebe, N.; Zhu, S. Theory of trotter error with commutator scaling. Phys. Rev. X 2021, 11, 011020. [Google Scholar] [CrossRef]
  23. Farhi, E.; Goldstone, J.; Gutmann, S. A quantum approximate optimization algorithm. arXiv 2014, arXiv:1411.4028. [Google Scholar]
  24. Peruzzo, A.; McClean, J.; Shadbolt, P.; Yung, M.H.; Zhou, X.Q.; Love, P.J.; Aspuru-Guzik, A.; O’brien, J.L. A variational eigenvalue solver on a photonic quantum processor. Nat. Commun. 2014, 5, 4213. [Google Scholar] [CrossRef] [Green Version]
  25. Nakaji, K.; Yamamoto, N. Expressibility of the alternating layered ansatz for quantum computation. arXiv 2020, arXiv:2005.12537. [Google Scholar] [CrossRef]
  26. Sim, S.; Johnson, P.D.; Aspuru-Guzik, A. Expressibility and Entangling Capability of Parameterized Quantum Circuits for Hybrid Quantum-Classical Algorithms. Adv. Quant. Technol. 2019, 2, 1900070. [Google Scholar] [CrossRef] [Green Version]
  27. Romero, J.; Babbush, R.; McClean, J.R.; Hempel, C.; Love, P.J.; Aspuru-Guzik, A. Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz. Quant. Sci. Technol. 2018, 4, 014008. [Google Scholar] [CrossRef] [Green Version]
  28. Hadfield, S.; Wang, Z.; O’Gorman, B.; Rieffel, E.G.; Venturelli, D.; Biswas, R. From the quantum approximate optimization algorithm to a quantum alternating operator ansatz. Algorithms 2019, 12, 34. [Google Scholar] [CrossRef] [Green Version]
  29. Wiersema, R.; Zhou, C.; de Sereville, Y.; Carrasquilla, J.F.; Kim, Y.B.; Yuen, H. Exploring entanglement and optimization within the hamiltonian variational ansatz. PRX Quant. 2020, 1, 020319. [Google Scholar] [CrossRef]
  30. Anand, A.; Schleich, P.; Alperin-Lea, S.; Jensen, P.W.; Sim, S.; Díaz-Tinoco, M.; Kottmann, J.S.; Degroote, M.; Izmaylov, A.F.; Aspuru-Guzik, A. A quantum computing view on unitary coupled cluster theory. Chem. Soc. Rev. 2022, 51, 1659–1684. [Google Scholar] [CrossRef]
  31. Lavrijsen, W.; Tudor, A.; Müller, J.; Iancu, C.; de Jong, W. Classical optimizers for noisy intermediate-scale quantum devices. In Proceedings of the 2020 IEEE International Conference on Quantum Computing and Engineering (QCE), Denver, CO, USA, 12–16 October 2020; pp. 267–277. [Google Scholar]
  32. Pan, F.; Zhang, P. Simulating the Sycamore quantum supremacy circuits. arXiv 2021, arXiv:2103.03074. [Google Scholar]
  33. Eleuch, H.; Rotter, I. Nearby states in non-Hermitian quantum systems I: Two states. Eur. Phys. J. D 2015, 69, 229. [Google Scholar] [CrossRef] [Green Version]
  34. Ruder, S. An overview of gradient descent optimization algorithms. arXiv 2016, arXiv:1609.04747. [Google Scholar]
  35. Gill, P.E.; Murray, W.; Wright, M.H. Practical Optimization; SIAM: Philadelphia, PA, USA, 2019. [Google Scholar]
  36. Homid, A.; Abdel-Aty, M.; Qasymeh, M.; Eleuch, H. Efficient quantum gates and algorithms in an engineered optical lattice. Sci. Rep. 2021, 11, 15402. [Google Scholar] [CrossRef]
  37. Vandersypen, L.M.; Chuang, I.L. NMR techniques for quantum control and computation. Rev. Mod. Phys. 2005, 76, 1037. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Conventional method to realize e i ρ t , where ρ is a hermitian matrix. (a) is via decomposition from Lie–Trotter products and (b) is the density matrix evolution using infinitesimal swaps e i S Δ t . Wherein, /denotes a register of multi qubits, and|means tracing out the corresponding register.
Figure 1. Conventional method to realize e i ρ t , where ρ is a hermitian matrix. (a) is via decomposition from Lie–Trotter products and (b) is the density matrix evolution using infinitesimal swaps e i S Δ t . Wherein, /denotes a register of multi qubits, and|means tracing out the corresponding register.
Entropy 24 00899 g001
Figure 2. Layered circuit and ancillary layered circuit. (a) gives an example of an m layered circuit, which is employed in (b) as U l ( θ ) . U i ( θ i ) are parametrized by tuneable θ i . (b) is ancillary layered circuit, which has three parts: encoding circuit (blue) U e , layered circuit (pink) U l I and decoding circuit (green) U d . Ω 0 , Ω 1 , Ω 2 , and Ω 3 represent temporary states.
Figure 2. Layered circuit and ancillary layered circuit. (a) gives an example of an m layered circuit, which is employed in (b) as U l ( θ ) . U i ( θ i ) are parametrized by tuneable θ i . (b) is ancillary layered circuit, which has three parts: encoding circuit (blue) U e , layered circuit (pink) U l I and decoding circuit (green) U d . Ω 0 , Ω 1 , Ω 2 , and Ω 3 represent temporary states.
Entropy 24 00899 g002
Figure 3. Ancillary layered circuit to implement method 2. U i + 1 ( θ ) is the circuit to be optimized with known U i . Notations are the same as Figure 2a. Ω 0 , Ω 1 , Ω 2 , and Ω 3 are temporary states.
Figure 3. Ancillary layered circuit to implement method 2. U i + 1 ( θ ) is the circuit to be optimized with known U i . Notations are the same as Figure 2a. Ω 0 , Ω 1 , Ω 2 , and Ω 3 are temporary states.
Entropy 24 00899 g003
Figure 4. Expressibility values computed for circuits specified in Appendix C and circuit used in our protocol. Circuit 1 (red circles), 2 (green diamonds) and 3 (blue squares) are repeatedly applied with up to 5, while the circuit employed in our strategy (cyan cross) keeps unchanged.
Figure 4. Expressibility values computed for circuits specified in Appendix C and circuit used in our protocol. Circuit 1 (red circles), 2 (green diamonds) and 3 (blue squares) are repeatedly applied with up to 5, while the circuit employed in our strategy (cyan cross) keeps unchanged.
Entropy 24 00899 g004
Figure 5. Results of Method 1 and Method 2 for simulating the evolutions by ρ for a period of t = 0.05 , 0.1 and 0.2 . The fidelities (vertical axis) vary with the iterations (horizon axis). ρ is chosen as (a) bell state; (b) GHZ state and (c) Hamiltonian of a molecular Crotonic acid system.
Figure 5. Results of Method 1 and Method 2 for simulating the evolutions by ρ for a period of t = 0.05 , 0.1 and 0.2 . The fidelities (vertical axis) vary with the iterations (horizon axis). ρ is chosen as (a) bell state; (b) GHZ state and (c) Hamiltonian of a molecular Crotonic acid system.
Entropy 24 00899 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, K.; Gao, P. A NISQ Method to Simulate Hermitian Matrix Evolution. Entropy 2022, 24, 899. https://doi.org/10.3390/e24070899

AMA Style

Li K, Gao P. A NISQ Method to Simulate Hermitian Matrix Evolution. Entropy. 2022; 24(7):899. https://doi.org/10.3390/e24070899

Chicago/Turabian Style

Li, Keren, and Pan Gao. 2022. "A NISQ Method to Simulate Hermitian Matrix Evolution" Entropy 24, no. 7: 899. https://doi.org/10.3390/e24070899

APA Style

Li, K., & Gao, P. (2022). A NISQ Method to Simulate Hermitian Matrix Evolution. Entropy, 24(7), 899. https://doi.org/10.3390/e24070899

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop