{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,30]],"date-time":"2026-04-30T02:41:01Z","timestamp":1777516861907,"version":"3.51.4"},"reference-count":31,"publisher":"Verein zur Forderung des Open Access Publizierens in den Quantenwissenschaften","license":[{"start":{"date-parts":[[2023,8,29]],"date-time":"2023-08-29T00:00:00Z","timestamp":1693267200000},"content-version":"unspecified","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"The Research Foundation \u2013 Flanders","award":["G006020N"],"award-info":[{"award-number":["G006020N"]}]},{"name":"Excellence of Science (EOS) research programme","award":["G0H1422N"],"award-info":[{"award-number":["G0H1422N"]}]}],"content-domain":{"domain":["quantum-journal.org"],"crossmark-restriction":false},"short-container-title":["Quantum"],"abstract":"<jats:p>Linear optical quantum circuits with photon number resolving (PNR) detectors are used for both Gaussian Boson Sampling (GBS) and for the preparation of non-Gaussian states such as Gottesman-Kitaev-Preskill (GKP), cat and NOON states. They are crucial in many schemes of quantum computing and quantum metrology. Classically optimizing circuits with PNR detectors is challenging due to their exponentially large Hilbert space, and quadratically more challenging in the presence of decoherence as state vectors are replaced by density matrices. To tackle this problem, we introduce a family of algorithms that calculate detection probabilities, conditional states (as well as their gradients with respect to circuit parametrizations) with a complexity that is comparable to the noiseless case. As a consequence we can simulate and optimize circuits with twice the number of modes as we could before, using the same resources. More precisely, for an <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>M<\/mml:mi><\/mml:math>-mode noisy circuit with detected modes <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>D<\/mml:mi><\/mml:math> and undetected modes <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>U<\/mml:mi><\/mml:math>, the complexity of our algorithm is <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>O<\/mml:mi><mml:mo stretchy=\"false\">(<\/mml:mo><mml:msup><mml:mi>M<\/mml:mi><mml:mn>2<\/mml:mn><\/mml:msup><mml:munder><mml:mo>&amp;#x220F;<\/mml:mo><mml:mrow class=\"MJX-TeXAtom-ORD\"><mml:mi>i<\/mml:mi><mml:mspace width=\"0.111em\"\/><mml:mo>&amp;#x2208;<\/mml:mo><mml:mspace width=\"0.111em\"\/><mml:mi>U<\/mml:mi><\/mml:mrow><\/mml:munder><mml:msubsup><mml:mi>C<\/mml:mi><mml:mi>i<\/mml:mi><mml:mn>2<\/mml:mn><\/mml:msubsup><mml:munder><mml:mo>&amp;#x220F;<\/mml:mo><mml:mrow class=\"MJX-TeXAtom-ORD\"><mml:mi>i<\/mml:mi><mml:mspace width=\"0.111em\"\/><mml:mo>&amp;#x2208;<\/mml:mo><mml:mspace width=\"0.111em\"\/><mml:mi>D<\/mml:mi><\/mml:mrow><\/mml:munder><mml:msub><mml:mi>C<\/mml:mi><mml:mi>i<\/mml:mi><\/mml:msub><mml:mo stretchy=\"false\">)<\/mml:mo><\/mml:math>, rather than <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>O<\/mml:mi><mml:mo stretchy=\"false\">(<\/mml:mo><mml:msup><mml:mi>M<\/mml:mi><mml:mn>2<\/mml:mn><\/mml:msup><mml:munder><mml:mo>&amp;#x220F;<\/mml:mo><mml:mrow class=\"MJX-TeXAtom-ORD\"><mml:mspace width=\"0.111em\"\/><mml:mi>i<\/mml:mi><mml:mspace width=\"0.111em\"\/><mml:mo>&amp;#x2208;<\/mml:mo><mml:mspace width=\"0.111em\"\/><mml:mi>D<\/mml:mi><mml:mspace width=\"0.167em\"\/><mml:mo>&amp;#x222A;<\/mml:mo><mml:mspace width=\"0.167em\"\/><mml:mi>U<\/mml:mi><\/mml:mrow><\/mml:munder><mml:msubsup><mml:mi>C<\/mml:mi><mml:mi>i<\/mml:mi><mml:mn>2<\/mml:mn><\/mml:msubsup><mml:mo stretchy=\"false\">)<\/mml:mo><\/mml:math>, where <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:msub><mml:mi>C<\/mml:mi><mml:mi>i<\/mml:mi><\/mml:msub><\/mml:math> is the Fock cutoff of mode <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\"><mml:mi>i<\/mml:mi><\/mml:math>. As a particular case, our approach offers a full quadratic speedup for calculating detection probabilities, as in that case all modes are detected. Finally, these algorithms are implemented and ready to use in the open-source photonic optimization library MrMustard.<\/jats:p>","DOI":"10.22331\/q-2023-08-29-1097","type":"journal-article","created":{"date-parts":[[2023,8,29]],"date-time":"2023-08-29T15:19:47Z","timestamp":1693322387000},"page":"1097","update-policy":"https:\/\/doi.org\/10.22331\/q-crossmark-policy-page","source":"Crossref","is-referenced-by-count":6,"title":["A Quadratic Speedup in the Optimization of Noisy Quantum Optical Circuits"],"prefix":"10.22331","volume":"7","author":[{"given":"Robbe","family":"De Prins","sequence":"first","affiliation":[{"name":"Photonics Research Group, INTEC, Ghent University \u2013 imec, Sint-Pietersnieuwstraat 41, 9000 Ghent, Belgium"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yuan","family":"Yao","sequence":"additional","affiliation":[{"name":"T\u00e9l\u00e9com Paris and Institut Polytechnique de Paris, LTCI, 20 Place Marguerite Perey, 91120 Palaiseau, France"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Anuj","family":"Apte","sequence":"additional","affiliation":[{"name":"Xanadu, Toronto, ON, M5G 2C8, Canada"},{"name":"Kadanoff Center for Theoretical Physics & Enrico Fermi Institute, Department of Physics, University of Chicago, Chicago, IL 60637"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Filippo M.","family":"Miatto","sequence":"additional","affiliation":[{"name":"T\u00e9l\u00e9com Paris and Institut Polytechnique de Paris, LTCI, 20 Place Marguerite Perey, 91120 Palaiseau, France"},{"name":"Xanadu, Toronto, ON, M5G 2C8, Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"9598","published-online":{"date-parts":[[2023,8,29]]},"reference":[{"key":"0","doi-asserted-by":"publisher","unstructured":"Juan Miguel Arrazola and Thomas R. Bromley. Using Gaussian boson sampling to find dense subgraphs. Physical Review Letters, 121 (3), July 2018. 10.1103\/physrevlett.121.030503.","DOI":"10.1103\/physrevlett.121.030503"},{"key":"1","doi-asserted-by":"publisher","unstructured":"Juan Miguel Arrazola, Thomas R. Bromley, and Patrick Rebentrost. Quantum approximate optimization with Gaussian boson sampling. Physical Review A, 98 (1), July 2018. 10.1103\/physreva.98.012322.","DOI":"10.1103\/physreva.98.012322"},{"key":"2","doi-asserted-by":"publisher","unstructured":"Leonardo Banchi, Mark Fingerhuth, Tomas Babej, Christopher Ing, and Juan Miguel Arrazola. Molecular docking with Gaussian boson sampling. Science Advances, 6 (23), June 2020a. 10.1126\/sciadv.aax1950.","DOI":"10.1126\/sciadv.aax1950"},{"key":"3","doi-asserted-by":"publisher","unstructured":"Leonardo Banchi, Nicol\u00e1s Quesada, and Juan Miguel Arrazola. Training Gaussian boson sampling distributions. Physical Review A, 102 (1): 012417, 2020b. 10.1103\/PhysRevA.102.012417.","DOI":"10.1103\/PhysRevA.102.012417"},{"key":"4","doi-asserted-by":"publisher","unstructured":"J. Eli Bourassa, Rafael N. Alexander, Michael Vasmer, Ashlesha Patil, Ilan Tzitrin, Takaya Matsuura, Daiqin Su, Ben Q. Baragiola, Saikat Guha, Guillaume Dauphinais, et al. Blueprint for a scalable photonic fault-tolerant quantum computer. Quantum, 5: 392, 2021. 10.22331\/q-2021-02-04-392.","DOI":"10.22331\/q-2021-02-04-392"},{"key":"5","doi-asserted-by":"publisher","unstructured":"Kamil Br\u00e1dler, Pierre-Luc Dallaire-Demers, Patrick Rebentrost, Daiqin Su, and Christian Weedbrook. Gaussian boson sampling for perfect matchings of arbitrary graphs. Physical Review A, 98 (3), September 2018. 10.1103\/physreva.98.032310.","DOI":"10.1103\/physreva.98.032310"},{"key":"6","doi-asserted-by":"publisher","unstructured":"Kamil Br\u00e1dler, Shmuel Friedland, Josh Izaac, Nathan Killoran, and Daiqin Su. Graph isomorphism and Gaussian boson sampling. Special Matrices, 9 (1): 166\u2013196, January 2021. 10.1515\/spma-2020-0132.","DOI":"10.1515\/spma-2020-0132"},{"key":"7","doi-asserted-by":"publisher","unstructured":"Thomas R. Bromley, Juan Miguel Arrazola, Soran Jahangiri, Josh Izaac, Nicol\u00e1s Quesada, Alain D. Gran, Maria Schuld, Jeremy Swinarton, Zeid Zabaneh, and Nathan Killoran. Applications of near-term photonic quantum computers: software and algorithms. Quantum Science and Technology, 5 (3): 034010, 2020. 10.1088\/2058-9565\/ab8504.","DOI":"10.1088\/2058-9565\/ab8504"},{"key":"8","doi-asserted-by":"publisher","unstructured":"Jacob F. F. Bulmer, Bryn A. Bell, Rachel S. Chadwick, Alex E. Jones, Diana Moise, Alessandro Rigazzi, Jan Thorbecke, Utz-Uwe Haus, Thomas Van Vaerenbergh, Raj B. Patel, et al. The boundary for quantum advantage in Gaussian boson sampling. Science advances, 8 (4): eabl9236, 2022. 10.1126\/sciadv.abl9236.","DOI":"10.1126\/sciadv.abl9236"},{"key":"9","doi-asserted-by":"publisher","unstructured":"Kevin E. Cahill and Roy J. Glauber. Density operators and quasiprobability distributions. Physical Review, 177 (5): 1882, 1969. 10.1103\/PhysRev.177.1882.","DOI":"10.1103\/PhysRev.177.1882"},{"key":"10","doi-asserted-by":"publisher","unstructured":"Kosuke Fukui, Shuntaro Takeda, Mamoru Endo, Warit Asavanant, Jun-ichi Yoshikawa, Peter van Loock, and Akira Furusawa. Efficient backcasting search for optical quantum state synthesis. Phys. Rev. Lett., 128: 240503, June 2022. 10.1103\/PhysRevLett.128.240503.","DOI":"10.1103\/PhysRevLett.128.240503"},{"key":"11","unstructured":"Christopher C. Gerry and Peter L. Knight. Introductory quantum optics. Cambridge university press, 2005."},{"key":"12","doi-asserted-by":"publisher","unstructured":"Daniel Gottesman, Alexei Kitaev, and John Preskill. Encoding a qubit in an oscillator. Phys. Rev. A, 64: 012310, June 2001. 10.1103\/PhysRevA.64.012310.","DOI":"10.1103\/PhysRevA.64.012310"},{"key":"13","doi-asserted-by":"publisher","unstructured":"Craig S. Hamilton, Regina Kruse, Linda Sansoni, Sonja Barkhofen, Christine Silberhorn, and Igor Jex. Gaussian boson sampling. Phys. Rev. Lett., 119: 170501, October 2017. 10.1103\/PhysRevLett.119.170501.","DOI":"10.1103\/PhysRevLett.119.170501"},{"key":"14","doi-asserted-by":"publisher","unstructured":"Joonsuk Huh and Man-Hong Yung. Vibronic boson sampling: Generalized Gaussian boson sampling for molecular vibronic spectra at finite temperature. Scientific Reports, 7 (1), August 2017. 10.1038\/s41598-017-07770-z.","DOI":"10.1038\/s41598-017-07770-z"},{"key":"15","doi-asserted-by":"publisher","unstructured":"Soran Jahangiri, Juan Miguel Arrazola, Nicol\u00e1s Quesada, and Nathan Killoran. Point processes with Gaussian boson sampling. Physical Review E, 101 (2), February 2020. 10.1103\/physreve.101.022134.","DOI":"10.1103\/physreve.101.022134"},{"key":"16","doi-asserted-by":"publisher","unstructured":"Regina Kruse, Craig S. Hamilton, Linda Sansoni, Sonja Barkhofen, Christine Silberhorn, and Igor Jex. Detailed study of Gaussian boson sampling. Phys. Rev. A, 100: 032326, September 2019. 10.1103\/PhysRevA.100.032326.","DOI":"10.1103\/PhysRevA.100.032326"},{"key":"17","doi-asserted-by":"publisher","unstructured":"Filippo M. Miatto and Nicol\u00e1s Quesada. Fast optimization of parametrized quantum optical circuits. Quantum, 4: 366, 2020. 10.22331\/q-2020-11-30-366.","DOI":"10.22331\/q-2020-11-30-366"},{"key":"18","doi-asserted-by":"publisher","unstructured":"Changhun Oh, Minzhao Liu, Yuri Alexeev, Bill Fefferman, and Liang Jiang. Tensor network algorithm for simulating experimental Gaussian boson sampling. arXiv preprint arXiv:2306.03709, 2023. 10.48550\/arXiv.2306.03709.","DOI":"10.48550\/arXiv.2306.03709"},{"key":"19","doi-asserted-by":"publisher","unstructured":"Nicol\u00e1s Quesada. Franck-Condon factors by counting perfect matchings of graphs with loops. The Journal of chemical physics, 150 (16): 164113, 2019. 10.1063\/1.5086387.","DOI":"10.1063\/1.5086387"},{"key":"20","doi-asserted-by":"publisher","unstructured":"Nicol\u00e1s Quesada, Luke G. Helt, Josh Izaac, Juan Miguel Arrazola, Reihaneh Shahrokhshahi, Casey R. Myers, and Krishna K. Sabapathy. Simulating realistic non-Gaussian state preparation. Phys. Rev. A, 100: 022341, August 2019. 10.1103\/PhysRevA.100.022341.","DOI":"10.1103\/PhysRevA.100.022341"},{"key":"21","doi-asserted-by":"publisher","unstructured":"Krishna K. Sabapathy, Haoyu Qi, Josh Izaac, and Christian Weedbrook. Production of photonic universal quantum gates enhanced by machine learning. Phys. Rev. A, 100: 012326, July 2019. 10.1103\/PhysRevA.100.012326.","DOI":"10.1103\/PhysRevA.100.012326"},{"key":"22","doi-asserted-by":"publisher","unstructured":"Maria Schuld, Ville Bergholm, Christian Gogolin, Josh Izaac, and Nathan Killoran. Evaluating analytic gradients on quantum hardware. Phys. Rev. A, 99 (3): 032331, 2019. 10.1103\/PhysRevA.99.032331.","DOI":"10.1103\/PhysRevA.99.032331"},{"key":"23","doi-asserted-by":"publisher","unstructured":"Maria Schuld, Kamil Br\u00e1dler, Robert Israel, Daiqin Su, and Brajesh Gupt. Measuring the similarity of graphs with a Gaussian boson sampler. Physical Review A, 101 (3), March 2020. 10.1103\/physreva.101.032314.","DOI":"10.1103\/physreva.101.032314"},{"key":"24","doi-asserted-by":"publisher","unstructured":"Daiqin Su, Casey R. Myers, and Krishna K. Sabapathy. Conversion of Gaussian states to non-Gaussian states using photon-number-resolving detectors. Phys. Rev. A, 100: 052301, November 2019a. 10.1103\/PhysRevA.100.052301.","DOI":"10.1103\/PhysRevA.100.052301"},{"key":"25","doi-asserted-by":"publisher","unstructured":"Daiqin Su, Casey R. Myers, and Krishna K. Sabapathy. Generation of photonic non-Gaussian states by measuring multimode Gaussian states. arXiv preprint arXiv:1902.02331, 2019b. 10.48550\/arXiv.1902.02331.","DOI":"10.48550\/arXiv.1902.02331"},{"key":"26","doi-asserted-by":"publisher","unstructured":"Kan Takase, Jun-ichi Yoshikawa, Warit Asavanant, Mamoru Endo, and Akira Furusawa. Generation of optical Schr\u00f6dinger cat states by generalized photon subtraction. Phys. Rev. A, 103: 013710, January 2021. 10.1103\/PhysRevA.103.013710.","DOI":"10.1103\/PhysRevA.103.013710"},{"key":"27","doi-asserted-by":"publisher","unstructured":"Kan Takase, Kosuke Fukui, Akito Kawasaki, Warit Asavanant, Mamoru Endo, Jun-ichi Yoshikawa, Peter van Loock, and Akira Furusawa. Gaussian breeding for encoding a qubit in propagating light. arXiv preprint arXiv:2212.05436, 2022. 10.48550\/arXiv.2212.05436.","DOI":"10.48550\/arXiv.2212.05436"},{"key":"28","unstructured":"Xanadu Quantum Technologies. MrMustard. https:\/\/github.com\/XanaduAI\/MrMustard, 2022."},{"key":"29","doi-asserted-by":"publisher","unstructured":"Ilan Tzitrin, J. Eli Bourassa, Nicolas C. Menicucci, and Krishna K. Sabapathy. Progress towards practical qubit computation using approximate Gottesman-Kitaev-Preskill codes. Phys. Rev. A, 101: 032315, March 2020. 10.1103\/PhysRevA.101.032315.","DOI":"10.1103\/PhysRevA.101.032315"},{"key":"30","doi-asserted-by":"publisher","unstructured":"Yuan Yao, Filippo M. Miatto, and Nicol\u00e1s Quesada. The recursive representation of Gaussian quantum mechanics. arXiv preprint arXiv:2209.06069, 2022. 10.48550\/arXiv.2209.06069.","DOI":"10.48550\/arXiv.2209.06069"}],"container-title":["Quantum"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/quantum-journal.org\/papers\/q-2023-08-29-1097\/pdf\/","content-type":"unspecified","content-version":"vor","intended-application":"text-mining"}],"deposited":{"date-parts":[[2023,9,27]],"date-time":"2023-09-27T16:12:19Z","timestamp":1695831139000},"score":1,"resource":{"primary":{"URL":"https:\/\/quantum-journal.org\/papers\/q-2023-08-29-1097\/"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,8,29]]},"references-count":31,"URL":"https:\/\/doi.org\/10.22331\/q-2023-08-29-1097","archive":["CLOCKSS"],"relation":{},"ISSN":["2521-327X"],"issn-type":[{"value":"2521-327X","type":"electronic"}],"subject":[],"published":{"date-parts":[[2023,8,29]]},"article-number":"1097"}}