---
title: 'Hamiltonian-based QAOA: Advances & Applications'
url: https://www.emergentmind.com/topics/hamiltonian-based-quantum-approximate-optimization-algorithm-qaoa
type: topic
---

# Hamiltonian-based QAOA: Advances & Applications

Hamiltonian-based Quantum Approximate Optimization Algorithm (QAOA) denotes the class of alternating-operator variational methods in which both the phase-separation family and the mixer family are generated by Hamiltonians. In the canonical formulation, a classical objective function \(f\) is encoded as a diagonal operator \(H_f\) through \(H_f\ket{\mathbf{x}}=f(\mathbf{x})\ket{\mathbf{x}}\), and a depth-\(p\) state is prepared by alternating the corresponding cost evolution with a mixer evolution, typically from the initial state \(\ket{+}^{\otimes n}\). Within the literature, this Hamiltonian picture has been extended well beyond the original MaxCut setting: to higher-order cost Hamiltonians, custom phase operators, feasible-subspace-preserving mixers, counterdiabatic corrections, local-Hamiltonian state preparation, continuous-variable optimization, and even computational universality [1709.03489, 1812.11075].

## 1. Formal definition and canonical operator structure

In the terminology of "From the Quantum Approximate Optimization Algorithm to a Quantum Alternating Operator Ansatz" [1709.03489], Hamiltonian-based QAOA is the subclass in which the phase separator and mixer are both Hamiltonian evolutions,
\[
\phaseUnitary(\gamma)=e^{-i\gamma H_P},\qquad \mixerUnitary(\beta)=e^{-i\beta H_M},
\]
with the output state
\[
\ket{\boldsymbol \beta,\boldsymbol \gamma}
=
e^{-i \beta_p H_M} e^{-i \gamma_p H_P}\cdots e^{-i \beta_1 H_M} e^{-i \gamma_1 H_P}\ket{s}.
\]
The same work also distinguishes local Hamiltonian-based QAOA, in which the mixer Hamiltonian is a sum of polynomially many local terms [1709.03489].

For unconstrained binary optimization, the standard mixer is the transverse field,
\[
H_M=-\sum_{i=1}^n X_i,
\]
whose ground state is \(\ket{+}^{\otimes n}\), prepared by Hadamards on \(\ket{0}^{\otimes n}\). For MaxCut, a canonical cost Hamiltonian is
\[
H_{MC}=\sum_{(v_i,v_j)\in E}\frac{1}{2}\left(I-Z_i Z_j\right),
\]
or equivalently the usual Ising/QUBO form with pairwise \(ZZ\) couplings and optional local \(Z\) fields [2511.12379, 2511.18377].

Hamiltonian construction is not restricted to quadratic objectives. The tutorial "An Introduction to the Quantum Approximate Optimization Algorithm" extends the mapping from QUBO to PUBO, yielding cost Hamiltonians of the form
\[
\mathcal{H}_{\text{f}}=\sum_{k=1}^d\sum_{i_1<\cdots<i_k} a_{k,i_1\ldots i_k}\,\sigma_z^{i_1}\cdots \sigma_z^{i_k},
\]
with corresponding \(R_{Z^k}\) gate decompositions for higher-order terms [2511.18377]. The graph-coloring study "QAOA of the Highest Order" makes the same point operationally: gate-model QAOA can directly exploit higher-order Pauli gadgets instead of forcing the objective into a quadratic annealing-style encoding [2111.12754].

## 2. Discretized annealing, integrated resources, and schedule physics

A standard interpretation of Hamiltonian-based QAOA is as a discretization of adiabatic evolution. "Quantum Optimization Algorithms" writes the continuous interpolation as
\[
H(t)=(1-s(t))H_{\text{initial}}+s(t)H_{\text{problem}},
\]
and approximates the time-ordered evolution by Suzuki–Trotter splitting,
\[
U(T)\approx \prod_{k=0}^{p-1}\exp\!\left(-i(1-s(t_k))H_{\text{initial}}\delta\right)\exp\!\left(-i\,s(t_k)H_{\text{problem}}\delta\right),
\]
with \(\beta_k=(1-s(t_k))\delta\) and \(\gamma_k=s(t_k)\delta\) giving the familiar QAOA angles [2511.12379].

