The Journal of Chemical Physics

Optimal control simulation of the Deutsch-Jozsa algorithm in a two-dimensional double well coupled to an environment

Mamadou Ndong, David Lauvergnat, Xavier Chapuisat, Michèle Desouter-Lecomte

Published Jun 26, 2007• 34 citations

Abstract

The quantum Deutsch-Jozsa algorithm is implemented by using vibrational modes of a two-dimensional double well. The laser fields realizing the different gates (NOT, CNOT, and HADAMARD) on the two-qubit space are computed by the multitarget optimal control theory. The stability of the performance index is checked by coupling the system to an environment. Firstly, the two-dimensional subspace is coupled to a small number Nb of oscillators in order to simulate intramolecular vibrational energy redistribution. The complete (2+Nb)D problem is solved by the coupled harmonic adiabatic channel method which allows including coupled modes up to Nb=5. Secondly, the computational subspace is coupled to a continuous bath of oscillators in order to simulate a confined environment expected to be favorable to achieve molecular computing, for instance, molecules confined in matrices or in a fullerene. The spectral density of the bath is approximated by an Ohmic law with a cutoff for some hundreds of cm(-1). The time scale of the bath dynamics (of the order of 10 fs) is then smaller than the relaxation time and the controlled dynamics (2 ps) so that Markovian dissipative dynamics is used.

DOI: 10.1063/1.2743429

Related papers

Vibrational molecular quantum computing: Basis set independence and theoretical realization of the Deutsch–Jozsa algorithm

Carmen M. Tesch, Regina de Vivie‐Riedle

The phase of quantum gates is one key issue for the implementation of quantum algorithms. In this paper we first investigate the phase evolution of global molecular quantum gates, which are realized by optimally shaped femtosecond laser pulses. The specific laser fields are calculated using the multitarget optimal control algorithm, our modification of the optimal control theory relevant for application in quantum computing. As qubit system we use vibrational modes of polyatomic molecules, here the two IR-active modes of acetylene. Exemplarily, we present our results for a Pi gate, which shows a strong dependence on the phase, leading to a significant decrease in quantum yield. To correct for this unwanted behavior we include pressure on the quantum phase in our multitarget approach. In addition the accuracy of these phase corrected global quantum gates is enhanced. Furthermore we could show that in our molecular approach phase corrected quantum gates and basis set independence are directly linked. Basis set independence is also another property highly required for the performance of quantum algorithms. By realizing the Deutsch-Jozsa algorithm in our two qubit molecular model system, we demonstrate the good performance of our phase corrected and basis set independent quantum gates.

The Journal of Chemical PhysicsDec 15, 2004106 citations

Quantum computing using molecular vibrational and rotational modes

Kaito Shioya, Kenji Mishima, Koichi Yamashita

Elementary quantum gates have been constructed using molecular vibrational and rotational modes to encode two qubits and implement the Deutsch–Jozsa algorithm. The molecular system chosen was 12C16O. Several phase-correct global quantum gates were constructed using optimal control theory (OCT). The gate fidelity was more than 90%. On implementing the Deutsch–Jozsa algorithm combining these elementary gates, a fidelity greater than 96.11% was obtained, which indicates that vibrational and rotational qubits are promising for performing quantum algorithms.

Molecular PhysicsMay 10, 200749 citations

Vibrational computing: Simulation of a full adder by optimal control

L. Bomble, David Lauvergnat, F. Remacle +1 more

Within the context of vibrational molecular quantum computing, we investigate the implementation of a full addition of two binary digits and a carry that provides the sum and the carry out. Four qubits are necessary and they are encoded into four different normal vibrational modes of a molecule. We choose the bromoacetyl chloride molecule because it possesses four bright infrared active modes. The ground and first excited states of each mode form the one-qubit computational basis set. Two approaches are proposed for the realization of the full addition. In the first one, we optimize a pulse that implements directly the entire addition by a single unitary transformation. In the second one, we decompose the full addition in elementary quantum gates, following a scheme proposed by Vedral et al. [Phys. Rev. A 54, 147 (1996)]. Four elementary quantum gates are necessary, two two-qubit CNOT gates (controlled NOT) and two three-qubit TOFFOLI gates (controlled-controlled NOT). All the logic operations consist in one-qubit flip. The logic implementation is therefore quasiclassical and the readout is based on a population analysis of the vibrational modes that does not take the phases into account. The fields are optimized by the multitarget extension of the optimal control theory involving all the transformations among the 2(4) qubit states. A single cycle of addition without considering the preparation or the measure or copy of the result can be carried out in a very competitive time, on a picosecond time scale.

The Journal of Chemical PhysicsFeb 14, 200838 citations

Relaxation dynamics in quantum dissipative systems: The microscopic effect of intramolecular vibrational energy redistribution

Llinersy Uranga-Piña, Jean Christophe Tremblay

We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.

The Journal of Chemical PhysicsAug 15, 201410 citations

Using quantum optimal control to drive intramolecular vibrational redistribution and to perform quantum computing

Ludovic Santos, Nathalie Vaeck

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.

Nov 28, 20170 citations

Complexity and simplicity of optimal control theory pulses shaped for controlling vibrational qubits

Dmytro Shyshlov, Dmitri Babikov

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.

The Journal of Chemical PhysicsNov 21, 20127 citations