Source author record

Rolando D. Somma

Rolando D. Somma appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Tensor networks for High Energy Physics: contribution to Snowmass 2021

Tensor network methods are becoming increasingly important for high-energy physics, condensed matter physics and quantum information science (QIS). We discuss the impact of tensor network methods on lattice field theory, quantum gravity and QIS in the context of High Energy Physics (HEP). These tools will target calculations for strongly interacting systems that are made difficult by sign problems when conventional Monte Carlo and other importance sampling methods are used. Further development of methods and software will be needed to make a significant impact in HEP. We discuss the roadmap to perform quantum chromodynamics (QCD) related calculations in the coming years. The research is labor intensive and requires state of the art computational science and computer science input for its development and validation. We briefly discuss the overlap with other science domains and industry.

preprint2021arXiv

Complexity of quantum state verification in the quantum linear systems problem

We analyze the complexity of quantum state verification in the context of solving systems of linear equations of the form $A \vec x = \vec b$. We show that any quantum operation that verifies whether a given quantum state is within a constant distance from the solution of the quantum linear systems problem requires $q=Ω(κ)$ uses of a unitary that prepares a quantum state $\left| b \right>$, proportional to $\vec b$, and its inverse in the worst case. Here, $κ$ is the condition number of the matrix $A$. For typical instances, we show that $q=Ω(\sqrt κ)$ with high probability. These lower bounds are almost achieved if quantum state verification is performed using known quantum algorithms for the quantum linear systems problem. We also analyze the number of copies of $\left| b \right>$ required by verification procedures of the prepare and measure type. In this case, the lower bounds are quadratically worse, being $Ω(κ^2)$ in the worst case and $Ω(κ)$ in typical instances with high probability. We discuss the implications of our results to known variational and related approaches to this problem, where state preparation, gate, and measurement errors will need to decrease rapidly with $κ$ for worst-case and typical instances if error correction is not used, and present some open problems.

preprint2020arXiv

Operator Sampling for Shot-frugal Optimization in Variational Algorithms

Quantum chemistry is a near-term application for quantum computers. This application may be facilitated by variational quantum-classical algorithms (VQCAs), although a concern for VQCAs is the large number of measurements needed for convergence, especially for chemical accuracy. Here we introduce a strategy for reducing the number of measurements (i.e., shots) by randomly sampling operators $h_i$ from the overall Hamiltonian $H = \sum_i c_i h_i$. In particular, we employ weighted sampling, which is important when the $c_i$'s are highly non-uniform, as is typical in chemistry. We integrate this strategy with an adaptive optimizer developed recently by our group to construct an improved optimizer called Rosalin (Random Operator Sampling for Adaptive Learning with Individual Number of shots). Rosalin implements stochastic gradient descent while adapting the shot noise for each partial derivative and randomly assigning the shots amongst the $h_i$ according to a weighted distribution. We implement this and other optimizers to find the ground states of molecules H$_2$, LiH, and BeH$_2$, without and with quantum hardware noise, and Rosalin outperforms other optimizers in most cases.

preprint2020arXiv

Quantum eigenvalue estimation via time series analysis

We present an efficient method for estimating the eigenvalues of a Hamiltonian $H$ from the expectation values of the evolution operator for various times. For a given quantum state $ρ$, our method outputs a list of eigenvalue estimates and approximate probabilities. Each probability depends on the support of $ρ$ in those eigenstates of $H$ associated with eigenvalues within an arbitrarily small range. The complexity of our method is polynomial in the inverse of a given precision parameter $ε$, which is the gap between eigenvalue estimates. Unlike the well-known quantum phase estimation algorithm that uses the quantum Fourier transform, our method does not require large ancillary systems, large sequences of controlled operations, or preserving coherence between experiments, and is therefore more attractive for near-term applications. The output of our method can be used to compute spectral properties of $H$ and other expectation values efficiently, within additive error proportional to $ε$.

preprint2018arXiv

