Method and system for solving problems with multiple conflicting objectives

Abstract

To provide a computer-implemented method for obtaining a set of Pareto-optimal solutions to a given multi-objective optimization problem with several different objectives.SOLUTION: The method comprises generating a set of parameters on a classical computation component using a classical optimization algorithm. Based on the set of parameters, a quantum component generates a quantum state by executing a parameterized quantum circuit. A set of basis states is selected from the one quantum state. The classical computation component calculates an objective values set consisting of objective values for each solution in the set of classical solutions, and a new estimate for the Pareto set of solutions and the corresponding Pareto front by the classical optimizer is generated. The internal state of the classical optimization algorithm is updated based on the current set of solutions and the Pareto set and their corresponding sets of objective values.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] The present invention relates to a computer-implemented method for obtaining a set of Pareto-optimal solutions for a given multi-objective optimization problem, and to a hybrid quantum-classical computer system for implementing such a method. [Background technology]

[0002] Quantum optimization is a very active area of ​​research. The current generation of quantum computing hardware devices is still limited in size and susceptible to noise. Currently, the usefulness of so-called NISQ (noisy intermedia scale quantum) devices for applied problems is being tested and evaluated against a number of real-world applications.

[0003] Current state-of-the-art approaches to solving optimization problems using quantum computers are hybrid quantum-classical algorithms, as described, for example, in "Noisy intermediate-scale quantum (NISQ) algorithms, Rev. Mod. Phys. 94, 015004." Examples are the variational quantum eigensolver (VQE) described in "The variational quantum eigensolver: a review of methods and best practices. Physics Reports, 986, 1-128"; "Quantum machine learning beyond kernel methods", arXiv preprint arXiv:2110.13162 (2021) as well as quantum approximate optimization algorithms (QAOA) such as "A Quantum Approximate Optimization Algorithm, arXiv:1411.4028 (2014)" or "From the quantum approximate optimization algorithm to a quantum alternating operator ansatz, Algorithms 12 (2019)". All of these methods consider problems where the goal is to find an optimal solution that minimizes a predetermined cost function. Technically, they aim to solve single-objective optimization tasks. What is still missing are equivalent quantum methods aimed at solving multi-objective optimization problems.

[0004] There are several methods for solving multi-objective problems using quantum hardware. In "Multiobjective variational quantum optimization for constrained problems: an application to Cash Management," arXiv:2302.04196, https: / / doi.org / 10.48550 / arXiv.2302.04196, a VQE quantum circuit and a combined classical NSGAII algorithm are applied. From the quantum state prepared by the quantum device, they classically calculate two objective functions: a cost function that, when minimized, achieves the actual goal, and the probability of finding a feasible solution at the final state. The goal is not to map the entire Pareto front, but to reach a single optimal solution point where the probability of finding a feasible solution is very close to 1 and the cost function is minimized. In essence, this method is a constraint-solving method that uses classical multi-objective optimization to efficiently explore a search space of feasible and infeasible solutions. The quantum aspects of quantum computing are not utilized to reach a solution to the multi-objective problem.

[0005] Although addressing these problems is naturally becoming more practical, one central feature of real-world problems remains largely underconsidered: Real-world problems often deal with multiple conflicting objectives, such as cost (lowest possible), quality (highest possible), user satisfaction (highest possible), CO2 footprint (lowest possible), etc. These objectives are usually in a trade-off relationship, where a solution with very good values ​​for one objective, e.g., high quality, is less favorable for another objective, i.e., cost, too. Thus, the goal of multi-objective optimization is to find a whole set of solutions (not just one) that gives a good approximation of the so-called Pareto front (PF), the best feasible set of trade-off solutions. Each solution in the PF is equal to or better than any other solution in the PF in at least one objective and worse than at least one other objective.

[0006] "Multi-Objective Routing Optimization for 6G Communication Networks Using a Quantum Approximate Optimization Algorithm" (Sensors 22(19), 7570 (2022) doi:10.3390 / s22197570) describes solving multi-objective problems by iteratively solving a large number of constrained single-objective problems. A quantum computer is used only to solve the single-objective problems. Each single-objective problem optimizes only one objective and imposes constraints on the other objectives. A quantum computer is not used to directly obtain a population of solutions.

