Also known as landscape theory · traps · critical points · kinematic landscape · saddles
The control landscape is the figure of merit viewed as a function of the controls,
\(u(\cdot) \mapsto J[u]\) — an extremely high-dimensional surface over the space of all pulses. Its shape decides how
hard an optimization is. If every local optimum is also a global one, any local method such as GRAPE or
Krotov will succeed from almost any starting point. If there are traps — local optima that are not
global — the optimizer can get stuck. A central, and at first surprising, result of the field is that for many
quantum control problems the landscape is trap-free. This explains why optimal control works so well in practice.
Kinematic and dynamical landscapes
The map from controls to the figure of merit can be split into two steps:
The kinematic landscape \(J(U)\) is a function on the unitary group, for example the transition probability
\(|\langle f|U|i\rangle|^2\) or an expectation value \(\mathrm{Tr}\big(O\,U\rho U^\dagger\big)\). Its critical points can be
classified completely: they are global maxima, global minima and saddles, with no local optima [1]. The
dynamical landscape \(J[u]\) inherits this structure if the first map is well behaved
[2, 3].
When is the landscape trap-free?
Three assumptions guarantee that there are no traps [4]:
the map from controls to \(U(T)\) is locally surjective — its Jacobian has full rank at the critical points;
the controls are unconstrained.
Traps: where the assumptions fail
Singular controls. Where the Jacobian loses rank — at singular controls, the singular extremals of
Pontryagin's maximum principle — the kinematic argument does not apply. A three-level example shows that such
points can indeed be traps [5], and the assumptions of landscape analysis have been examined critically
[6]. A later geometric analysis proved that almost all closed, finite-dimensional systems
are nonetheless trap-free in a broad class of situations, including qubit evolutions in quantum computing. Being
singular is not by itself enough for a control to stop an optimization [7].
Constraints. The practically most important source of traps is the third assumption. Real optimizations restrict
the controls:
few parameters — a truncated basis, as in CRAB, or a coarse time grid;
short durations — close to the minimal time;
bounded amplitudes or bandwidth.
Simulations show that such constraints can create artificial traps, severe enough to prevent successful
optimization [4]. The minimum number of controls needed for a high yield depends on the number of
accessible states [8]. Even a two-level system develops a rich landscape structure when the total
time approaches the speed limit [9].
Near the minimal time. In many-body state preparation, the landscape changes character as the protocol duration
decreases. A glass-like transition appears: an exponential number of near-optimal solutions separated by barriers,
with the true optimum exponentially hard to find [10, 11]. See Quantum speed limits and
Many-body optimal control.
What it means for optimization
Why local methods usually work. With enough control freedom and enough time, a gradient method from a random
start reaches a global optimum. This is the common experience with GRAPE and Krotov.
When to expect trouble. Few parameters, short times, strong bounds: use restarts, continuation in time or
parameters, or methods that escape false traps. dCRAB changes its basis to reach the trap-free landscape
[12]; global methods find the best solution of the constrained problem directly
(QCPOP, Gradient-free & global optimizers).
Beyond quantum control. The same questions — how controllability shapes the landscape — decide the trainability of
variational quantum algorithms, where the dimension of the dynamical Lie algebra determines whether gradients vanish
(barren plateaus) [13]. See Variational algorithms & QAOA as control.
Key papers
Herschel A. Rabitz, Michael M. Hsieh and Carey M. RosenthalQuantum Optimally Controlled Transition LandscapesScience303, 1998-2001 (2004)doi
The landscape conjecture. Under mild assumptions, the optimal-control landscape for transition probabilities has no local traps — only global optima and saddles.
Raj Chakrabarti and Herschel RabitzQuantum control landscapesInternational Reviews in Physical Chemistry26, 671-735 (2007)doi
A comprehensive review of quantum control landscapes — kinematic and dynamical pictures, critical-point structure, consequences for optimization.
Alexander N. Pechen and David J. TannorAre there Traps in Quantum Control Landscapes?Phys. Rev. Lett.106, 120402 (2011)doi
A counterexample. Special (singular) critical points can create traps in the dynamical landscape, shown for a three-level system.
Benjamin Russell, Herschel Rabitz and Re-Bing WuControl landscapes are almost always trap free: a geometric assessmentJ. Phys. A: Math. Theor.50, 205302 (2017)doi
Proves that almost all closed, finite-dimensional systems have trap-free landscapes, and clarifies when singular controls can and cannot cause traps.
Gregory Riviello, Katharine Moore Tibbetts, Constantin Brif, Ruixing Long, Re-Bing Wu, et al.Searching for quantum optimal controls under severe constraintsPhys. Rev. A91, 043401 (2015)doi
How constraints on the number of controls, the duration and the field strength create artificial traps in practice.
References
Tak-San Ho and Herschel RabitzWhy do effective quantum controls appear easy to find?Journal of Photochemistry and Photobiology A: Chemistry180, 226-240 (2006)doi
Herschel A. Rabitz, Michael M. Hsieh and Carey M. RosenthalQuantum Optimally Controlled Transition LandscapesScience303, 1998-2001 (2004)doi
Raj Chakrabarti and Herschel RabitzQuantum control landscapesInternational Reviews in Physical Chemistry26, 671-735 (2007)doi
Gregory Riviello, Katharine Moore Tibbetts, Constantin Brif, Ruixing Long, Re-Bing Wu, et al.Searching for quantum optimal controls under severe constraintsPhys. Rev. A91, 043401 (2015)doi
Alexander N. Pechen and David J. TannorAre there Traps in Quantum Control Landscapes?Phys. Rev. Lett.106, 120402 (2011)doi
PIERRE DE FOUQUIERES and SOPHIE G. SCHIRMERA CLOSER LOOK AT QUANTUM CONTROL LANDSCAPES AND THEIR IMPLICATION FOR CONTROL OPTIMIZATIONInfin. Dimens. Anal. Quantum. Probab. Relat. Top.16, 1350021 (2013)doi
Benjamin Russell, Herschel Rabitz and Re-Bing WuControl landscapes are almost always trap free: a geometric assessmentJ. Phys. A: Math. Theor.50, 205302 (2017)doi
Katharine W. Moore and Herschel RabitzExploring constrained quantum control landscapesThe Journal of Chemical Physics137, 134113 (2012)doi
Martín Larocca, Pablo M Poggi and Diego A WisniackiQuantum control landscape for a two-level system near the quantum speed limitJ. Phys. A: Math. Theor.51, 385305 (2018)doi
Marin Bukov, Alexandre G. R. Day, Dries Sels, Phillip Weinberg, Anatoli Polkovnikov and Pankaj MehtaReinforcement Learning in Different Phases of Quantum ControlPhys. Rev. X8, 031086 (2018)doi
Alexandre G. R. Day, Marin Bukov, Phillip Weinberg, Pankaj Mehta and Dries SelsGlassy Phase of Optimal Quantum ControlPhys. Rev. Lett.122, 020601 (2019)doi
N. Rach, M. M. Müller, T. Calarco and S. MontangeroDressing the chopped-random-basis optimization: A bandwidth-limited access to the trap-free landscapePhys. Rev. A92, 062343 (2015)doi
Martin Larocca, Piotr Czarnik, Kunal Sharma, Gopikrishnan Muraleedharan, Patrick J. Coles and M. CerezoDiagnosing Barren Plateaus with Tools from Quantum Optimal ControlQuantum6, 824 (2022)doi
We revisit an established Quantum Control objective, characterize its exact local metric geometry, and reposition it as a rigorous benchmark for Search and Optimization. The trap-free topology of the landscape, established two decades ago, guarantees unhindered optimal pathways, but global topology alone does not govern convergence speed. Although the so-called Quantum Sphere is defined by the fourth power of the control field, we show that its yield gap admits an exact representation as a nonnegative linear combination of $\sin^2$ terms. From this single identity, we derive a two-sided quadratic sandwich with ratio $4/π^2$ on a certified, worst-case quarter-wave region, proving the optimal peak nondegenerate modulo its affine gauge and well-conditioned. For white-box optimization, on a halved (eighth-wave) region, the objective satisfies a Polyak-Łojasiewicz inequality with an explicit constant, granting Gradient Descent a linear rate of $1 - 16/(π^4κ_r)$, where $κ_r$ is the condition number of the relevant Hessian block. For black-box optimization, we prove that the elitist single-parent Evolution Strategy converges linearly, almost surely and in expected hitting time, from a certified sublevel set, despite the landscape's continuous gauge symmetry and its lattice of periodic optima; the certified level is explicit whenever the landscape is trap-free, and near the optimum the strategy provably sees only the quadratic tangent of the peak. Systematic numerical simulations corroborate the theory. Finally, we show that additive input (actuator) noise reweights each term of the representation by a Debye-Waller factor without disrupting the landscape's structure, and outline open challenges: output (signal) noise, far-field plateau escape, angular-synchronization certificates, and another physically scalable benchmark (SHG with Glass).
Engineering dissipation has emerged as a powerful strategy for controlling quantum many-body dynamics, complementing purely coherent control for state preparation. However, finite-time driving across a quantum critical point generates nonadiabatic excitations and coherences, and existing dissipation-engineering approaches have largely characterized how a fixed environment modifies critical dynamics rather than actively steering a system toward a prescribed thermal trajectory. Here we develop a Markovian dissipative framework for finite-time thermal-state tracking based on engineered population currents with non-negative transition rates, and apply it to the transverse-field Ising model driven across its quantum critical point. Compared with unitary dynamics and a natural thermal bath, the engineered dynamics strongly suppresses thermal excitation deviations and provides additional control over coherence damping through a free common current $Q_k$, while exhibiting a substantially steeper residual scaling, $δn_{\mathrm{ex}}^{\mathrm{th}}\sim(v/v_{\mathrm{ch}})^2$, than the Kibble--Zurek behavior recovered for unitary evolution. These results point toward engineered dissipation as a route for controlling thermal-state preparation in driven quantum critical systems.
Variational quantum algorithms based on qudits have attracted significant attention in recent years. However, as in their qubit-based counterparts, challenges such as barren plateaus and the design of efficient ansatz remain major obstacles. In this work, we propose to address these issues through a qudit implementation of the ADAPT-VQE algorithm, which constructs the ansatz iteratively. Specifically, we introduce an operator pool inspired by adiabatic evolution enhanced with counterdiabatic driving for ansatz construction and employ it to solve Max 3-Cut. We show that the warm-start strategy inherent to ADAPT-VQE, together with an ansatz construction based on counterdiabatic operators, achieves higher accuracy and lower native gates implementation than approaches on fixed ansatz. Furthermore, we show that, in qudit-based quantum computing, ADAPT-VQE with a counterdiabatic operator pool can navigate rough optimization landscapes with local traps through the burrowing mechanism, suggesting robustness against barren plateau effects and providing a scalable framework for variational quantum algorithms with qudits.
Optimal control problems arise in a wide range of scientific disciplines, but the corresponding optimization algorithms often display a strong dependence on hyperparameters that significantly influence performance and convergence. For the optimal implementation of quantum algorithms, these challenges are further amplified by high-dimensional control landscapes and the need for high-fidelity operations. Here, we propose an algorithm for optimal control problems in quantum computing to efficiently generate high-fidelity control protocols for multi-qubit systems in a hyperparameter and gradient free manner. The method, referred to as Pulse Engineering via Projection of response functions at infinite nonlinear order (PEPRino), leverages the framework of response theory to navigate the control landscape to find high-fidelity implementations. This is achieved by determining the control landscape via response functions to infinite order, efficiently evaluated by resummation in terms of the first and second order response function. To demonstrate the approach, we apply it to quantum systems consisting of two and three qubits for the optimal implementation of the Quantum Fourier Transform (QFT). We benchmark the proposed algorithm against the Chopped Random Basis (CRAB) algorithm utilizing the Nelder-Mead method, focusing on the 2-qubit scenario. The results indicate faster convergence regarding iteration steps and computational time, highlighting the advantages of our approach.
Quantum optimal control methods are widely used to design experimental control pulses such as laser amplitudes, phases, or detunings, that implement a target unitary evolution. In practice, what makes a pulse "good" depends not only on its fidelity, but also on the experimental setting and the relevant hardware constraints. Here, we introduce geometric quantum control with kernel optimisation (GECKO), a model-agnostic method for improving control pulses after a high-fidelity solution has been found. GECKO uses the Riemannian geometry of the special unitary group to identify directions in pulse space that leave the implemented unitary unchanged to first order, allowing one to traverse level sets of the control landscape while optimising a chosen differentiable pulse-quality function. We demonstrate GECKO on a transverse-field Ising Hamiltonian implementing CZ and CNOT gates, optimising pulse properties including spectral filtering, smoothness, robustness to parameter deviations, and pulse duration. In all cases, GECKO finds substantially improved pulse solutions.