Here are some of the papers and work I’ve been involved with

Quantum Computing

Classical Shadows for Quantum Process Tomography on Near-term Quantum Computers

Ryan Levy, Di Luo, Bryan K. Clark

Quantum process tomography is a powerful tool for understanding quantum channels and characterizing properties of quantum devices. Inspired by recent advances using classical shadows in quantum state tomography[1], we have developed a classical shadow method, ShadowQPT, for quantum process tomography. ShadowQPT allows for the reconstruction of the Choi matrix for unitary and non-unitary processes including an efficient reconstruction of fixed-sized reduced processes; we also show how to predict the overlap between any arbitrary state and the output of the quantum channel on a different arbitrary state. We introduce both a scheme using ancilla qubits as well as a two-sided scheme with unitaries before and after the channel. A number of additional approximations and improvements are developed including the use of a pair-factorized Clifford shadow and a series of post-processing techniques which significantly enhance the accuracy for recovering the quantum channel. Both the theoretical scaling for large systems and the practicality of using shadow tomography on NISQ-era hardware are considered. Our algorithms have been implemented with both Pauli and Clifford measurements on the IonQ trapped ion quantum computer for quantum processes up to n=4 qubits (equivalent to the experimental complexity of n=8 qubits for quantum state tomography) and achieved good performance.


Tensor Network Methods

Entanglement Entropy Transitions with Random Tensor Networks

Ryan Levy, Bryan K. Clark

Entanglement is a key quantum phenomena and understanding transitions between phases of matter with different entanglement properties are an interesting probe of quantum mechanics. We numerically study a model of a 2D tensor network proposed to have an entanglement entropy transition first considered by Vasseur et al.[Phys. Rev. B 100, 134203 (2019)]. We find that by varying the bond dimension of the tensors in the network we can observe a transition between an area and volume phase with a logarithmic critical point around \(D\approx 2\). We further characterize the critical behavior measuring a critical exponent using entanglement entropy and the tripartite quantum mutual information, observe a crossover from a 'nearly pure' to entangled area law phase using the the distributions of the entanglement entropy and find a cubic decay of the pairwise mutual information at the transition. We further consider the dependence of these observables for different Rényi entropy. This work helps further validate and characterize random tensor networks as a paradigmatic examples of an entanglement transition.


Distributed-Memory DMRG via Sparse and Dense Parallel Tensor Contractions

Ryan Levy, Edgar Solomonik, Bryan K. Clark

The Density Matrix Renormalization Group (DMRG) algorithm is a powerful tool for solving eigenvalue problems to model quantum systems. DMRG relies on tensor contractions and dense linear algebra to compute properties of condensed matter physics systems. However, its efficient parallel implementation is challenging due to limited concurrency, large memory footprint, and tensor sparsity. We mitigate these problems by implementing two new parallel approaches that handle block sparsity arising in DMRG, via Cyclops, a distributed memory tensor contraction library. We benchmark their performance on two physical systems using the Blue Waters and Stampede2 supercomputers. Our DMRG performance is improved by up to 5.9X in runtime and 99X in processing rate over ITensor, at roughly comparable computational resource use. This enables higher accuracy calculations via larger tensors for quantum state approximation. We demonstrate that despite having limited concurrency, DMRG is weakly scalable with the use of efficient parallel tensor contraction mechanisms.

SC’20 Proceedings arXiv:2007.05540

Topology and the one-dimensional Kondo-Heisenberg model

Julian May-Mann, Ryan Levy, Rodrigo Soto-Garrido, Gil Young Cho, Bryan K. Clark, Eduardo Fradkin

The Kondo-Heinsberg chain is an interesting model of a strongly correlated system which has a broad superconducting state with pair-density wave (PDW) order. Some of us have recently proposed that this PDW state is a symmetry-protected topological (SPT) state, and the gapped spin sector of the model supports Majorana zero modes. In this work, we reexamine this problem using a combination of numeric and analytic methods. In extensive density matrix renormalization group calculations, we find no evidence of a topological ground state degeneracy or the previously proposed Majorana zero modes in the PDW phase of this model. This result motivated us to reexamine the original arguments for the existence of the Majorana zero modes. A careful analysis of the effective continuum field theory of the model shows that the Hilbert space of the spin sector of the theory does not contain any single Majorana fermion excitations. This analysis shows that the PDW state of the doped 1D Kondo-Heisenberg model is not an SPT with Majorana zero modes.

