Simulation of the elementary evolution operator with the motional states of an ion in an anharmonic trap

TL;DR

This paper demonstrates how the motional states of a trapped ion in an anharmonic potential can be used to simulate the time-dependent Schrödinger equation, enabling quantum simulation of wave packet dynamics.

Contribution

It introduces a method to simulate quantum wave packet evolution using ion motional states and optimal control theory for the driving fields.

Findings

Successful discretization of wave packets on ion motional states

Implementation of quantum gates via optimal control theory

Simulation stability tested against decoherence effects

Abstract

Following a recent proposal of L. Wang and D. Babikov, J. Chem. Phys. 137, 064301 (2012), we theoretically illustrate the possibility of using the motional states of a $Cd^+$ ion trapped in a slightly anharmonic potential to simulate the single-particle time-dependent Schr\"odinger equation. The simulated wave packet is discretized on a spatial grid and the grid points are mapped on the ion motional states which define the qubit network. The localization probability at each grid point is obtained from the population in the corresponding motional state. The quantum gate is the elementary evolution operator corresponding to the time-dependent Schr\"odinger equation of the simulated system. The corresponding matrix can be estimated by any numerical algorithm. The radio-frequency field able to drive this unitary transformation among the qubit states of the ion is obtained by multi-target optimal control theory. The ion is assumed to be cooled in the ground motional state and the preliminary step consists in initializing the qubits with the amplitudes of the initial simulated wave packet. The time evolution of the localization probability at the grids points is then obtained by successive applications of the gate and reading out the motional state population. The gate field is always identical for a given simulated potential, only the field preparing the initial wave packet has to be optimized for different simulations. We check the stability of the simulation against decoherence due to fluctuating electric fields in the trap electrodes by applying dissipative Lindblad dynamics.