[0007] There are several papers that utilize so-called non-dominated quantum optimization (NDQO) and its variants, such as multi-objective decomposition quantum optimization (MODQO) and non-dominated quantum iterative optimization (NDQIO), such as "Quantum-Assisted Routing Optimization for Self-Organizing Networks," IEEE Access 2, 614-632 (2014) doi: 10.1109 / ACCESS.2014.2327596, "Non-Dominated Quantum Iterative Routing Optimization for Wireless Multihop Networks," IEEE Access 3, 1704-1728 (2015) doi: 10.1109 / ACCESS.2015.2478793, or "Quantum-Assisted Joint Multi-Objective Routing and Load Balancing for Socially-Aware Networks," IEEE Access 4, 9993-10028 (2016) doi: 10.1109 / ACCESS.2016.2629671). All of these methods aim to find a Pareto-optimal solution set for a multi-objective problem. They do so by utilizing quantum computing devices. Specifically, they design quantum circuits that use auxiliary qubits to evaluate the Pareto-dominance relationships between individual solutions. A variant / extension of Grover search can then be used to search for new suitable candidate solutions, which can then be inserted into the current approximation of the Pareto front. These algorithms are specialized to solve one single-objective problem and evaluate the Pareto-dominance relationships between sets of solutions using complex circuits.They cannot be directly applied to general multi-objective optimization problems, and in particular cannot be run on NISQ devices, requiring fully fault-tolerant and error-corrected quantum hardware that will not be available in the near to medium future. Summary of the Invention [Problem to be solved by the invention]

[0008] However, there remains a need to directly address the optimization of multiple, conflicting objectives to enable the efficient computation of multiple optimized solutions to the problem. Current technology does not allow multi-objective optimization problems to be processed natively on quantum hardware devices.

[0009] It should be noted that in the following description, a "classical computing component" may be a classical computer, and a "quantum computing component" or "quantum component" may be a general-purpose or special-purpose quantum computing device. [Means for solving the problem]

[0010] The present invention addresses the above-mentioned problems and finds solutions to practical real-world problems known in particular from the area of ​​energy management using quantum-specific multi-objective quantum optimization techniques according to the methods and computer systems defined in the independent claims on current and future quantum hardware devices.

[0011] The dependent claims define advantageous aspects and features.

[0012] The present invention involves utilizing quantum superpositions, which are intrinsic to quantum states, as a resource for encoding the complete population of solutions for a multi-objective optimization algorithm. Quantum circuits are designed to increase the diversity of the population and improve all objectives so that the objectives provide a good approximation of the Pareto front.

[0013] The method and apparatus provide for, for a given multi-objective optimization problem with M distinct objectives:

[0014]

number

[0015]

number

[0016]

number

[0017]

number

[0018] At each iteration t, the following steps are performed: a. The classical computation component uses a classical optimization algorithm to parameterize γ t and sending it to the quantum component; b. In the quantum computing component, the parameter γ t By performing PQC using the set of

[0019]

number

[0020]

number

[0021]

number

[0022]

number

[0023]

number

[0024]

number

[0025]

number

[0026]

number

[0027]

number

[0028]

number

[0029]

number

[0030]

number

[0031]

number

[0032]

number

[0033] After completing the iterations, the final result is a good set of Pareto optimal solutions.

[0034]

number

[0035]

number

[0036] The present invention uses quantum superposition to encode the entire population of solutions into one quantum state. The quantum computing components are designed such that the layout of the quantum circuit collectively achieves a multi-objective evolution of Pareto-optimal solutions using true quantum encoding.

[0037] According to a preferred embodiment, the quantum component comprises N quantum components each having a different set of selectable parameters.L A unitary quantum operation U(γ) consisting of layers is

[0038]

number

[0039]

number

[0040]

number

[0041]

number

[0042]

number

[0043]

number

[0044] According to another preferred embodiment, the quantum component comprises N quantum components each having a different set of selectable parameters. L A unitary quantum operation U(γ) consisting of layers is

[0045]

number

[0046]

number

[0047]

number

[0048]

number

[0049]

number

[0050] The coefficients in these sums

[0051]

number

[0052]

number

[0053]

number

[0054] Each cost Hamiltonian H α is a function of one classical objective function C in the sense that the measurement basis state |x〉 is a simultaneous eigenstate of all cost Hamiltonians and the eigenvalues ​​of the problem Hamiltonian are given by the values ​​of the corresponding classical objective function for the corresponding classical solution x. α represents, i.e., H α |x〉=C α (x)|x〉 and the mixed unitary

[0055]

number

[0056]

number

[0057]

number

[0058] According to another preferred embodiment, the quantum component comprises N quantum components each having a different set of selectable parameters. L A unitary quantum operation U(γ) consisting of layers is

[0059]

number

[0060]

number

[0061]

number

[0062]

number

[0063]

number

[0064]

number

[0065]

number

[0066]

number

[0067]

number

[0068]

number

[0069]

number

[0070]

number

[0071]

number

[0072] More preferably, the mixing operator for each layer is given by a sum over all x angular momentum operators for each qudit variable, and the coefficient of each individual angular momentum operator is given by the freely selectable parameter

[0073]

number

[0074]

number

[0075] In an alternative embodiment, the mixing operator for each layer is given by a sum over all generalized Pauli X operators for each qudit variable, as described, for example, in Y. Wang, Z. Hu, B.C. Sanders, and S. Kais, "Quudits and High-Dimensional Quantum Computing, Frontiers in Physics 8 (2020)," with one coefficient for each physical qudit, which is given by the freely selectable parameter

[0076]

number

[0077]

number

[0078]

number

