Exploring the top and bottom of the quantum control landscape

A controlled quantum system possesses a search landscape defined by the target physical objective as a function of the controls. This paper focuses on the landscape for the transition probability P(i → f) between the states of a finite level quantum system. Traditionally, the controls are applied fields; here, we extend the notion of control to also include the Hamiltonian structure, in the form of time independent matrix elements. Level sets of controls that produce the same transition probability value are shown to exist at the bottom P(i → f)=0.0 and top P(i → f)=1.0 of the landscape with the field and/or Hamiltonian structure as controls. We present an algorithm to continuously explore these level sets starting from an initial point residing at either extreme value of P(i → f). The technique can also identify control solutions that exhibit the desirable properties of (a) robustness at the top and (b) the ability to rapidly rise towards an optimal control from the bottom. Numerical simulations are presented to illustrate the varied control behavior at the top and bottom of the landscape for several simple model systems.     

Local topology at limited resource induced suboptimal traps on the quantum control landscape

In a quantum optimal control experiment a system is driven towards a target observable value with a tailored external field. The underlying quantum control landscape, defined by the observable as a function of the control variables, lacks suboptimal extrema upon satisfaction of certain physical assumptions. This favorable topology implies that upon climbing the landscape to seek an optimal control field, a steepest ascent algorithm should not halt prematurely at suboptimal critical points, or traps. One of the important aforementioned assumptions is that no limitations are imposed on the control resources. Constraints on the control restricts access to certain regions of the landscape, potentially preventing optimal performance through convergence to limited resource induced suboptimal traps. This work develops mathematical tools to explore the local landscape structure around suboptimal critical points. The landscape structure may be favorably altered by systematically relaxing the control resources. In this fashion, isolated suboptimal critical points may be transformed into extensive level sets and then to saddle points permitting further landscape ascent. Time-independent kinematic controls are employed as stand-ins for traditional dynamic controls to allow for performing a simpler constrained resource landscape analysis. The kinematic controls can be directly transferred to their dynamic counterparts at any juncture of the kinematic analysis. The numerical simulations employ a family of landscape exploration algorithms while imposing constraints on the kinematic controls. Particular algorithms are introduced to meet the goals of either climbing the landscape or seeking specific changes in the topology of the landscape by relaxing the control resources.     

Optimal control landscapes for quantum observables

The optimal control of quantum systems provides the means to achieve the best outcome from redirecting dynamical behavior. Quantum systems for optimal control are characterized by an evolving density matrix and a Hermitian operator associated with the observable of interest. The optimal control landscape is the observable as a functional of the control field. The features of interest over this control landscape consist of the extremum values and their topological character. For controllable finite dimensional quantum systems with no constraints placed on the controls, it is shown that there is only a finite number of distinct values for the extrema, dependent on the spectral degeneracy of the initial and target density matrices. The consequences of these findings for the practical discovery of effective quantum controls in the laboratory is discussed.     

Exploiting time-independent Hamiltonian structure as controls for manipulating quantum dynamics

The opportunities offered by utilizing time-independent Hamiltonian structure as controls are explored for manipulating quantum dynamics. Two scenarios are investigated using different manifestations of Hamiltonian structure to illustrate the generality of the concept. In scenario I, optimally shaped electrostatic potentials are generated to flexibly control electron scattering in a two-dimensional subsurface plane of a semiconductor. A simulation is performed showing the utility of optimally setting the individual voltages applied to a multi-pixel surface gate array in order to produce a spatially inhomogeneous potential within the subsurface scattering plane. The coherent constructive and destructive electron wave interferences are manipulated by optimally adjusting the potential shapes to alter the scattering patterns. In scenario II, molecular vibrational wave packets are controlled by means of optimally selecting the Hamiltonian structure in cooperation with an applied field. As an illustration of the concept, a collection (i.e., a level set) of dipole functions is identified where each member serves with the same applied electric field to produce the desired final transition probability. The level set algorithm additionally found Hamiltonian structure controls exhibiting desirable physical properties. The prospects of utilizing the applied field and Hamiltonian structure simultaneously as controls is also explored. The control scenarios I and II indicate the gains offered by algorithmically guided molecular or material discovery for manipulating quantum dynamics phenomenon.     

