Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
QUANTUM COMPUTATION ON A SPECIFIC QPU TOPOLOGY
Document Type and Number:
WIPO Patent Application WO/2023/170334
Kind Code:
A1
Abstract:
The invention relates to the field of quantum computing. The invention includes a method for determining a control sequence for performing digital-analog quantum computation on a quantum processing unit with a star topology. The use of a QPU with this topology enables a significant reduction in the number of single qubit gates that must be applied to the system, thereby reducing circuit depth and reducing error.The invention also include a method for performing digital-analog quantum computation on a QPU with a star topology.

Inventors:
PINA CANELLES VICENTE (FI)
HEIMONEN HERMANNI (FI)
LAMPRICH JULIA (FI)
WURZ NICOLA (FI)
KU HSIANG-SHENG (FI)
Application Number:
PCT/FI2023/050086
Publication Date:
September 14, 2023
Filing Date:
February 13, 2023
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
IQM FINLAND OY (FI)
International Classes:
G06N10/20
Other References:
ADRIAN PARRA-RODRIGUEZ ET AL: "Digital-Analog Quantum Computation", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 16 July 2020 (2020-07-16), XP081702413, DOI: 10.1103/PHYSREVA.101.022305
FRANCESCO TACCHINO ET AL: "Quantum computers as universal quantum simulators: state-of-art and perspectives", ARXIV.ORG, CORNELL UNIVERSITY LIBRARY, 201 OLIN LIBRARY CORNELL UNIVERSITY ITHACA, NY 14853, 8 July 2019 (2019-07-08), XP081439186
Attorney, Agent or Firm:
BERGGREN OY (FI)
Download PDF:
Claims:
Claims 1. A computer-implemented method for determining a control sequence for performing digital- analog quantum computation, DAQC, on a quantum processing unit, the quantum processing unit comprising a central qubit or resonator q0 and a plurality of external qubits q1 to qN−1, each external qubit being connected directly to the central qubit, the method comprising: - receiving coupling coefficients of a target Hamiltonian and a time tF after which the quantum state of a system described by the target Hamiltonian is to be calculated; - determining a sequence of analog blocks to be applied to the central and external qubits, during which the quantum states of the central qubit and external qubits evolve according to a resource Hamiltonian; and - determining a sequence of single-qubit gates to be applied to the external qubits. 2. The method of claim 1, wherein determining the sequence of analog blocks to be applied to the central and external qubits comprises determining the durations t1 to tN−1 of the analog blocks A1 to AN−1 based on the coupling coefficients of the target Hamiltonian and the resource Hamiltonian. 3. The method of claim 2, wherein determining the sequence of analog blocks Ak further comprises placing the analog blocks in sequence based on their indices k and omitting analog blocks Ak with duration tk = 0. 4. The method of any preceding claim, wherein determining the sequence of single qubit gates comprises placing the single qubit gates in the sequence such that single-qubit gates are applied to qubits qj+1 to qk−1 at the end of or after each analog block Aj with non-zero duration tj and j < N− 1, and where Ak is the next analog block with non-zero duration tk. 5. The method of claim any preceding claim, wherein the analog blocks are applied continuously, and the single-qubit gates are applied simultaneously with the analog blocks at the end of each analog block. 6. The method of any of claims 1 to 4, wherein the analog blocks are applied discretely, and the single-qubit gates are applied serially in between the analog blocks. 7. The method of any preceding claim, wherein determining a sequence of analog blocks comprises determining the time durations t1 to tN−1 of the analog blocks according to the following equations: where coupling coefficient of the target Hamiltonian and is th h e jt coupling coefficient of the resource Hamiltonian. 8. The method of claim 7, wherein determining a sequence of analog blocks comprises determining a sequence of analog blocks with a non-zero duration tk. 9. The method of any preceding claim, wherein determining a sequence of single-qubit gates comprises placing single-qubit gates applied to qubits qj+1 to qN−1 at the beginning and end, or before and after each analog block Aj with non-zero duration tj and j < N − 1, and removing consecutive pairs of single-qubit gates. 10. The method of any preceding claim, wherein the sequence of single-qubit gates further comprises applying single-qubit gates to all qubits q2 to qN−1 at the start of or before the first analog block. 11. The method of any preceding claim, wherein the method further comprises outputting the determined control sequence. 12. A method for performing digital-analog quantum computation, DAQC, on a quantum processing unit, the quantum processing unit comprising a central qubit or resonator q0 and a plurality of external qubits q1 to qN−1, each external qubit being connected directly to the central qubit, the method comprising: - preparing initial quantum states in the central qubit and external qubits; - turning on interactions between the central qubit and external qubits such that the quantum states of the central qubit and external qubits evolve according to a resource Hamiltonian; - applying a sequence of single-qubit gates to the external qubits to modify the evolution of the quantum states of the central qubit and external qubits such that the quantum state of the central qubit and external qubits after a time period tF corresponds to evolution of the quantum states according to a target Hamiltonian after the time period tF. 13. The method of claim 12, wherein the sequence of single-qubit gates and time intervals between the application of single-qubit gates are controlled in order to modify the evolution of the state of the central qubit and external qubits to correspond to the target Hamiltonian.