[0079] In a preferred embodiment, the quantum state |ψ γ 〉, i.e., the probability of the Born rule being realized by repeatedly measuring the quantum state, P(x,γ)=|〈x|ψ γ 〉| 2 Sample x from each measured state |x m 〉, the corresponding classical solution x m is a set of selected states

[0080]

number

[0081]

number

[0082] The method and apparatus may further be applied to solving a multi-objective bidirectional electric vehicle charging problem, where the goal is a discrete integer variable σ that encodes a full charge and service schedule for a fleet of electric vehicles.

[0083]

number

[0084]

number

[0085]

number

[0086] The method and apparatus according to one of the above-described configurations may be implemented by using an integer variable x={ω n,s}, this efficient frontier is spanned by the portfolio configuration that has the optimal trade-off between maximizing return and minimizing risk, thus presenting a two-objective optimization problem.

[0087] The invention will now be further described with reference to the accompanying drawings. [Brief explanation of the drawings]

[0088] [Figure 1] 1 is a schematic diagram of a hybrid quantum-classical computer system illustrating method steps according to a preferred implementation of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0089] The proposed computer-implemented method and hybrid quantum-classical computer system aim to solve combinatorial integer multi-objective optimization problems, where there are two or more objectives to minimize. Typically, these objectives cannot be simultaneously minimized by several optimal solutions because they are in trade-off relationships with each other. Therefore, the goal of multi-objective optimization methods is to find the entire set of distinct solutions, called the Pareto set, where each solution achieves one optimal trade-off between the objectives. Therefore, the set of solutions must map all possible trade-offs between the objectives, called the Pareto front.

[0090] The proposed method and apparatus of the present invention allows such problems to be addressed natively on current NISQ quantum hardware. A crucial aspect of the present invention is provided by representing the entire population of candidate solutions to an optimization problem as a single quantum state and extracting the entire set of Pareto-optimal solutions from this single quantum state. In accordance with the present invention, the proposed computer-implemented method and quantum circuit of the hybrid quantum-classical computer system increase the diversity of the population by utilizing all cost function Hamiltonians corresponding to all classical objectives.

[0091] According to this description, the term "classical" always refers to known computations that were performed in the past using only digital computer systems.

[0092] The proposed computer-implemented method and hybrid computer system comprises a classical computing component and a quantum computing component, both configured to send and receive information to and from each other so as to achieve a shared computation for solving a multi-objective integer optimization problem. N ) (where a i∈[0,...,L-1] is a finite integer) is represented by the degrees of freedom of a qubit (for L=2) or qudit (for L>2) in a quantum hardware device, i.e., |x〉=|α1,α2,...,α N 〉 is.

[0093] Each of these states is a basis state of the measurement basis and satisfies the orthonormality condition:

[0094]

number

[0095]

number

[0096] This encoding means that all possible solutions to the original multi-objective integer optimization problem are represented by the measurement basis states of the quantum computation component. One advantage of using a quantum computer to tackle this problem is that a general quantum state is naturally a superposition of all basis states, i.e., a superposition of all possible solutions to the problem.

[0097]

number

[0098]

number

[0099] The present invention utilizes parameterized quantum circuits in quantum computing components, which use unitary gate operations applied to a predetermined fixed initial state |ψ0〉. Such application of unitary gate operations results in a final quantum state |ψ〉=U(γ)|ψ0〉 This leads to:

[0100] According to the present invention, this state is S Encode the entire population of candidate solutions;

[0101]

number

[0102]

number

[0103]

number

[0104]

number

[0105]

number

[0106] State selection mechanisms, such as full-state tomography, classical shadowing, or measurement-based sampling, allow the entire set of solutions to be extracted from the quantum state.

[0107]

number

[0108] Suitable techniques are described in "Efficient quantum state tomography" (Nat. Commun 1, 149 (2010) doi 10.1038 / ncomms114), "Predicting many properties of a quantum system from very few measurements" (Nat. Phys. 16, 1050-1057 (2020) doi 10.1038 / s41567-020-0932-7), or "Quantum Computation and Quantum Information" (Cambridge: Cambridge University Press. (2010) doi 10.1017 / CBO9780511976667). The unitary of the present invention is specifically designed to increase the contributions to the superposition from different states, which gives an optimal trade-off solution. That is, the population of candidate solutions represented by a quantum state is divided into a Pareto set of solutions.

[0109]

number

[0110]

number

[0111] In the proposed invention, the set of solutions is sent from the quantum computing component to the classical computing component, where for each solution a vector of objective values ​​is calculated, which is then used to find the Pareto frontier.

[0112]

number

[0113]

number

[0114]

number

[0115] The final result after running the proposed method using the hybrid quantum-classical computing system of the present invention is a Pareto set of solutions.

[0116]

number

[0117]

number

[0118] In the following, the invention is described in two different implementation examples that are typical for multi-objective problems:

[0119] The operation of a fleet of electric vehicles (EVs) constitutes a typical optimization problem with several challenges: a central aspect of this operation is to serve customers in a way that satisfies customer demand for all vehicles, while at the same time reducing the electricity costs for charging all EVs as much as possible.