Quantum circuit synthesis for generalized coherent states

We present a method that outputs a sequence of simple unitary operations to prepare a given quantum state that is a generalized coherent state. Our method takes as inputs the expectation values of some relevant observables on the state to be prepared. Such expectation values can be estimated by performing projective measurements on $O(M^3 \log(M/δ)/ε^2)$ copies of the state, where $M$ is the dimension of an associated Lie algebra, $ε$ is a precision parameter, and $1-δ$ is the required confidence level. The method can be implemented on a classical computer and runs in time $O(M^4 \log(M/ε))$. It provides $O(M \log(M/ε))$ simple unitaries that form the sequence. The number of all computational resources is then polynomial in $M$, making the whole procedure very efficient in those cases where $M$ is significantly smaller than the Hilbert space dimension. When the algebra of relevant observables is determined by some Pauli matrices, each simple unitary may be easily decomposed into two-qubit gates. We discuss applications to quantum state tomography and classical simulations of quantum circuits.

preprint2016arXiv

Quantum simulations of one dimensional quantum systems

We present quantum algorithms for the simulation of quantum systems in one spatial dimension, which result in quantum speedups that range from superpolynomial to polynomial. We first describe a method to simulate the evolution of the quantum harmonic oscillator (QHO) based on a refined analysis of the Trotter-Suzuki formula that exploits the Lie algebra structure. For total evolution time $t$ and precision $ε>0$, the complexity of our method is $ O(\exp(γ\sqrt{\log(N/ε)}))$, where $γ>0$ is a constant and $N$ is the quantum number associated with an "energy cutoff" of the initial state. Remarkably, this complexity is subpolynomial in $N/ε$. We also provide a method to prepare discrete versions of the eigenstates of the QHO of complexity polynomial in $\log(N)/ε$, where $N$ is the dimension or number of points in the discretization. This method may be of independent interest as it provides a way to prepare, e.g., quantum states with Gaussian-like amplitudes. Next, we consider a system with a quartic potential. Our numerical simulations suggest a method for simulating the evolution of sublinear complexity $\tilde O(N^{1/3+o(1)})$, for constant $t$ and $ε$. We also analyze complex one-dimensional systems and prove a complexity bound $\tilde O(N)$, under fairly general assumptions. Our quantum algorithms may find applications in other problems. As an example, we discuss the fractional Fourier transform, a generalization of the Fourier transform that is useful for signal analysis and can be formulated in terms of the evolution of the QHO.

preprint2015arXiv

A Trotter-Suzuki approximation for Lie groups with applications to Hamiltonian simulation

We present a product formula to approximate the exponential of a skew-Hermitian operator that is a sum of generators of a Lie algebra. The number of terms in the product depends on the structure factors. When the generators have large norm with respect to the dimension of the Lie algebra, or when the norm of the effective operator resulting from nested commutators is less than the product of the norms, the number of terms in the product is significantly less than that obtained from well-known results. We apply our results to construct product formulas useful for the quantum simulation of some continuous-variable and bosonic physical systems, including systems whose potential is not quadratic. For many of these systems, we show that the number of terms in the product can be sublinear or subpolynomial in the dimension of the relevant local Hilbert spaces, where such a dimension is usually determined by an energy scale of the problem. Our results emphasize the power of quantum computers for the simulation of various quantum systems.

preprint2015arXiv

Quantum algorithms for simulated annealing

This paper summarizes a quantum algorithm of [R.D. Somma, et.al., Phys. Rev. Lett. 101, 130504 (2008)] that simulates a classical annealing process for solving discrete optimization problems. The complexity of the quantum algorithm scales with the inverse square root of the spectral gap of an associated stochastic matrix. This represents a quadratic quantum speedup, in terms of the gap, with respect to classical simulated annealing.

preprint2014arXiv

Exponential improvement in precision for simulating sparse Hamiltonians

