Riemannian optimization of photonic quantum circuits in phase and Fock space

Abstract

We propose a framework to design and optimize generic photonic quantum circuits composed of Gaussian objects (pure and mixed Gaussian states, Gaussian unitaries, Gaussian channels, Gaussian measurements) as well as non-Gaussian effects such as photon-number-resolving measurements. In this framework, we parametrize a phase space representation of Gaussian objects using elements of the symplectic group (or the unitary or orthogonal group in special cases), and then we transform it into the Fock representation using a single linear recurrence relation that computes the Fock amplitudes of any Gaussian object recursively. We also compute the gradient of the Fock amplitudes with respect to phase space parameters by differentiating through the recurrence relation. We can then use Riemannian optimization on the symplectic group to optimize M-mode Gaussian objects, avoiding the need to commit to particular realizations in terms of fundamental gates. This allows us to "mod out" all the different gate-level implementations of the same circuit, which now can be chosen after the optimization has completed. This can be especially useful when looking to answer general questions, such as bounding the value of a property over a class of states or transformations, or when one would like to worry about hardware constraints separately from the circuit optimization step. Finally, we make our framework extendable to non-Gaussian objects that can be written as linear combinations of Gaussian ones, by explicitly computing the change in global phase when states undergo Gaussian transformations. We implemented all of these methods in the freely available open-source library MrMustard, which we use in three examples to optimize the 216-mode interferometer in Borealis, and 2- and 3-modes circuits (with Fock measurements) to produce cat states and cubic phase states.

References in corpus (21)

• QuTiP 2: A Python framework for the dynamics of open quantum systems
• Graphical calculus for Gaussian pure states
• Quantum Computational Advantage via High-Dimensional Gaussian Boson Sampling
• The Boundary for Quantum Advantage in Gaussian Boson Sampling
• Perceval: A Software Platform for Discrete Variable Photonic Quantum Computing
• Franck-Condon factors by counting perfect matchings of graphs with loops
• The Complexity of Bipartite Gaussian Boson Sampling
• Further Compactifying Linear Optical Unitaries
• Fast simulation of bosonic qubits via Gaussian functions in phase space
• Training Gaussian Boson Sampling Distributions
• Matrix decompositions in Quantum Optics: Takagi/Autonne, Bloch-Messiah/Euler, Iwasawa, and Williamson
• Efficient backcasting search for optical quantum state synthesis
• Digital Discovery of 100 diverse Quantum Experiments with PyTheus
• Strong Simulation of Linear Optical Processes
• Threshold detection statistics of bosonic states
• Simulation of quantum optics by coherent state decomposition
• Automated design of quantum optical experiments for device-independent quantum key distribution
• Noise thresholds for classical simulability of non-linear Boson sampling
• Taming numerical errors in simulations of continuous variable non-Gaussian state preparation
• A Quadratic Speedup in the Optimization of Noisy Quantum Optical Circuits
• Linear cross-entropy certification of quantum computational advantage in Gaussian Boson Sampling

Cited by in corpus (11)

• Classification of dynamical Lie algebras for translation-invariant 2-local spin systems in one dimension
• Quantifying the imaginarity of quantum states via Tsallis relative entropy
• Classical simulation of non-Gaussian bosonic circuits
• Piquasso: A Photonic Quantum Computer Simulation Software Platform
• Photon-number moments and cumulants of Gaussian states
• Classical simulation and quantum resource theory of non-Gaussian optics
• Accurate Unsupervised Photon Counting from Transition Edge Sensor Signals
• The stellar decomposition of Gaussian quantum states
• What's my phase again? Computing the vacuum-to-vacuum amplitude of quadratic bosonic evolution
• Linear-optical protocols for mitigating and suppressing noise in bosonic systems
• The symplectic rank of non-Gaussian quantum states