[0120] In one formulation of this problem, the goal is to T N over time frames EV The goal is to find an optimal charging schedule for a fleet of electric vehicles (EVs). The search variable x defines the charging schedule for the entire fleet of EVs. The schedule is EV N T consists of discrete integer variables,

[0121]

number

[0122]

number

[0123]

number

[0124]

number

[0125]

number

[0126] The first objective (OBJ1) represents the total electricity cost for the entire charging schedule. PV(s) represents the excess energy produced by the photovoltaic (PV) system in time window s that is available for charging the EV, and the function P(s) gives the price for buying or selling energy from or to the grid in time window s. P max represents the maximum unit cost of electricity. If the total electricity demand is greater than the PV surplus, energy must be purchased from the grid; if it is less than the surplus, the sign of the term in brackets is reversed and energy is sold to the grid, resulting in a negative cost, i.e., income.

[0127] The second objective (OBJ2) is to calculate the energy required by each vehicle, E tot,n This represents a measure of fairness between different EVs / customers by minimizing the difference between the charging energy provided to each EVn and its required value E tot,n This is because the difference between the variances is proportional to the variance of the difference between the variances, which strongly favors equal variances.

[0128] The third objective (OBJ3) seeks to use only excess PV energy to charge EVs in each time slot s, so that this objective function seeks to minimize the difference between the excess PV energy used in each time slot and the current total charging energy. Minimizing this objective is a means of making the service independent from the power grid.

[0129] The fourth objective (OBJ4) aims at an equal distribution of charging power over time frames, and in particular tries to reduce the peaks of charging power over all time frames. The rationale behind this is two-fold: (i) from a company's perspective, the peak power over a billing period has a direct impact on the electricity price charged, and therefore strong demand peaks need to be avoided; and (ii) strong fluctuations in electricity demand generally put a lot of pressure on energy grid providers, which they need to balance. Both aspects are met by reducing the overall variation, or variance, of the total energy demand over time, as formulated in this objective.

[0130] The problem characterized by the example objectives (OBJ1) through (OBJ4) described above defines a combinatorial integer multi-objective optimization problem with M=4 objectives, all of which must be minimized. However, the present invention is not limited to four objectives; fewer or more objectives may define a multi-objective optimization problem. These objectives typically trade off against one another and cannot be simultaneously minimized by several optimal solutions. Therefore, the goal of the multi-objective optimization technique of the present invention is to provide a computationally efficient means for finding the entire set of different solutions, called the Pareto set, where each solution achieves one optimal tradeoff between the objectives. Therefore, the set of solutions must map all possible tradeoffs between the objectives, called the Pareto front.

[0131] For the EV charging problem, the discrete search variables are the variables

[0132]

number

[0133] For the EV charging problem, the quantum formulation is as follows: each qudit in the quantum computation component is assigned to one time window s∈[1,...,N T ], n∈[1,...,N EV ] one charging power level l n,s [0,...,L-1]. The quantum state of the measurement basis is

[0134]

number

[0135]

number

[0136]

number

[0137] The M classical objective functions C1(x) through C4(x) defined above are represented by M cost Hamiltonians whose functional forms are given by quadratic polynomials in the search variables. This takes the classical objective functions C1(x) through C4(x) and calculates the search variables l that specify the charge level of EVn in time window s. n,s the corresponding z angular momentum operator L of the qudit that encodes the corresponding variables z;n,sby replacing

[0138]

number

[0139] The measurement basis state is given by the z angular momentum operator L z;n,s and the eigenvalues ​​are the charge levels, i.e.,

[0140]

number

[0141] This gives us an explicit form for the corresponding cost Hamiltonian representing the objective function:

[0142]

number

[0143] The measurement ground state is also a simultaneous eigenstate of these cost Hamiltonians (the latter being L z;n,s(since they are polynomials in ), the eigenvalues ​​are given by the values ​​of the objective function of the corresponding configuration.

[0144]

number

[0145] A quantum computation component that runs a parameterized quantum circuit (PQC) generates a state

[0146]

number

[0147]

number

[0148]

number

[0149]

number

[0150]

number

[0151]

number

[0152] The mixed unitary is N P N EV N T freely selectable parameters

[0153]

number

[0154]

number

[0155]

number

[0156] N in each layer P Each of the phase terms is a freely selectable parameter

[0157]

number

[0158]

number

[0159] Coefficient a k,a ≧0 is the reference vector

[0160]

number

[0161]

number

[0162]

number

[0163]

number

[0164]

number

[0165] In a simple but typical case, a quantum circuit has only one layer and N L =1, in this case the quantum circuit prepares the following state:

[0166]

number

[0167]

number

[0168] The claimed method executed on the claimed hybrid quantum-classical computer system includes the following steps: A classical optimization algorithm COA, such as stochastic gradient descent (SGD), is launched on the classical computation component. The search parameters of this classical optimizer are the 10 parameters of PQC, namely:

[0169]

number

[0170] Additionally, an archive A may be initiated in the classical computation component, where the search parameters, candidate solutions, and corresponding values ​​of the objective function are stored. Furthermore, the quantum device is configured such that the measurement basis represents the classical solution as shown in (ENC), and the quantum circuit generates a parameterized quantum state (PQCEX) using the corresponding gate operations (MIXEX), (PQCEX), and (HAM1)-(HAM4). Then, an iteration counter is initialized to t=0, and the following steps are repeated: i) based on its internal state, a classical optimizer algorithm generates a set of search parameter values;