We provide a quantum algorithm for simulating the dynamics of sparse Hamiltonians with complexity sublogarithmic in the inverse error, an exponential improvement over previous methods. Specifically, we show that a $d$-sparse Hamiltonian $H$ acting on $n$ qubits can be simulated for time $t$ with precision $ε$ using $O\big(τ\frac{\log(τ/ε)}{\log\log(τ/ε)}\big)$ queries and $O\big(τ\frac{\log^2(τ/ε)}{\log\log(τ/ε)}n\big)$ additional 2-qubit gates, where $τ= d^2 \|{H}\|_{\max} t$. Unlike previous approaches based on product formulas, the query complexity is independent of the number of qubits acted on, and for time-varying Hamiltonians, the gate complexity is logarithmic in the norm of the derivative of the Hamiltonian. Our algorithm is based on a significantly improved simulation of the continuous- and fractional-query models using discrete quantum queries, showing that the former models are not much more powerful than the discrete model even for very small error. We also simplify the analysis of this conversion, avoiding the need for a complex fault correction procedure. Our simplification relies on a new form of "oblivious amplitude amplification" that can be applied even though the reflection about the input state is unavailable. Finally, we prove new lower bounds showing that our algorithms are optimal as a function of the error.

preprint2014arXiv

Simulating Hamiltonian dynamics with a truncated Taylor series

We describe a simple, efficient method for simulating Hamiltonian dynamics on a quantum computer by approximating the truncated Taylor series of the evolution operator. Our method can simulate the time evolution of a wide variety of physical systems. As in another recent algorithm, the cost of our method depends only logarithmically on the inverse of the desired precision, which is optimal. However, we simplify the algorithm and its analysis by using a method for implementing linear combinations of unitary operations to directly apply the truncated Taylor series.

preprint2013arXiv

An exact real-space renormalization method and applications

We present a numerical method based on real-space renormalization that outputs the exact ground space of "frustration-free" Hamiltonians. The complexity of our method is polynomial in the degeneracy of the ground spaces of the Hamiltonians involved in the renormalization steps. We apply the method to obtain the full ground spaces of two spin systems. The first system is a spin-1/2 Heisenberg model with four-spin cyclic-exchange interactions defined on a square lattice. In this case, we study finite lattices of up to 160 spins and find a triplet ground state that differs from the singlet ground states obtained in C.D. Batista and S. Trugman, Phys. Rev. Lett. 93, 217202 (2004). We characterize such a triplet state as consisting of a triplon that propagates in a background of fluctuating singlet dimers. The second system is a family of spin-1/2 Heisenberg chains with uniaxial exchange anisotropy and next-nearest neighbor interactions. In this case, the method finds a ground-space degeneracy that scales quadratically with the system size and outputs the full ground space efficiently. Our method can substantially outperform methods based on exact diagonalization and is more efficient than other renormalization methods when the ground-space degeneracy is large.

preprint2013arXiv

Exponential improvement in precision for Hamiltonian-evolution simulation

We provide a quantum method for simulating Hamiltonian evolution with complexity polynomial in the logarithm of the inverse error. This is an exponential improvement over existing methods for Hamiltonian simulation. In addition, its scaling with respect to time is close to linear, and its scaling with respect to the time derivative of the Hamiltonian is logarithmic. These scalings improve upon most existing methods. Our method is to use a compressed Lie-Trotter formula, based on recent ideas for efficient discrete-time simulations of continuous-time quantum query algorithms.

preprint2013arXiv

Improved Bounds for Eigenpath Traversal

We present a bound on the length of the path defined by the ground states of a continuous family of Hamiltonians in terms of the spectral gap G. We use this bound to obtain a significant improvement over the cost of recently proposed methods for quantum adiabatic state transformations and eigenpath traversal. In particular, we prove that a method based on evolution randomization, which is a simple extension of adiabatic quantum computation, has an average cost of order 1/G^2, and a method based on fixed-point search, has a maximum cost of order 1/G^(3/2). Additionally, if the Hamiltonians satisfy a frustration-free property, such costs can be further improved to order 1/G^(3/2) and 1/G, respectively. Our methods offer an important advantage over adiabatic quantum computation when the gap is small, where the cost is of order 1/G^3.

