Dionisis Stefanatos, Heinz Schaettler, Jr-Shin Li
Frictionless atom cooling in harmonic traps is formulated as a time-optimal control problem and a synthesis of optimal controlled trajectories is obtained.
Justin Ruths, Jr-Shin Li
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.
Anatoly Zlotnik, Jr-Shin Li
An emerging and challenging area in mathematical control theory called Ensemble Control encompasses a class of problems that involves the guidance of an uncountably infinite collection of structurally identical dynamical systems, which are indexed by a parameter set, by applying the same open-loop control. The subject originates from the study of complex spin dynamics in Nuclear Magnetic Resonance (NMR) spectroscopy and imaging (MRI). A fundamental question concerns ensemble controllability, which determines the existence of controls that transfer the system between desired initial and target states. For ensembles of finite-dimensional time-varying linear systems, the necessary and sufficient controllability conditions and analytical optimal control laws have been shown to depend on the singular system of the operator characterizing the system dynamics. Because analytical solutions are available only in the simplest cases, there is a need to develop numerical methods for synthesizing these controls. We introduce a direct, accurate, and computationally efficient algorithm based on the singular value decomposition (SVD) that approximates ensemble controls of minimum norm for such systems. This method enables the application of ensemble control to engineering problems involving complex, time-varying, and high-dimensional linear dynamic systems.
Dionisis Stefanatos, Jr-Shin Li
We formulate the problem of efficient transport of a quantum particle trapped in a harmonic potential which can move with a bounded velocity, as a minimum-time problem on a linear system with bounded input. We completely solve the corresponding optimal control problem and obtain an interesting bang-bang solution. These results are expected to find applications in quantum information processing, where quantum transport between the storage and processing units of a quantum computer is an essential step. They can also be extended to the efficient transport of Bose-Einstein condensates, where the ability to control them is crucial for their potential use as interferometric sensors.
Isuru Dasanayake, Jr-Shin Li
In this paper, we study the optimal control of phase models for spiking neuron oscillators. We focus on the design of minimum-power current stimuli that elicit spikes in neurons at desired times. We furthermore take the charge-balanced constraint into account because in practice undesirable side effects may occur due to the accumulation of electric charge resulting from external stimuli. Charge-balanced minimum-power controls are derived for a general phase model using the maximum principle, where the cases with unbounded and bounded control amplitude are examined. The latter is of practical importance since phase models are more accurate for weak forcing. The developed optimal control strategies are then applied to both mathematically ideal and experimentally observed phase models to demonstrate their applicability, including the phase model for the widely studied Hodgkin-Huxley equations.
Anatoly Zlotnik, Jr-Shin Li
We derive optimal periodic controls for entrainment of a self-driven oscillator to a desired frequency. The alternative objectives of minimizing power and maximizing frequency range of entrainment are considered. A state space representation of the oscillator is reduced to a linearized phase model, and the optimal periodic control is computed from the phase response curve using formal averaging and the calculus of variations. Computational methods are used to calculate the periodic orbit and the phase response curve, and a numerical method for approximating the optimal controls is introduced. Our method is applied to asymptotically control the period of spiking neural oscillators modeled using the Hodgkin-Huxley equations. This example illustrates the optimality of entrainment controls derived using phase models when applied to the original state space system.
Dionisis Stefanatos, Jr-Shin Li
In this article we study the frictionless cooling of atoms trapped in a harmonic potential, while minimizing the transient energy of the system. We show that in the case of unbounded control, this goal is achieved by a singular control, which is also the time-minimal solution for a "dual" problem, where the energy is held fixed. In addition, we examine briefly how the solution is modified when there are bounds on the control. The results presented here have a broad range of applications, from the cooling of a Bose-Einstein condensate confined in a harmonic trap to adiabatic quantum computing and finite time thermodynamic processes.
Wei Miao, Gong Cheng, Jr-Shin Li
In this paper, we study the problem of learning dynamical properties of ensemble systems from their collective behaviors using statistical approaches in reproducing kernel Hilbert space (RKHS). Specifically, we provide a framework to identify and cluster multiple ensemble systems through computing the maximum mean discrepancy (MMD) between their aggregated measurements in an RKHS, without any prior knowledge of the system dynamics of ensembles. Then, leveraging the gradient flow of the newly proposed notion of aggregated Markov parameters, we present a systematic framework to recognize and identify an ensemble systems using their linear approximations. Finally, we demonstrate that the proposed approaches can be extended to cluster multiple unknown ensembles in RKHS using their aggregated measurements. Numerical experiments show that our approach is reliable and robust to ensembles with different types of system dynamics.
Wei Miao, Vignesh Narayanan, Jr-Shin Li
The reservoir computing networks (RCNs) have been successfully employed as a tool in learning and complex decision-making tasks. Despite their efficiency and low training cost, practical applications of RCNs rely heavily on empirical design. In this paper, we develop an algorithm to design RCNs using the realization theory of linear dynamical systems. In particular, we introduce the notion of $α$-stable realization, and provide an efficient approach to prune the size of a linear RCN without deteriorating the training accuracy. Furthermore, we derive a necessary and sufficient condition on the irreducibility of number of hidden nodes in linear RCNs based on the concepts of controllability and observability matrices. Leveraging the linear RCN design, we provide a tractable procedure to realize RCNs with nonlinear activation functions. Finally, we present numerical experiments on forecasting time-delay systems and chaotic systems to validate the proposed RCN design methods and demonstrate their efficacy.