[0171]

number

[0172] ii) The quantum device converts the parameter set from the classical computation component

[0173]

number

[0174]

number

[0175] The quantum device then executes the quantum circuit, and the resulting quantum state |Ψ(γ t )〉 is N S Encode the entire population of candidate solutions as a superposition,

[0176]

number

[0177]

number

[0178]

number

[0179]

number

[0180]

number

[0181]

number

[0182] One way to practically extract these solutions from a quantum state can be by repeated sampling, where a quantum state is prepared and the same parameter γ t measured multiple times (i.e., N measure >N S ), where each measurement leads to one ground state, i.e., one solution

[0183]

number

[0184] Another method is to perform a full-state tomography protocol or classical shadow tomography (see, for example, “Efficient quantum state tomography”, Nat Commun 1, 149 (2010), doi 10.1038 / ncomms1147, or “Predicting many properties of a quantum system from very few measurements”, Nat. Phys. 16, 1050-1057 (2020). https: / / doi.org / 10.1038 / s41567-020-0932-7) and explicitly select the ground state with the highest probability.

[0185] In either case, the procedure is to use the present invention to allow a quantum computing component to generate a quantum state |Ψ(γ t )〉 to the entire population of solutions

[0186]

number

[0187] iii) The classical computation component receives a population of solutions from the quantum computation component.

[0188]

number

[0189]

number

[0190]

number

[0191]

number

[0192] From this set of non-dominated solutions, the hypervolume index

[0193]

number

[0194]

number

[0195] Fitness score f t is then fed back into a classical optimization algorithm, COA←f(γ t ) The optimizer updates its internal state accordingly.

[0196] In addition, the search parameter γ t , the corresponding fitness value f(γ t ), the current estimate of the Pareto front

[0197]

number

[0198]

number

[0199]

number

[0200] iv) A convergence criterion is then checked in the classical computation component, e.g., whether the maximum allowed number of cost function evaluations has been reached, or the fitness value f(γ t ) was not enough compared to a given threshold. a. If the convergence criterion is met, the optimization is stopped and as a final result, the final Pareto front and corresponding Pareto set are obtained by selecting all non-dominated solutions from all approximations of the Pareto front at each iteration.

[0201]

number

[0202] The final result is a Pareto set of solutions

[0203]

number

[0204]

number

[0205] Because all possible solutions typically exist in a superposition of the initial states of the quantum device, the proposed method and hybrid computer system utilizing PQC are effective in generating a good approximation of the Pareto front. The PQC phase operator assigns a different phase to each individual basis state based on the value of the reference vector used in the phase operator in equation (PQCEX). In this way, a classical optimizer can adjust the PQC parameters to provide the greatest diversity and include mostly non-dominated solutions in the quantum state, thereby maximizing the hypervolume index.

[0206] The set of Pareto-optimal solutions serves as the basis for realizing efficient control of EV charging and scheduling services. A solution selection routine may be implemented, in which the best-fit solution is automatically determined from the entire set of Pareto-optimal solutions. A natural solution selection routine utilizes current-art knee point detection methods, for example, as described in "Performance Analysis on Knee Point Selection Methods for Multi-Objective Sparse Optimization Problems" (IEEE Congress on Evolutionary Computation, 1-8 (2018) doi 10.1109 / CEC.2018.8477915) or "A Multi-objective Evolutionary Algorithm for "Finding Knee Regions Using Two Localized Dominance Relationships" (IEEE Transactions on Evolutionary Computation 25(1), 145-158 (2021) doi 10.1109 / TEVC.2020.3008877). The knee point naturally provides the best trade-off solution that is preferable.

[0207] A second application example concerns Markowitz portfolio optimization, whose goal is to manage a portfolio of financial assets over time so as to maximize expected returns while simultaneously minimizing risk (see, for example, "Computational aspects of alternative portfolio selection models in the presence of discrete asset choice constraints," Quantitative finance 1(5), 489-501 (2001), "Solving the Optimal Trading Trajectory Problem Using a Quantum Annealer," IEEE Journal of Selected Topics in Signal Processing 10(6), 1053-1060, (2016) doi 10.1109 / JSTSP.2016.2574703, or "Dynamic portfolio optimization with real datasets using quantum processors and quantum-inspired tensor networks," Phys. Rev. Research 4, 013006 (2022) for recent introductions). Since there is usually a trade-off between these two objectives, it is desirable to map the efficient frontier, i.e., the Pareto front of the best possible trade-off solutions.

