Ramil Nigmatullin

Spectral functions measured with angle-resolved photoemission spectroscopy (ARPES) provide key insight to elucidate the band structure of materials. Comparison with theory requires computing dynamical one-point functions in some equilibrium state, which can be classically challenging. Their measurement on quantum computers poses multiple problems and comes with a large sampling overhead when standard techniques are used. We introduce an efficient way of measuring spectral functions on a quantum computer by directly modeling the interaction of the system with the environment involved in ARPES experiments. We develop quantum circuits whose local expectation values are proportional to the spectral function $A(k,ω)$ for all momentum $k$ and a specific chosen frequency $ω$. Although coming with a qubit and two-qubit gate overhead, our approach requires $O(N)$ times less sampling than previous approaches, translating into a factor $O(N)$ faster in runtime, and is particularly adapted to ion-trap quantum computers. The algorithm requires to implement a fermionic Fourier transform (FFT). We write out an efficient gate decomposition for generic radix-$n$ FFT and benchmark it on hardware for radix-$3$ on $27$ qubits. We finally demonstrate our algorithm on a Quantinuum System Model H2 ion-trap system, computing the spectral function on a one-dimensional system of $27$ sites, using $54$ qubits.

The Fermi-Hubbard model is the starting point for the simulation of many strongly correlated materials, including high-temperature superconductors, whose modelling is a key motivation for the construction of quantum simulation and computing devices. However, the detection of superconducting pairing correlations has so far remained out of reach, both because of their off-diagonal character - which makes them inaccessible to local density measurements - and because of the difficulty of preparing superconducting states. Here, we report measurement of significant pairing correlations in three different regimes of Fermi-Hubbard models simulated on Quantinuum's Helios trapped-ion quantum computer. Specifically, we measure non-equilibrium pairing induced by an electromagnetic field in the half-filled square lattice model, d-wave pairing in an approximate ground state of the checkerboard Hubbard model at $1/6$-doping, and s-wave pairing in a bilayer model relevant to nickelate superconductors. These results show that a quantum computer can reliably create and probe physically relevant states with superconducting pairing correlations, opening a path to the exploration of superconductivity with quantum computers.

The utility of solving the Fermi-Hubbard model has been estimated in the billions of dollars. Digital quantum computers can in principle address this task, but have so far been limited to quasi one-dimensional models. This is because of exponential overheads caused by the interplay of noise and the non-locality of the mapping between fermions and qubits. Here, we show experimentally that a recently developed local encoding can overcome this problem. We develop a new compilation scheme, called "corner hopping", that reduces the cost of simulating fermionic hopping by 42% which allows us to conduct the largest digital quantum simulations of a fermionic model to date, using a trapped ion quantum computer to prepare adiabatically the ground state of a 6 x 6 spinless Fermi-Hubbard model encoded in 48 physical qubits. We also develop two new error mitigation schemes for systems with conserved quantities, one based on local postselection and one on extrapolation of local observables. Our results suggest that Fermi-Hubbard models beyond classical simulability can be addressed by digital quantum computers without large increases in gate fidelity.

Calculating the equilibrium properties of condensed matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude $\langle ψ| e^{-i \hat H t}|ψ\rangle$ from initial states $|ψ\rangle$ and a time evolution under the Hamiltonian $\hat H$ up to short times $t$. In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a $16$-site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm.

We study the non-equilibrium dynamics of two dimensional planar ion Coulomb crystals undergoing a structural buckling transition to a three plane configuration, driven by a reduction of the transverse confining frequency. This phase transition can be theoretically modeled using a mapping to a two dimensional Ginzburg-Landau theory with complex order parameter field. We demonstrate that finite rate quenches result in creation of stable topological vortices, which are localized point regions around which the phase of the order parameter field winds a multiple of 2π. The density of the defects as a function of quench rate is investigated using molecular dynamics simulations and its scaling is shown to be consistent with Kibble-Zurek theory of defect formation. Following the quench, the annihilation of vortex and anti-vortex pairs results in the relaxation of defect density that follows a diffusive scaling with a logarithmic correction. This work highlights the potential for investigating complex non-equilibrium statistical physics of topological defects in an experimentally accessible ion trap setting.