Phys Rev B arXiv:2002.01483

Quantum Monte Carlo

Mitigating the Sign Problem Through Basis Rotations

Ryan Levy and Bryan K. Clark

Quantum Monte Carlo simulations of quantum many body systems are plagued by the Fermion sign problem. The computational complexity of simulating Fermions scales exponentially in the projection time \(\beta\) and system size. The sign problem is basis dependent and an improved basis, for fixed errors, lead to exponentially quicker simulations. We show how to use sign-free quantum Monte Carlo simulations to optimize over the choice of basis on large two-dimensional systems. We numerically illustrate these techniques decreasing the 'badness' of the sign problem by optimizing over single-particle basis rotations on one and two-dimensional Hubbard systems. We find a generic rotation which improves the average sign of the Hubbard model for a wide range of \(U\) and densities for \(4\times L\) systems. In one example improvement, the average sign (and hence simulation cost at fixed accuracy) for the \(16\times 4\) Hubbard model at \(U/t=4\) and \(n=0.75\) increases by \(\exp\left[8.64(6)\beta\right]\). For typical projection times of \(\beta\gtrapprox 100\), this accelerates such simulation by many orders of magnitude.

Phys Rev Lett arXiv:1907.02076

Implementation of the Maximum Entropy Method for Analytic Continuation

Ryan Levy J.P.F. LeBlanc, and Emanuel Gull

We present `Maxent`, a tool for performing analytic continuation of spectral functions using the maximum entropy method. The code operates on discrete imaginary axis datasets (values with uncertainties) and transforms this input to the real axis. The code works for imaginary time and Matsubara frequency data and implements the 'Legendre' representation of finite temperature Green's functions. It implements a variety of kernels, default models, and grids for continuing bosonic, fermionic, anomalous, and other data. Our implementation is licensed under GPLv2 and extensively documented. This paper shows the use of the programs in detail.

Comp Phys Comm arXiv:1606.00368

Magnetic susceptibility and simulated neutron signal in the two-dimensional Hubbard model

J. P. F. LeBlanc, Shaozhi Li, Xi Chen, Ryan Levy, A. E. Antipov, Andrew J. Millis, and Emanuel Gull

We compute dynamic spin susceptibilities in the two-dimensional Hubbard model usingthe method of dual fermions, and we provide a comparison to lattice Monte Carlo and cluster dynamical mean-field theory. We examine the energy dispersion identified by peaks in Imχ(ω,q), which define spin modes, and we compare the exchange scale and magnon dispersion to neutron experiments on the parent La2CuO4 cuprate. We present the evolution of the spin excitations as a function of Hubbard interaction strengths and doping, and we explore the particle-hole asymmetry of the spin excitations. We also study the correlation lengths and the spin excitation dispersion peak structure, and we find a Y-shaped dispersion similar to neutron results on doped HgBa2CuO4+δ.

Phys Rev B arXiv:1904.10782

Experimental Work

In another life, I worked on analyzing quantum dot formation.

Mechanisms of InAs/GaAs quantum dot formation during annealing of In island

S. Huang, S.J. Kim, Ryan Levy X.Q. Pan, and R.S. Goldman

We have examined the formation mechanisms of InAs quantum dots (QDs) via annealing In islands under As flux. We report two distinct mechanisms, droplet epitaxy (DE) and solid phase epitaxy (SPE), which depend on As surface coverage. On c(4 × 4) GaAs surfaces, QDs form by DE. For c(4 × 4)α, one-to-one conversion from In islands to InAs QDs is observed. For c(4 × 4)β, lower densities of larger QDs are observed, presumably due to enhanced In surface diffusion in the absence of metastable Ga-As dimers. For the As capped surface, In deposition leads to an amorphous film, from which QDs nucleate by SPE during annealing.

Appl Phys Lett