"Universal Resources for QAOA and Quantum Annealing" refines this correspondence by introducing integrated coordinates
\[
\Theta(t)=\int_0^t -A(\tau)\,d\tau,\qquad
\Gamma(\Theta)=\int_0^\Theta B(\tau(\Theta_1)) \left(\frac{d\Theta}{d\tau}\right)^{-1}d\Theta_1,
\]
so that quantum annealing becomes a path in Hamiltonian space,
\[
i\partial_\Theta\ket{\psi(\Theta)}
=
\left[
\hat H_x+\frac{d\Gamma}{d\Theta}\hat H_{\rm QSNet}
\right]\ket{\psi(\Theta)}.
\]
In this language, QAOA is a first-order Trotter approximation to the same path, with cumulative angles \(\Theta_n\) and \(\Gamma_n\) identified with the continuous resources [2506.03241].

That work also reports that optimized QAOA angles collapse onto universal annealing-like trajectories in the \((\Theta,\Gamma)\) plane and interprets QAOA and QA outputs through a bimodal pseudo-Boltzmann distribution
\[
B(E)=\frac{1}{\mathcal Z_B}\left(c_{\rm high}e^{-\beta_{\rm high}E}+c_{\rm low}e^{-\beta_{\rm low}E}\right).
\]
Within that phenomenology, the coldest temperature scales as \(T\sim 1/p\), while \(\beta_{\rm high}\propto p\) and \(\Gamma_{\max}\propto p\) [2506.03241].

A related but distinct schedule-centered analysis appears in quantum chemistry. "Quantum Alternating Operator Ansatz (QAOA) Phase Diagrams and Applications for Quantum Chemistry" uses low-parameter linear ramps,
\[
\gamma(f)=\Delta f,\qquad \beta(f)=\Delta(1-f),\qquad f_j=\frac{j}{p+1},
\]
and studies performance as a function of \((\Delta,p)\) via squared overlap with the target ground state. The resulting phase diagrams show a regime in which non-adiabatic schedules perform better than the adiabatic limit while employing lower quantum circuit depth [2108.13056].

## 3. Hamiltonian design: custom phase operators, higher-order costs, and shortcut terms

A major branch of Hamiltonian-based QAOA research treats the Hamiltonians themselves as design variables rather than fixed problem encodings. In "Quantum approximate optimization algorithm with random and subgraph phase operators", the phase Hamiltonian \(C\) is allowed to differ from the original cost Hamiltonian \(C'\), with QAOA still optimizing the original objective. For MaxCut, the paper derives a closed-form \(p=1\) expression for each edge expectation \(\langle C'_{uv}\rangle\) in terms of: \(\chi_{uv}\), indicating whether \(Z_uZ_v\) is present in the custom phase operator; \(d\) and \(e\), counting incident phase terms; and \(f\), counting triangle-condition pairs [2402.18412].

The same study evaluates random, subgraph, triangle-removed, and maximal-degree-edge-removed phase operators on all non-isomorphic 8-vertex graphs. At \(p=1\), the percentage of tested graphs with at least one custom phase operator outperforming standard QAOA was \(0.036\%\) for random phase operators, \(75.9\%\) for subgraph phase operators, \(95.1\%\) for triangle-removed phase operators, and \(93.9\%\) for maximal-degree-edge-removed phase operators [2402.18412]. The paper’s interpretation is structural: triangle terms and high-degree vertices can worsen the one-layer expression, so Hamiltonian pruning can improve approximation ratio while also reducing circuit complexity.

Higher-order Hamiltonian design leads to a different conclusion: some objectives should not be quadratized at all. In the four-corners graph-coloring example of "QAOA of the Highest Order", the native higher-order binary encoding yields a quartic Ising Hamiltonian
\[
H_1=(a_{0}a_{1} + c_{0}c_{1})(b_{0}b_{1} + d_{0}d_{1}) + (a_{0} + c_{0})(b_{0} + d_{0}) + (a_{1} + c_{1})(b_{1} + d_{1}) + 4,
\]
implemented with \(40\) two-qubit gates, compared with \(96\) for the order-reduced binary formulation and \(80\) for the unary formulation. The reported QAOA simulations favor the native higher-order encoding in both gate cost and optimization behavior [2111.12754].

A third design axis adds explicit auxiliary Hamiltonians. "Shortcuts to Quantum Approximate Optimization Algorithm" introduces S-QAOA, where the standard alternating structure is enriched by an extra two-body interaction
\[
H_M=\sum_{(i,j)\in E}\frac{P_iQ_j + Q_iP_j}{2},\qquad PQ\in\{YZ,YY,XX,XZ,XY\}.
\]
For the MaxCut and SK instances studied there, the \(YY\) interaction performs best numerically. The paper attributes this to counterdiabatic structure visible in the BCH expansion and reports, for weighted 3-regular MaxCut, that at \(p=10\) S-QAOA still gives about a \(38\%\) fidelity improvement over QAOA, while for the SK model it reaches about \(80\%\) fidelity at \(p=1\) [2112.10943].