Photonic reagent control of dynamically homologous quantum systems

The general objective of quantum control is the manipulation of atomic scale physical and chemical phenomena through the application of external control fields. These tailored fields, or photonic reagents, exhibit systematic properties analogous to those of ordinary laboratory reagents. This analogous behavior is explored further here by considering the controlled response of a family of homologous quantum systems to a single common photonic reagent. A level set of dynamically homologous quantum systems is defined as the family that produces the same value(s) for a target physical observable(s) when controlled by a common photonic reagent. This paper investigates the scope of homologous quantum system control using the level set exploration technique (L-SET). L-SET enables the identification of continuous families of dynamically homologous quantum systems. Each quantum system is specified by a point in a hypercube whose edges are labeled by Hamiltonian matrix elements. Numerical examples are presented with simple finite level systems to illustrate the L-SET concepts. Both connected and disconnected families of dynamically homologous systems are shown to exist.     

Quantum optimal control: Hessian analysis of the control landscape

Seeking an effective quantum control entails searching over a landscape defined as the objective as a functional of the control field. This paper considers the problem of driving a state-to-state transition in a finite level quantum system, and analyzes the local topology of the landscape of the final transition probability in terms of the variables specifying the control field. Numerical calculation of the eigenvalues of the Hessian of the transition probability with respect to the control field variables reveals systematic structure in the spectra reflecting the existence of a generic and simple control landscape topology. An illustration shows that the number of nonzero Hessian eigenvalues is determined by the number of quantum states in the system. The Hessian eigenvectors associated with its nonzero eigenvalues are shown to give insight into the cooperative roles of the control variables. The practical consequences of these findings for quantum control are discussed.     

Quantum observable homotopy tracking control

This paper presents a new tracking method where the target observable O(s,T) at the final dynamical time T follows a predefined track P(s) with respect to a homotopy tracking variable s>or=0. The procedure calculates the series of control fields E(s,t) required to accomplish observable homotopy tracking by solving a first-order differential equation in s for the evolution of the control field. Controls produced by this technique render the desired track for all s without encountering field singularities. This paper also extends the technique to the case where the field-free Hamiltonian and dipole moment operator change with s in order to explore the control of new physical systems along the track. Several simulations are presented illustrating the various uses for this quantum tracking control technique.     

Modeling the influence of a laser pulse on the potential energy surface in optimal molecular control theory

Understanding and modeling the interaction between light and matter is essential to the theory of optical molecular control. While the effect of the electric field on a molecule's electronic structure is often not included in control theory, it can be modeled in an optimal control algorithm by a set or toolkit of potential energy surfaces indexed by discrete values of the electric field strength where the surfaces are generated by Born-Oppenheimer electronic structure calculations that directly include the electric field. Using a new optimal control algorithm with a trigonometric mapping to limit the maximum field strength explicitly, we apply the surface-toolkit method to control the hydrogen fluoride molecule. Potential energy surfaces in the presence and absence of the electric field are created with two-electron reduced-density-matrix techniques. The population dynamics show that adjusting for changes in the electronic structure of the molecule beyond the static dipole approximation can be significant for designing a field that drives a realistic quantum system to its target observable.     

Optimal control of multiple two-photon transitions

Optimal strategies to maximize the two-photon transition amplitude have been well studied. However, for a system with multiple intermediate states linking the same initial and final states, the question that how to achieve an optimal population transfer remains nontrivial. In this work, we propose a systematic block scheme to maximize the transition amplitude by explicitly considering the interferences among different transition pathways. The scheme can probably provide a quasi-optimal solution, even considering the uncertainties and noises in experiments. Cases with infinitesimal and finite spectral resolution are both investigated. A special example with all first transition frequencies being larger than half of the two-photon transition frequencies is employed to demonstrate our scheme. The analysis provides valuable insights on how to manipulate the interferences in control of quantum systems.     

Some Mathematical and Algorithmic Challenges in the Control of Quantum Dynamics Phenomena