14. The method of claim 12, wherein the method further comprises: - before each single-qubit gate is applied to one of the external qubits, turning off interactions between the central qubit and the remaining external qubits; and - after each single-qubit gate is applied to one of the external qubits, turning interactions between the central qubit and the remaining external qubits back on. 15. The method of claim 14, wherein the sequence of single-qubit gates and time periods between the application of single-qubit gates during which interactions between the central qubit and the external qubits are turned on are controlled in order to modify the evolution of the state of the central qubit and external qubits to correspond to the target Hamiltonian. 16. The method of any of claims 12 to 15, wherein a single-qubit gate is applied to each of the qubits q2 to qN−1 before turning on the interactions between the central qubit and external qubits, and wherein only one further single-qubit gate is applied to each of the qubits q2 to qN−1. 17. The method of any preceding claim, wherein the single-qubit gates are involutory. 18. The method of claim 17, wherein the single-qubit gates are Pauli X, Y or Z gates. 19. The method of any preceding claim, wherein the resource Hamiltonian and target Hamiltonian are Ising Hamiltonians of the form where the operator Aq is one of the Pauli matrices X, Y or Z operating on qubit q and where the resource and target Hamiltonians have different coupling coefficients gk. 20. The method of any preceding claim, wherein the coupling coefficients of the connections between the central qubit or resonator q0 and external qubits q1 to qN−1 are tuneable, and wherein the method further comprises tuning the coupling coefficients in order to shorten or reduce the number of analog blocks. 21. A computer system configured to perform the method of any preceding claim. 22. A computer system comprising a classical processing unit configured to perform the method of any of claim 1 to 11 and a quantum processing unit configured to perform the method of any of claims 12 to 20.

23. A computer program product comprising instructions which, when executed by a computer, cause the computer to perform the method of any of claims 1 to 20. 24. A computer-readable medium comprising instructions which, when executed by a computer, cause the computer to perform the method of any of claims 1 to 20.

