Kick and fix: the roots of quantum control

TL;DR

This paper surveys the mathematical techniques for approximating the exponential of sums of non-commuting operators, crucial for quantum control and related fields, highlighting their historical development and interconnections.

Contribution

It provides a comprehensive historical overview and unifies pulsed and continuous approaches in quantum operator exponentiation.

Findings

Highlights the link between operator exponentiation and the strong coupling regime.

Shows that pulsed and continuous formulations are fundamentally connected.

Discusses the importance of these techniques in quantum control and technologies.

Abstract

When two operators $A$ and $B$ do not commute, the calculation of the exponential operator $e^{A+B}$ is a difficult and crucial problem. The applications are vast and diversified: to name but a few examples, quantum evolutions, product formulas, quantum control, Zeno effect. The latter are of great interest in quantum applications and quantum technologies. We present here a historical survey of results and techniques, and discuss differences and similarities. We also highlight the link with the strong coupling regime, via the adiabatic theorem, and contend that the "pulsed" and "continuous" formulations differ only in the order by which two limits are taken, and are but two faces of the same coin.

Equations

$$U_N(t)$$

$$V_N(t)$$

$$U_N(t)$$

$${\mathrm{e}}^{-\gamma Q}=P+{\mathrm{e}}^{-\gamma }Q,$$

$$U_N(t)$$

$$H_{\mathrm{Z}}=\sum _n{P}_{n}H{P}_n$$

$$U_{\mathrm{kick}}=\sum _n{\mathrm{e}}^{-\mathrm{i}{\lambda }_n}{P}_n.({\mathrm{e}}^{-\mathrm{i}{\lambda }_n}\ne {\mathrm{e}}^{-\mathrm{i}{\lambda }_l},\text{mbox}forn\ne l.)$$

$$\underset{N\mathrm{times}}{\underbrace{\left({\mathrm{e}}^{-\mathrm{i}\tau T}{\mathrm{e}}^{-\mathrm{i}{\tau }_{0}V}\right)\left({\mathrm{e}}^{-\mathrm{i}\tau T}{\mathrm{e}}^{-\mathrm{i}{\tau }_{0}V}\right)\cdots \left({\mathrm{e}}^{-\mathrm{i}\tau T}{\mathrm{e}}^{-\mathrm{i}{\tau }_{0}V}\right)}}={\left({\mathrm{e}}^{-\mathrm{i}\tau T}{\mathrm{e}}^{-\mathrm{i}{\tau }_{0}V}\right)}^N,\text{as }N\to \infty ,$$

$$H_K=H+K{H}_{\mathrm{c}},$$

$$U_K(t)={\mathrm{e}}^{-\mathrm{i}t{H}_K}\sim {\mathrm{e}}^{-\mathrm{i}t{H}_{\mathrm{Z}}}{\mathrm{e}}^{-\mathrm{i}tK{H}_{\mathrm{c}}},\text{as }K\to \infty ,$$

$$H_{\mathrm{c}}=\sum _n{\eta }_n{P}_n,({\eta }_n\ne {\eta }_m,\text{mbox}forn\ne m).$$

$$H_{\mathrm{kick}}=H+{\tau }_0\sum _n\delta (t-n\tau )H_{\mathrm{c}},$$

$$H(\tau,K)=H+K\sum_{n}g\big{(}t-n(\tau+\tau_{0}/K)\big{)}H_{\mathrm{c}},$$

$$\underset{K\to \infty }{\lim }\underset{\tau \to 0}{\lim }{U}_{\tau ,K}(t)=\underset{\tau \to 0}{\lim }\underset{K\to \infty }{\lim }{U}_{\tau ,K}(t),$$

Full text

\pubvolume

xx \issuenum1 \articlenumber1

\history \Title Kick and fix: the roots of quantum control

\AuthorPaolo Facchi 1,2*\orcidA and Saverio Pascazio 1,2,3\orcidB \AuthorNamesPaolo Facchi and Saverio Pascazio

\corresCorrespondence: paolo.facchi@ba.infn.it

1 Introduction

Quantum control is a relatively novel discipline, in which one seeks to control the evolution or the final state of a quantum system or to determine the output of a quantum process. The main objective is to counter the detrimental effects of (uncontrollable) noise and environment. The inspiration goes back to classical control, where one would make some measurement of the output and then use the results of such measurement to obtain the desired outcome (feedback mechanism). Similar ideas (but different techniques) motivated the development, in the 1970s and 1980s, of nuclear magnetic resonance imaging, towards very successful applications in medical diagnosis. The objective in that case was to control (spin) coherence.

In quantum mechanics and quantum applications things become more involved, because a quantum measurement collapses the wave function, thereby strongly affecting the quantum evolution. The question becomes: what can one actually do to a quantum system (from carefully tailoring measurements to accurately designing Hamiltonians or even Lindbladians) that will enable some control of the output? In an even more recent twist, some researchers have also considered the option of engineering the environment itself, towards specific quantum goals.