## 4. Constraint-preserving mixers and hard-feasibility Hamiltonians

Constrained optimization makes the mixer Hamiltonian a primary design object. The general design criteria in [1709.03489] are that the mixer should preserve the feasible subspace and connect all feasible states. This is the Hamiltonian analogue of replacing soft penalty enforcement by hard subspace-preserving dynamics.

Several concrete constructions instantiate that principle. One family uses mixers tailored to a feasible superposition \(\ket{F}=\frac{1}{\sqrt{|F|}}\sum_{x\in F}\ket{x}\). In the tutorial treatment of constrained QAOA, the Grover mixer satisfies
\[
e^{-i\beta H_M}=I-\left(1-e^{-i\beta}\right)|F\rangle\langle F|,
\]
and at \(\beta=\pi\) becomes the Grover-like reflection \(I-2|F\rangle\langle F|\) [2511.12379]. Another family uses excitation-preserving Hamiltonians. For fixed-cardinality portfolio optimization, "Constrained Counterdiabatic QAOA for Portfolio Optimization" employs the Hamming-weight-preserving XY mixer
\[
H_M^{XY}=\sum_{i<j}(X_iX_j+Y_iY_j),
\]
together with the Dicke initial state \(\ket{D^N_B}\), so that the baseline alternating dynamics remain inside the feasible \(|z|=B\) sector [2605.06858].

A different approach is to learn the mixer itself. "Quantum constraint learning for quantum approximate optimization algorithm" parameterizes a learned constrained unitary as
\[
V_{lk}(\boldsymbol{\alpha}, t)=W_l(\boldsymbol{\theta}) D_k(\boldsymbol{\gamma}, t) W_l(\boldsymbol{\theta})^{\dagger}=e^{-i H_{qcl}(\boldsymbol{\theta},\boldsymbol{\gamma}) t},
\]
with explicit conditions that \(e^{-iH_{qcl}t}\ket{\psi_i}\in\mathcal H_c\) and that feasible states remain connected under the evolution. The paper also introduces the projection operator \(\hat P_c=\sum_{x_i\in\mathcal H_c}\ket{x_i}\bra{x_i}\), the notion of quantum leakage when the learned mixer is imperfect, and the Wasserstein-based metric
\[
\zeta\left(\ket{\psi^*}\right)=1-W_1\left(\ket{\tilde{0}},\ket{\psi^*}\right)
\]
to assess constrained performance [2105.06770].

Hard-feasibility can also be encoded algebraically in the driver Hamiltonian. Choco-Q defines a commute Hamiltonian \(H^d\) satisfying \([\hat C,H^d]=0\), where \(\hat C\) is the operator form of the linear constraint. Its universal formulation is
\[
H^d=\sum_{\vec{u}\in \Delta} \left(\sigma_1^{u_1}\cdots \sigma_n^{u_n}+ \sigma_1^{-u_1}\cdots \sigma_n^{-u_n}\right),\qquad C\vec{u}=\vec{0},
\]
and the paper claims \(100\%\) in-constraints rate, more than \(235\times\) algorithmic improvement in successfully finding the optimal solution, and \(4.69\times\) end-to-end acceleration compared to prior QAOA designs [2503.23941].

Constraint preservation can be combined with counterdiabatic augmentation. In CCD-QAOA, approximate adiabatic gauge potentials derived from nested commutators of the Ising portfolio Hamiltonian and the XY mixer are inserted into each layer,
\[
\ket{\psi_{\mathrm{CD}}}
=
\prod_{k=1}^{p} e^{-i \eta_k H_{CD}} e^{-i \beta_k H_M^{XY}} e^{-i \gamma_k H_C} \ket{D^N_B},
\]
and the paper reports consistently higher approximation ratios than standard XY-mixer QAOA, Grover-mixer QAOA, and penalty-based QAOA at fixed depth, while also noting leakage and circuit-overhead tradeoffs introduced by the CD terms [2605.06858].

## 5. Extensions to chemistry, local Hamiltonians, continuous variables, and other nonclassical objectives