preprint2013arXiv

Network-Centric Quantum Communications with Application to Critical Infrastructure Protection

Network-centric quantum communications (NQC) - a new, scalable instantiation of quantum cryptography providing key management with forward security for lightweight encryption, authentication and digital signatures in optical networks - is briefly described. Results from a multi-node experimental test-bed utilizing integrated photonics quantum communications components, known as QKarDs, include: quantum identification; verifiable quantum secret sharing; multi-party authenticated key establishment, including group keying; and single-fiber quantum-secured communications that can be applied as a security retrofit/upgrade to existing optical fiber installations. A demonstration that NQC meets the challenging simultaneous latency and security requirements of electric grid control communications, which cannot be met without compromises using conventional cryptography, is described.

preprint2013arXiv

Security of Decoy-State Protocols for General Photon-Number-Splitting Attacks

Decoy-state protocols provide a way to defeat photon-number splitting attacks in quantum cryptography implemented with weak coherent pulses. We point out that previous security analyses of such protocols relied on assumptions about eavesdropping attacks that considered treating each pulse equally and independently. We give an example to demonstrate that, without such assumptions, the security parameters of previous decoy-state implementations could be worse than the ones claimed. Next we consider more general photon-number splitting attacks, which correlate different pulses, and give an estimation procedure for the number of single photon signals with rigorous security statements. The impact of our result is that previous analyses of the number of times a decoy-state quantum cryptographic system can be reused before it makes a weak key must be revised.

preprint2012arXiv

Condensation of Anyons in Frustrated Quantum Magnets

We derive the exact ground space of a family of spin-1/2 Heisenberg chains with uniaxial exchange anisotropy (XXZ) and interactions between nearest and next-nearest-neighbor spins. The Hamiltonian family, H(Q), is parametrized by a single variable Q. By using a generalized Jordan-Wigner transformation that maps spins into anyons, we show that the exact ground states of H(Q) correspond to a condensation of anyons with statistical phase phi=-4Q. We also provide matrix-product state representations of some ground states that allow for the efficient computation of spin-spin correlation functions.

preprint2012arXiv

Quantum Speedup by Quantum Annealing

We study the glued-trees problem of Childs et. al. in the adiabatic model of quantum computing and provide an annealing schedule to solve an oracular problem exponentially faster than classically possible. The Hamiltonians involved in the quantum annealing do not suffer from the so-called sign problem. Unlike the typical scenario, our schedule is efficient even though the minimum energy gap of the Hamiltonians is exponentially small in the problem size. We discuss generalizations based on initial-state randomization to avoid some slowdowns in adiabatic quantum computing due to small gaps.

preprint2012arXiv

Spectral Gap Amplification

A large number of problems in science can be solved by preparing a specific eigenstate of some Hamiltonian H. The generic cost of quantum algorithms for these problems is determined by the inverse spectral gap of H for that eigenstate and the cost of evolving with H for some fixed time. The goal of spectral gap amplification is to construct a Hamiltonian H' with the same eigenstate as H but a bigger spectral gap, requiring that constant-time evolutions with H' and H are implemented with nearly the same cost. We show that a quadratic spectral gap amplification is possible when H satisfies a frustration-free property and give H' for these cases. This results in quantum speedups for optimization problems. It also yields improved constructions for adiabatic simulations of quantum circuits and for the preparation of projected entangled pair states (PEPS), which play an important role in quantum many-body physics. Defining a suitable black-box model, we establish that the quadratic amplification is optimal for frustration-free Hamiltonians and that no spectral gap amplification is possible, in general, if the frustration-free property is removed. A corollary is that finding a similarity transformation between a stoquastic Hamiltonian and the corresponding stochastic matrix is hard in the black-box model, setting limits to the power of some classical methods that simulate quantum adiabatic evolutions.