The theory and practice of control over quantum mechanical phenomena is receiving increasing attention, underscored by striking experimental successes. Nevertheless, many questions of fundamental and practical relevance to the field remain unresolved. With the aim of stimulating further development, this paper formulates a number of theoretical questions, divided into three categories. First, questions related to control law design are discussed, with an emphasis on controllability and optimal control theory. This leads to the second category of open problems relevant to closed loop laboratory implementation of quantum control, including learning and feedback methods. The sensitive dependence of control on basic quantum mechanical interactions motivates the third section, which treats coherent dynamical techniques for identifying the system Hamiltonian. An open issue overarching all of these directions is the need to discover general rules for the control of quantum systems. Although the list of issues raised in this paper is extensive, it should be viewed not as a complete menu for exploration, but rather as a springboard to new challenges as the field evolves.     

Attaining persistent field-free control of open and closed quantum systems

Persistent quantum control (PQC) aims to maintain an observable objective value over a period of time following the action of an applied field. This paper assesses the feasibility of achieving PQC for arbitrary finite-level systems and observables. The analysis is carried out independent of the particular method used for state preparation. The PQC behavior is optimized over the set of physically accessible prepared states for both open and closed systems. The quality of observable value persistence in the postcontrol period was found to vary with the required duration of persistence, the system temperature, the chosen observable operator, and the energy levels of the system. The alignment of a rigid diatomic rotor is studied as a model system. The theoretical estimates of PQC behavior are encouraging and suggest feasible exploration in the laboratory using currently available technology.     

On the relationship between quantum control landscape structure and optimization complexity

It has been widely observed in optimal control simulations and experiments that state preparation is surprisingly easy to achieve, regardless of the dimension N of the system Hilbert space. In contrast, simulations for the generation of targeted unitary transformations indicate that the effort increases exponentially with N. In order to understand such behavior, the concept of quantum control landscapes was recently introduced, where the landscape is defined as the physical objective, as a function of the control variables. The present work explores how the local structure of the control landscape influences the effectiveness and efficiency of quantum optimal control search efforts. Optimizations of state and unitary transformation preparation using kinematic control variables (i.e., the elements of the action matrix) are performed with gradient, genetic, and simplex algorithms. The results indicate that the search effort scales weakly, or possibly independently, with N for state preparation, while the search effort for the unitary transformation objective increases exponentially with N. Analysis of the mean path length traversed during a search trajectory through the space of action matrices and the local structure along this trajectory provides a basis to explain the difference in the scaling of the search effort with N for these control objectives. Much more favorable scaling for unitary transformation preparation arises upon specifying an initial action matrix based on state preparation results. The consequences of choosing a reduced number of control variables for state preparation is also investigated, showing a significant reduction in performance for using fewer than 2N-2 variables, which is consistent with the topological analysis of the associated landscape.     

Hamiltonian theory for vibrational line shapes of atoms adsorbed on surfaces

The vibrational motions of atomic adsorbates on surfaces can be probed by helium atom scattering. The experimental observable is the dynamic structure factor, which shows an inelastic peak around the vibrational frequency of the isolated adsorbates known as the frustrated translational or T-mode peak. In this paper we develop a theory for the line shape of this peak, as well as for its temperature-dependent shift and broadening, based on a Hamiltonian equivalent of the generalized Langevin equation. The theory can be used to infer physical parameters of the adatom-surface interaction, such as the friction coefficient, the barrier height to diffusion, and the anharmonicity parameter. Numerical simulations are used to ascertain the range of validity of the theory, which is also generalized to describe multidimensional systems and to include quantum corrections. We compare the theoretical predictions for the shift and broadening with experimental results for the Na/Cu(001) system, showing quantitative agreement within experimental resolution.     

Nuclear magnetic resonance proton dipolar order relaxation in thermotropic liquid crystals: a quantum theoretical approach

