The number of topological defects created in a system driven through a quantum phase transition exhibits a power-law scaling with the driving time. This universal scaling law is the key prediction of the Kibble-Zurek mechanism (KZM), and testing it using a hardware-based quantum simulator is a coveted goal of quantum information science. Here we provide such a test using quantum annealing. Specifically, we report on extensive experimental tests of topological defect formation via the one-dimensional transverse-field Ising model on two different D-Wave quantum annealing devices. We find that the quantum simulator results can indeed be explained by the KZM for open-system quantum dynamics with phase-flip errors, with certain quantitative deviations from the theory likely caused by factors such as random control errors and transient effects. In addition, we probe physics beyond the KZM by identifying signatures of universality in the distribution and cumulants of the number of kinks and their decay, and again find agreement with the quantum simulator results. This implies that the theoretical predictions of the generalized KZM theory, which assumes isolation from the environment, applies beyond its original scope to an open system. We support this result by extensive numerical computations. To check whether an alternative, classical interpretation of these results is possible, we used the spin-vector Monte Carlo model, a candidate classical description of the D-Wave device. We find that the degree of agreement with the experimental data from the D-Wave annealing devices is better for the KZM, a quantum theory, than for the classical spin-vector Monte Carlo model, thus favoring a quantum description of the device. Our work provides an experimental test of quantum critical dynamics in an open quantum system, and paves the way to new directions in quantum simulation experiments.

The Kibble–Zurek mechanism (KZM) describes the dynamics across a phase transition leading to the formation of topological defects, such as vortices in superfluids and domain walls in spin systems. Here, we experimentally probe the distribution of kink pairs in a one-dimensional quantum Ising chain driven across the paramagnet-ferromagnet quantum phase transition, using a single trapped ion as a quantum simulator in momentum space. The number of kink pairs after the transition follows a Poisson binomial distribution, in which all cumulants scale with a universal power law as a function of the quench time in which the transition is crossed. We experimentally verified this scaling for the first cumulants and report deviations due to noise-induced dephasing of the trapped ion. Our results establish that the universal character of the critical dynamics can be extended beyond KZM, which accounts for the mean kink number, to characterize the full probability distribution of topological defects.

We consider the number distribution of topological defects resulting from the finite-time crossing of a continuous phase transition and identify signatures of universality beyond the mean value, predicted by the Kibble-Zurek mechanism. Statistics of defects follows a binomial distribution with $\mathcal{N}$ Bernouilli trials associated with the probability of forming a topological defect at the locations where multiple domains merge. All cumulants of the distribution are predicted to exhibit a common universal power-law scaling with the quench time in which the transition is crossed. Knowledge of the distribution is used to discuss the onset of adiabatic dynamics and bound rare events associated with large deviations.

In the nonadiabatic dynamics across a quantum phase transition, the Kibble-Zurek mechanism predicts that the formation of topological defects is suppressed as a universal power law with the quench time. In inhomogeneous systems, the critical point is reached locally and causality reduces the effective system size for defect formation to regions where the velocity of the critical front is slower than the sound velocity, favoring adiabatic dynamics. The reduced density of excitations exhibits a much steeper dependence on the quench rate and is also described by a universal power law that we demonstrated in a quantum Ising chain.

Traversing a continuous phase transition at a finite rate leads to the breakdown of adiabatic dynamics and the formation of topological defects, as predicted by the celebrated Kibble-Zurek mechanism (KZM). We investigate universal signatures beyond the KZM, by characterizing the distribution of vortices generated in a thermal quench leading to the formation of a holographic superconductor. The full counting statistics of vortices is described by a binomial distribution, in which the mean value is dictated by the KZM and higher-order cumulants share the universal power-law scaling with the quench time. Extreme events associated with large fluctuations no longer exhibit a power-law behavior with the quench time and are characterized by a universal form of the Weibull distribution for different quench rates.