We study the non-equilibrium dynamics of second-order phase transitions in a simplified Ginzburg-Landau model using the Fokker-Planck formalism. In particular, we focus on deriving the Kibble-Zurek scaling laws that dictate the dependence of spatial correlations on the quench rate. In the limiting cases of overdamped and underdamped dynamics, the Fokker-Planck method confirms the theoretical predictions of the Kibble-Zurek scaling theory. The developed framework is computationally efficient, enables the prediction of finite-size scaling functions and is applicable to microscopic models as well as their hydrodynamic approximations. We demonstrate this extended range of applicability by analyzing the non-equilibrium linear to zigzag structural phase transition in ion Coulomb crystals confined in a trap with periodic boundary conditions.

Optically linked ion traps are promising as components of network-based quantum technologies, including communication systems and modular computers. Experimental results achieved to date indicate that the fidelity of operations within each ion trap module will be far higher than the fidelity of operations involving the links; fortunately internal storage and processing can effectively upgrade the links through the process of purification. Here we perform the most detailed analysis to date on this purification task, using a protocol which is balanced to maximise fidelity while minimising the device complexity and the time cost of the process. Moreover we 'compile down' the quantum circuit to device-level operations including cooling and shutting events. We find that a linear trap with only five ions (two of one species, three of another) can support our protocol while incorporating desirable features such as 'global control', i.e. laser control pulses need only target an entire zone rather than differentiating one ion from its neighbour. To evaluate the capabilities of such a module we consider its use both as a universal communications node for quantum key distribution, and as the basic repeating unit of a quantum computer. For the latter case we evaluate the threshold for fault tolerant quantum computing using the surface code, finding acceptable fidelities for the 'raw' entangling link as low as 83% (or under 75% if an additional ion is available).

We use laser-cooled ion Coulomb crystals in the well-controlled environment of a harmonic radiofrequency ion trap to investigate phase transitions and defect formation. Topological defects in ion Coulomb crystals (kinks) have been recently proposed for studies of nonlinear physics with solitons and as carriers of quantum information. Defects form when a symmetry breaking phase transition is crossed nonadiabatically. For a second order phase transition, the Kibble-Zurek mechanism predicts that the formation of these defects follows a power law scaling in the rate of the transition. We demonstrate a scaling of defect density and describe kink dynamics and stability. We further discuss the implementation of mass defects and electric fields as first steps toward controlled kink preparation and manipulation.

Symmetry breaking phase transitions play an important role in nature. When a system traverses such a transition at a finite rate, its causally disconnected regions choose the new broken symmetry state independently. Where such local choices are incompatible, defects will form with densities predicted to follow a power law scaling in the rate of the transition. The importance of this Kibble-Zurek mechanism (KZM) ranges from cosmology to condensed matter [1-4]. In previous tests in homogeneous systems, defect formation was seen, but weak dependence on the transition rate and limited control of external parameters so far prevented tests of KZM scaling. As recently predicted [5-9], in inhomogeneous systems propagation of the critical front enhances the role of causality and steepens scaling of defect density with the transition rate. We use ion Coulomb crystals in a harmonic trap to demonstrate, for the first time, scaling of the number of topological defects with the transition rate - the central prediction of KZM - in a well-controlled environment.

We study the nonequilibrium dynamics of the linear to zigzag structural phase transition exhibited by an ion chain confined in a trap with periodic boundary conditions. The transition is driven by reducing the transverse confinement at a finite quench rate, which can be accurately controlled. This results in the formation of zigzag domains oriented along different transverse planes. The twists between different domains can be stabilized by the topology of the trap and under laser cooling the system has a chance to relax to a helical chain with nonzero winding number. Molecular dynamics simulations are used to obtain a large sample of possible trajectories for different quench rates. The scaling of the average winding number with different quench rates is compared to the prediction of the Kibble-Zurek theory, and a good quantitative agreement is found.