Description:
QUANTUM COMPUTATION ON A SPECIFIC QPU TOPOLOGY Technical Field The present invention relates to the field of quantum computing. Background Digital-analog quantum computation (DAQC) is a paradigm for quantum computation in which the “analog” Hamiltonian corresponding to the natural evolution of qubits can be used, along with “digital” single-qubit gates, to perform universal quantum computations by effectively engineering an arbitrary Hamiltonian. DAQC can be implemented in two ways: stepwise DAQC (sDAQC) and banged DAQC (bDAQC). In sDAQC, entangling operations, i.e. analogue Hamiltonian evolution of the qubits, are turned on and off repeatedly throughout the computation. In bDAQC, single-qubit gates are applied simultaneously with the entangling operations, which are constantly on throughout the computation. A significant source of errors in stepwise DAQC arises from turning multi-qubit interactions on and off. Banged DAQC was introduced in order to mitigate these errors by leaving the analog multi- qubit interactions turned on throughout the computation; however, the singe qubit gates applied to each qubit do not commute with the terms relevant to that qubit in the Hamiltonian describing the analog multi-qubit interactions, therefore the evolution of the state of the qubit when the single qubit gate is applied is not exactly the same as in the idealised stepwise DAQC model. Summary of the Invention A first aspect of the invention relates to a computer-implemented method for determining a control sequence for performing digital-analog quantum computation, DAQC, on a quantum processing unit. At a high level, the control sequence may be a sequence of analog blocks and single qubit gates to be applied to the quantum processing unit. Alternatively, the control sequence may be an intermediate representation such as OpenQASM, at a hardware interaction level such as defining control signals to be applied to the qubits, or any suitable description of the parameters needed to perform the DAQC. The quantum processing unit for which the control sequence is determined comprises a central qubit q 0 and a plurality of external qubits q 1 to q N−1 , each external qubit being connected directly to the central qubit. The method comprises receiving coupling coefficients of a target Hamiltonian and a time t F after which the quantum state of a system described by the target Hamiltonian is to be calculate, determining a sequence of analog blocks to be applied to the central and external qubits, during which the quantum states of the central qubit and external qubits evolve according to a resource Hamiltonian, and determining a sequence of single-qubit gates to be applied to the external qubits.

Determining the sequence of analog blocks to be applied to the central and external qubits may comprise determining the durations t 1 to t N-1 of the analog blocks to N _ 1 based on the coupling coefficients of the target Hamiltonian and the resource Hamiltonian. Determining the sequence of analog blocks A k may further comprise placing the analog blocks in sequence based on their indices k and omitting analog blocks A k with duration t k = 0.

Determining the sequence of single qubit gates may comprise placing the single qubit gates in the sequence such that single-qubit gates are applied to qubits q j+1 to q k-1 at the end of or after each analog block with non-zero duration t j and j < N - 1, and where A k is the next analog block with non-zero duration t k .

The analog blocks may be applied continuously, and the single-qubit gates are applied simultaneously with the analog blocks at the end of each analog block according to the banged DAQC paradigm. Alternatively, the analog blocks may be applied discretely, and the single-qubit gates are applied serially in between the analog blocks according to the stepwise DAQC paradigm.

Determining the sequence of analog blocks may comprise determining the time durations t t to t N-r of the analog blocks according to the following equations: where is the j th coupling coefficient of the target Hamiltonian and g j is the j th coupling coefficient of the resource Hamiltonian. Determining the sequence of analog blocks may comprises determining a sequence of analog blocks with a non-zero duration t k .

Determining the sequence of single-qubit gates may comprise placing single-qubit gates applied to qubits q j+1 to q N-1 at the beginning and end, or before and after, each analog block A j with non- zero duration t j and j < N - 1, and removing consecutive pairs of single-qubit gates.