Hamiltonian-based QAOA is not limited to finding good classical bit strings. In quantum chemistry, the target is a many-body quantum ground state, not a computational-basis optimum. "Quantum Alternating Operator Ansatz (QAOA) Phase Diagrams and Applications for Quantum Chemistry" therefore chooses
\[
H_C=H_e,\qquad H_B=H_{\text{HF}},
\]
with the Hartree–Fock ground state as the initial state, so that the ansatz becomes
\[
\ket{\psi_p(\vec{\gamma},\vec{\beta})}
=
e^{-i\beta_pH_B}e^{-i\gamma_pH_C}\cdots e^{-i\beta_1H_B}e^{-i\gamma_1H_C}\ket{\psi_0}.
\]
The paper studies \(\mathrm{P}_2\), \(\mathrm{CO}_2\), \(\mathrm{Cl}_2\), and \(\mathrm{CH}_2\), with \(N\) values around \(114\)–\(128\) qubits, and uses ASCI to simulate the relevant subspace [2108.13056].

A more direct many-body generalization appears in "A Quantum Approximate Optimization Algorithm for Local Hamiltonian Problems", which introduces a specialized algorithm called HamQAOA for Local Hamiltonian Problems and Quantum MaxCut. For a general 2-local Hamiltonian, it uses four drivers,
\[
A=\sum_{u\sim v} Z_uZ_v,\qquad
B=\sum_v X_v,\qquad
C=\sum_v Z_v,\qquad
D=\sum_v n_v\cdot \vec\sigma_v,
\]
and the depth-\(p\) state
\[
\ket{\mathrm{HQS}_p}
=
e^{-i\delta_p D}e^{-i\gamma_p C}e^{-i\beta_p B}e^{-i\alpha_p A}\cdots e^{-i\delta_1 D}e^{-i\gamma_1 C}e^{-i\beta_1 B}e^{-i\alpha_1 A}\bigotimes_v \ket{m_v}.
\]
For Quantum MaxCut on Heisenberg systems, the paper reports rigorous high-girth guarantees, improved ground-energy-density bounds with depth, and numerical evidence that linear-depth HamQAOA can deterministically prepare exact ground states of 1-dimensional antiferromagnetic Heisenberg spin chains; it explicitly finds exact ground states for \(N=4\) at \(p=4\) and \(N=6\) at \(p=7\) [2412.09221].

Continuous-variable Hamiltonian-based QAOA goes further by replacing bit strings with positions \(\mathbf{x}\in\mathbb R^N\). In CV-QAOA, the cost Hamiltonian is
\[
\hat H_C=f(\hat{\mathbf{x}}),
\]
and the default mixer is kinetic,
\[
\hat H_M=\frac{1}{2}\hat{\mathbf{p}}^2.
\]
The key Heisenberg-picture update derived in [1902.00409] is
\[
\hat{\mathbf{x}}\rightarrow \hat{\mathbf{x}}+\gamma \hat{\mathbf{p}}-\eta\gamma\nabla f(\hat{\mathbf{x}}),
\]
so each layer acts like gradient descent with momentum, and for broad initial superpositions the algorithm becomes "gradient descent in superposition". The same framework incorporates equality and inequality constraints through penalty potentials and was numerically tested on the Styblinski–Tang function [1902.00409].

Problem-specific Hamiltonian synthesis has also been developed in domains such as digital communications. For ML detection with Gray-labelled constellations, the objective is transformed into a pseudo-Boolean polynomial and then into a diagonal Hamiltonian
\[
H_f=\sum_{\mathcal S\subseteq \mathcal N} g_{\mathcal S}\prod_{n\in\mathcal S}\sigma_z^{(n)},
\]
with the paper showing that for Gray-labelled rectangular MQAM the in-phase and quadrature qubits are independent in the Hamiltonian [2204.05126].

## 6. Universality and the limits of computational expressivity

Hamiltonian-based QAOA is not merely an optimization heuristic. "Quantum approximate optimization is computationally universal" proves that the alternating-Hamiltonian template can implement universal quantum computation when the layer times are treated as control parameters rather than variational angles [1812.11075].

The construction uses
\[
H_X=\sum_j X_j
\]
and a specially engineered one-dimensional \(ZZ\)-type Hamiltonian
\[
H_Z=\sum_j \Big( \omega_A Z_{2j} +\omega_B Z_{2j+1} +\gamma_{AB} Z_{2j}Z_{2j+1} +\gamma_{BA} Z_{2j+1}Z_{2j+2} \Big),
\]
with coefficients \(\omega_A,\omega_B,\gamma_{AB},\gamma_{BA}\) chosen to be not rationally related. By choosing a time \(t\) so that three unwanted phases nearly wrap around the circle while one desired term accumulates the target angle, the paper obtains effective evolutions such as
\[
\bigl\|e^{-itH_Z} - e^{-i\theta_A H_A}\bigr\| < \varepsilon,
\]
and similarly for the other components [1812.11075].

