Energy-Selective Quantum Search

• Energy-selective quantum search is a framework that reformulates search as dynamic evolution toward energetically distinguished states via spectral bias, phase, or relaxation.
• It employs tunable energy ratios, non-Hermitian adiabatic evolution, and dissipative mechanisms to achieve faster search times with controlled fidelity trade-offs.
• Models such as Hamiltonian phase oracles and exciton trapping leverage engineered energy spectra and system–bath couplings for optimized search performance.

Searching arXiv for relevant papers on energy-selective quantum search and closely related formulations. In the literature surveyed here, energy-selective quantum search encompasses search procedures in which the search target is singled out through energetic bias, spectral phase, or relaxation dynamics rather than through a purely Boolean marked-set oracle. The relevant constructions include a generalized analog search Hamiltonian with a tunable energy ratio $r=E'/E$, non-Hermitian adiabatic search with a finite complex gap, dissipative relaxation of an $N$-level system into an unknown ground state, Hamiltonian phase oracles based on direct Ising evolution $\exp(-iTH)$, and exciton trapping at an impurity site whose energy defect marks the “winner” configuration (Cafaro et al., 2019, Nesterov et al., 2012, Allahverdyan et al., 2021, Plyashechnik et al., 2 Jun 2026, Thilagam, 2010).

1. Conceptual setting and relation to oracle search

A recurring feature of these models is that search is reformulated as dynamics toward energetically distinguished configurations. In the generalized analog Hamiltonian search, the target $\lvert w\rangle$ and source $\lvert s\rangle$ are assigned distinct energy weights through

$$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$

with $r=E'/E$ acting as an additional control parameter. In the non-Hermitian adiabatic construction, the target state is the ground state of $H_0=-|m\rangle\langle m|$, while an auxiliary term proportional to $(g+i\delta)(1-s)H_1$ controls the instantaneous gap. In the dissipative database search, the unique marked element is made the ground state of an $N$-level system and the task is recast as relaxation under a master equation. In the Ising Hamiltonian phase-oracle construction, the native evolution $N$0 marks configurations continuously by phases and selects a finite resonance band rather than a preassigned marked set. In the Frenkel-exciton model, a shallow isotopic impurity introduces a site-energy shift $N$1 that singles out $N$2 energetically (Cafaro et al., 2019, Nesterov et al., 2012, Allahverdyan et al., 2021, Plyashechnik et al., 2 Jun 2026, Thilagam, 2010).

These formulations differ materially from the standard Boolean-oracle model. Grover’s algorithm and its unitary Hamiltonian analogue accomplish unstructured search in $N$3 time steps, while the dissipative model is described as incoherent Markov dynamics in $N$4 levels with no “oracle calls,” and the Ising phase-oracle model uses a continuous spectral oracle rather than a Boolean marking rule (Allahverdyan et al., 2021, Plyashechnik et al., 2 Jun 2026). This suggests that complexity comparisons must be read together with the assumed oracle, spectral, and control resources.

2. Generalized analog search and the energy-ratio parameter

The generalized analog construction operates in an $N$5-dimensional Hilbert space with two special normalized states, the target $N$6 and the source $N$7, with overlap

$N$8

Because the Hamiltonian acts only in $N$9, the dynamics reduces to a two-dimensional problem after introducing

$\exp(-iTH)$0

so that $\exp(-iTH)$1. The exact propagator is $\exp(-iTH)$2, with eigenfrequency

$\exp(-iTH)$3

The success probability is

$\exp(-iTH)$4

and its peak value is

$\exp(-iTH)$5

while the first time at which the maximum is attained is

$\exp(-iTH)$6

For $\exp(-iTH)$7, one recovers the original Farhi–Gutmann result, with $\exp(-iTH)$8 and $\exp(-iTH)$9 (Cafaro et al., 2019).

The central feature of this model is a speed-versus-fidelity trade-off. When $\lvert w\rangle$0, the maximum fidelity drops below $\lvert w\rangle$1 but $\lvert w\rangle$2 can be shortened because the denominator $\lvert w\rangle$3 becomes larger than $\lvert w\rangle$4. The fidelity deficit is parametrized by

$\lvert w\rangle$5

with

$\lvert w\rangle$6

The analysis imposes the condition that the search infidelity not exceed the minimum-error probability for ambiguous discrimination of two pure states,

$\lvert w\rangle$7

For symmetric priors $\lvert w\rangle$8, this reproduces a Fuchs–van-de-Graaf-type bound. Numerically, the advantage region