The sequence of single-qubit gates may further comprise applying single-qubit gates to all qubits q 2 to q N-1 at the start of or before the first analog block. The method further comprises outputting the determined control sequence. Outputting the control sequence may comprise providing the output in some machine- or human-readable format such as a quantum circuit or intermediate representation such as OpenQASM. Alternatively, the control sequence may be an intermediate representation such as OpenQASM, at a hardware interaction level such as defining control signals to be applied to the qubits, or any suitable description of the parameters needed to perform the DAQC. A second aspect of the invention relates to method for performing digital-analog quantum computation, DAQC, on a quantum processing unit. The quantum processing unit on which the DAQC is performed comprised a central qubit q 0 and a plurality of external qubits q 1 to q N-1 , each external qubit being connected directly to the central qubit. The method comprises: - preparing initial quantum states in the central qubit and external qubits; - turning on interactions between the central qubit and external qubits such that the quantum states of the central qubit and external qubits evolve according to a resource Hamiltonian; and - applying a sequence of single-qubit gates to the external qubits to modify the evolution of the quantum states of the central qubit and external qubits such that the quantum state of the central qubit and external qubits after a time period t F corresponds to evolution of the quantum states according to a target Hamiltonian after the time period t F . The sequence of single-qubit gates and time intervals between the application of single-qubit gates may be controlled in order to modify the evolution of the state of the central qubit and external qubits to correspond to the target Hamiltonian. When the stepwise DAQC paradigm is used, the method may further comprise: - before each single-qubit gate is applied to the central qubit and one of the external qubits, turning off interactions between the central qubit and the remaining external qubits; and - after each single-qubit gate is applied to the central qubit and one of the external qubits, turning interactions between the central qubit and the remaining external qubits back on. When the banged DAQC paradigm is used, the sequence of single-qubit gates and time periods between the application of single-qubit gates during which interactions between the central qubit and the external qubits are turned on may be controlled in order to modify the evolution of the state of the central qubit and external qubits to correspond to the target Hamiltonian. A single-qubit gate may be applied to each of the qubits q 2 to q N−1 before turning on the interactions between the central qubit and external qubits, and only one further single-qubit gate may be applied to each of the qubits q 2 to q N−1 . In both of the first and second aspects of the invention, the single-qubit gates may be involutory. The single-qubit gates may be Pauli X, Y or Z gates. The resource Hamiltonian and target Hamiltonian may be Ising Hamiltonians of the form where the operator A q is one of the Pauli matrices X, Y or Z operating on qubit q and where the resource and target Hamiltonians have different coupling coefficients g k . The coupling coefficients of the connections between the central qubit or resonator q 0 and external qubits q 1 to q N−1 may be tuneable, in which case the method may further comprise tuning the coupling coefficients in order to shorten or reduce the number of analog blocks. A third aspect of the invention relates to a computer system configured to perform the method of the first or second aspect described above. A fourth aspect of the invention relates to a computer system comprising a classical processing unit configured to perform the method of the first aspect and a quantum processing unit configured to perform the method of the second aspect. A fifth aspect of the invention relates to a computer program product comprising instructions which, when executed by a computer, cause the computer to perform the method of the first or second aspect. A sixth aspect relates to a computer-readable medium comprising instructions which, when executed by a computer, cause the computer to perform the method of the first or second aspect described above. Brief Description of the Drawings Figure 1 A “star” topology of qubits in a quantum processing unit according to the present invention. Figure 2 A quantum circuit diagram showing the placement of X gates in a DAQC algorithm on a star topology QPU. Figure 3 A quantum circuit diagram showing a reduced number of X gates in a DAQC algorithm on a star topology QPU. Figure 4 A quantum circuit diagram showing the banged DAQC algorithm corresponding to the stepwise DAQC algorithm of Figure 3. Figure 5 A flow diagram depicting a method for determining a control sequence for performing DAQC on a QPU with a star topology according to the present invention. Figure 6 A flow diagram depicting a method for performing DAQC on a QPU with a star topology according to the present invention. Detailed Description of the Invention In digital-analog quantum computation (DAQC) the “analog” evolution of the quantum states of qubits in a quantum processing unit (QPU) according to a “resource Hamiltonian” is modified by applying “digital” single-qubit gates (SQGs) to the qubits of the QPU such that the quantum states of the qubits after a time ^^ ^^ is the same, or at least approximates, that following evolution of the states under a different Hamiltonian, referred to as the “target Hamiltonian”. In this context, a resource Hamiltonian is an entangling Hamiltonian, denoted according to which qubits of the QPU that are used in the computation evolve naturally, when all qubit-qubit interactions between the qubits used in the computation are turned on. It will be appreciated that the present invention may be performed on a dedicated QPU including only the qubits required for the computation, or it may be performed on a QPU with additional qubits and qubit-qubit connections beyond those strictly required for the computation. In the latter case, the additional qubits remain unused and qubit-qubit interactions between pairs of unused qubits and between unused qubits and qubits used in the computation remain turned off. When, in the present disclosure, the phrase “turning on qubit-qubit interactions” or similar is used, it is intended to mean only turning on interactions between qubits used in the computation, unless otherwise specified. A target Hamiltonian is an entangling Hamiltonian, denoted H, that implements a specific, desired computation. Its coupling coefficients can be chosen arbitrarily, depending on the computation to be implemented. Under the stepwise DAQC (sDAQC) paradigm, a number of “analog blocks” are applied to the qubits, each analog block corresponding to a controllable time period t during which qubit-qubit interactions are turned on and the quantum states of the qubits are allowed to evolve under the resource Hamiltonian. Such analog blocks are therefore described by the unitary operator U analog : According to the present invention, the DAQC is performed on a QPU having qubits arranged in a “star” topology, which includes a “central” qubit, labelled q 0 and “external” qubits q 1 to q N −1 , where ^^ is the number of qubits used in the computation. An example of such a star topology is shown in Figure 1. The qubits are referred to as “central” and “external” due to their logical relationship and no inference should be made regarding their physical arrangement beyond the physical constraints implied by the requirement that each of the external qubits be connected to the central qubit. The central qubit may indeed be alternative referred to as a “hub” qubit, or a “common” qubit, and the external qubits may be referred to as “spoke” qubits, for example. The qubits depicted in Figure 1 may be physical qubits or logical qubits, where logical qubits are made up of one or more physical qubits or qubit-like elements, such as a resonator, that are capable of storing a quantum state. Where a logical qubit is made up of more than one physical qubit or qubit-like elements, there exist connections between the physical qubits or qubit like elements within the logical qubit, and the connections between the logical qubit and other logical or physical qubits are distributed between the physical qubits or qubit like elements of the logical qubit. Coupling large numbers of qubits in a quantum processing unit is technically challenging as the direct coupling of one qubit to many requires increasing capacitance or inductance, resulting in larger qubit sizes. Such constraints limits both the number of qubits that can be coupled and the number of qubits that can be formed on a single die or wafer. Thus, where a large number of external qubits are employed, it may be necessary to use multiple physical qubits or qubit-like elements as a logical central qubit. The central qubit ^^ 0 in particular may be a qubit-like element such as a resonator, since in the method of the present invention no single qubit gates need be applied to the central qubit. Throughout this disclosure, where the term “qubit” is used, it should be considered to encompass both physical and logical qubits, including those formed of qubit-like elements such as resonators and waveguides. For a star topology, such as that shown in Figure 1 a typical the resource Hamiltonian is an Ising Hamiltonian, such as a ZZ Ising Hamiltonian of the form: where g k is the coupling coefficient of the k th qubit with the central qubit in the resource Hamiltonian. The index 0 in Z 0 indicates that the Z operator acts on the central qubit q 0 and Z k operates on the k th qubit q k . Such a topology is particularly suited to target Hamiltonians that are also Ising Hamiltonians, such as a ZZ Ising Hamiltonian of the form: where g k is the coupling coefficient of the k th qubit with the central qubit in the target Hamiltonian. Again, the index 0 in Z 0 indicates that the Z operator acts on the central qubit q 0 and Z k operates on the k th qubit q k . The coupling coefficients of the resource Hamiltonian are determined by the physical properties of the qubits and qubit-qubit connections in the QPU and are therefore essentially fixed for QPUs without tuneable couplings. For QPUs with tunable couplings between the qubits, the accessible range of the coupling coefficients of the resource Hamiltonian is nevertheless determined by the physical properties of the qubits and qubit-qubit connections, including the properties of the tunable couplers. By implementing a number of single-qubit gates and time-controlled analog blocks, in the case of sDAQC, or time periods between SQGs, in the case of bDAQC, the coupling coefficients of the target Hamiltonian can be simulated. By interleaving suitable SQGs between the analog blocks, such that one SQG is performed before the analog block and one SQG is performed after the analog block, the sign of the coupling coefficient, i.e. positive or negative, can be modified. For example, the non-commutativity of X andZ gates allows X gates to be used as the SQGs for a ZZ Ising Hamiltonian, where the identity X a Z b X a = (−1) δ ab Z b and property Re it H R = e it RH R † leads to a unitary operator describing the evolution of a modified resource Hamiltonian over time t of the following form, for a ≠ 0: The unitary operator U DA QC describing evolution of the state of the system in which SQGs are interleaved between the analog blocks is: The effect of the four X gates in each term of the product in equation (5) is to flip the sign +/- of the term depending on the index m. The X gates in equation (5) can therefore be re-expressed as(−1) 1+ δ m +δ k+ δ mk , which can be viewed as a set of elements M mk , where M mk = −1 for m < k and M mk = +1 for m ≤ k: This set of elements M mk can be interpreted as a matrix M with elements M mk . Since the unitary operator corresponding to the target Hamiltonian can be written as: it can be seen that equations (6) and (7) are equal if wheret is the vectorized timest m and G is the vectorized ratio of coupling coefficients, such that each element of the vector G j ≡ g j ⁄ g . The run times t of the analog blocks required to implement the target Hamiltonian can be calculated by the inversion of matrix M: The individual times t k for the analog blocks A k can therefore be calculated according to the following equations: These equations always give real times t k ≥ 0 as long as G k ≥ 0 and G k ≥ G k+1 . To meet the first condition, it may be necessary to shift the phase of the evolution in equation (7), Φ k → Φ k ± 2 π, in order to change the signs of the coupling coefficients of the target Hamiltonian while keeping the evolution of the target Hamiltonian exactly the same. To meet the second condition, the elementsG k of G may need to be re-ordered. The placement of SQGs for the star topology shown in Figure 1 is given by a matrix M of the form: where each column in the matrix corresponds to an analog block, and each row corresponds to a coupling of the target Hamiltonian. If element M ab = +1, the coupling coefficient for the interactionZ 0 Z a for the b th analog block should be positive, and if M ab = −1, the coupling coefficient should be negative. For example, based on the matrix ^^ shown in equation (11), the coupling of q 0 and q 1 is always positive, therefore no X gates need to be performed on qubit q 1 . The coupling of q 0 and q 2 is negative for the first analog block, so X gates are performed before and after the first analog block on qubit q 2 , but the remaining analog blocks are performed without performing any more X gates on q 2 . This is shown in more detail in Figure 2, which shows the placement of X gates according to the matrix M shown in equation (11) on a quantum circuit diagram. Since X gates are involutory, i.e. XX = II, where two X gates are applied consecutively in Figure 2, they do not need to be performed. As a result, the circuit diagram of Figure 2 can be simplified, resulting in the quantum circuit diagram of Figure 3. While the above has been described in context of analog blocks where SQGs are applied before and after each block, i.e. in the context of stepwise DAQC, it will be appreciated that the same derivation is equally applicable to a context in which there are no gaps between analog blocks, i.e. without turning the interactions between qubits on or off, and SQGs are applied simultaneously with analog blocks, i.e. banged DAQC. Where bDAQC is used, a SQG performed “before” the analog block in the description above is performed at the beginning of the analog block, and a gate performed “after” an analog block is performed at the end of an analog block. The quantum circuit diagram of Figure 4 depicts a quantum circuit for bDAQC, where the times t k for each analog block referred to above are the time periods that elapse between the X gates. Since real-world quantum gates cannot be applied instantaneously, there is some length of time At required to perform the X gates. As a result, the time period between the start of the continuous analog block and the start of the first X gate in Figure 4 is — At/2, the time period between subsequent X gates is t k - Δt, and the time period between the final X gate and the end of the continuous analog block is t N-1 - Δt/2.. Alternatively to the circuit shown in Figure 4, the first X gates on qubits q 2 to q N-1 may also be applied simultaneously with the analog block.