With these effective interactions, the alternating evolution synthesizes pairwise operations
\[
U_{AB}=e^{-i\theta_A H_A}e^{-i\theta_B H_B}e^{-i\theta_{AB}H_{AB}},\qquad
U_{BA}=e^{-i\theta_A H_A}e^{-i\theta_B H_B}e^{-i\theta_{BA}H_{BA}},
\]
acting in parallel on alternating nearest-neighbor pairs. The same paper further notes that appropriate global \(X\)-rotations give access to a \(Y\)-type Hamiltonian through
\[
e^{-i\frac{3\pi}{4}H_X}\, e^{-itH_Z}\, e^{-i\frac{5\pi}{4}H_X} = e^{-itH_Y}.
\]
This combination is then identified with a broadcast quantum cellular automaton architecture, implying that the QAOA-generated dynamics can simulate arbitrary quantum circuits to arbitrary accuracy, given enough layers and sufficiently precise timing [1812.11075].

The significance of this result is conceptual as much as constructive. The original variational interpretation remains intact, but the same alternating template
\[
e^{-i T_p H_X}e^{-i t_p H_Z}\cdots e^{-i T_1 H_X}e^{-i t_1 H_Z}
\]
can serve either as a heuristic optimizer or as a universal quantum processor, depending on how the times are chosen [1812.11075].

## 7. Parameter landscapes, optimizer behavior, and implementation strategies

The optimization problem induced by Hamiltonian-based QAOA can itself be studied analytically. For MaxCut at \(p=1\), "Quantum Approximate Optimization Algorithm for MaxCut: A Fermionic View" derives an explicit edgewise formula in terms of the degrees of the endpoints and the number of triangles containing the edge, showing that the one-layer contribution depends only on local graph structure [1706.02998]. For the 1D antiferromagnetic ring, the paper maps QAOA to control of an ensemble of independent pseudospins via a Jordan–Wigner transform, identifies symmetry-reduced critical manifolds such as
\[
\gamma_i+\beta_{p+1-i}=0,
\]
and reports that no local optima were observed numerically on the relevant reduced landscape for the ring case [1706.02998].

Practical optimization under hardware constraints has motivated a separate literature. "The QAOA with Few Measurements" studies the standard Hamiltonian formulation
\[
U(\vec{\gamma},\vec{\beta})=\prod_{j=1}^{p}e^{-i\beta_j B}e^{-i\gamma_j H}
\]
for MaxCut and shows that classical optimization is possible even with \(N=1\) shot per objective evaluation on a \(p=1\) problem with \(n=16\) active qubits. The reported optimizers are dual annealing, which used \(326\) parameter evaluations, and natural evolution strategies, which used \(M=10\), \(N=1\), \(30\) generations, and \(300\) total quantum evaluations [2205.06845].

Other implementation strategies modify the Hamiltonian between shallow runs. Loop-QAOA keeps the circuit at very small depth, typically \(p=1\), but updates the MaxCut Hamiltonian weights using the measured shallow-circuit output distribution. The paper reports that under bit-flip, phase-flip, and depolarizing noise, loop-QAOA continues improving as the number of loops increases, whereas conventional QAOA improves only up to a small depth and may then degrade [2109.11350].

Symmetry reduction offers another route. Automorphism-assisted QAOA identifies edge-equivalence classes under graph automorphisms and replaces many repeated Pauli terms by weighted representatives in a reduced Hamiltonian. On tree-structured graphs, the paper reports nearly unchanged approximation ratios together with large runtime and memory reductions; for example, on a \((31,30)\) binary tree at \(p=1\), optimization time drops from \(2067.36\) s for the full Hamiltonian to \(13.34\) s for the reduced Hamiltonian, with peak memory dropping from \(37.4\) GB to \(5.4\) GB [2410.22247]. The same work notes a caveat: if the reverse causal cone at the chosen depth already covers all relevant gates, Hamiltonian term reduction alone may not produce a comparable speedup [2410.22247].

Taken together, these results present Hamiltonian-based QAOA as a family of alternating Hamiltonian control schemes whose performance depends on far more than the nominal depth \(p\). Hamiltonian choice, locality, symmetry, commutator structure, schedule design, measurement budget, and feasible-subspace engineering all directly shape the attainable state manifold and the trainability of the resulting variational problem.

Source: https://www.emergentmind.com/topics/hamiltonian-based-quantum-approximate-optimization-algorithm-qaoa