We introduce a method for efficient, in situ characterization of linear-optical networks (LONs) in randomized boson-sampling (RBS) experiments. We formulate RBS as a distributed task between two parties, Alice and Bob, who share two-mode squeezed-vacuum states. In this protocol, Alice performs local measurements on her modes, either photon counting or heterodyne. Bob implements and applies to his modes the LON requested by Alice; at the output of the LON, Bob performs photon counting, the results of which he sends to Alice via classical channels. In the ideal situation, when Alice does photon counting, she obtains from Bob samples from the probability distribution of the RBS problem, a task that is believed to be classically hard to simulate. When Alice performs heterodyne measurements, she converts the experiment to a problem that is classically efficiently simulable, but more importantly enables her to characterize a lossy LON on the fly, without Bob's knowledge and without changing anything at Bob's end (this is what we mean by in situ). We introduce and calculate the fidelity between the joint states shared by Alice and Bob after the ideal and lossy LONs as a measure of distance between the two LONs. Using this measure, we obtain an upper bound on the total variation distance between the ideal probability distribution for the RBS problem and the probability distribution achieved by a lossy LON. Our method displays the power of the entanglement of the two-mode squeezed-vacuum states: The entanglement allows Alice to choose for each run of the experiment between RBS and a simple characterization protocol based on first-order coherence between complex amplitudes.
We introduce a method for efficient, in situ characterization of linear-optical networks (LONs) in randomized boson-sampling (RBS) experiments. We formulate RBS as a distributed task between two parties, Alice and Bob, who share two-mode squeezed-vacuum states. In this protocol, Alice performs local measurements on her modes, either photon counting or heterodyne. Bob implements and applies to his modes the LON requested by Alice; at the output of the LON, Bob performs photon counting, the results of which he sends to Alice via classical channels. In the ideal situation, when Alice does photon counting, she obtains from Bob samples from the probability distribution of the RBS problem, a task that is believed to be classically hard to simulate. When Alice performs heterodyne measurements, she converts the experiment to a problem that is classically efficiently simulable, but more importantly enables her to characterize a lossy LON on the fly, without Bob's knowledge and without changing anything at Bob's end (this is what we mean by in situ). We introduce and calculate the fidelity between the joint states shared by Alice and Bob after the ideal and lossy LONs as a measure of distance between the two LONs. Using this measure, we obtain an upper bound on the total variation distance between the ideal probability distribution for the RBS problem and the probability distribution achieved by a lossy LON. Our method displays the power of the entanglement of the two-mode squeezed-vacuum states: The entanglement allows Alice to choose for each run of the experiment between RBS and a simple characterization protocol based on first-order coherence between complex amplitudes.
This thesis presents results under the topics of quantum correlations and quantum information processing using linear-optical networks. The first part of this thesis concerns with quantum correlations, specially in Gaussian states, which play a central role as resource in quantum protocols. The second part is about quantum information processing using linear-optical networks that are of great practical interest due to their simple physical realization. In the first part, we present a measurement-based method for verifying nonzero quantum discord, a measure of quantum correlations, in bipartite systems that can be used for both continuous and discrete-variable systems. Based on this method, we show that the only bipartite Gaussian states with zero quantum discord are product states. We also present a simple and efficient experimental method for verifying quantum correlations in Gaussian states using homodyne measurements. Moreover, we introduce an operational measure for quantifying quantum correlations in Gaussian states that is purely based on Gaussian measurements. We then illustrate the operational significance of the measure in terms of a Gaussian quantum protocol. In the second part, we present results on linear quantum-optics experiments, specifically boson sampling, an intermediate model of quantum computation that can be implemented using single-photon states, linear-optical networks, and photodetectors. We first present an efficient experimental method for characterization of linear-optical networks, which is essential for practical implementation of boson sampling. We show that using coherent states one can directly measure all the nontrivial phases and moduli of the matrix describing the network. Next, we consider the problem of boson sampling with Gaussian states. We derive a general formula for calculating the output photon-counting probabilities. We show that for input thermal states, sampling from the output probability distribution can be efficiently classically simulated. Using this, we find that permanents of positive-semidefinite Hermitian matrices can be approximated to within a multiplicative error in BPPNP. We also discuss boson sampling with squeezed vacuum states. Finally, we provide a sufficient condition for the efficient classical simulation of the general sampling problem using linear-quantum optics, which involves a multimode input quantum state to a linear-optical network and set of measurements at the output. This condition is based on the phase-space quasiprobability distributions, and suggests that negativity is a quantum resource. We apply this condition to implementations of boson sampling using single-photon or spontaneous parametric down-conversion sources. We show that above some threshold for loss and noise in the experiment, boson sampling is classically simulatable.
Boson-Sampling is a classically computationally hard problem that can - in principle - be efficiently solved with quantum linear optical networks. Very recently, a rush of experimental activity has ignited with the aim of developing such devices as feasible instances of quantum simulators. Even approximate Boson-Sampling is believed to be hard with high probability if the unitary describing the optical network is drawn from the Haar measure. In this work we show that in this setup, with probability exponentially close to one in the number of bosons, no symmetric algorithm can distinguish the Boson-Sampling distribution from the uniform one from fewer than exponentially many samples. This means that the two distributions are operationally indistinguishable without detailed a priori knowledge. We carefully discuss the prospects of efficiently using knowledge about the implemented unitary for devising non-symmetric algorithms that could potentially improve upon this. We conclude that due to the very fact that Boson-Sampling is believed to be hard, efficient classical certification of Boson-Sampling devices seems to be out of reach.
BosonSampling, which we proposed three years ago, is a scheme for using linear-optical networks to solve sampling problems that appear to be intractable for a classical computer. \ In a recent manuscript, Gogolin et al.\ claimed that even an ideal BosonSampling device's output would be operationally indistinguishable\textquotedblright\ from a uniform random outcome, at least \textquotedblleft without detailed a priori knowledge; or at any rate, that telling the two apart might itself be a hard problem. We first answer these claims---explaining why the first is based on a definition of a priori knowledge such that, were it adopted, almost no quantum algorithm could be distinguished from a pure random-number source; while the second is neither new nor a practical obstacle to interesting BosonSampling experiments.However, we then go further, and address some interesting research questions inspired by Gogolin et al.'s arguments. We prove that, with high probability over a Haar-random matrix A, the BosonSampling distribution induced by A is far from the uniform distribution in total variation distance. More surprisingly, and counter to Gogolin et al., we give an efficient algorithm that distinguishes these two distributions with constant bias. Finally, we offer three bonus results about BosonSampling. First, we report an observation of Fernando Brandao: that one can efficiently sample a distribution that has large entropy and that's indistinguishable from a BosonSampling distribution by any circuit of fixed polynomial size. Second, we show that BosonSampling distributions can be efficiently distinguished from uniform even with photon losses and for general initial states. Third, we offer the simplest known proof that Fermion Sampling is solvable in classical polynomial time, and we reuse techniques from our Boson Sampling analysis to characterize random FermionSampling distributions.
Sampling from probability distributions of quantum circuits is a fundamentally and practically important task which can be used to demonstrate quantum supremacy using noisy intermediate-scale quantum devices. In the present work, we examine classical simulability of sampling from the output photon-number distribution of linear-optical circuits composed of random beam splitters with equally distributed squeezed vacuum states and single-photon states input. We provide efficient classical algorithms to simulate linear-optical random circuits and show that the algorithms' error is exponentially small up to a depth less than quadratic in the distance between sources using a classical random walk behavior of random linear-optical circuits. Notably, the average-case depth allowing an efficient classical simulation is larger than the worst-case depth limit, which is linear in the distance. Besides, our results together with the hardness of boson sampling give a lower-bound on the depth for constituting global Haar-random unitary circuits.
Boson sampling is a candidate for quantum protocols to truly realize the quantum computation advantage and to be used in advanced fields where complex computations are needed, such as quantum chemistry. However, this proposal is hard to achieve due to the existence of noise sources such as photon losses. In order to quantificationally analyze the influences of photon losses in optical networks, boson sampling is classically simulated based on the equivalent beam splitter mechanism, where the photon loss happening in optical units is equivalent to the photon transmission into the environmental paths through a virtual beam splitter. In our simulation, networks corresponding to random unitary matrices are made up, considering both the Reck structure and the Clements structure. The photon loss probability in an optical unit is well controlled by adjusting the parameters of the virtual beam splitter. Therefore, to simulate boson sampling with photon losses in optical networks is actually to simulate ideal boson sampling with more modes. It is found that when the photon loss probability is constant, boson sampling with Clements structures distinctly performs much better than that with Reck structures. Furthermore, the photon loss probability is also set to follow the normal distribution, which is thought to be closer to the situation in reality. It is found that when the mean value of photon loss probability is constant, for both network structures, errors of outputs become more obvious with the increase of standard deviation. It can be inferred that the increase of error rate can be explained by the network depth and the conclusion is suitable for larger-scale boson sampling. Finally, the number of output photons is taken into consideration, which is directly related to the classical computation complexity. It is found that with the photon loss probability, the ratio of output combinations without photon losses decreases sharply, implying that photon losses can obviously affect the quantum computation advantage of boson sampling. Our results indicate that photon losses can result in serious errors for boson sampling, even with a stable network structure such as that of Clements. This work is helpful for boson sampling experiments in reality and it is desired to develop a better protocol, for example, a well-designed network or excellent optical units, to well suppress photon losses.
BosonSampling, which we proposed three years ago, is a scheme for using linear-optical networks to solve sampling problems that appear to be intractable for a classical computer. In arXiv:1306.3995, Gogolin et al. claimed that even an ideal BosonSampling device's output would be "operationally indistinguishable" from a uniform random outcome, at least "without detailed a priori knowledge"; or at any rate, that telling the two apart might itself be a hard problem. We first answer these claims---explaining why the first is based on a definition of "a priori knowledge" so strange that, were it adopted, almost no quantum algorithm could be distinguished from a pure random-number source; while the second is neither new nor a practical obstacle to interesting BosonSampling experiments. However, we then go further, and address some interesting research questions inspired by Gogolin et al.'s mistaken arguments. We prove that, with high probability over a Haar-random matrix A, the BosonSampling distribution induced by A is far from the uniform distribution in total variation distance. More surprisingly, and directly counter to Gogolin et al., we give an efficient algorithm that distinguishes these two distributions with constant bias. Finally, we offer three "bonus" results about BosonSampling. First, we report an observation of Fernando Brandao: that one can efficiently sample a distribution that has large entropy and that's indistinguishable from a BosonSampling distribution by any circuit of fixed polynomial size. Second, we show that BosonSampling distributions can be efficiently distinguished from uniform even with photon losses and for general initial states. Third, we offer the simplest known proof that FermionSampling is solvable in classical polynomial time, and we reuse techniques from our BosonSampling analysis to characterize random FermionSampling distributions.