However, some of the ZZ terms in the resource Hamiltonian act on qubit k and thus “overlap” with the X k term. Due to the non-com mutativity of the operators [Z a ,X a ] = 2iY a , this operation is not equal to the evolution of the analog block and the X gate independently. Each term Z j Z k thus introduces a mis-rotation that depends on the population of qubits j and k. The error introduced by the non-com mutativity of the operators, which is a necessary characteristic of the operators in order to flip the signs of the coupling coefficients, is proportional to the number of ZZ terms that overlap with the X a gate. In the circuit shown in Figure 4, the number of X gates applied simultaneously with the analog block is significantly reduced compared to a more general circuit required to perform the method on other qubit topologies. The error is therefore significantly reduced in the circuit shown in Figure 4 as a result of the star topology of the QPU compared to other QPU topologies.

The above has been described in the context of a ZZ Ising Hamiltonian comprising a sum of ZZ terms and single-qubit X gates; however, the invention is also applicable to other two-body resource and target Hamiltonians made up of sums of XX, in which case the SQGs may be Z gates, or YY terms, in which case the SQGs gates may be X gates, for example, since the gates follow the relationship Z a X b Z a = (-1) δab X b and X a Y b X a = (-1) δab Y b . The use of ZZ, XX or YY Hamiltonians therefore does not significantly change the derivation presented above and ultimately still leads to equation (9), which fully encodes the problems of determining the analog block run times and identifying a sequence of single-qubit gates corresponding to matrix M.