The quantum control problem translates to rotating a vector state from a specified initial direction to a desired final direction in a certain time. In a more general setting, this means evolving a generic state (density matrix) from a given initial state to a desired final state in a certain time. This is a difficult task, which essentially consists in the construction of a suitable transformation, in agreement with the guiding principles of quantum mechanics. We highlight in this article two main techniques that are studied nowadays.

2 Divide et impera (Divide and conquer)

2.1 History

In 1948 Feynman, working on his path-integral formulation of quantum mechanics ref:Feynman(1948) ; ref:FeynmanHibbs(1965) ; Schulman1981 noticed that one can obtain the quantum evolution by applying the following formula

[TABLE]

where $H=T+V$ is the Hamiltonian of the system, $T$ and $V$ being respectively the kinetic and potential energy. Feynman strived to bypass a crucial difficulty that plagues (exponentiated) quantum operators. Formula (1) is valid although $T$ and $V$ do not commute: the key observation is that the error one makes in approximating the evolution with the $N$-product is of order $1/N$, and disappears in taking the limit.

Feynman built on a crucial observation that Dirac had made 15 years earlier ref:Dirac1933 , concerning the infinitesimal expression of the quantum propagator. Feynman made no attempt at rigor (beyond the reasonable level of theoretical physics), unlike (presumably) Dirac, who did not dare to go beyond infinitesimal evolutions.

In fact, mathematicians independently worked on what they called product formulae. Trotter in 1959 ref:Trotter and Kato in 1978 ref:KatoTrotter investigated the limits of products of the following type

[TABLE]

and proved convergence for certain unbounded (linear) operators $A$ and $B$. Both Trotter and Kato were extending results by Lie ref:Lie , who had proved the validity of Eq. (2) for $n\times n$ matrices.

2.2 Zeno

Misra and Sudarshan’s formulation of the quantum Zeno effect ref:QZEMisraSudarshan ; ref:QZEreview-JPA makes use of the following formula

[TABLE]

where $P$ is a projection operator and $H_{\textrm{Z}}=PHP$ the "Zeno" Hamiltonian ref:QZS . Notice that in general $U_{N}(t)$ is not unitary on the range of $P$, while the limit is, under suitable conditions and in particular cases, such as bounded $H$ or finite-dimensional $P$ ref:artzeno . A general proof of the above formula is still missing, together with a rigorous definition of the Zeno Hamiltonian for infinite-dimensional projections and unbounded Hamiltonian.

Interestingly, the above formula can be re-written in terms of an absorbing "optical potential". Let $P+Q=1$, with both $P$ and $Q$ projections. Then, by making use of the identity (valid in general, but we take $\gamma>0$)

[TABLE]

one gets

[TABLE]

which is essentially the same as (3). In a few words, the above formula works because the imaginary optical potential $\mathrm{i}\gamma Q$ quickly (on a timescale $1/N\gamma$) absorbs away the unwanted component of the wave function, acting as a projection $P$.

Observe also that formula (3) is at the basis of the Faddeev-Popov method ref:FaddeevPopov to quantize gauge quantum field theories within the framework of the path integral formulation. In this case, the projection imposes the gauge condition.

The afore-mentioned formulas are at the basis of the "pulsed" formulation of quantum control. Physically, the validity of the constraint is preserved during the evolution, but the quantum system is free to move within the (Zeno) subspaces of the (multidimensional constraint) ref:QZS . In applications, $P$ must be designed in such a way as to reduce (or hinder) decoherence.

2.3 Kicks

In a similar spirit, physically equivalent dynamics can be obtained for a quantum system undergoing repeated unitary "kicks", the physical duration of the kick being the shortest timescale of the problem (notice the analogy with a projection, that is also supposed to take place instantaneously). The evolution reads

[TABLE]

where

[TABLE]

is the Zeno Hamiltonian and $P_{n}$ the spectral projections of the kick (taken with a discrete spectrum)

[TABLE]

This is again a Zeno dynamics. The Zeno subspaces are now a consequence of rapidly oscillating phases between different eigenspaces of the kick, yielding a superselection rule. This phenomenon was discovered in chemical physics, where it is known as NMR anderson ; ernst ; freeman ; levitt , and (re)baptized "bang-bang" control in the quantum-information literature viola98 . A very nice review is lidarrev .

Finally, we mention that the same formulas are of crucial importance in the study of quantum chaos ref:CasatiChaos ; ref:BerryChaos ; ref:Gutzwiller . Set $t=N\tau$, keep $\tau$ fixed, and let $U_{\mathrm{kick}}=\mathrm{e}^{-\mathrm{i}\tau_{0}{V}}$, where $\tau_{0}$ is needed for dimensional purposes. The evolution is

[TABLE]

