We develop monotonically convergent free-time and fixed end-point optimal control theory (OCT) in the density-matrix representation to deal with quantum systems showing dissipation. Our theory is more general and flexible for tailoring optimal laser pulses in order to control quantum dynamics with dissipation than the conventional fixed-time and fixed end-point OCT in that the optimal temporal duration of laser pulses can also be optimized exactly. To show the usefulness of our theory, it is applied to the generation and maintenance of the vibrational entanglement of carbon monoxide adsorbed on the copper (100) surface, CO/Cu(100). We demonstrate the numerical results and clarify how to combat vibrational decoherence as much as possible by the tailored shapes of the optimal laser pulses. It is expected that our theory will be general enough to be applied to a variety of dissipative quantum dynamics systems because the decoherence is one of the quantum phenomena sensitive to the temporal duration of the quantum dynamics.
We develop monotonically convergent free-time and fixed end-point optimal control theory (OCT) in the density-matrix representation to deal with quantum systems showing dissipation. Our theory is more general and flexible for tailoring optimal laser pulses in order to control quantum dynamics with dissipation than the conventional fixed-time and fixed end-point OCT in that the optimal temporal duration of laser pulses can also be optimized exactly. To show the usefulness of our theory, it is applied to the generation and maintenance of the vibrational entanglement of carbon monoxide adsorbed on the copper (100) surface, CO/Cu(100). We demonstrate the numerical results and clarify how to combat vibrational decoherence as much as possible by the tailored shapes of the optimal laser pulses. It is expected that our theory will be general enough to be applied to a variety of dissipative quantum dynamics systems because the decoherence is one of the quantum phenomena sensitive to the temporal duration of the quantum dynamics.
We have constructed free-time and fixed end-point optimal control theory for quantum systems and applied it to entanglement generation between rotational modes of two polar molecules coupled by dipole-dipole interaction. The motivation of the present work is to solve optimal control problems more flexibly by extending the popular fixed time and fixed end-point optimal control theory for quantum systems to free-time and fixed end-point optimal control theory. As a demonstration, the theory that we have constructed in this paper will be applied to entanglement generation in rotational modes of NaCl-NaBr polar molecular systems that are sensitive to the strength of entangling interactions. Our method will significantly be useful for the quantum control of nonlocal interaction such as entangling interaction, which depends crucially on the strength of the interaction or the distance between the two molecules, and other general quantum dynamics, chemical reactions, and so on.
Open-system quantum optimal control theory for optical control of the dynamics of a quantum system in contact with a dissipative bath is used here for explicitly time-dependent target operators, O(t). Global and local control strategies are combined in a novel algorithm by defining a set of time slices, into which the total control time is subdivided. The optimization then proceeds locally forward in time from subinterval to subinterval, while within each subinterval global control theory is used with iterative forward-backward propagation. The subintervals are connected by appropriate boundary conditions. In the present paper, all operators are represented in the basis of the eigenstates of the field-free system Hamiltonian. The algorithm is first applied to and its computational performance tested for a two-level system with energy and phase relaxation, and later extended to a many-level model. Model parameters are chosen to represent the IR pulse excitation of the adsorbate-surface stretch mode of vibrationally relaxing CO on a Cu(100) surface. Various time-dependent targets are formulated to achieve (i) population inversion, (ii) the creation of a wavepacket, and (iii) overtone excitation by "ladder climbing."
Quantum optimal control theory is applied to find optimal pulses for controlling the motion of an ion and a molecule for two different applications. Those optimal pulses enable the control of the dynamics of the system by driving the atom or the molecule from an initial state to desired states.The evolution equations obtained by means of the quantum optimal control theory are resolved iteratively using a monotonic convergent algorithm. A number of simulation parameters are varied in order to get the optimal pulses including the duration of the pulses, the time step of the time grid, a penalty factor that limits the maximal intensity of the fields, and a guess pulse which is used to start the optimal control.The optimal pulses obtained for each application are analyzed by Fourier transform, and also by looking at the time evolution of the populations that they generate in the system.The first application is the preparation of specific vibrational states of acetylene that are usually not reachable from the ground state, and that would remain unpopulated by usual spectroscopy. Relevant state energies and transition dipole moments are extracted from the experimental literature and especially from the global acetylene Hamiltonian conferring an uncommon precision to the control simulation. The control starts from the ground state. The target states belongs to the polyad Ns=1, Nr=5 of acetylene which includes two vibrational dark states and one vibrational bright state. First, the simulation is performed with the Schrödinger equation and in a second step, with the Liouville--von Neumann equation, as mixed states are prepared. Indeed, the control starts from a Boltzmann distribution of population in the rotational levels of the vibrational ground state chosen in order to simulate an experimental condition. But the distribution is truncated to limit the computational effort. One of the dark states appears to be a potential target for a realistic experimental investigation because the average population of the Rabi oscillation remains high and decoherence is expected to be weak. The optimal pulses obtained have a high fidelity, have a spectrum with well-resolved peak frequencies, and their experimental feasibility seems achievable within the current abilities of experimental laboratories.The second application is to propose an experimental realization of a microscopic physical device able to simulate quantum dynamics. The idea is to use the motional states of a Cd^+ ion trapped in an anharmonic potential to realize a quantum dynamics simulator of a single-particle Schrödinger equation. In this way, the motional states store the information and the optimal pulse manipulates this information to realize operations. In the present case, the simulated dynamics was the propagation of a wave packet in a harmonic potential. Starting from an initial quantum state, the pulse acts on the system to modify the motional states of the ion in such a way that the final superposition of motional states corresponds to the results of the dynamics. This simulation is performed with the Liouville--von Neumann equation and also with the Lindblad equation as dissipation is included to test the robustness of the pulse against perturbations of the potential. The optimal pulses that are obtained have a high fidelity which shows that the ion trap system has correctly realized the quantum dynamics simulation. The optimal pulses are valid for any initial condition if the potential of the simulation or the mass of the propagated wave packet is unchanged.
In the context of molecular quantum computation the optimal control theory (OCT) is used to obtain shaped laser pulses for high-fidelity control of vibrational qubits. Optimization is done in time domain and the OCT algorithm varies values of electric field in each time step independently, tuning hundreds of thousands of parameters to find one optimal solution. Such flexibility is not available in experiments, where pulse shaping is done in frequency domain and the number of "tuning knobs" is much smaller. The question of possible experimental interpretations of theoretically found OCT solutions arises. In this work we analyze very accurate optimal pulse that we obtained for implementing quantum gate CNOT for the two-qubit system encoded into the exited vibrational states of thiophosgene molecule. Next, we try to alter this pulse by reducing the number of available frequency channels and intentionally introducing systematic and random errors (in frequency domain, by modifying the values of amplitudes and phases of different frequency components). We conclude that a very limited number of frequency components (only 32 in the model of thiophosgene) are really necessary for accurate control of the vibrational two-qubit system, and such pulses can be readily constructed using OCT. If the amplitude and phase errors of different frequency components do not exceed ±3% of the optimal values, one can still achieve accurate transformations of the vibrational two-qubit system, with gate fidelity of CNOT exceeding 0.99.
Optimal control theory (OCT) has wide applicability in the quest for seeking optimal laser pulses in coherent control of quantum phenomena. We have revised OCT making it possible to distill robust electric fields realizable via femtosecond pulse-shaping technology. The simple pulse structures offer direct physical insight into the control mechanism and selectively isolate the most stable pathways involving only few frequencies. Moreover, a generalized functional makes a vital extension in including the laser system used in control experiments.
Optimal control within the density matrix formalism is applied to the creation of a specified superposition state in condensed phases. The primary system modeled by a displaced harmonic oscillator is surrounded by a boson heat bath, the dynamics of which is described by a non-Markovian master equation. A newly developed monotonically convergent algorithm is used to solve the pulse design equations. The control mechanisms are strongly dependent on the bath correlation time that is characterized by the inverse of an exponential decay constant, γ. If the correlation time is shorter than the temporal width of a typical subpulse involved in an optimal pulse, the solution is reduced to that in the Markovian case. If we assume a longer correlation time, by weighing less physical significance on the penalty due to pulse fluence, an optimal pulse with high intensity is obtained, the temporal width of which approaches ∼1/γ. We also see considerable changes in the shape of the optimal pulse with increasing intensity, suggesting that strong fields open up control mechanisms that are qualitatively different from those in weak fields.