Furthermore, by selecting different SQGs to be applied to different qubits, resource and target Hamiltonians that are sums of mixed terms, i.e. sums of XT, YX, XZ, ZX, YZ and ZY terms, or even resource and target Hamiltonians that have different terms for different couplings between the central qubit and external qubits can be used with the method of the present invention, as long as the SQGs S performed before and after each analog block follow the relationship S a D b S a = (-1) δ a b D b with the operator D in the coupling of the central and external qubits. The large reduction in the number of SQGs shown in Figure 3 and Figure 4 can be achieved as long as the SQGs S are involutory. Based on the above, the present invention includes a method for determining a control sequence for performing DAQC on a QPU with a star topology. The control sequence may be a quantum circuit or equivalent, an intermediate representation such as OpenQASM, Alternatively, the control sequence may be an intermediate representation such as OpenQASM, at a hardware interaction level such as defining control signals to be applied to the qubits, or any suitable description of the parameters needed to perform the DAQC, such as defining the couplings to be turned on/off, timings and locations of single qubit gates. The method is depicted in the flow diagram of Figure 5. The method of Figure 5 is typically performed on a classical computer in order to determine the control sequence to be applied to a quantum computer. The method may therefore be performed in a hybrid classical-quantum computer system. At step 501, the target Hamiltonian and evolution time t F are received. The target Hamiltonian and evolution time t F may be input directly by a user, may be received from another computer program running on the same computer or on a networked computer, or may be retrieved from memory or from a network location. The target Hamiltonian may be received in equation form, similar to that identified above in equation (3) or (7). Alternatively, the characteristics of the target Hamiltonian, such as the coupling coefficients g k may be received. The coupling coefficients may be received in an array or list, for example [g 0 , g 1 , g 2 , …, g N-1 ], or any other suitable format. The resource Hamiltonian and/or its coupling coefficients g k may already be known, for example in a hybrid-classical computer system where the classic computer component already has access to key information defining the quantum system. Alternatively, the resource Hamiltonian and/or its coupling coefficients may also be received or retrieved by the method at step 501. At step 502, the method determines the times t 1 to t N−1 of the sequence of analog blocks A 1 toA N−1 . The times are calculated according to equations (10) above based on ratios of the coupling coefficients of the resource and target Hamiltonians. The determined times t k are used to determine the sequence of analog blocks A k to be applied to the qubits q 0 to q N−1 : the analog blocks A k are applied in sequence according to their index k, except for any analog blocks with time t k = 0, which do not need to be applied. At step 503, the sequence of SQGs to be applied in order to flip the signs of the coupling coefficients in the resource Hamiltonian is determined based on the determined sequence of analog blocks. The SQGs are applied to qubits q 2 to q N−1 such that single-qubit gates are applied to qubits q j +1 to q k −1 at after each analog block A j with non-zero duration t j and j < N − 1, and where A k is the next analog block with non-zero duration t k . This definition results in a placement of SQGs equivalent to that described above, following the process of applying SQGs to every qubitq j +1 to q N −1 before and after every analog block A j and subsequently removing pairs of consecutively applied SQGs. Such a process may be followed in order to determine the placement of SQGs. As described above, in order for consecutive pairs of SQGs to cancel, the SQGs should be involutory. Where the resource and target Hamiltonians comprise sums of two-qubit terms of Pauli X, Y and/or Z operators, the SQGs are themselves Pauli X , Y or Z operators. Step 503 also includes placing SQGs on all qubits q 2 to q N −1 before the first analog block. Following step 503, the full control sequence required to perform DAQC on the star topology QPU is known. At an optional set 504, the control sequence may be output. Outputting the control sequence may comprise providing the output in some machine- or human-readable format such as a quantum circuit or intermediate representation such as OpenQASM. ,Alternatively, the control sequence may be an intermediate representation such as OpenQASM, at a hardware interaction level such as defining control signals to be applied to the qubits, or any suitable description of the parameters needed to perform the DAQC, such as defining the couplings to be turned on/off, timings and locations of single qubit gates. Again, while the method of Figure 5 has been described in the context of stepwise DAQC, it will be appreciated that it applies equally to banged DAQC, in which case the determined sequence of analog blocks A 1 to A N −1 are applied consecutively, without turning the interactions between qubits on or off, a SQG performed “before” the analog block in the description above is performed at the beginning of the analog block, and a gate performed “after” an analog block is performed at the end of an analog block. Figure 6 is a flow chart depicting a method of performing a digital-analog quantum computation algorithm on a star topology QPU according to the present invention. At step 601, initial quantum states are prepared in the qubits q 0 to q N −1 . Preparing the initial quantum states may comprise preparing |0 states in each qubit and manipulating the qubit states by the application of quantum gates, such as Pauli X, Y and Z gates and the Hadamard gate H. Alternatively, preparing the initial quantum states may comprise performing another quantum computation algorithm, such as a quantum Fourier transform, such that the resulting quantum states of the previous quantum algorithm are stored in the qubits q 0 to q N −1 . At step 602, interactions between the central qubit q 0 and each external qubit q 1 to q N −1 are turned on such that the quantum states of the central qubit and external qubits evolve according to a resource Hamiltonian. At step 603, a sequence of SQGs is applied to the external qubits to modify the evolution of the quantum states of the central qubit and external qubits such that the quantum state of the central qubit and external qubits after a time period t F corresponds to evolution of the quantum states according to a target Hamiltonian after the time period t F . In the case of bDAQC, the sequence of single-qubit gates and the time intervals between the application of single-qubit gates are controlled in order to modify the evolution of the state of the central qubit and external qubits to correspond to the target Hamiltonian. The time intervals correspond to the analog blocks described above. In the case of sDAQC, interactions between the central qubit and the remaining external qubits are turned off before each single-qubit gate is applied to the central qubit and one of the external qubits and turned back on after each single-qubit gate is applied to the central qubit and one of the external qubits. The periods during which interactions are turned on correspond to the analog blocks described above. Where there are fewer analog blocks than external qubits q 1 to q N −1 , i.e. fewer than N − 2 analog blocks, some of the SQGs are applied parallel. This arises when the ratio between the coupling coefficients g k and g of the target and resource Hamiltonians for coupling k is equal to the ratio of coupling coefficients g k+1 and g k̅+1 for the subsequent coupling k + 1, resulting in an analog block with t k = 0. As described above, it will be understood that SQGs may be applied to qubits q 2 to q N−1 before the interactions between the central qubit and each external qubit are turned on. In the case of bDAQC, these SQGs may be applied before the interactions are turned on, or may be applied at the start of the time period for which the interactions are turned on. In all circumstances and for both sDAQC and bDAQC, only one further SQG is applied to each qubit q 2 to q N−1 due to the use of a QPU with a star topology. The description of the system and method above has generally been presented in the context of a quantum processing unit with fixed couplings between qubits; however, it will be appreciated that the method is equally applicable to the quantum processing unit with tunable couplings. Indeed, the use of tunable couplings can even provide advantages, such as allowing for shorter or fewer analog blocks with suitable adjustments to the coupling coefficients within the accessible range of values for those coefficients.