[0208] Technically, the problem is solved by solving a set of problems over a fixed time window, s (s=1,...,N T ) in asset n(n=1,...,N Asset ) ratio ω n,sThe objective of the present invention is to determine the composition of a portfolio of financial assets, denoted by ω. Typically, rebalancing, or changing, the portfolio composition is done over a fixed, predetermined number of times, resulting in these discrete time steps. In addition, funds typically trade assets in discrete chunks, resulting in the ratio ω n,s is expressed in terms of integer variables, i.e., ω n,s ∈[0,...,L-1].

[0209] The two objectives of this problem are

[0210]

number

[0211] Thus, the first term gives the cumulative total expected return for the entire portfolio over the time horizon considered.

[0212] Among the second objectives, Σ nm,s represents the covariance of assets n and m at time s, and for this purpose represents the covariance of expected portfolio returns summed over all times. This provides an indication of risk because (a) a large covariance suggests that true returns may differ significantly from expected returns, whereas a small covariance suggests that true returns are more similar to expected returns, and (b) a large covariance between two assets suggests correlation between the assets, which in turn reduces the robustness of the portfolio and therefore makes it riskier.

[0213] This dual-objective portfolio optimization problem is solved by dividing the portfolio by the proportion ω n,s The integer variables that describe the measurement basis states of the quantum device are

[0214]

number

[0215]

number

[0216]

number

[0217] This results in a cost Hamiltonian of

[0218]

number

[0219] The measurement basis states are also simultaneous eigenstates of these cost Hamiltonians, with eigenvalues ​​given by the values ​​of the objective functions of the corresponding configurations.

[0220]

number

[0221] In this example, the quantum computing component running the parameterized quantum circuit PQC uses two layers to prepare the state: |ψ〉=U(γ2)U(γ1)|ψ0〉 (PQC) Here, the unitary operation U can be freely parameterized by a set of parameters γ = {γ1, γ2}. The initial state of the circuit is given by an equal superposition of all possible states,

[0222]

number

[0223]

number

[0224]

number

[0225]

number

[0226]

number

[0227] The mixed unitary is a freely selectable parameter

[0228]

number

[0229]

number

[0230]

number

[0231] Each of the two phase terms in each layer is determined by one freely selectable parameter

[0232]

number

[0233] Putting this together, the quantum circuit prepares the following state:

[0234]

number

[0235]

number

[0236] The proposed method and apparatus performs the following steps:

[0237] First, a classical optimization algorithm, COA, such as stochastic gradient descent (SGD), is initiated in the classical computation component. The search parameters of this classical optimizer are the eight parameters of PQC, i.e., D=8.

[0238]

number

[0239] Second, the quantum computation component is configured such that the measurement basis represents the classical solution as shown in (ENC), and the quantum computation component generates a parameterized quantum state (PQCEX) using (HAM1) and (HAM2) along with the corresponding gate operations (MIXEX) and (PQCEX). Then, an iteration counter is initialized to t=0, and the following steps are repeated: i) based on its internal state, a classical optimizer algorithm generates a set of search parameter values;

[0240]

number

[0241] ii) The quantum computation component is a parameter set

[0242]

number

[0243]

number

[0244] The quantum device then executes the quantum circuit, and the resulting quantum state |Ψ(γ t )〉 is N S The entire population of candidate solutions is encoded as a superposition.

[0245]

number

[0246] For this state, a full state tomography (see, for example, "Efficient quantum state tomography", Nat Commun 1, 149 (2010) doi 10.1038 / ncomms1147) is performed on the quantum device, thereby obtaining a list of all weights contributing to that state, T={|ξ t (x c )| 2 c is 1,...d H}, and select the solution with the largest amplitude from T to N at iteration t. S A set of solutions

[0247]

number

[0248]

number

[0249] iii) The classical computation component receives the population of solutions from the quantum computation component.

[0250]

number

[0251]

number

[0252] From this set of objective vectors, a set of non-dominated solutions is

[0253]

number

[0254]

number

[0255]

number

[0256]

number

[0257]

number

[0258] The fitness values ​​are then fed back into a classical optimization algorithm, COA←f(γt ) The optimizer updates its internal state accordingly.

[0259] In addition, the search parameter γ t , the corresponding fitness value f(γ t ), the current estimate of the Pareto front

[0260]

number

[0261]

number

[0262]

number

[0263] iv) A convergence criterion is then checked in the classical computing device, e.g., whether the maximum allowed number of cost function evaluations has been reached, or whether the fitness value f over the last few iterations has been reached. t It is possible that the improvements were not sufficient. a. If the convergence criterion is met, the optimization is stopped and as a final result, a final Pareto front and corresponding Pareto set are obtained by selecting all non-dominated solutions from all approximations of the Pareto front at each iteration.

[0264]

number

[0265] The final result is a Pareto set of solutions

[0266]

number

[0267]

number

