Jr-Shin Li

Ensemble Control of Bloch Equations

In this article, we study a class of control problems which involves controlling a continuum of dynamical systems with different values of parameters characterizing the system dynamics by using the same control signal. We call such problems control of ensembles. The motivation for looking into these problems comes from the manipulation of an ensemble of nuclear spins in nuclear magnetic resonance (NMR) spectroscopy and imaging (MRI) with dispersions in natural frequencies and the strengths of the applied radio frequency (rf) field. From the standpoint of mathematical control theory, the challenge is to simultaneously steer a continuum of systems between points of interest with the same control input. This raises some new and unexplored questions about controllability of such systems. We show that controllability of an ensemble can be understood by the study of the algebra of polynomials defined by the noncommuting vector fields that govern the system dynamics. A systematic study of these systems has immediate applications to broad areas of control of ensembles of quantum systems as arising in coherent spectroscopy and quantum information processing. The new mathematical structures appearing in such problems are excellent motivation for new developments in control theory.

Optimal pulse design in quantum control: A unified computational method

Many key aspects of control of quantum systems involve manipulating a large quantum ensemble exhibiting variation in the value of parameters characterizing the system dynamics. Developing electromagnetic pulses to produce a desired evolution in the presence of such variation is a fundamental and challenging problem in this research area. We present such robust pulse designs as an optimal control problem of a continuum of bilinear systems with a common control function. We map this control problem of infinite dimension to a problem of polynomial approximation employing tools from geometric control theory. We then adopt this new notion and develop a unified computational method for optimal pulse design using ideas from pseudospectral approximations, by which a continuous-time optimal control problem of pulse design can be discretized to a constrained optimization problem with spectral accuracy. Furthermore, this is a highly flexible and efficient numerical method that requires low order of discretization and yields inherently smooth solutions. We demonstrate this method by designing effective broadband π/2 and π pulses with reduced rf energy and pulse duration, which show significant sensitivity enhancement at the edge of the spectrum over conventional pulses in 1D and 2D NMR spectroscopy experiments.

Control and Synchronization of Neuron Ensembles

Synchronization of oscillations is a phenomenon prevalent in natural, social, and engineering systems. Controlling synchronization of oscillating systems is motivated by a wide range of applications from surgical treatment of neurological diseases to the design of neurocomputers. In this paper, we study the control of an ensemble of uncoupled neuron oscillators described by phase models. We examine controllability of such a neuron ensemble for various phase models and, furthermore, study the related optimal control problems. In particular, by employing Pontryagins maximum principle, we analytically derive optimal controls for spiking single- and two-neuron systems, and analyze the applicability of the latter to an ensemble system. Finally, we present a robust computational method for optimal control of spiking neurons based on pseudospectral approximations. The methodology developed here is universal to the control of general nonlinear phase oscillators.

Optimal Control of Inhomogeneous Ensembles

Inhomogeneity, in its many forms, appears frequently in practical physical systems. Readily apparent in quantum systems, inhomogeneity is caused by hardware imperfections, measurement inaccuracies, and environmental variations, and subsequently limits the performance and efficiency achievable in current experiments. In this paper, we provide a systematic methodology to mathematically characterize and optimally manipulate inhomogeneous ensembles with concepts taken from ensemble control. In particular, we develop a computational method to solve practical quantum pulse design problems cast as optimal ensemble control problems, based on multidimensional pseudospectral approximations. We motivate the utility of this method by designing pulses for both standard and novel applications. We also show the convergence of the pseudospectral method for optimal ensemble control. The concepts developed here are applicable beyond quantum control, such as to neuron systems, and furthermore to systems with by parameter uncertainty, which pervade all areas of science and engineering.

Optimal entrainment of neural oscillator ensembles

