Random samples of quantum channels have many applications in quantum information processing tasks. Due to the Choi-Jamio\l{}kowski isomorphism, there is a well-known correspondence between channels and states, and one can imagine adapting state sampling methods to sample quantum channels. Here, we discuss such an adaptation, using the Hamiltonian Monte Carlo method, a well-known classical method capable of producing high-quality samples from arbitrary, user-specified distributions. Its implementation requires an exact parametrization of the space of quantum channels, with no superfluous parameters and no constraints. We construct such a parametrization, and demonstrate its use in three common channel sampling applications.
High-quality random samples of quantum states are needed for a variety of tasks in quantum information and quantum computation. Searching the high-dimensional quantum state space for a global maximum of an objective function with many local maxima or evaluating an integral over a region in the quantum state space are but two exemplary applications of many. These tasks can only be performed reliably and efficiently with Monte Carlo methods, which involve good samplings of the parameter space in accordance with the relevant target distribution. We show how the Markov-chain Monte Carlo method known as Hamiltonian Monte Carlo, or hybrid Monte Carlo, can be adapted to this context. It is applicable when an efficient parameterization of the state space is available. The resulting random walk is entirely inside the physical parameter space, and the Hamiltonian dynamics enable us to take big steps, thereby avoiding strong correlations between successive sample points while enjoying a high acceptance rate. We use examples of single and double qubit measurements for illustration.
High-quality random samples of quantum states are needed for a variety of tasks in quantum information and quantum computation. Searching the high-dimensional quantum state space for a global maximum of an objective function with many local maxima or evaluating an integral over a region in the quantum state space are but two exemplary applications of many. These tasks can only be performed reliably and efficiently with Monte Carlo methods, which involve good samplings of the parameter space in accordance with the relevant target distribution. We show how the standard strategies of rejection sampling, importance sampling, and Markov-chain sampling can be adapted to this context, where the samples must obey the constraints imposed by the positivity of the statistical operator. For illustration, we generate sample points in the probability space of qubits, qutrits, and qubit pairs, both for tomographically complete and incomplete measurements. We use these samples for various purposes: establish the marginal distribution of the purity; compute the fractional volume of separable two-qubit states; and calculate the size of regions with bounded likelihood.
This is the documentation for generating random samples from the quantum state space in accordance with a specified distribution, associated with this webpage: http://tinyurl.com/QSampling . Ready-made samples (each with at least a million points) from various distributions are available for download, or one can generate one's own samples from a chosen distribution using the provided source codes. The sampling relies on the Hamiltonian Monte Carlo algorithm as described in New J. Phys. 17, 043018 (2015). The random samples are reposited in the hope that they would be useful for a variety of tasks in quantum information and quantum computation. Constructing credible regions for tomographic data, optimizing a function over the quantum state space with a complicated landscape, testing the typicality of entanglement among states from a multipartite quantum system, or computing the average of some quantity of interest over a subset of quantum states are but some exemplary applications among many.
Provided a maximally entangled state vector in
H⊗H, the Choi-Jamiolkowski isomorphism of quantum information theory (or channel/state duality) produces a mapping of
CP(H,K) (completely positive maps of density operators on
H tothose on
K) onto the density operators on
K⊗H, which extends to a unitary map. Here,
H and
K are {finite-dimensional} Hilbert spaces.The nature of this construction, the form of the associated formulas, and the possibility of versions which do not require an entangled vector seed, are illuminated from the perspective of {naturality}, in the sense of Eilenberg and MacLane.There is a unique such correspondence (``Choi map'') natural in both
H and
K which produces density operators on either
L(H,K) (Hilbert-Schmidt inner product) or
K⊗H∗ instead of
K⊗H. In particular, no entangled state input is required. There are no such binatural maps carried on
K⊗H, but those natural in only
K are canonically associated with unitary maps
H→H∗ and derive from the binatural maps. The construction recovers the standard Choi-Jamiolkowski formulae and shows how they express naturality.
Random samples of quantum states with specific properties are useful for various applications, such as Monte Carlo integration over the state space. In the high-dimensional situations that one already encounters when working with a few qubits, the quantum state space has a very complicated boundary, and it is challenging to incorporate the specific properties into the sampling algorithm. In this paper, we present the Sequentially Constrained Monte Carlo (SCMC) algorithm as a practical and versatile method for sampling quantum states in accordance with properties that can be stated as inequalities. We apply the SCMC algorithm to the generation of samples of bound entangled states; for example, we obtain nearly ten thousand bound, entangled, two-qutrit states in a few minutes, compared with less than ten such states per day from independence sampling in our implementation. In the second application, we draw samples of high-dimensional quantum states from a narrowly peaked target distribution and observe, for the system sizes investigated, that SCMC sampling remains computationally manageable as the dimensions grow. In yet another application, the SCMC algorithm produces uniformly distributed quantum states in regions bounded by values of the problem-specific target distribution; such samples are needed when estimating parameters from the probabilistic data acquired in quantum experiments.
One of the open challenges in quantum computing is to find meaningful and practical methods to leverage quantum computation to accelerate classical machine-learning workflows. A ubiquitous problem in machine-learning workflows is sampling from probability distributions that we only have access to via their log probability. To this end, we extend the well-known Hamiltonian Monte Carlo (HMC) method for Markov chain Monte Carlo (MCMC) sampling to leverage quantum computation in a hybrid manner as a proposal function. Our new algorithm, Quantum Dynamical Hamiltonian Monte Carlo (QD-HMC), replaces the classical symplectic integration proposal step with simulations of quantum-coherent continuous-space dynamics on digital or analog quantum computers. We show that QD-HMC maintains key characteristics of HMC, such as maintaining the detailed balanced condition with momentum inversion, while also having the potential for polynomial speedups over its classical counterpart in certain scenarios. As sampling is a core subroutine in many forms of probabilistic inference, and MCMC in continuously parametrized spaces covers a large class of potential applications, this work widens the areas of applicability of quantum devices. Published by the American Physical Society 2024