[0268] Because all possible solutions typically exist in a superposition of the initial states of the quantum device, the proposed method utilizing PQC and a hybrid quantum-classical computer system is effective in generating a good approximation of the Pareto front. The PQC phase operator assigns a different phase to each individual basis state based on the value of the scalarization function used in the phase operator in equation (PQCEX). In this way, the classical optimizer can adjust the PQC parameters to give the greatest diversity, allowing the quantum state to contain mostly non-dominated solutions, thereby maximizing the hypervolume index.

Claims

1. The set of Pareto-optimal solutions for a given multi-objective optimization problem with M distinct objectives [Equation 1] , wherein x represents a classical solution to the multi-objective optimization problem, the method comprising performing an iterative procedure, wherein each iteration t a. The classical computation component uses a classical optimization algorithm to find the parameter γ ｔ and sending it to the quantum component; b. In the quantum component, the initial quantum state |ψ ０ 〉, the parameter γ ｔ Quantum state [Equation 2] generating the quantum state [Equation 3] N Ｓ a step, which represents a set of distinct solutions; c. The one quantum state [Equation 4] From N Ｓ ≧N ＰＦ a set of ground states [Equation 5] and sending information about which ground state was selected to the classical computation component; d. computing, with the classical computation component, a set of objective values ​​comprising an objective value for each solution of a set of classical solutions; e) taking the set of solutions and the corresponding set of objectives with a classical optimizer to generate a new estimate of a Pareto set of solutions and a corresponding Pareto front; f. updating the internal state of the classical optimization algorithm based on the current set of solutions and the Pareto set and their corresponding sets of objective values; g. checking a convergence criterion, and if said convergence criterion is not met, repeating the process from step a onwards until said convergence criterion is met; A computer-implemented method in which the results after the most recent iteration is completed are taken as a set of Pareto-optimal solutions, and for this set of Pareto-optimal solutions, the corresponding set of objective value vectors is an approximation of the true Pareto front.

2. 2. The method of claim 1, wherein, in addition to updating the internal state of the classical optimization algorithm, an archive is updated based on the current solution set and the Pareto set and their corresponding sets of objective values, and the step of generating new estimates further incorporates values ​​from a previous iteration obtained from an archive.

