By leveraging intrinsic properties of quantum mechanics, quantum computers may provide a disruptive computational tool for solving problems such as optimization (travelling salesman, machine learning), prime factorization, simulation of matter (materials, molecules and proteins), and financial market modeling.
The rise of increasingly powerful new quantum devices drives research to find algorithms increasingly able to leverage the computational power of small quantum processors. Hybrid quantum-classical algorithms have been introduced to this end. Known as variational quantum algorithms (VQA), such methods make use of a quantum computer only for the most performance-critical part of a program, while a more robust classical co-processor serves as an optimizer which operates in a loop with the quantum device to find the quantum circuit that transforms some state into the ground state of a given Hamiltonian.
Variational quantum algorithms leverage the use of a variational circuit, the Ansatz, that is to be executed on the quantum computer. Existing Ansatze typically aim to preserve symmetries of the problem Hamiltonian with the goal of restricting the variational search to a small subspace of the Hilbert space, and ultimately limit the expected amount of iterations to solve a problem. While existing approaches were satisfactory to a certain degree, there remained room for improvement.
It was found that using, in parameterized quantum circuits (Ansatze), Hamiltonians which break the symmetry (i.e. which does not commute with the problem Hamiltonian), could counter-intuitively lead to efficiency gains. Depending on the exact scenario, greater gains may be achievable when the symmetry-breaking Hamiltonians, while structured on the problem, allow to explore to a certain extent, but not so much as to become “lost”.
In accordance with one aspect, there is provided a method comprising performing a quantum computation including, in sequence, initializing a plurality of qubits of a quantum processor, applying a sequence of quantum logic gates onto the qubits in accordance with a parameterized quantum circuit which is based on a problem Hamiltonian Ĥproblem and at least one additional Hamiltonian Ĥk which does not commute with Ĥproblem, and measuring the expectation value of Ĥproblem in the final state of the qubits.
In practice, a large number of iterations will likely be performed wherein, after each iteration, the measured expectation value is provided to a classical computer, which calculates variational parameters for the next iteration, until the process converges onto the ground state.
In accordance with another aspect, there is provided a quantum circuit in the form of computer readable instructions stored in a non-transitory memory which, when executed by a classical processor is operable to drive the application of quantum gates onto qubits of a quantum processor, wherein the quantum circuit is based on a problem Hamiltonian Ĥproblem and at least one additional Hamiltonian Ĥk which does not commute with Ĥproblem.
In accordance with another aspect, there is provided a quantum circuit in the form of computer readable instructions stored in a non-transitory memory which, when executed by a processor is operable to simulate a quantum computation, wherein the quantum circuit is based on a problem Hamiltonian Ĥproblem and at least one additional Hamiltonian Ĥk which does not commute with Ĥproblem.
In accordance with another aspect, there is provided a method comprising simulating a quantum computation including, in sequence, applying a sequence of quantum logic gates onto a plurality of qubits in accordance with a parameterized quantum circuit which is based on a problem Hamiltonian Ĥproblem and at least one additional Hamiltonian Ĥk which does not commute with Ĥproblem, and measuring the expectation value of Ĥproblem in the final state of the qubits.
Many further features and combinations thereof concerning the present improvements will appear to those skilled in the art following a reading of the instant disclosure.
It will be understood that the expression “computer”, or “controller”, as used herein is not to be interpreted in a limiting manner. “computer” is rather used in a broad sense to generally refer to the combination of some form of one or more processing units and some form of memory system accessible by the processing unit(s). “controller” is used in a broad sense to generally refer to a device which performs a function of controlling, and may be a computer or another type of device. The memory system if a computer can be of the non-transitory type. The use of the expression “computer” in its singular form as used herein includes within its scope the combination of a two or more computers working collaboratively to perform a given function, independently of whether these two or more computers are local, remote, or distributed. Moreover, the expression “computer” as used herein includes within its scope the use of partial capabilities of a given processing unit.
A processing unit can be embodied in the form of a general-purpose micro-processor or microcontroller, a digital signal processing (DSP) processor, an integrated circuit, a field programmable gate array (FPGA), a reconfigurable processor, a programmable read-only memory (PROM), to name a few examples.
The memory system can include a suitable combination of any suitable type of computer-readable memory located either internally, externally, and accessible by the processor in a wired or wireless manner, either directly or over a network such as the Internet. A computer-readable memory can be embodied in the form of random-access memory (RAM), read-only memory (ROM), compact disc read-only memory (CDROM), electro-optical memory, magneto-optical memory, erasable programmable read-only memory (EPROM), and electrically-erasable programmable read-only memory (EEPROM), Ferroelectric RAM (FRAM) to name a few examples.
A computer can have one or more input/output (I/O) interface to allow communication with a human user and/or with another computer via an associated input, output, or input/output device such as a keyboard, a mouse, a touchscreen, an antenna, a port, etc. Each I/O interface can enable the computer to communicate and/or exchange data with other components, to access and connect to network resources, to serve applications, and/or perform other computing applications by connecting to a network (or multiple networks) capable of carrying data including the Internet, Ethernet, plain old telephone service (POTS) line, public switch telephone network (PSTN), integrated services digital network (ISDN), digital subscriber line (DSL), coaxial cable, fiber optics, satellite, mobile, wireless (e.g. Wi-Fi, Bluetooth, WiMAX), SS7 signaling network, fixed line, local area network, wide area network, to name a few examples.
It will be understood that a computer can perform functions or processes via hardware or a combination of both hardware and software. For example, hardware can include logic gates included as part of a silicon chip of a processor. Software (e.g. application, process) can be in the form of data such as computer-readable instructions stored in a non-transitory computer-readable memory accessible by one or more processing units. With respect to a computer or a processing unit, the expression “configured to” relates to the presence of hardware or a combination of hardware and software which is operable to perform the associated functions.
In the figures,
The controller can be expected to include a memory which can include functions in the form of computer readable instructions driving the operation of the controller, and data. The functions can include a “driving program”, for instance, which, in the case of gate-based quantum computing, can include a sequence of gates, typically referred to as a quantum circuit, stored as data in the memory. When the state of the qubits are read, the measured values can be stored in the form of data, for instance. As we will soon see, the data can also include variables in the form of values for parameters according to which the operations are to be performed.
During operation of the quantum computer 14, the sequence of quantum gates is performed in accordance with a quantum algorithm typically referred to as a quantum circuit in the art. In computing operations, a goal is to find the quantum circuit which will transform some state of qubits into the ground state of the problem. The typical process is iterative and essentially consists of making a sequence of educated guesses for the quantum circuit in a manner for the process to converge to the ground state. The quantum algorithm can be part of a classical/quantum algorithm where other processing steps are performed by the classical co-processor. In a family of approaches commonly referred to as “variational”, the quantum circuit is parameterized, meaning that the structure of the circuit (the configuration of the boxes along the lines in the upper right portion of
An example hybrid algorithm is presented in greater detail in
It will be noted here that in some alternate embodiments, instead of actually being performed on a quantum computer, the quantum computing iterations can be simulated by a quantum computing simulation program which can run on a classical computer. In this latter case, the algorithm can remain the same with the exception that the quantum computing iterations are simulated by a classical computer instead of being actually performed on quantum processing hardware.
The exact set of quantum gates allowable on a given quantum processor will vary based on the quantum processor’s build, but software solutions are available to translate (compile) a given Ansatz into terms of quantum gate elements which are available on a given quantum processor, as a function of the specific application. Some specific types of quantum processors are more efficient at running specific types of quantum circuits.
More specifically, in an example scheme of Variational Quantum Algorithms (VQAs), one uses a quantum computer to prepare a quantum state |ψ(θ)〉, which is parametrized by a collection of a manageable number of classical variational parameters θ. This state is typically generated from a reference state, |ψ0〉, to which a parametrized quantum circuit is applied to find |ψ(θ)〉 = U(θ)|ψ0〉. This parametrized quantum circuit is often referred to as the variational Ansatz or variational form. The value of the variational parameters is then classically tuned in order to minimize the result of a cost function, which in the case of chemistry and physics problems is often the energy of the state given by the expectation value of the problem Hamiltonian
VQAs thus involve four potential optimization areas:
The following text concentrates on point 2, designing an efficient parameterized quantum circuit. The following questions are of interest in this task:
Because the Hilbert space is large, one may be motivated to investigate Ansâtze that preserves symmetries of the problem Hamiltonian, to constrain the system to a much smaller Hilbert space which, in principle, would increase the chances of converging to the target state in the case where it is contained in the subspace, even though it is not guaranteed that this approach allows for an efficient search through parameter space.
Counter-intuitively, it was found that introducing a symmetry-breaking component into the otherwise symmetry-preserving Variational Hamiltonian Ansatz (VHA) could allow to find shortcuts in the parameter landscape and enable convergence from a greater number of reference states. Indeed, as will be presented below, it was found that adding symmetry-breaking terms could increase the performance of VQAs in solving the Fermi-Hubbard model for 4 or more sites.
The flexibility of VQAs comes from the freedom one has for the design of the variational form U(θ), with a plethora of choices applicable to different problems.
In the following text, we will explore the Variational Hamiltonian Ansatz (VHA), but also the Quantum Control Theory Ansatz, and an Ansatz based on the combination of VHA and Quantum Control Theory, which will be referred to as Quantum Optimal Control Inspired Ansatz (QOCA), before discussing an example specific application to the Fermi-Hubbard model.
Ansâtze that emulate the physics of the problem can be used to favor scalability. To construct physics-inspired Ansâtze, one uses prior knowledge of the problem’s structure when designing the variational form. The Variational Hamiltonian Ansatz (VHA) is one such Ansatz, which consists of a parametrized version of the quantum circuit implementing time evolution (Trotterization) under the problem Hamiltonian. The VHA is the translation of the Quantum Approximate Optimization Algorithm (QAOA) from combinatorial optimization problems to chemistry and physics problems. In the original QAOA, the variational Ansatz is constructed by alternating the application of a diagonal cost-function based Hamiltonian with a mixing Hamiltonian (typically
This idea can be extended to include a more general class of mixing Hamiltonians. In the VHA framework, the cost-function and mixing Hamiltonians are given by the different non-commuting terms of the problem Hamiltonian Ĥ= Σj Ĥj. The state-preparation unitary therefore reads :
where θ = {θj,d} are the variational parameters. Typically, one can arrange the terms of the Hamiltonian into a manageable number of groups of commuting terms, minimizing the Trotter error. In the analogy of time evolution, the depth d is associated with each time increment of the Trotterization.
This approach can be implemented using few variational parameters, therefore easing the classical optimization, but the circuit depth, which depends on the complexity of the problem, can be very large.
Because it is built around the problem Hamiltonian, the VHA implements operations that conserve the symmetries of the problem. For example, if no term of the Hamiltonian allows the number of particles to change, this quantity will be conserved in the variational state. This restricts the variational search to a small subspace of the Hilbert space which, in the case where it comprises the target state, can increase the performance of the VQA.
Quantum control theory can be applied to the quantum control of chemical reactions, spins in nuclear magnetic resonance experiments, and to pulse shaping in superconducting qubits, to name a few examples.
To manipulate a quantum system of Hamiltonian Ĥ0, it is coupled to one or many control apparatuses. The interactions that follow are given by a set of time-independent control Hamiltonians {Ĥk} which are parametrized by the time-dependent coefficients {ck(t)} ∈ ℝ. In general, the control Hamiltonians do not commute with Ĥ0. The total Hamiltonian is
leading to a time evolution described by
The solution to equation 4 is the unitary , which can propagate pure states through time as
. Taking this into equation 4, the differential equation for
can be obtained:
If the system is fully controllable, quantum control theory states that, for any initial state |ψ(0)〉, there exists a set of admissible controls {ck(t)} and a time T > 0 for which the time-evolved state |ψ(T)〉 is a desired target state |ψtarget〉. Or equivalently, that the resulting time propagator Û(T) implements a desired unitary operation Ûideal.
Quantum optimal control refers to methods by which the control pulses are designed for achieving a desired state preparation. This can be realized by seeking the set of controls and time T that optimize a cost functional characterizing the state-preparation fidelity C[{ck(t)},T], which may include constraints such as the control time and the maximum pulse strength.
In QOC methods time can be discretized into N increments, or pixels, of duration Δt such that the total evolution is accomplished in a time T = NΔt. This discretization may correspond to the one used for the numerical integration of the Schrödinger equation if the system is simulated or to the sample rate of the control electronics in the case of an experiment. It can also be set arbitrarily.
Using this discretization, the continuous control fields ck(t) are parametrized by the discrete control fields uk as
where ⊓j (t,Δt) ≡ Θ(t - jΔt) - Θ(t - (j + 1)Δt) with Θ the Heaviside function. The time evolution operator for the jth pixel is
and thefore the propagator for a time T is
Because this time propagator incorporates the symmetry-breaking terms {Ĥk}, it can potentially allow for quantum operations that are much faster than adiabatic evolution governed by Ĥ0 alone, at least when the control fields are optimized. This method is thus particularly well suited to cases where the target state (or operation) is known beforehand. In such a context, these fast quantum operations are possible in view of the introduction of control, or drive, Hamiltonians {Ĥk} in the system.
Let us now considering applying the drive terms {Ĥk} to the VHA Ansatz to create a new Ansatz referred to as Quantum-optimal-control-inspired Ansatz (QOCA). More specifically, the VHA Ansatz, which resembles time evolution under the problem Hamiltonian Ĥproblem, is modified to incorporate a set of drive terms {Ĥk}. This can yield a time evolution under the new Hamiltonian
The variational state preparation circuit for QOCA can now be by parametrizing the time evolution operator that follows from the procedure of the previous section as
Here, θ = {θj,d} and δ = {δk,d} are respectfully the VHA VHA parameters and the drive parameters. As for VHA, the d parameter is the number of layers of the Ansatz. δk,d can be used here instead of uk,j to denote the discrete control fields, and thus create a notation gap between QOC and VQAs.
In principle, one has the freedom to select drive Hamiltonians that do not commute with Ĥproblem, and therefore several alternate drive Hamiltonians may be used in alternate embodiments. To facilitate the identification of drive Hamiltonians which will have the most positive impact on the outcome of the VQA, search strategies such as adapt-vqe may be used.
One of numerous possible potential implementations of QOCA is the Fermi-Hubbard model (FHM). Indeed, since FHM is considered to be notoriously hard to solve using classical computers, it can make a good candidate for a proof of principle experiment, and it is thus believed that achieving performance with FHM is a good indication that adaptations of the Ansatz would perform on other problems. Indeed, FHM offers one of the simplest ways to investigate systems with strong interactions between fermions (strongly-correlated systems). Despite its apparent simplicity, it can give rise to metallic, insulating, magnetic, and even superconducting phases of matter which motivated its study to understand systems ranging from heavy fermions to high-temperature superconductors.
A particularly interesting instance of the FHM is the so-called half-filling regime at intermediate coupling. In this regime, even some of the biggest classical supercomputers have failed to solve the FHM for more than 20 lattice sites.
In terms of fermionic ladder operators , the FHM Hamiltonian for L lattice sites is given by
where i,j are the lattice site indices and σ={↑,↓} is the spin index. (i,j) denotes a sum over nearest-neighbor sites and is the occupation operator of spin-orbital iσ. The first term of eq. 11 represents the nearest-neighbor hopping of amplitude t and will often be denoted T. The second term is the on-site Coulomb repulsion of strength U, and the last term is the chemical potential. We often use
to denote the last two terms. The value of µ sets the filling of the ground state, with µ=U/2 ensuring half-filling. Because the Pauli exclusion principle only allows for occupations of up to two fermions per lattice site, half-filling occurs when there is on average one electron per site,
Here, a Hybrid quantum-classical methods is used as a way to investigate the thermodynamical properties of the FHM. The challenge associated to preparing the ground state of the model on a quantum processor is tackled by a VQA.
More specifically, the Jordan-Wigner (JW) transformation can be used to encode the FHM into a qubit Hamiltonian. In the JW framework, each fermionic site is encoded into the state of two qubits as (0, ↑, ↓, ↑↓) ↦ (00,01,10,11) and the fermionic ladder operators are given by where and the greek subscripts denote the spin-orbital indices as v=jo. For a system of L sites, the N=2L spin-orbitals are arranged as
but the choice is purely conventional.
Hopping terms between spin-orbitals v and µ with v<µ transforms as where the product of %operators, referred to as the JW string, vanishes when µ=v+1. The number operator on orbital v and therefore the onsite Coulomb interaction between orbitals v and µ become, using the convention
At half-filling, i.e. setting µ=U/2 in eq. 11, the single coming from the onsite interaction terms are cancelled by the chemical potential, leaving a simpler Hamiltonian composed solely of X.X·Y.Y (plus JW strings when fermionic orbitals are not adjacent) and ZZ qubit interactions.
The performance of VQAs strongly depends on the choice of initial state and variational parameters. The initial state acts as an educated guess to the target state and is often one that is easily computed classically. However, its preparation on a quantum computer is generally not an easy task and may require a large overhead in terms of circuit depth.
To compare all Ansâtze on the same footing, we first used the equal superposition of all basis states as initial state. It is prepared by simply applying Hadamard gates on all qubits . One advantage of this state is that it is easy to prepare. Also, for simulation done at half-filling, this state already has the desired number of particles and total spin since it is akin to placing half a fermion on every fermionic orbital.
This choice of initial state allows for demonstration of the robustness of the Ansâtze to unstructured initial conditions. The performance can be increased at the expense of deeper circuits by using the ground state of the non-interacting FHM i.e. fixing U = µ = 0 in eq 13, as initial state.
In the application of drive terms to a VHA, as opposed to QOC, the drive terms do not need to describe actual physical interactions present in a control hardware, providing for more freedom to engineer desired interactions.
In a first step, the drive was adapted from QOC experiments where control is done through RF voltage pulses. In this case, the X and Y quadratures of the electromagnetic field are driven generating the interaction Hamiltonians where
and
are bosonic operators.
Expressing these terms by means of fermionic operators and applying a drive on every lattice site, the following interaction Hamiltonians are achieved: where are now fermionic operators. This formulation no longer describes the X and Y quadratures of a bosonic field but rather describes a coupling of each lattice site to a bath of Majorana fermions
and
.
The drive equations for a spinless system can be derived, and the resulting circuit can be applied to both species of spins. Doing the JW transformations of eq. 12 on equations 15 and 16 et find:
To incorporate equations 17 and 18 into the QOCA variational form of equation 10, the drive terms are exponentiated. By performing a first-order Trotter-Suzuki decomposition, we obtain the drive circuit equation for one layer of the Ansatz, where {δk,d} are the variational parameters associated with each drive term. Note that we intertwined the application of Ĥ1 and Ĥ2 to better emulate a simultaneous drive of all sites. A schematic of the circuit implementing equation 19 for four qubits is illustrated in
More specifically,
is the Hadamard gate which changes between the X̂ and Ẑ bases and
is the equivalent transformation between the Ŷ and Ẑ bases. The angles of the Ra(θ) = exp[-iθâtσa/2] rotations are the variational parameters, where âtσa is a Pauli matrix.
The circuit designer also has the freedom to decide where to incorporate the variational parameters within the Ansatze. It is possible to interpolate between two parametrization strategies to achieve this, for instance. In the first one of these strategies, the Ansatz is fully parametrized and the second one presents a number of parameters that is independent of the number of qubits used for the computation. It will be noted that in various embodiments, various variational parameter incorporation strategies can be used, such as focussing on only the first one, focussing only on the second one, using a different variational parameter incorporation strategy and interpolating between any two variational parameter incorporation strategies, to name some possible examples.
The first strategy corresponds to defining all (or almost all) gate angles as variational parameters. This can give the classical optimizer a lot of freedom to explore the Hilbert space spanned by the Ansatz, but comes at the cost of longer optimization time. Moreover, this strategy might not be suitable for larger instances of problems which could require an larger number of parameters.
Using this strategy for VHA can consist of assigning one parameter to every âdiσâjσ + h. c. hopping terms and duplicating the parameter for both flavors of spins. This is because at half-filling and zero total spin, there is a spin-inversion symmetry which removes the need to treat spins up and down differently. Additionally, every term of the on-site interaction is associated with a variational parameter.
QOCA can be implemented similar to VHA, to which we add the drive terms where all
and
become variational parameter in equation 19. The fully parameterized version of QOCA will be referred to as QOCAf hereinbelow.
In the second strategy, scalable parametrization, we employ a number of variational parameters that is independent of the system size. Because there are fewer parameters, we expect the optimization to be faster, but larger circuit depths might be necessary to achieve the same accuracy a full parametrization of the Ansatz would provide.
In physics-inspired Ansatze, a scalable parametrization can be performed by collecting terms of the Hamiltonian into a constant number of groups and affiliating a single parameter to each group. For example, a common way of grouping the terms of the 2D FHM is
where the first four Hamiltonians group the even and odd, vertical and horizontal hopping terms and ĤU collects the on-site interaction terms. Note that for the 3D FHM, two additional groups of hopping terms covering the third dimension may be necessary. This partitioning can be chosen in order to minimize the Trotter error when implementing the Ansatze since every term within each group commute with all the others.
The scalable parametrization version of QOCA will be referred to as QOCAs hereinbelow.
One drawback of QOCA as presented above is that the drive circuits can be deep depending on the nature of the drive. It can be preferred to shorten the total Ansatz circuit without compromising the effect of the drive and therefore the overall performance.
Because we apply the drive D of
In at least some embodiments, the kinetic part of the VHA can be removed with an acceptable tradeoff in performance. The kinetic part is more costly in circuit depth than the interaction part, which can lead to efficiencies. The result will be referred to as short QOCA (sQOCA) Ansatz. It can be expressed as
where V̂ is the on-site interaction part of eq. 11 and v = {vj,d} are the associated variational parameters. sQOCA is presented schematically at
More specifically,
Numerical simulations were conducted using the VQA tools provided by Qiskit Aqua™. A unitary statevector simulator was used since no noise was considered and we assume an all-to-all coupling graph. The COBYLA method was used as the classical optimizer with a maximum number of function evaluation of ~ 105. This number can be justified using experimentally realistic arguments.
So far, our discussion of variational Ansatze remains hardware-agnostic in the sense that any digital quantum computer can be used for our purpose. The list of available technology includes: superconducting qubits, trapped ions, spin qubits, photonic qubits, Rydberg atoms, among others. Each one of these architectures operates differently and allow for different quantum operations to be executed natively on the hardware. Moreover, different approaches within a particular technology might also lead to different quantum operations.
Before executing a quantum circuit on a real quantum computer, it must be compiled into native gates, that is the set of quantum logical gates the hardware is calibrated to perform. This native gates set is typically composed of arbitrary single-qubit rotations and a two-qubit interaction. Furthermore, the compilation should take into account the qubit connectivity layout of the quantum computer. Because of these constraints, a particular architecture might be more suitable than others to run a quantum circuit.
For example, Google’s approach to superconducting qubits seems suitable for fermionic problems. In particular, their choice of two-qubit gate, the so-called fermionic simulation gate
can natively perform parametrized X̂X̂+ ŶŶ and ẐẐ interactions, which are common in fermionic problems since they are particle conserving. Other architectures might have advantages for other quantum or classical problems.
A quick sidenote will be made here that quantum computers are typically noisy, but the nature of the noise depends strongly on the exact nature of the hardware of the quantum computer on which the quantum circuit will ultimately be run. Accordingly, in the simulations which were conducted, noise was not simulated. It is believed that simulation of noise may be preferably conducted at a later time, in view of establishing whether, say, a specific build of quantum computer is well adapted to run a specific quantum circuit, for instance.
The results of simulations using VHA, QOCAs, QOCAf, and sQOCAf are presented in
The results of the simulations using VHA, full parametrization QOCA (QOCAf), Scalable parametrization QOCA (QOCAs), and another, fully parametrizable, Ansatz HEA, are presented in
All variational parameters can be initialized to 0 at every opportunity, so that Hamiltonian-based Ansatze implement the identity operator at first. This choice can be made to ensure that the variational search truly starts from the unaltered initial state. In some cases, however, this strategy resulted in premature convergence of the optimizer into states close to the initial guess, which could be circumvented by using a random initialization of the parameters. Finally, all layers of the Ansatze where optimized simultaneously. In alternate embodiments, improvement may be achievable by adopting a layer-by-layer optimization strategy.
It is interesting to see in
Variational quantum algorithms are based on heuristics. Variational forms can be built using some physical intuition, but theory may not allow to go so far as to explain why a particular construction might or might not work. A first reasonable approach would be to build a variational Ansatz only using terms that preserve symmetries of the problem Hamiltonian Hproblem. This way, we ensure that the variational landscape is completely structured with respect to Hproblem in the sense that every possible variational state, corresponding to every point in this variational landscape, will conserve symmetries of Hproblem. This approach is typically good as it prevents the VQA from searching an exponentially large Hilbert space. Indeed, the search is restricted only the small sub-space of the full Hilbert space corresponding to the symmetry sectors of Hproblem.
However, allowing some symmetries to be broken can lead to faster convergence to the target state or unitary operation. One way to achieve this is by controlling the system using drive terms that do not commute with the system’s Hamiltonian. This allows for a richer set of states to be reached, which often leads to faster quantum operations. In the case where the control is operated through laser pulses, these drive terms typically take the form (b̂† + b̂) and i(b̂† - b̂), where b̂† and b̂ are the bosonic creation and annihilation operators of the controlled system.
As a first attempt to build drive terms for fermionic systems, we textually translated these bosonic drives into fermionic operators to obtain
where
and âj are now fermionic ladder operatorsoperators respecting the anti-commutation relations
and
Intuitively, one can expect these linear terms to perform well since they do not commute with any physical fermionic Hamiltonian, which are composed of quadratic and quartic fermionic terms only. Indeed, these latter equations do not conserve neither the number of particles nor the parity symmetries. Moreover, when transformed to a qubit form via the Jordan-Wigner transformation,
we remark that these drive terms decompose into a series of multi-qubit interactions which produce entanglement amongst the qubits. This is desired since the target state is likely to be highly entangled. As it will be explained later, this choice of drive terms produced promising results for the Fermi-Hubbard model and the water molecule, two hard problems of condensed matter physics and quantum chemistry respectively.
In what follows, we first formalize why non-commuting terms can help produce a richer set of states and then give some intuition as to why the selected drive terms of equations 23 and 24 are a good choice for fermionic problems.
To know if adding a non-commuting drive term to the Hamiltonian help reaching more states, we must first know what states can be reached without the addition of the drive term. For this purpose, we are interested in the notion of controllability of a given system Ĥ0 + Σk ck(t)Ĥk.
We start by writing the Schrödinger equation of the control problem described above
The solution of which is the unitary Û(t), which can propagate pure states through time as |ψ(t)〉 = Û(t)|ψ(0)〉.
The system Ĥ0 + Σk ck(t)Ĥk is said to be controllable is the set of possible matrices that can be obtained from solving equation 27 is the set of all the unitary matrices. Let us begin by detailing what a reachable set is.
The reachable set at time T > 0 for system of equation 27, ℜ(T), is the set of all unitary matrices
This relates to the system through a theorem which states that the set of reachable states for system Ĥ0 + Σk ck(t)Ĥk is the connected Lie group associated with the Lie algebra L generated by span
where the Lie group eL is defined as
The system is controllable if dim(L) = n2 = dim(u(n)), where n is the size of the Hilbert space and u(n) is the Lie algebra of skew-Hermitian n × n matrices. This is equivalent to L = u(n) or eL = U(n), where U(n) is the Lie group of unitary matrices of dimension n. The system is also said to be controllable when dim(L) = n2 - 1 = dim(su(n)), where su(n) is the Lie algebra of matrices of u(n) with zero trace, or equivalently when L = su(n) or eL = SU(n), where SU(n) is the set of matrices of U(n) with determinant equal to one.
The dynamical Lie algebra of the system L can be generated iteratively by taking all (possibly repeated) commutators of elements of L and appending L when a new matrix is found. In other words, dim(L), therefore the set of reachable states ℜ, increases for every non-commuting control Hamiltonian Ĥk that is added to the set and generates a new term.
This shows that adding drive terms {Hk} that don’t commute with Ĥproblem in the QOCA variational form (equation 10) can increase the set of reachable variational states and therefore can possibly lead to shortcuts in the parameter landscape.
In the QOCA variational form, the goal is to add drive terms that could expand the set of reachable states. A good starting point would be to sample these terms from a fully controllable fermionic Hamiltonian since we know from the theorem that the set of reachable states is given by the span of its terms. The fermionic Hamiltonian
is universal in the sense that any unitary matrix can be generated by solving its Schrödinger equation as in equation 27. This control Hamiltonian (eq. 30) is unphysical since it contains linear terms that break the parity symmetry, which is not allowed for physical Hamiltonians. Indeed, physical Hamiltonians must only contain terms that are at least quadratic in fermionic operators and are an even power thereof. For example, the Fermi-Hubbard model only contains quadratic â†â and quartic â†ââ†â terms which are captured by the last two sums of eq. 30.
Because physical Hamiltonians never contain linear terms, using drives of the form of α(t)â + a∗(t)↠guaranties that they will not commute with Ĥproblem and therefore expanding the set of reachable states.
In the context of VQAs, one could build the most general fermionic variational form by performing a Trotterization of the time evolution operator when the control functions α, β, γ are constant piecewise functions. The variational circuit therefore reads
where α = {αj,d}, β = {βi,j,d} and γ = {γi,j,d} are the variational parameters. Because this is a Trotterization, the universality is guaranteed only in the limit of large d. The variational form can be seen as a more general version of QOCA (eq 10).
Although the Ansatz of quation 31 is appealing because of its universality, it will lead to very long quantum circuits and is completely unstructured with respect to Ĥproblem. To this end, with the QOCA variational form, we propose to only keep the quadratic and quartic terms of equation 30 that are also present in Ĥproblem alongside the linear drive terms. This way, the variational landscape is structured around the symmetries of Ĥproblem and the drive terms allow for exploration outside of that sub-space of the Hilbert space.
If some insights are known about which symmetries are broken in the target state, one can use this information to construct drive terms that are tailored to the problem. For example, if we know that we are looking for a state that has an antiferromagnetic order and no term of describes such an order, we can think of adding a drive that has the shape of an alternating magnetic field
where
is the total spin of site j and h is the strength of the magnetic field. Moreover, if one wishes to enforce singlet superconductivity, a drive of the form
or more generally,
might be useful since they create and annihilate Cooper pairs.
A non-exhaustive list of potential drive terms is presented in table 1, below, for the case of bosonic and fermionic systems. Note that not all drive terms are physical interactions.
The QOCA method was tested on the Fermi-Hubbard model (ref eq 13 above) with the drive terms 1 and 2. The VHA method was applied on the same problems. To put all Ansatze on an equal footing and to test their robustness to simple initial states, we initialize the circuit in the equal superposition of all basis states |ψ0〉 = H⊗N|0〉⊗N, which is easy to prepare.
To quantify the quality of the variational state |⊗(θ)〉 that is produced, we calculate its fidelity with the target state |Ω〉, which is the ground state of the FHM. We define the fidelity cost-function as
The number of particles is always constant for VHA since this Ansatz does not allow this quantity to change. This Ansatz presents bad performance, which can be seen from the infidelity remaining high until the maximum number of iteration is reached.
For QOCA with the full and scalable parametrizations, we observe that the number of particles deviates from one as it is enabled by the drive terms which are particles non-conserving. The interesting fact is that peaks in the number of particles are often associated with the triggering of abrupt descents in the infidelity plot.
To corroborate this fact we also present results obtained with the hardware-efficient Ansatz, another Ansatz that does not preserve symmetries of Ĥproblem. These results also show peaks in particle number associated with drops of infidelity.
These results give insight that allowing the breaking of symmetries, such as the number of particles, helps the classical optimizer to find regions of steep gradients that are otherwise impossible to reach. This leads to faster convergence to the target state.
As a very preliminary test of QOCA for quantum chemistry problems, we tackled the H 2O molecule with 12 qubits. The Ansatz used was one based on the Hamiltonian of a periodic 1 × 6 Fermi-Hubbard chain with the drive terms 1 and 2 above. A full parametrization strategy us used and the initial state is the easily prepared Hartree-Fock approximation to the ground state.
Convergence process results for the H 2O molecule are presented in
The same reasoning that was done for fermions to build the QOCA method can be applied to bosonic and spin -
(qubit) problems. The differences are that the universal control Hamiltonians equivalent to equation 30 differ in both cases and so are the associated drive terms. The type of problems and how they are mapped to a quantum computer is also different. Here, we briefly state how one would address bosonic and
problems using a QOCA-like method.
Bosons are described by the Lie algebra of bosonic ladder operators
and
which obey the commutation relations:
A universal control Hamiltonian for bosonic systems is of the form
where the α, β, γ and δ are controllable. Note that contrary to the case of fermions, every term here describes physical processes. This Hamiltonian, and the corresponding Lie algebra generated by {-iĤb}, can be sampled to constitute the symmetry-breaking terms {Ĥk} of the QOCA variational form (ref eq 10) applied to bosonic problems.
To closely relate to the discussion on fermions, we consider the Bose-Hubbard model as an example of bosonic problem
This Hamiltonian describes interacting bosons on a lattice. n̂i = b̂i†b̂i is the number operator, t is the hopping amplitude, U is the interaction strength and µ is the chemical potential.
There are several ways to map bosonic modes to qubits which have varying ressource requirements in terms of number of qubits and allowed operations.
As an example, one could use the Single-Excitation-Subspace encoding in which we use nmax qubits to encode the first nmax energy level of a bosonic mode. Such a truncation is necessary as the Hilbert space associated with bosonic modes is infinite-dimensional. We map the state containing n excitation as
and the associated ladder operators as
As for the Fermi-Hubbard model, the variational circuit could take the form of
where Ĥproblem is now the Bose-Hubbard model (equation 39). Because there are no linear terms in equation 39, potentially good drive Hamiltonians are of the form
as it was the case for the Fermi-Hubbard model. Other possible bosonic drive terms are presented in Table 1. All of these term are elements of the Lie algebra generated by the universal control Hamiltonian for bosonic systems presented above.
Qubits are described by the Lie algebra of Pauli matrices σ̂x, σ̂y and σ̂z which satisfy the commutation and anti-commutation relations:
where a, b, c ∈, x, y, z and εabc is the Levi-Civita symbol.
A theorem states that generic single-qubit control and any interaction between two qubits is sufficient for building any unitary matrices. A universal control Hamiltonian for qubits systems is therefore of the form
where
denotes σa, applied on qubit i. Again, this Hamiltonian, and the corresponding Lie algebra generated by {-iĤq}, can be sampled to constitute the symmetry-breaking terms {Ĥk} of the QOCA variational form applied to
problems.
An example of a
problem would be the Heisenberg XYZ chain with periodic boundary conditions defined as
where we associate the indices N + 1 = 1. The Ja are coupling constants and h is the energy of every identical spins.
As opposed to fermions and bosons,
problems are trivially mapped to qubits. Indeed, they share the same Lie algebra su(2).
Again, as for the last two examples, the variational circuit could take the form of (eq 12) where Ĥproblem is now the Heisenberg XYZ chain with periodic boundary conditions defined as above. The drive terms could be of the form ξab(t)σ̂aσ̂b with a ≠ b, since these terms do not commute with ĤXYZ, or a generic single-site drive of the form
The last two sections illustrate how one would tackle quantum mechanical problems of different nature using a QOCA-like method. However, it is possible the apply the same technique to a broader class of problems e.g. a class that would include classical Hamiltonians. The main steps of how one could proceed to build the QOCA variational form may be broken down in the following way:
As can be understood in view of the above, the examples described above and illustrated are intended to be exemplary only. It will be further understood that, for instance, it is possible to specifically dictate parameters to be used in a parameterized quantum circuit, rather than have these parameters calculated iteratively using classical processing. Moreover, while typically a same circuit structure is used for multiple quantum computation iterations, some meta-algorithms such as adapt-vqe can allow changing the circuit structure itself during the processing. This would typically be done following a plurality of quantum processing iterations, for example. Various initialization schemes can be used in different embodiments. For instance, the initialization routine of the quantum computation can include a simple preparation containing single qubit gates such as H⊗n|0〉⊗n , the Hartree-Fock state of quantum chemistry, or a more complex initialization routine containing one or more single-qubit gates and one or more multi-qubit gates. A suitable choice can be selected and tested as a function of the details of a specific embodiment. The scope is indicated by the appended claims.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CA2021/050468 | 4/8/2021 | WO |
Number | Date | Country | |
---|---|---|---|
63007573 | Apr 2020 | US |