By means of the Jeener-Broekaert nuclear magnetic resonance pulse sequence, the proton spin system of a liquid crystal can be prepared in quasiequilibrium states of high dipolar order, which relax to thermal equilibrium with the molecular environment with a characteristic time (T1D). Previous studies of the Larmor frequency and temperature dependence of T1D in thermotropic liquid crystals, that included field cycling and conventional high-field experiments, showed that the slow hydrodynamic modes dominate the behavior of T1D, even at high Larmor frequencies. This noticeable predominance of the cooperative fluctuations (known as order fluctuations of the director, OFD) could not be explained by standard models based on the spin-lattice relaxation theory in the limit of high temperature (weak order). This fact points out the necessity of investigating the role of the quantum terms neglected in the usual high temperature theory of dipolar order relaxation. In this work, we present a generalization of the proton dipolar order relaxation theory for highly correlated systems, which considers all the spins belonging to correlated domains as an open quantum system interacting with quantum bath. As starting point, we deduce a formulation of the Markovian master equation of relaxation for the statistical spin operator, valid for all temperatures, which is suitable for introducing a dipolar spin temperature in the quantum regime, without further assumptions about the form of the spin-lattice Hamiltonian. In order to reflect the slow dynamics occurring in correlated systems, we lift the usual short-correlation-time assumption by including the average over the motion of the dipolar Hamiltonian together with the Zeeman Hamiltonian into the time evolution operator. In this way, we calculate the time dependence of the spin operators in the interaction picture in a closed form, valid for high magnetic fields, bringing into play the spin-spin interactions within the microscopic time scale. Then, by adopting the spin-temperature density operator to represent the collective state of the spin system, and removing the traditional hypothesis of high temperature, we deduce an expression for the first order quantum contribution to T1D (-1), in terms of spectral densities, with coefficients in form of spin traces. The properties that distinguish our result from the high-temperature T1D (-1) are as follows. (a) It is exclusively associated to cooperative fluctuations. (b) Because of its quantum character, it relies on both considering the lattice degrees of freedom quantum mechanically and including the spin-spin interactions in the microscopic time scale. With regard to the average dipolar Hamiltonian, only the nonsecular part plays a relevant role. (c) Associated with the structure of the spin operator involved in the quantum contribution, a term arises which is proportional to the number of spins in the correlated molecular domains, showing that the quantum contribution may be of macroscopic size in highly correlated systems. When applied to nematic liquid crystals, the new term exhibits the typical nu(-1/2) Larmor frequency dependence through the spectral density of the OFD, in consistence with the experimental results.     

Quantum control mechanism analysis through field based Hamiltonian encoding

Optimal control of quantum dynamics in the laboratory is proving to be increasingly successful. The control fields can be complex, and the mechanisms by which they operate have often remained obscure. Hamiltonian encoding (HE) has been proposed as a method for understanding mechanisms in quantum dynamics. In this context mechanism is defined in terms of the dominant quantum pathways leading to the final state of the controlled system. HE operates by encoding a special modulation into the Hamiltonian and decoding its signature in the dynamics to determine the dominant pathway amplitudes. Earlier work encoded the modulation directly into the Hamiltonian operators. This present work introduces the alternative scheme of field based HE, where the modulation is encoded into the control field and not directly into the Hamiltonian operators. This distinct form of modulation yields a new perspective on mechanism and is computationally faster than the earlier approach. Field based encoding is also an important step towards a laboratory based algorithm for HE as it is the only form of encoding that may be experimentally executed. HE is also extended to cover systems with noise and uncertainty and finally, a hierarchical algorithm is introduced to reveal mechanism in a stepwise fashion of ever increasing detail as desired. This new hierarchical algorithm is an improvement over earlier approaches to HE where the entire mechanism was determined in one stroke. The improvement comes from the use of less complex modulation schemes, which leads to fewer evaluations of Schrodinger's equation. A number of simulations are presented on simple systems to illustrate the new field based encoding technique for mechanism assessment.     

Effect of surface modes on the six-dimensional molecule-surface scattering dynamics of H2-Cu(100) and D2-Cu(111) systems