3. The parameterized quantum circuit PQC is a set of N quantum circuits each having a different set of selectable parameters. Ｌ A unitary quantum operation U(γ) consisting of layers is [Equation 6] Implement it as Each layer l has a separately selectable parameter set γ ｌ = {γ Ｘ；ｌ , γ Ｐ；ｌ }, and a topological operator, each of which has a unitary [Equation 7] and Topological unitary of each layer [Equation 8] is a freely selectable parameter [Equation 9] , M, and the cost Hamiltonian H α and each cost Hamiltonian H α is a function of one classical objective function C in the sense that the measurement basis state |x〉 is a simultaneous eigenstate of all the problem Hamiltonians, and the eigenvalues ​​of the problem Hamiltonians are given by the values ​​of the corresponding classical objective function of the corresponding classical solution x. α represents Mixed unitary for each layer [Equation 10] is a freely selectable parameter γ Ｘ；ｌ，α , M parameterized mixing operators X α (γ Ｘ；ｌ，α 3. The method of claim 1, wherein the cost Hamiltonian is generated by:

4. The parameterized quantum circuit PQC is a set of N quantum circuits each having a different set of selectable parameters. Ｌ A unitary quantum operation U(γ) consisting of layers is [0011] Implement it as Each layer l has a separately selectable parameter set γ ｌ ∈{γ Ｘ；ｌ , γ Ｐ；ｌ }, and a plurality of mixing operators and topological operators, [0012] and Topological unitary of each layer [0013] is a freely selectable parameter [0014] and the cost Hamiltonian H α A specific sum over [Equation 15] and each cost Hamiltonian H α is a function of one classical objective function C in the sense that the measurement basis state |x〉 is a simultaneous eigenstate of all cost Hamiltonians and the eigenvalues ​​of the problem Hamiltonian are given by the values ​​of the corresponding classical objective function of the corresponding classical solution x. α represents Mixed unitary for each layer [0016] is a freely selectable parameter γ Ｘ；ｌ，ｋ , k=1,...,N Ｐ M parameterized mixing operators X ｋ (γ Ｘ；ｌ，ｋ 3. The method of claim 1, wherein the cost Hamiltonian is generated by:

5. The coefficients in these sums [Equation 17] is the reference vector [Equation 18] and the reference vector is [Equation 19] and α ｋ，ｎ ≥ 0, and for each k, define a specific direction in the space of objective values, thereby reflecting one user preference. First, define the problem: [Equation 20] 5. The method of claim 4, wherein the value is determined as part of

6. The parameterized quantum circuit PQC is a set of N quantum circuits each having a different set of selectable parameters. Ｌ A unitary quantum operation U(γ) consisting of layers is [Equation 21] Implement it as Each layer l has a separately selectable parameter set γ ｌ ∈{γ Ｘ；ｌ , γ Ｖ；ｌ }, and a plurality of mixing operators and topological operators with [Equation 22] and Topological unitary of each layer [Equation 23] is a freely selectable parameter [0000] and generated by an operator V that estimates a hypervolume in the object space spanned by the population of states when encoded by a quantum superposition of states, function [Equation 25] is the set of solutions [Equation 26] The objective value vector of [0000] Given the hypervolume spanned by the function [0000] has a function value between zero and one and allows downscaling of the hypervolume of a particular set of solutions if the individual amplitudes of some of the contributing states are smaller than a threshold or the heterogeneity of the amplitude distribution exceeds a threshold; Each cost Hamiltonian H α is a function of one classical objective function C in the sense that the measurement basis state |x〉 is a simultaneous eigenstate of all cost Hamiltonians and the eigenvalues ​​of the problem Hamiltonian are given by the values ​​of the corresponding classical objective function of the corresponding classical solution x. α represents Mixed unitary for each layer [0000] is at least one freely selectable parameter γ Ｘ；ｌ and a parameterized mixing operator X ｌ (γ Ｘ；ｌ 3. The method of claim 1, wherein the cost Hamiltonian is generated by:

7. The choice of the ground state is the quantum state [Equation 30] The sampling is realized by repeatedly measuring ｍ 〉, the corresponding classical solution x ｍ is a set of selected states [Equation 31] To, N Ｓ Measuring different states of [Equation 32] Until Ｓ 3. The method according to claim 1, wherein the method is carried out by adding different states.

8. The choice of the ground state is the quantum state [Equation 33] , where the set of weights for each basis state in the full-state tomography is [Equation 34] is determined, and N Ｓ The basis states corresponding to the subsets with the largest weights are selected, and these are called the set of selected states [Equation 35] 3. The method according to claim 1 or 2, characterized in that:

9. The mixing operator for each layer is given by a sum over all x angular momentum operators or generalized Pauli X operators for each physical qudit, and the coefficient of each individual angular momentum operator for each layer is given by one freely selectable parameter [Equation 36] 4. The method of claim 3, wherein k / α enumerates the mixing operators, l denotes the sheave, and i enumerates all physical qudits.

10. the mixing operator of each layer is given by a sum over all x angular momentum operators or generalized Pauli X operators of each physical qudit, and there is only one global coefficient for the complete sum over all individual angular momentum or generalized Pauli X operators that apply to each physical qudit of each layer, and the global coefficient is determined by the freely selectable parameters [Equation 37] 4. The method of claim 3, wherein k / α denotes a blending operator and l denotes a layer.

11. Reference vector of the phase Hamiltonian [Equation 38] Coefficient α of ｋ，α is determined, and said coefficient α is adjusted so as to ensure that a different set of reference vectors is used in each iteration, thereby enabling a speedup of the optimization process and a diversification of the Pareto solutions obtained. ｋ，α 5. The method of claim 4, wherein at each iteration step t, is adapted according to state-of-the-art procedures for reference vector adaptation.

12. The problem is a multi-objective bidirectional electric vehicle charging problem, where the goal is to find a set of Pareto-optimal configurations of discrete integer variables that encode full charging and service schedules for a fleet of electric vehicles, and the Pareto front [Number 39] 3. The method according to claim 1 or 2, characterized in that, to achieve an approximation of (i) above, multiple objectives are taken into account, being at least two of the following: cost of electricity to meet all charging demands, customer satisfaction, maximizing the use of renewable energy, minimizing electricity leaving the power grid, and minimizing peak power consumption.

13. 3. The method of claim 1, wherein the method is applied to obtain an efficient frontier of possible portfolio configurations, as formulated by the Markowitz model of modern portfolio theory, by finding the configuration of a portfolio over time, specified by an integer variable specifying the amount of investment in asset n at time s, wherein the efficient frontier is spanned by portfolio configurations that have an optimal trade-off between maximizing return and minimizing risk, and thus present a two-objective optimization problem.

14. A hybrid quantum-classical computer system, Quantum components configured to realize a parameterized quantum circuit PQC by implementing a unitary quantum operation U(γ), the unitary quantum operation U(γ) being parameterized by a freely selectable parameter γ, the unitary quantum operation U(γ) being parameterized by a predetermined initial state |ψ of N qubits or qudits. ０ By applying a unitary operator to |ψ, the quantum state |ψ of N qubits or qudits γ 〉, i.e., |ψ γ 〉 = U(γ)|ψ ０ 〉, which can be decomposed into a measurement basis for the quantum device. [Equation 40] , d Ｈ is the dimension of the Hilbert space, and each measurement basis state |x Ｃ 〉 is one classical solution x Ｃ a quantum component encoding the a classical computing component that provides a new parameter value γ to the quantum computing component and executes the PQC and determines the state |ψ γ and a classical computation component configured to perform a classical optimization algorithm by initiating the preparation of a classical solution set from the quantum module and then receiving a set of classical solutions as a result and computing all cost functions for each solution; A hybrid quantum-classical computer system in which both components are used together to solve classical multi-objective optimization problems by iteration as defined in the method of claim 1 or 2.

Images