where $T$ is a kinetic energy operator and $V$ a potential. In the above formula $V$ is not divided by $N$, so that the evolution is governed by a singular time-dependent Hamiltonian $H=T+\tau_{0}\sum_{n}\delta(t-n\tau)V$. Notice that unlike in formula (6), since $\tau$ is kept fixed, this is a large-time limit $t=N\tau\to\infty$. Alternatively, the evolution (9) can be thought as a sequence of purely kinetic and purely potential motions in which alternately the force is switched off (for a time $\tau$) and the mass made infinite (for a time $\tau_{0}$) ref:BerryChaos .

3 Persuade et rege (Persuade and rule)

Short timescales are physically associated with strong couplings. It must therefore be possible to rephrase the above evolutions in terms of a strong (continuous) coupling. Let

[TABLE]

where $H$ is the Hamiltonian of the system, $H_{\mathrm{c}}$ an interaction Hamiltonian performing a "continuous measurement" and $K$ a coupling constant.

We are interested in the $K\to\infty$ limit of the evolution operator

[TABLE]

where $H_{\mathrm{Z}}$ is given in Eq. (7), $P_{n}$ being the spectral projection of $H_{\mathrm{c}}$ (taken with a discrete spectrum)

[TABLE]

These results are consequences of the adiabatic theorem ref:KatoAdiabatic . The Zeno subspaces are again a consequence of the wildly oscillating phases between different eigenspaces. One can interpret the above results by saying that the external field/potential takes a steady, "persuading gaze" at the system.

4 Equivalence between continuous and pulsed formulations

The similarity between the "kicked" and "continuous" formulations, outlined in Secs. 2.3 and 3, is in fact even more profound. The two procedures are physically equivalent, and only differ in the order in which two limits are computed ref:BBZeno .

The continuous case considers Eq. (10) in the strong coupling limit $K\rightarrow\infty$, while the kicked dynamics is generated by the time-dependent Hamiltonian

[TABLE]

where $\tau$ is the period between two kicks and $U_{\mathrm{kick}}=\exp(-\mathrm{i}\tau_{0}H_{\mathrm{c}})$. The $N\rightarrow\infty$ limit in (6) corresponds to $\tau\rightarrow 0$. The Hamiltonians (10) and (13) are both limiting cases of the following one

[TABLE]

where $g(t)=\chi_{[-\tau_{0}/2K,\tau_{0}/2K]}(t)$, with a fixed $\tau_{0}>0$. In Eq. (14) the period between two kicks is $\tau_{0}/K+\tau$, while the kick lasts for a time $\tau_{0}/K$. By taking the limit $\tau\rightarrow 0$ in Eq. (14) (a sequence of pulses of finite duration $\tau_{0}/K$ without any time interval among them), one recovers the continuous case (10): $H(\tau,K)\to H_{K}$, as $\tau\to 0$. One takes afterwards the strong coupling limit $K\rightarrow\infty$ and gets the Zeno dynamics and subspaces. Let us now invert the order of the limits: first take the $K\rightarrow\infty$ limit (short pulses, but with the same global—integral—effect), to obtain the kicked case (13): $H(\tau,K)\to H_{\mathrm{kick}}$, as $K\to\infty$. Then, take the vanishing time interval limit $\tau\rightarrow 0$ to get again the Zeno dynamics and subspaces.

In short, by denoting with $U_{\tau,K}(t)$ the unitary evolution generated by $H(\tau,K)$, the formal equivalence between the two protocols is expressed by the relation

[TABLE]

with the left (right) side expressing the continuous (pulsed) case. This equivalence is only valid in the limit. It is physically legitimate if the inverse Zeno regime ref:InverseZeno is avoided. In practice, for finite $N$ and $K$, there can be important differences ref:ControlDecoZeno .

5 Comments and conclusions

The above discussion hinged upon Schrödinger equations, unitary dynamics, and projections à la von Neumann. This is only a part of the whole picture. Most results can be generalized to master (GKLS ref:GKLS-DariuszSaverio ) equations, quantum semigroups, and sequences of generic quantum operations ref:ZanardiDFS-PRL2014 ; ref:ZanardiDFS-PRA2017 ; ref:unity1 ; ref:unity2 . Moreover, the analysis can be extended to multidimensional spaces, highlighting interesting relations with geometry ref:berry ; ref:geo and complexity ref:plato .

American composer Harold Budd once said: "But that’s fine, because I like to have control of the ambience." Lucky him. Physicists and mathematicians can only endeavor to control the system, out of the (often detrimental) effects of the environment.

Acknowledgements.

We thank many colleagues with whom we discussed the quantum Zeno effect and quantum control over the last 20 years. A special thank to our students G. Gramegna, D. Lonigro, D. Pomarico and G. Scala for their constant interest and for motivating our research. This work is partially supported by Istituto Nazionale di Fisica Nucleare (INFN) through the project "QUANTUM". PF is partially supported by the Italian National Group of Mathematical Physics (GNFM-INdAM). \reftitleReferences