$\lvert w\rangle$9

is non-empty, so the modified algorithm can attain a prescribed near-optimal fidelity in strictly less time than the $\lvert s\rangle$0 search (Cafaro et al., 2019).

3. Non-Hermitian adiabatic search and engineered finite gaps

The non-Hermitian adiabatic algorithm for Grover-type search begins with

$\lvert s\rangle$1

where

$\lvert s\rangle$2

Using the schedule parameter $\lvert s\rangle$3 and the complex coupling

$\lvert s\rangle$4

the full Hamiltonian is

$\lvert s\rangle$5

Because both $\lvert s\rangle$6 and $\lvert s\rangle$7 have degeneracy $\lvert s\rangle$8 at zero energy, the dynamics reduces exactly to the two-dimensional subspace spanned by $\lvert s\rangle$9. The instantaneous eigenvalues are

$$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$0

with

$$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$1

The complex gap is $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$2, and for large $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$3 its minimum magnitude is finite:

$$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$4

The minimum occurs near the avoided crossing $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$5 (Nesterov et al., 2012).

Because the minimum gap is controlled by $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$6 rather than by $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$7, the adiabatic run time required for adiabatic following is

$$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$8

In the Hermitian limit $$H = E\,|w\rangle\langle w| + E'\,|s\rangle\langle s| ,$$9, the minimum gap is $r=E'/E$0, giving $r=E'/E$1 for a global schedule. The success probability tends to $r=E'/E$2 already for $r=E'/E$3, whereas the overall norm decays because of the imaginary component:

$r=E'/E$4

The trade-off is therefore explicit: adding a small imaginary part $r=E'/E$5 opens a finite minimum gap and allows an adiabatic run time only logarithmic in $r=E'/E$6, but the system norm decays as $r=E'/E$7, so $r=E'/E$8 must be chosen small enough to keep $r=E'/E$9 over the evolution interval (Nesterov et al., 2012).

4. Dissipative relaxation search in logarithmic spectra

The dissipative search paradigm maps an unstructured database onto an $H_0=-|m\rangle\langle m|$0-level system weakly coupled to a thermal bath. The marked element corresponds to an unknown index $H_0=-|m\rangle\langle m|$1 and is made the ground state by

$H_0=-|m\rangle\langle m|$2

A commuting splitting Hamiltonian $H_0=-|m\rangle\langle m|$3 with nondegenerate spectrum $H_0=-|m\rangle\langle m|$4 removes the degeneracy of the remaining $H_0=-|m\rangle\langle m|$5 states, so that

$H_0=-|m\rangle\langle m|$6

A central construction chooses

$H_0=-|m\rangle\langle m|$7

and

$H_0=-|m\rangle\langle m|$8

which yields the gap

$H_0=-|m\rangle\langle m|$9

The weak-coupling, Born–Markov limit produces a classical master equation for the populations,

$(g+i\delta)(1-s)H_1$0

with detailed balance

$(g+i\delta)(1-s)H_1$1

An equivalent population-only Lindblad description uses jump operators

$(g+i\delta)(1-s)H_1$2

The generalized Glauber rates

$(g+i\delta)(1-s)H_1$3

satisfy detailed balance together with the bounded-rate condition

$(g+i\delta)(1-s)H_1$4

(Allahverdyan et al., 2021).

The relaxation time is governed by the slowest nonzero eigenvalue $(g+i\delta)(1-s)H_1$5 of the rate matrix $(g+i\delta)(1-s)H_1$6:

$(g+i\delta)(1-s)H_1$7

Numerically, for $(g+i\delta)(1-s)H_1$8 and $(g+i\delta)(1-s)H_1$9,

$N$0

so $N$1. The mechanism is attributed to three simultaneous conditions: at low temperature $N$2 the Gibbs weight overwhelmingly concentrates on level $N$3; the long-range rates allow direct jumps between any two levels; and the logarithmic spacing keeps upward and downward hops moderately biased while each site’s total exit rate remains bounded. In the trivial case $N$4, by contrast, one finds a purely classical $N$5 search. The model is therefore presented as an exponentially better nominal scaling in system size, but with explicit limitations: it requires fine engineering of the many-level energy spectrum $N$6 and of system–bath couplings, plus cooling to low $N$7 so that $N$8, and the resource cost of constructing a $N$9-level system with log-spaced levels and all-to-all couplings may offset the algorithmic gain (Allahverdyan et al., 2021).

