Control of Quantum Systems


There's currently much interests among certain control engineers, physicists and physical chemists in trying to integrate ideas from control design with those of quantum dynamics. The applications range from laser chemistry (molecular vibrations control) and NMR (nuclear spin control) to quantum computing.
As a control problem, the dynamics are governed by the Schroedinger equation

where u(t) is an external (control) signal (e.g. an electromagnetic field strength in the laser chemistry case or a magnetic field in the NMR case), Ho is the nominal (internal) Hamiltonian, and V is the "control" Hamiltonian representing the coupling of the external control to the internal dynamics (e.g. dipole moment coupling). Notice that this is a so-called "bilinear system". These equations look deceptively simple. However, this hides the fact that infinite-dimensional bilinear systems are a universal class representing effectively all dynamical systems with an input (see "Myths, ...." in the talks page).
One of the simplest problems to pose is that of optimal population transfer: find the optimal control u(t) to transfer the system from one given state to another while minimizing control energy. Actually, this problem turns out to be very hard to solve for large quantum systems (read: molecules with two or more bonds). Except for very simple cases, it can only be solved numerically, and even then, one is only guaranteed to converge to local minima. None the less, numerical algorithms provide some insight. The picture below is that of a candidate control signal for a quantum Morse oscillator system with dipole moment coupling. Such signals look very complex, but if one looks at its spectrogram (i.e. a time-frequency decomposition), one sees that it is composed of weighted combinations of Bohr frequencies with time varying weights.
In fact one can be much more precise. We can show that in the limit of large transfer times T, the optimal control has the following form:

where the exponential term represents the Bohr frequency between the k and l states of the system. Thus, up to first order in 1/T, the optimal control is of the form of Bohr frequencies modulated by slowly varying (realtive to the time scale 1/T) envelopes. This makes precise the notion of time scale separation for such problems, and reduces the optimal control problem (for O(1/T) terms) to that of designing the envelopes. The slow envelopes obey optimal-control-type differential equations (which are derived via an averaging procedure) that are much less stiff than the original ones, and thus much easier to solve numerically.

The picture here shows a cartoon spectrogram of the control signal and the respective state populations for a 10 state representation of the quantum Morse oscillator. This particular control transitions from the ground state (1st state) to the 10th state. The spectrogram (bottom) is labeled by the Bohr frequencies (labeled by the corresponding transitions rather than the actual frequency). Notice that there's no component for the (1,10) transition Bohr frequency, but a collection of intermediate transitions. One can intuitively think of the optimal control algorithm as providing an optimally timed "schedule" of intermediate transitions.







Related Publications
  1. S. Grivopoulos and B. Bamieh. Optimal population transfers in a quantum system for large transfer time. IEEE Transactions on Automatic Control, 53(4):980-992, May 2008.




[BACK TO Pictures PAGE]