We include the phonon modes originating from the three layers of Cu(100)/Cu(111) surface atoms on the dynamics of molecular [H(2)(v,j)/D(2)(v,j)] degrees of freedom (DOFs) through a mean field approach, where the surface temperature is incorporated into the effective Hamiltonian (potential) either by considering Boltzmann probability (BP) or by including the Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes. The formulation of effective potential has been carried out by invoking the expression of transition probabilities for phonon modes known from the "stochastic" treatment of linearly forced harmonic oscillator (LFHO). We perform four-dimensional (4D?2D) as well as six-dimensional (6D) quantum dynamics on a parametrically time and temperature-dependent effective Hamiltonian to calculate elastic/inelastic scattering cross-section of the scattered molecule for the H(2)(v,j)-Cu(100) system, and dissociative chemisorption-physisorption for both H(2)(v,j)-Cu(100) and D(2)(v,j)-Cu(111) systems. Calculated sticking probabilities by either 4D?2D or 6D quantum dynamics on an effective potential constructed by using BP factor for the initial state distribution of the phonon modes could not show any surface temperature dependence. In the BEP case, (a) both 4D?2D and 6D quantum dynamics demonstrate that the phonon modes of the Cu(100) surface affect the state-to-state transition probabilities of the scattered H(2) molecule substantially, and (b) the sticking probabilities due to the collision of H(2) on Cu(100) and D(2) on Cu(111) surfaces show noticeable and substantial change, respectively, as function of surface temperature only when the quantum dynamics of all six molecular DOFs are treated in a fully correlated manner (6D).     

Quantum control mechanism analysis through field based Hamiltonian encoding: a laboratory implementable algorithm

While closed-loop control of quantum dynamics in the laboratory is proving to be broadly successful, the control mechanisms induced by the fields are often left obscure. Hamiltonian encoding (HE) was originally introduced as a method for understanding mechanisms in quantum dynamics in the context of computational simulations, based on access to the system wavefunction. As a step towards laboratory implementation of HE, this paper addresses the issues raised by the use of observables rather than the wavefunction in HE. The goal of laboratory based HE is to obtain an understanding of control mechanism through a sequence of systematic control experiments, whose collective information can identify the underlying control mechanism defined as the set of significant amplitudes connecting the initial and final states. Mechanism is determined by means of observing the dynamics of special sequences of system Hamiltonians encoded through the control field. The proposed algorithm can handle complex systems, operates with no recourse to dynamical simulations, and functions with limited understanding of the system Hamiltonian. As with the closed-loop control experiments, the HE control mechanism identification algorithm performs a new experiment each time the dynamical outcome from an encoded Hamiltonian is called for. This paper presents the basic HE algorithm in the context of physical systems described by a finite dimensional Hilbert space. The method is simulated with simple models, and the extension to more complex systems is discussed.     

Dimension reduction by balanced truncation: application to light-induced control of open quantum systems

In linear control, balanced truncation is known as a powerful technique to reduce the state-space dimension of a system. Its basic principle is to identify a subspace of jointly easily controllable and observable states and then to restrict the dynamics to this subspace without changing the overall response of the system. This work deals with a first application of balanced truncation to the control of open quantum systems which are modeled by the Liouville-von Neumann equation within the Lindblad formalism. Generalization of the linear theory has been proposed to cope with the bilinear terms arising from the coupling between the control field and the quantum system. As an example we choose the dissipative quantum dynamics of a particle in an asymmetric double well potential driven by an external control field, monitoring population transfer between the potential wells as a control target. The accuracy of dimension reduction is investigated by comparing the populations obtained for the truncated system versus those for the original system. The dimension of the model system can be reduced very efficiently where the degree of reduction depends on temperature and relaxation rate.     

A quantum theory of chemical processes and reaction rates based on diabatic electronic functions coupled in an external field

The method discussed in this work provides a theoretical framework where simple chemical reactions resemble any other standard quantum process, i.e., a transition in quantum state mediated by the electromagnetic field. In our approach, quantum states are represented as a superposition of electronic diabatic basis functions, whose amplitudes can be modulated by the field and by the external control of nuclear configurations. Using a one-dimensional three-state model system, we show how chemical structure and dynamics can be represented in terms of these control parameters, and propose an algorithm to compute the reaction probabilities. Our analysis of effective energy barriers generalizes previous ideas on structural similarity between reactant, and product, and transition states using the geometry of conventional reaction paths. In the present context, exceptions to empirical rules such as the Hammond postulate appear as effects induced by the environment that supplies the external field acting on the quantum system.