5. Hamiltonian phase oracles and resonance-band amplification in Ising systems

The Ising Hamiltonian phase-oracle primitive uses the native evolution of an $N$00-spin Ising Hamiltonian directly:

$N$01

so that on a basis state $N$02,

$N$03

The oracle is

$N$04

and one iteration alternates it with the Grover diffusion operator:

$N$05

starting from the uniform superposition

$N$06

Configurations with $N$07 are only partially marked at first, but repeated reflections build a narrow resonance band whose amplitudes grow as in Grover-type amplification. The exact dynamics is organized through the spectral measure

$N$08

the characteristic-function samples

$N$09

and generating functions $N$10 and $N$11, with the algebraic identity

$N$12

The denominator

$N$13

acts as the spectral response; its minima in modulus determine resonance center, width, and height (Plyashechnik et al., 2 Jun 2026).

For an annealed Gaussian density of states with variance

$N$14

the characteristic samples are

$N$15

Targeting a prescribed energy $N$16 uses

$N$17

so that the principal resonance satisfies $N$18. In the high-density tail of the spectrum, the resonance contains $N$19 levels and the peak iteration count satisfies

$N$20

recovering the Grover-type square-root speedup; the abstract states the same scaling as $N$21 when the resonance contains $N$22 configurations. In a true random-Ising ensemble, overlap-induced correlations shift and distort the peak. The resonance condition is shifted by a small $N$23 with mean zero and variance

$N$24

so $N$25 in r.m.s. Although this is negligible on the $N$26 spectral scale, it is large compared to the exponentially narrow resonance width $N$27, so precise energy targeting requires correction. Two correction methods are given. Spectral symmetrization replaces $N$28 by $N$29, making the relevant characteristic function real and restoring an exact zero at $N$30. Iterative calibration instead updates $N$31 using measured output energies and is a contraction mapping; achieving final phase accuracy $N$32 requires only $N$33 calibration steps. Each iterate costs $N$34 gates for a dense two-body Ising Hamiltonian, and the overall gate count in the Gaussian tail regime is

$N$35

(Plyashechnik et al., 2 Jun 2026).

6. Physical realization by exciton trapping and broader interpretive issues

A concrete physical implementation of energy-selective search is the Frenkel-exciton trapping model in a one-dimensional molecular crystal with one impurity site $N$36. The exciton Hamiltonian is

$N$37

with long-range couplings

$N$38

The impurity is marked by

$N$39

and the oracle-like exciton–phonon interaction is

$N$40

with

$N$41

Because $N$42 contains the projector $N$43, it couples phonons only to the impurity and acts as the Grover oracle. The selectivity is explicitly energy-based: resonant exchange requires

$N$44

so only at site $N$45 can the exciton emit or absorb an acoustic phonon at low temperature (Thilagam, 2010).

The search time is extracted from the damping rate of the Green’s-function exponent and has the form

$N$46

with

$N$47

for $N$48. For $N$49, one obtains an $N$50-independent coherent search time; at the marginal $N$51, the $N$52-term recovers a Grover-like scaling; and for dipole-type coupling $N$53 one finds $N$54 as $N$55. Using naphthalene parameters $N$56, $N$57, $N$58, $N$59, with $N$60 and $N$61, the model gives $N$62 and $N$63–$N$64, while competing decoherence times are much longer at sufficiently low temperature (Thilagam, 2010).

Taken together, these proposals motivate several interpretive cautions. The dissipative $N$65 construction uses incoherent Markov dynamics in $N$66 levels and no “oracle calls,” the Ising primitive uses a Hamiltonian phase oracle that marks configurations continuously by their phases and selects a finite resonance band rather than a preassigned marked set, and the non-Hermitian adiabatic construction exchanges a finite gap for norm loss $N$67 (Allahverdyan et al., 2021, Plyashechnik et al., 2 Jun 2026, Nesterov et al., 2012). This suggests that the apparent improvements over $N$68 do not operate within a single resource model. The common pattern is not the replacement of one oracle by a faster oracle in the same formal setting, but the introduction of additional spectral, dissipative, or calibration structure: logarithmic level design and long-range bounded rates, fine control of the evolution time $N$69, ancilla-based symmetrization with three-body terms, or impurity-specific exciton–phonon coupling. Within those assumptions, energy selectivity becomes the operative mechanism that directs amplitude, probability, or excitation transport toward the sought configuration.