Real-time simulation of flux qubits used for quantum annealing

The real-time dynamics of systems with up to three SQUIDs is studied by numerically solving the time-dependent Schrödinger equation. The numerical results are used to scrutinize the mapping of the flux degrees of freedom onto two-level systems (the qubits) as well as the performance of the intermediate SQUID as a tunable coupling element. It is shown that the two-level representation yields a good description of the flux dynamics during quantum annealing, and the presence of the tunable coupling element does not have negative effects on the overall performance. Additionally, data obtained from a two-level spin dynamics simulation of quantum annealing is compared to experimental data produced by the D-Wave 2000Q quantum annealer. The effects of finite temperature are incorporated in the simulation by coupling the qubit-system to a bath of spin-1/2 particles. It is shown that including an environment modeled as non-interacting two-level systems that couple only to the qubits can produce data which matches the experimental data much better than the simulation data of the isolated qubits, and better than data obtained from a simulation including an environment modeled as interacting two-level systems coupling to the qubits.