In this paper, we derive the minimum-energy periodic control that entrains an ensemble of structurally similar neural oscillators to a desired frequency. The state-space representation of a nominal oscillator is reduced to a phase model by computing its limit cycle and phase response curve, from which the optimal control is derived by using formal averaging and the calculus of variations. We focus on the case of a 1:1 entrainment ratio and suggest a simple numerical method for approximating the optimal controls. The method is applied to asymptotically control the spiking frequency of neural oscillators modeled using the Hodgkin-Huxley equations. Simulations are used to illustrate the optimality of entrainment controls derived using phase models when applied to the original state-space system, which is crucial for using phase models in control synthesis for practical applications. This work addresses a fundamental problem in the field of neural dynamics and provides a theoretical contribution to the optimal frequency control of uncertain oscillating systems.

A pseudospectral method for optimal control of open quantum systems

In this paper, we present a unified computational method based on pseudospectral approximations for the design of optimal pulse sequences in open quantum systems. The proposed method transforms the problem of optimal pulse design, which is formulated as a continuous-time optimal control problem, to a finite-dimensional constrained nonlinear programming problem. This resulting optimization problem can then be solved using existing numerical optimization suites. We apply the Legendre pseudospectral method to a series of optimal control problems on open quantum systems that arise in nuclear magnetic resonance spectroscopy in liquids. These problems have been well studied in previous literature and analytical optimal controls have been found. We find an excellent agreement between the maximum transfer efficiency produced by our computational method and the analytical expressions. Moreover, our method permits us to extend the analysis and address practical concerns, including smoothing discontinuous controls as well as deriving minimum-energy and time-optimal controls. The method is not restricted to the systems studied in this article and is applicable to optimal manipulation of both closed and open quantum systems.

A multidimensional pseudospectral method for optimal control of quantum ensembles

In our previous work, we have shown that the pseudospectral method is an effective and flexible computation scheme for deriving pulses for optimal control of quantum systems. In practice, however, quantum systems often exhibit variation in the parameters that characterize the system dynamics. This leads us to consider the control of an ensemble (or continuum) of quantum systems indexed by the system parameters that show variation. We cast the design of pulses as an optimal ensemble control problem and demonstrate a multidimensional pseudospectral method with several challenging examples of both closed and open quantum systems from nuclear magnetic resonance spectroscopy in liquid. We give particular attention to the ability to derive experimentally viable pulses of minimum energy or duration.

Minimum-Time Frictionless Atom Cooling in Harmonic Traps

Frictionless atom cooling in harmonic traps is formulated as a time-optimal control problem and a synthesis of optimal controlled trajectories is obtained.

Ensemble Controllability of the Bloch Equations

Finding control fields (pulse sequences) that can compensate for the dispersion in the parameters governing the evolution of a quantum system is an important problem in coherent spectroscopy and quantum information processing. The use of composite pulses for compensating dispersions in system dynamics is widely known and applied. In this paper, we make explicit the key aspects of the dynamics that make such compensation possible. We highlight the role of Lie algebras and non-commutativity in the design of a compensating pulse sequence. Finally, we investigate three common dispersions in NMR spectroscopy, the Larmor dispersion, rf inhomogeneity and the strength of couplings between the spins

PHASE-SELECTIVE ENTRAINMENT OF NONLINEAR OSCILLATOR ENSEMBLES

The ability to organize and finely manipulate the hierarchy and timing of dynamic processes is important for understanding and influencing brain functions, sleep and metabolic cycles, and many other natural phenomena. However, establishing spatiotemporal structures in biological oscillator ensembles is a challenging task that requires controlling large collections of complex nonlinear dynamical units. In this report, we present a method to design entrainment signals that create stable phase patterns in ensembles of heterogeneous nonlinear oscillators without using state feedback information. We demonstrate the approach using experiments with electrochemical reactions on multielectrode arrays, in which we selectively assign ensemble subgroups into spatiotemporal patterns with multiple phase clusters. The experimentally confirmed mechanism elucidates the connection between the phases and natural frequencies of a collection of dynamical elements, the spatial and temporal information that is encoded within this ensemble, and how external signals can be used to retrieve this information.