{"id":2513,"date":"2018-07-19T15:00:42","date_gmt":"2018-07-19T14:00:42","guid":{"rendered":"https:\/\/www.quantum-bits.org\/?p=2513"},"modified":"2022-08-12T17:15:11","modified_gmt":"2022-08-12T16:15:11","slug":"on-photonic-quantum-computing","status":"publish","type":"post","link":"https:\/\/www.quantum-bits.org\/?p=2513","title":{"rendered":"On photonic quantum computing"},"content":{"rendered":"<p>The worldwide quest to build practical quantum computers is undergoing a critical period. Fault-tolerant quantum computers will soon provide significant computational speedups for problems like factoring, search, or linear algebra, etc. During the next few years, a number of different quantum devices will become available to the public.&nbsp;<\/p>\n<p>Over the past decade, a wide variety of physical architectures have been investigated for their suitability for quantum technologies, including (among others) <a href=\"https:\/\/en.wikipedia.org\/wiki\/Trapped_ion_quantum_computer\" target=\"_blank\" rel=\"noopener\">trapped atoms and ions<\/a> quantum computers, <a href=\"https:\/\/www.quantum-bits.org\/?p=2226\" target=\"_blank\" rel=\"noopener\">topological quantum computers<\/a>, or photonic quantum computers.<\/p>\n<p>In 2000 by<a href=\"https:\/\/arxiv.org\/abs\/quant-ph\/0006088\" target=\"_blank\" rel=\"noopener\"> E. Knill, R. Laflamme and G. Milburn<\/a> proposed a protocol (now named <strong>KLM<\/strong> scheme) using <a href=\"https:\/\/en.wikipedia.org\/wiki\/Photon\" target=\"_blank\" rel=\"noopener\">photons&nbsp;<\/a>as information carriers to implement linear optical quantum computing. This protocol makes it possible to create <a href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_Turing_machine\" target=\"_blank\" rel=\"noopener\">universal&nbsp;quantum computers<\/a>&nbsp;solely with&nbsp;linear optical&nbsp;tools.<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-1612\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/03\/quantum-physics-formulas-over-blackboard.jpg\" alt=\"\" width=\"850\" height=\"332\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/03\/quantum-physics-formulas-over-blackboard.jpg 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/03\/quantum-physics-formulas-over-blackboard-300x117.jpg 300w\" sizes=\"(max-width: 850px) 100vw, 850px\" \/><\/p>\n<p><strong>Photons<\/strong><\/p>\n<p>Let&#8217;s start first with a general recap on photons, which are <a title=\"Elementary particle\" href=\"https:\/\/en.wikipedia.org\/wiki\/Elementary_particle\" target=\"_blank\" rel=\"noopener\">elementary<\/a> spin-1 particles. By virtue of their integer spin, they are <a title=\"Boson\" href=\"https:\/\/en.wikipedia.org\/wiki\/Boson\" target=\"_blank\" rel=\"noopener\">bosons<\/a>&nbsp;(as opposed to&nbsp;<a title=\"Fermion\" href=\"https:\/\/en.wikipedia.org\/wiki\/Fermion\" target=\"_blank\" rel=\"noopener\">fermions<\/a>&nbsp;with half-integer spin). By the&nbsp;<a class=\"mw-redirect\" title=\"Spin-statistics theorem\" href=\"https:\/\/en.wikipedia.org\/wiki\/Spin-statistics_theorem\" target=\"_blank\" rel=\"noopener\">spin-statistics theorem<\/a>, all bosons obey <a href=\"https:\/\/en.wikipedia.org\/wiki\/Bose%E2%80%93Einstein_statistics\" target=\"_blank\" rel=\"noopener\">Bose\u2013Einstein statistics<\/a> (whereas all fermions obey&nbsp;<a title=\"Fermi\u2013Dirac statistics\" href=\"https:\/\/en.wikipedia.org\/wiki\/Fermi%E2%80%93Dirac_statistics\" target=\"_blank\" rel=\"noopener\">Fermi\u2013Dirac statistics<\/a>). As such, they do not obey the&nbsp;<a title=\"Pauli exclusion principle\" href=\"https:\/\/en.wikipedia.org\/wiki\/Pauli_exclusion_principle\" target=\"_blank\" rel=\"noopener\">Pauli exclusion principle<\/a>&nbsp;restrictions (no two identical fermions&nbsp; may occupy the same quantum state simultaneously).<\/p>\n<p>The following image sketches the difference between bosons and fermions:<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2533\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/bosons-and-fermions.png\" alt=\"\" width=\"400\" height=\"171\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/bosons-and-fermions.png 817w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/bosons-and-fermions-300x129.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/bosons-and-fermions-768x329.png 768w\" sizes=\"(max-width: 400px) 100vw, 400px\" \/><\/p>\n<p>Imagine two apparati, one containing a gas of bosons, while the other one contains a gas of (non interacting) fermions. As the temperature drops near&nbsp;<a href=\"http:\/\/en.wikipedia.org\/wiki\/Absolute_zero\" target=\"_blank\" rel=\"noopener\">absolute zero<\/a>, the gas of bosons collapses, forming a&nbsp;<a href=\"http:\/\/en.wikipedia.org\/wiki\/Bose-Einstein_condensate\" target=\"_blank\" rel=\"noopener\">Bose-Einstein Condensate<\/a>.&nbsp; Fermions can\u2019t reach this state, since they can\u2019t occupy the same quantum state. The total energy of the&nbsp;<a href=\"http:\/\/en.wikipedia.org\/wiki\/Fermi_gas\" target=\"_blank\" rel=\"noopener\">Fermi gas<\/a>&nbsp;is then larger than the sum of the single-particle ground state. The energy of the highest occupied quantum state is called the&nbsp;<a href=\"http:\/\/en.wikipedia.org\/wiki\/Fermi_energy\" target=\"_blank\" rel=\"noopener\">Fermi energy<\/a>.<\/p>\n<p>At the most elementary level, one can see&nbsp;fermions&nbsp;as the perfect candidate for&nbsp;building blocks of matter, while&nbsp;bosons&nbsp;are candidates for&nbsp;interactions. As such, photons are the carriers of the electromagnetic interaction.<\/p>\n<p>Photons are indeed the&nbsp;quantum&nbsp;of the&nbsp;<a title=\"Electromagnetic field\" href=\"https:\/\/en.wikipedia.org\/wiki\/Electromagnetic_field\" target=\"_blank\" rel=\"noopener\">electromagnetic field<\/a> which is understood as a&nbsp;<a class=\"mw-redirect\" title=\"Gauge field\" href=\"https:\/\/en.wikipedia.org\/wiki\/Gauge_field\" target=\"_blank\" rel=\"noopener\">gauge field<\/a>, i.e., as a field that results from requiring that a gauge symmetry holds independently at every position in&nbsp;spacetime. For the&nbsp;electromagnetic field, this gauge symmetry is the&nbsp;<a title=\"Abelian group\" href=\"https:\/\/en.wikipedia.org\/wiki\/Abelian_group\">Abelian<\/a>&nbsp;<a title=\"Unitary group\" href=\"https:\/\/en.wikipedia.org\/wiki\/Unitary_group\" target=\"_blank\" rel=\"noopener\">U(1) symmetry<\/a>, which reflects the ability to vary the&nbsp;<a title=\"Complex geometry\" href=\"https:\/\/en.wikipedia.org\/wiki\/Complex_geometry\" target=\"_blank\" rel=\"noopener\">phase<\/a>&nbsp;of a complex field without affecting&nbsp;<a title=\"Observable\" href=\"https:\/\/en.wikipedia.org\/wiki\/Observable\">observables<\/a> such as the&nbsp;energy&nbsp;or the&nbsp;<a title=\"Lagrangian (field theory)\" href=\"https:\/\/en.wikipedia.org\/wiki\/Lagrangian_(field_theory)\" target=\"_blank\" rel=\"noopener\">Lagrangian<\/a>.<\/p>\n<p>The quanta of an&nbsp;<a title=\"Gauge theory\" href=\"https:\/\/en.wikipedia.org\/wiki\/Gauge_theory\" target=\"_blank\" rel=\"noopener\">Abelian gauge field<\/a>&nbsp;must be massless, uncharged bosons, as long as the symmetry is not broken. The photon thus is massless, moves at the&nbsp;<a title=\"Speed of light\" href=\"https:\/\/en.wikipedia.org\/wiki\/Speed_of_light\">speed<\/a> <a title=\"Speed of light\" href=\"https:\/\/en.wikipedia.org\/wiki\/Speed_of_light\">of light<\/a> in vacuum, has zero&nbsp;<a title=\"Electric charge\" href=\"https:\/\/en.wikipedia.org\/wiki\/Electric_charge\">electric charge<\/a>.&nbsp;<\/p>\n<p>Each photon carries one quantum of energy, equal to <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chbar%5Comega&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hbar\\omega' title='\\hbar\\omega' class='latex' \/> (where <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Comega&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\omega' title='\\omega' class='latex' \/> is its <a href=\"https:\/\/en.wikipedia.org\/wiki\/Angular_frequency\" target=\"_blank\" rel=\"noopener\">angular frequency<\/a>). Photons are emitted in many natural processes. For example:<\/p>\n<ul>\n<li>When a charge is&nbsp;accelerated&nbsp;it emits&nbsp;<a title=\"Synchrotron radiation\" href=\"https:\/\/en.wikipedia.org\/wiki\/Synchrotron_radiation\" target=\"_blank\" rel=\"noopener\">synchrotron radiation<\/a>:<\/li>\n<\/ul>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2521\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/synchrotron.png\" alt=\"\" width=\"200\" height=\"105\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/synchrotron.png 524w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/synchrotron-300x157.png 300w\" sizes=\"(max-width: 200px) 100vw, 200px\" \/><\/p>\n<ul>\n<li>During an atomic&nbsp;or&nbsp;nuclear transition to a lower&nbsp;energy level, photons of various energy will be emitted, ranging from&nbsp;radio waves&nbsp;to&nbsp;gamma rays. The energy possessed by a single photon corresponds exactly to the transition between the discrete energy levels that emitted the photon:<\/li>\n<\/ul>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2520\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/spontaneous-emission.png\" alt=\"\" width=\"350\" height=\"180\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/spontaneous-emission.png 1029w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/spontaneous-emission-300x154.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/spontaneous-emission-768x395.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/spontaneous-emission-1024x526.png 1024w\" sizes=\"(max-width: 350px) 100vw, 350px\" \/><\/p>\n<ul>\n<li>Photons can also be emitted when a particle and its corresponding&nbsp;<a title=\"Antiparticle\" href=\"https:\/\/en.wikipedia.org\/wiki\/Antiparticle\" target=\"_blank\" rel=\"noopener\">antiparticle<\/a>&nbsp;are&nbsp;<a title=\"Annihilation\" href=\"https:\/\/en.wikipedia.org\/wiki\/Annihilation\" target=\"_blank\" rel=\"noopener\">annihilated<\/a>, like in the following <a title=\"Electron\u2013positron annihilation\" href=\"https:\/\/en.wikipedia.org\/wiki\/Electron%E2%80%93positron_annihilation\" target=\"_blank\" rel=\"noopener\">annihilation<\/a> diagrams:<\/li>\n<\/ul>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2517\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/annihilation.png\" alt=\"\" width=\"500\" height=\"175\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/annihilation.png 1536w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/annihilation-300x105.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/annihilation-768x269.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/annihilation-1024x358.png 1024w\" sizes=\"(max-width: 500px) 100vw, 500px\" \/><\/p>\n<p>Quantum computing with linear quantum optics has the advantage that its smallest unit of quantum information (the photon) is potentially free from decoherence: the quantum information stored in a photon tends to stay there.<\/p>\n<p>The downside is that photons do not naturally interact with each other, and in order to apply 2-qubit quantum gates such interactions are essential.<\/p>\n<p><strong>Linear Optical Quantum Computing<\/strong><\/p>\n<p>There are many implementations for&nbsp;quantum information processing and quantum computation. <a title=\"Quantum optics\" href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_optics\" target=\"_blank\" rel=\"noopener\">Optical quantum systems<\/a>&nbsp;are interesting candidates because they link quantum computation and&nbsp;<a class=\"mw-redirect\" title=\"Quantum communication\" href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_communication\" target=\"_blank\" rel=\"noopener\">quantum communication<\/a> within a same framework.&nbsp;<\/p>\n<p>Linear Optical Quantum Computing (LOQC) uses (mostly) linear optical&nbsp;elements to process&nbsp;quantum information, photon detectors and&nbsp;quantum memories&nbsp;to detect and store quantum information. These optical elements uses combinations of beam splitters,&nbsp;phase shifters, and&nbsp;mirrors:<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2535\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic.png\" alt=\"\" width=\"400\" height=\"351\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic.png 694w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-300x263.png 300w\" sizes=\"(max-width: 400px) 100vw, 400px\" \/><\/p>\n<p>Each linear optical element equivalently applies a&nbsp;<a title=\"Unitary transformation\" href=\"https:\/\/en.wikipedia.org\/wiki\/Unitary_transformation\">unitary transformation<\/a>&nbsp;on a finite number of qubits. This network of finite linear optics can realize any&nbsp;<a title=\"Quantum circuit\" href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_circuit\">quantum circuit<\/a>.&nbsp;<\/p>\n<p>Operations via linear optical elements preserve the photon statistics of input light. For example<\/p>\n<ul>\n<li>a&nbsp;<a class=\"mw-redirect\" title=\"Coherent light\" href=\"https:\/\/en.wikipedia.org\/wiki\/Coherent_light\" target=\"_blank\" rel=\"noopener\">coherent<\/a>&nbsp;(classical) light input produces a coherent light output<\/li>\n<li>a superposition of quantum states input yields a&nbsp;<a title=\"Nonclassical light\" href=\"https:\/\/en.wikipedia.org\/wiki\/Nonclassical_light\" target=\"_blank\" rel=\"noopener\">quantum light state<\/a>&nbsp;output.<\/li>\n<\/ul>\n<p>An intrinsic problem in using photons as information carriers is that photons hardly interact with each other. This potentially causes a scalability problem for LOQC, since nonlinear operations are hard to implement, which can increase the complexity of operators and hence can increase the resources required to realize a given computational function.<\/p>\n<p>One way to solve this problem is to bring <strong>nonlinear<\/strong> devices into the quantum network. For instance, the&nbsp;<a title=\"Kerr effect\" href=\"https:\/\/en.wikipedia.org\/wiki\/Kerr_effect\" target=\"_blank\" rel=\"noopener\">Kerr effect<\/a>&nbsp;can be applied into LOQC to make a single-photon <a title=\"Controlled NOT gate\" href=\"https:\/\/en.wikipedia.org\/wiki\/Controlled_NOT_gate\" target=\"_blank\" rel=\"noopener\">CNOT<\/a> gate and other operations.<\/p>\n<p>Quantum computing with continuous variables is also possible under the linear optics scheme.<\/p>\n<p><strong>KLM<br \/>\n<\/strong><\/p>\n<p>Knill, Laflamme and Milburn proved that it is possible to create universal quantum computers solely with linear optical tools. Ignoring error correction and other issues, one can implement elementary quantum gates using only mirrors, beam splitters and phase shifters using 1-qubit unitary operations.<\/p>\n<p>The unitary matrix associated with a beam splitter <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cmathbf%7BB%7D_%7B%5Ctheta%2C%5Cphi%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\mathbf{B}_{\\theta,\\phi}' title='\\mathbf{B}_{\\theta,\\phi}' class='latex' \/> is:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+U%28%5Cmathbf%7BB%7D_%7B%5Ctheta%2C%5Cphi%7D%29%3D+%5Cbegin%7Bbmatrix%7D+%5Cmathrm%7Bcos%7D%5C%3A%5Ctheta+%26+-e%5E%7Bi%5Cphi%7D%5C%3A%5Cmathrm%7Bsin%7D%5C%3A%5Ctheta+%5C%5C+e%5E%7B-i%5Cphi%7D%5C%3A%5Cmathrm%7Bsin%7D%5C%3A%5Ctheta+%26+%5Cmathrm%7Bcos%7D%5C%3A%5Ctheta+%5Cend%7Bbmatrix%7D+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle U(\\mathbf{B}_{\\theta,\\phi})= \\begin{bmatrix} \\mathrm{cos}\\:\\theta &amp; -e^{i\\phi}\\:\\mathrm{sin}\\:\\theta \\\\ e^{-i\\phi}\\:\\mathrm{sin}\\:\\theta &amp; \\mathrm{cos}\\:\\theta \\end{bmatrix} ' title='\\displaystyle U(\\mathbf{B}_{\\theta,\\phi})= \\begin{bmatrix} \\mathrm{cos}\\:\\theta &amp; -e^{i\\phi}\\:\\mathrm{sin}\\:\\theta \\\\ e^{-i\\phi}\\:\\mathrm{sin}\\:\\theta &amp; \\mathrm{cos}\\:\\theta \\end{bmatrix} ' class='latex' \/><\/p>\n<p>For a mirror the corresponding unitary operator is a rotation matrix given by:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Cmathbf%7BR%7D%28%5Ctheta%29+%3D+%5Cbegin%7Bbmatrix%7D+%5Cmathrm%7Bcos%7D%5C%3A%5Ctheta+%26+-%5Cmathrm%7Bsin%7D%5C%3A%5Ctheta+%5C%5C+%5Cmathrm%7Bsin%7D%5C%3A%5Ctheta+%26+%5Cmathrm%7Bcos%7D%5C%3A%5Ctheta+%5Cend%7Bbmatrix%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\mathbf{R}(\\theta) = \\begin{bmatrix} \\mathrm{cos}\\:\\theta &amp; -\\mathrm{sin}\\:\\theta \\\\ \\mathrm{sin}\\:\\theta &amp; \\mathrm{cos}\\:\\theta \\end{bmatrix}' title='\\displaystyle \\mathbf{R}(\\theta) = \\begin{bmatrix} \\mathrm{cos}\\:\\theta &amp; -\\mathrm{sin}\\:\\theta \\\\ \\mathrm{sin}\\:\\theta &amp; \\mathrm{cos}\\:\\theta \\end{bmatrix}' class='latex' \/><\/p>\n<p>Similarly, a phase shifter operator <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cmathbf%7BP%7D_%5Cphi&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\mathbf{P}_\\phi' title='\\mathbf{P}_\\phi' class='latex' \/> is associated with a unitary operator<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+U%28%5Cmathbf%7BP%7D_%5Cphi%29+%3D+%5Cbegin%7Bbmatrix%7D+e%5E%7Bi%5Cphi%7D+%26+0+%5C%5C+0+%26+1+%5Cend%7Bbmatrix%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle U(\\mathbf{P}_\\phi) = \\begin{bmatrix} e^{i\\phi} &amp; 0 \\\\ 0 &amp; 1 \\end{bmatrix}' title='\\displaystyle U(\\mathbf{P}_\\phi) = \\begin{bmatrix} e^{i\\phi} &amp; 0 \\\\ 0 &amp; 1 \\end{bmatrix}' class='latex' \/><\/p>\n<p>Since any two SU(2) rotations&nbsp;along orthogonal rotating axes can generate arbitrary rotations in the <a href=\"https:\/\/en.wikipedia.org\/wiki\/Bloch_sphere\" target=\"_blank\" rel=\"noopener\">Bloch sphere<\/a>, one can use a set of symmetric beam splitters and mirrors to realize an arbitrary SU(2) operators for QIP.<\/p>\n<p>The following figures are examples of implementing a&nbsp;<a href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_gate#Hadamard_gate\" target=\"_blank\" rel=\"noopener\">Hadamard gate<\/a>&nbsp;and a&nbsp;<a href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_gate#Pauli-X_gate_(=_NOT_gate)\" target=\"_blank\" rel=\"noopener\">Pauli-X gate<\/a>&nbsp;(NOT gate) by using beam splitters and mirrors:<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2543\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-gates.png\" alt=\"\" width=\"700\" height=\"235\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-gates.png 1638w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-gates-300x101.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-gates-768x257.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/optic-gates-1024x343.png 1024w\" sizes=\"(max-width: 700px) 100vw, 700px\" \/><\/p>\n<p>For further reading on LQOC &amp; KLM, I would recommend these two papers:<\/p>\n<ul>\n<li><a href=\"https:\/\/arxiv.org\/pdf\/quant-ph\/0006088.pdf\" target=\"_blank\" rel=\"noopener\">Efficient Linear Optics Quantum Computation<\/a> (E. Knill, R. Laflamme and G. Milburn)<\/li>\n<li><a href=\"https:\/\/arxiv.org\/pdf\/quant-ph\/0512071.pdf\" target=\"_blank\" rel=\"noopener\">Linear optical quantum computing<\/a> (Pieter Kok, W.J. Munro, Kae Nemoto,&nbsp; T.C. Ralph, Jonathan P. Dowling and G.J. Milburn)<\/li>\n<\/ul>\n<p><strong>Continuous-Variable model<\/strong><\/p>\n<p>Many physical systems are intrinsically continuous, with light being the prototypical example. Such systems reside in an infinite-dimensional Hilbert space, <strong>offering a paradigm for quantum computation which is distinct from the qubit model<\/strong>.<\/p>\n<p>This&nbsp;<strong>continuous-variable model<\/strong> (CV) takes its name from the fact that the quantum operators underlying the model have continuous spectra. The CV model is a natural fit for simulating bosonic systems (electromagnetic field, <a href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_harmonic_oscillator\" target=\"_blank\" rel=\"noopener\">harmonic oscillators<\/a>, <a href=\"https:\/\/en.wikipedia.org\/wiki\/Phonon\" target=\"_blank\" rel=\"noopener\">phonons<\/a>, Bose-Einstein condensates, &#8230;) and for settings where continuous quantum operators (such as <a href=\"https:\/\/en.wikipedia.org\/wiki\/Position_operator\" target=\"_blank\" rel=\"noopener\">position<\/a> or <a href=\"https:\/\/en.wikipedia.org\/wiki\/Momentum_operator\" target=\"_blank\" rel=\"noopener\">momentum<\/a>) are present.<\/p>\n<p>In <a title=\"Quantum optics\" href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_optics\" target=\"_blank\" rel=\"noopener\">quantum optics<\/a>, an <b>optical phase space<\/b> is a <a title=\"Phase space\" href=\"https:\/\/en.wikipedia.org\/wiki\/Phase_space\" target=\"_blank\" rel=\"noopener\">phase space<\/a> in which all <a title=\"Quantum state\" href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_state\">quantum states<\/a> of an optical system are described. Each point in the optical phase space corresponds to a unique state. For any such system, the plot of the quadratures (i.e position and momentum) against each other is called a <a title=\"Phase diagram\" href=\"https:\/\/en.wikipedia.org\/wiki\/Phase_diagram\">phase diagram<\/a>:<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2557\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/opd.png\" alt=\"\" width=\"450\" height=\"208\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/opd.png 1158w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/opd-300x138.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/opd-768x354.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/opd-1024x472.png 1024w\" sizes=\"(max-width: 450px) 100vw, 450px\" \/><\/p>\n<p>If the quadratures are functions of time then the optical phase diagram can show the evolution of a quantum optical system with time.<\/p>\n<p>When discussing the quantum theory of light, it is very common to use an electromagnetic oscillator as a model. An electromagnetic oscillator describes an oscillation of the electric field. And since the magnetic field is proportional to the rate of change of the electric field, this too oscillates. Systems composed of such oscillators can be described by an optical phase space.<\/p>\n<p>When quantized, such oscillators are described by <a href=\"https:\/\/en.wikipedia.org\/wiki\/Quantum_harmonic_oscillator\" target=\"_blank\" rel=\"noopener\">quantum harmonic oscillators<\/a>. which can be written in terms of <a href=\"https:\/\/en.wikipedia.org\/wiki\/Creation_and_annihilation_operators\" target=\"_blank\" rel=\"noopener\">creation and annihilation operators<\/a>.&nbsp;<\/p>\n<p>The <a href=\"https:\/\/en.wikipedia.org\/wiki\/Fock_space\" target=\"_blank\" rel=\"noopener\">Fock space<\/a> is commonly used in quantum mechanics to construct the quantum states space of a variable number of identical particles from a single particle Hilbert space <img src='https:\/\/s0.wp.com\/latex.php?latex=H&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='H' title='H' class='latex' \/>. It is the direct sum of the symmetric or antisymmetric tensors (depending on the bosonic or fermionic nature of the particles) in the tensor powers of a single-particle Hilbert:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+F_v+%3D+%5Cbigoplus_%7Bn%3D0%7D%5E%5Cinfty+S_v+H%5E%7B%5Cotimes+n%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle F_v = \\bigoplus_{n=0}^\\infty S_v H^{\\otimes n}' title='\\displaystyle F_v = \\bigoplus_{n=0}^\\infty S_v H^{\\otimes n}' class='latex' \/><\/p>\n<p>The creation operator <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Ba%7D%5E%5Cdagger+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{a}^\\dagger ' title='\\hat{a}^\\dagger ' class='latex' \/> acts on the Fock space, changing a <img src='https:\/\/s0.wp.com\/latex.php?latex=N&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N' title='N' class='latex' \/> particles state into a <img src='https:\/\/s0.wp.com\/latex.php?latex=N%2B1&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N+1' title='N+1' class='latex' \/> particles state. In the case of bosons, the <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Ba%7D%5E%5Cdagger+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{a}^\\dagger ' title='\\hat{a}^\\dagger ' class='latex' \/> operator creates a particle in state <img src='https:\/\/s0.wp.com\/latex.php?latex=i&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='i' title='i' class='latex' \/> such as<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Chat%7Ba%7D_i%5E%5Cdagger+%7C+N_1%2CN_2%2C+%5Ccdots%2C+N_i%2C+%5Ccdots+%5Crangle+%3D+%5Csqrt%7BN_i+%2B+1%7D+%7CN_1%2C+N_2%2C+%5Ccdots%2C+N_%7Bi%2B1%7D%2C%5Ccdots%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\hat{a}_i^\\dagger | N_1,N_2, \\cdots, N_i, \\cdots \\rangle = \\sqrt{N_i + 1} |N_1, N_2, \\cdots, N_{i+1},\\cdots\\rangle' title='\\displaystyle \\hat{a}_i^\\dagger | N_1,N_2, \\cdots, N_i, \\cdots \\rangle = \\sqrt{N_i + 1} |N_1, N_2, \\cdots, N_{i+1},\\cdots\\rangle' class='latex' \/><\/p>\n<p>Conversely, the annihilation operator <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Chat%7Ba%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\hat{a}' title='\\displaystyle \\hat{a}' class='latex' \/> acts on the Fock space, changing a <img src='https:\/\/s0.wp.com\/latex.php?latex=N&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N' title='N' class='latex' \/> particles state into a <img src='https:\/\/s0.wp.com\/latex.php?latex=N-1&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N-1' title='N-1' class='latex' \/> particles state. For bosonic states:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Chat%7Ba%7D_i%5E%7C+N_1%2CN_2%2C+%5Ccdots%2C+N_i%2C+%5Ccdots+%5Crangle+%3D+%5Csqrt%7BN_i%7D+%7CN_1%2C+N_2%2C+%5Ccdots%2C+N_%7Bi-1%7D%2C%5Ccdots%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\hat{a}_i^| N_1,N_2, \\cdots, N_i, \\cdots \\rangle = \\sqrt{N_i} |N_1, N_2, \\cdots, N_{i-1},\\cdots\\rangle' title='\\displaystyle \\hat{a}_i^| N_1,N_2, \\cdots, N_i, \\cdots \\rangle = \\sqrt{N_i} |N_1, N_2, \\cdots, N_{i-1},\\cdots\\rangle' class='latex' \/><\/p>\n<p>The number operator <img src='https:\/\/s0.wp.com\/latex.php?latex=N&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N' title='N' class='latex' \/> is defined as <img src='https:\/\/s0.wp.com\/latex.php?latex=N+%3D+%5Chat%7Ba%7D%5E%5Cdagger+%5Chat%7Ba%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N = \\hat{a}^\\dagger \\hat{a}' title='N = \\hat{a}^\\dagger \\hat{a}' class='latex' \/>, leading to the following canonical commutation relations:<\/p>\n<ul>\n<li><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Ba%2C%5Chat%7Ba%7D%5E%5Cdagger%5D+%3D+1+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='[a,\\hat{a}^\\dagger] = 1 ' title='[a,\\hat{a}^\\dagger] = 1 ' class='latex' \/><\/li>\n<li><img src='https:\/\/s0.wp.com\/latex.php?latex=%5BN%2C%5Chat%7Ba%7D%5E%5Cdagger%5D+%3D+%5Chat%7Ba%7D%5E%5Cdagger+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='[N,\\hat{a}^\\dagger] = \\hat{a}^\\dagger ' title='[N,\\hat{a}^\\dagger] = \\hat{a}^\\dagger ' class='latex' \/><\/li>\n<li><img src='https:\/\/s0.wp.com\/latex.php?latex=%5BN%2C%5Chat%7Ba%7D%5D+%3D+-%5Chat%7Ba%7D+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='[N,\\hat{a}] = -\\hat{a} ' title='[N,\\hat{a}] = -\\hat{a} ' class='latex' \/><\/li>\n<\/ul>\n<p>The <em>quadrature<\/em> operators can be defined in terms of creation an annihilation operators:<\/p>\n<ul>\n<li><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Chat%7Bx%7D+%5Cpropto+%5Csqrt%7B%5Cfrac%7B%5Chbar%7D%7B2%7D%7D+%28%5Chat%7Ba%7D+%2B+%5Chat%7Ba%7D%5E%5Cdagger%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\hat{x} \\propto \\sqrt{\\frac{\\hbar}{2}} (\\hat{a} + \\hat{a}^\\dagger)' title='\\displaystyle \\hat{x} \\propto \\sqrt{\\frac{\\hbar}{2}} (\\hat{a} + \\hat{a}^\\dagger)' class='latex' \/><\/li>\n<li><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Chat%7Bp%7D+%5Cpropto+-i+%5Csqrt%7B%5Cfrac%7B%5Chbar%7D%7B2%7D%7D+%28%5Chat%7Ba%7D+-+%5Chat%7Ba%7D%5E%5Cdagger%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\hat{p} \\propto -i \\sqrt{\\frac{\\hbar}{2}} (\\hat{a} - \\hat{a}^\\dagger)' title='\\displaystyle \\hat{p} \\propto -i \\sqrt{\\frac{\\hbar}{2}} (\\hat{a} - \\hat{a}^\\dagger)' class='latex' \/><\/li>\n<li>with <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5B%5Chat%7Bx%7D+%2C+%5Chat%7Bp%7D%5D+%3D+i+%5Chbar&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle [\\hat{x} , \\hat{p}] = i \\hbar' title='\\displaystyle [\\hat{x} , \\hat{p}] = i \\hbar' class='latex' \/><\/li>\n<\/ul>\n<p>The dichotomy between qubit and CV systems is perhaps most evident in the basis expansions of quantum states. For qubits, a discrete set of coefficients is used whereas CV systems have a&nbsp;<em>continuum<\/em>:<\/p>\n<ul>\n<li><strong>qubit<\/strong>: <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%7C%5Cpsi%5Crangle+%3D+%5Calpha+%7C0%5Crangle+%2B+%5Cbeta+%7C1%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle |\\psi\\rangle = \\alpha |0\\rangle + \\beta |1\\rangle' title='\\displaystyle |\\psi\\rangle = \\alpha |0\\rangle + \\beta |1\\rangle' class='latex' \/><\/li>\n<li><strong>qumode<\/strong>: <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%7C%5Cpsi%5Crangle+%3D+%5Cint+dx%5C%3B+%5Cpsi%28x%29+%7Cx%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle |\\psi\\rangle = \\int dx\\; \\psi(x) |x\\rangle' title='\\displaystyle |\\psi\\rangle = \\int dx\\; \\psi(x) |x\\rangle' class='latex' \/><\/li>\n<\/ul>\n<p>The following table draws a comparison of CV quantum computation with the qubit model:<\/p>\n<table width=\"450\" cellspacing=\"0\" cellpadding=\"0\" border=\"1\">\n<tbody>\n<tr style=\"border-top: none !important;\">\n<td style=\"text-align: center; border-top: none !important;\">&nbsp;<\/td>\n<td style=\"text-align: left; border-top: none !important;\"><b>qubit model<\/b><\/td>\n<td style=\"text-align: left; border-top: none !important;\"><b>Continuous-Variable model<\/b><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: left;\"><strong>Basic element<\/strong><\/td>\n<td style=\"text-align: left;\">qubits<\/td>\n<td style=\"text-align: left;\">qmodes<\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left; vertical-align: middle;\"><strong>Relevant operators<\/strong><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Pauli operators <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Csigma_x%2C+%5Csigma_y%2C+%5Csigma_z&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\sigma_x, \\sigma_y, \\sigma_z' title='\\sigma_x, \\sigma_y, \\sigma_z' class='latex' \/><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Quadrature operators <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Bx%7D%2C%5Chat%7Bp%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{x},\\hat{p}' title='\\hat{x},\\hat{p}' class='latex' \/> and creation-annihilation operators <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Ba%7D%5E%5Cdagger%2C+%5Chat%7Ba%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{a}^\\dagger, \\hat{a}' title='\\hat{a}^\\dagger, \\hat{a}' class='latex' \/><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left; vertical-align: middle;\"><strong>Common states<\/strong><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Pauli eigenstates <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C0%5Crangle%2C+%7C1%5Crangle%2C+%5Ccdots&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|0\\rangle, |1\\rangle, \\cdots' title='|0\\rangle, |1\\rangle, \\cdots' class='latex' \/><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Coherent states <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C%5Calpha%7C%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|\\alpha|\\rangle' title='|\\alpha|\\rangle' class='latex' \/>, Squeezed states <img src='https:\/\/s0.wp.com\/latex.php?latex=%7Cz%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|z\\rangle' title='|z\\rangle' class='latex' \/>,<br \/>\nNumber states <img src='https:\/\/s0.wp.com\/latex.php?latex=%7Cn%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|n\\rangle' title='|n\\rangle' class='latex' \/><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left; vertical-align: middle;\"><strong>Common gates<\/strong><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Phase Shift, Hadamard, CNOT, Pauli, &#8230;<\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Rotation, Displacement, Squeezing,<br \/>\nBeam Splitter, &#8230;<\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left; vertical-align: middle;\"><strong>Common measurements<\/strong><\/td>\n<td style=\"text-align: left; vertical-align: middle;\">Pauli-basis measurements <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C0%5Crangle+%5Clangle+0%7C%2C+%7C1%5Crangle%5Clangle+1%7C%2C+%5Ccdots&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|0\\rangle \\langle 0|, |1\\rangle\\langle 1|, \\cdots' title='|0\\rangle \\langle 0|, |1\\rangle\\langle 1|, \\cdots' class='latex' \/><\/td>\n<td style=\"text-align: left; vertical-align: middle;\"><a href=\"https:\/\/en.wikipedia.org\/wiki\/Homodyne_detection\" target=\"_blank\" rel=\"noopener\">Homodyne<\/a>, <a href=\"https:\/\/en.wikipedia.org\/wiki\/Heterodyne\" target=\"_blank\" rel=\"noopener\">Heterodyne<\/a>, Photon-counting <img src='https:\/\/s0.wp.com\/latex.php?latex=%7Cn%5Crangle+%5Clangle+n%7C&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|n\\rangle \\langle n|' title='|n\\rangle \\langle n|' class='latex' \/><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Note: qubit-based computations can be embedded into the CV picture, e.g., by using the Gottesman-Knill-Preskill embedding, so the CV model is as computationally powerful as its qubit counterparts.<\/p>\n<p>The starting point is the vacuum state <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C0%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|0\\rangle' title='|0\\rangle' class='latex' \/>. Other states can be created by evolving <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C0%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|0\\rangle' title='|0\\rangle' class='latex' \/> according to:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%7C%5Cpsi%5Crangle+%3D+%5Cmathrm%7Bexp%7D%28-itH%29+%7C0%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle |\\psi\\rangle = \\mathrm{exp}(-itH) |0\\rangle' title='\\displaystyle |\\psi\\rangle = \\mathrm{exp}(-itH) |0\\rangle' class='latex' \/><\/p>\n<p style=\"text-align: left;\">where <img src='https:\/\/s0.wp.com\/latex.php?latex=H&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='H' title='H' class='latex' \/> is a bosonic Hamiltonian. States where the Hamiltonian is at most quadratic in the operators <img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Bx%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{x}' title='\\hat{x}' class='latex' \/> and<img src='https:\/\/s0.wp.com\/latex.php?latex=%5Chat%7Bp%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\hat{p}' title='\\hat{p}' class='latex' \/> are called <strong>Gaussian<\/strong>. Typically, for a one-dimensional harmonic oscillator, the Hamiltonian of the particle is given by:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+H+%3D+%5Cfrac%7Bp%5E2%7D%7B2m%7D+%2B+%5Cfrac%7B1%7D%7B2%7Dm%5Comega%5E2+x%5E2&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle H = \\frac{p^2}{2m} + \\frac{1}{2}m\\omega^2 x^2' title='\\displaystyle H = \\frac{p^2}{2m} + \\frac{1}{2}m\\omega^2 x^2' class='latex' \/><\/p>\n<p>Solving <a href=\"https:\/\/en.wikipedia.org\/wiki\/Schr%C3%B6dinger_equation\" target=\"_blank\" rel=\"noopener\">Schr\u00f6dinger&#8217;s equation<\/a> <img src='https:\/\/s0.wp.com\/latex.php?latex=H%7C%5Cpsi%5Crangle+%3D+E%7C%5Cpsi%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='H|\\psi\\rangle = E|\\psi\\rangle' title='H|\\psi\\rangle = E|\\psi\\rangle' class='latex' \/> leads indeed to a Gaussian ground state<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%5Cpsi_0+%28x%29+%3D+%5CBig%28+%5Cfrac%7Bm%5Comega%7D%7B%5Cpi%5Chbar%7D%5CBig%29%5E%7B%5Cfrac%7B1%7D%7B4%7D%7D%5C%3B+e%5E%7B-m%5Comega+x%5E2%2F+2%5Chbar%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle \\psi_0 (x) = \\Big( \\frac{m\\omega}{\\pi\\hbar}\\Big)^{\\frac{1}{4}}\\; e^{-m\\omega x^2\/ 2\\hbar}' title='\\displaystyle \\psi_0 (x) = \\Big( \\frac{m\\omega}{\\pi\\hbar}\\Big)^{\\frac{1}{4}}\\; e^{-m\\omega x^2\/ 2\\hbar}' class='latex' \/><\/p>\n<p style=\"text-align: left;\">The <em>displacement<\/em> gives the center of the Gaussian, while the <em>squeezing<\/em> determines the variance and rotation of the distribution.<\/p>\n<p>Of course, complementary to (continuous) Gaussian states are the discrete Fock states. Let&#8217;s introduce the eigenstate <img src='https:\/\/s0.wp.com\/latex.php?latex=%7Cn%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|n\\rangle' title='|n\\rangle' class='latex' \/> of the number operator <img src='https:\/\/s0.wp.com\/latex.php?latex=N+%3D+%5Chat%7Ba%7D%5E%5Cdagger+%5Chat%7Ba%7D&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='N = \\hat{a}^\\dagger \\hat{a}' title='N = \\hat{a}^\\dagger \\hat{a}' class='latex' \/>.&nbsp; Each Gaussian state can be expanded in the number state basis. For example, coherent states <img src='https:\/\/s0.wp.com\/latex.php?latex=%7C%5Calpha%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='|\\alpha\\rangle' title='|\\alpha\\rangle' class='latex' \/> have the form:<\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+%7C%5Calpha%5Crangle+%3D+%5Cmathrm%7Bexp%7D%5CBig%28+-%5Cfrac%7B%7C%5Calpha%7C%5E2%7D%7B2%7D%5Csum_%7Bn%3D0%7D%5E%7B%5Cinfty%7D%5Cfrac%7B%5Calpha%5E2%7D%7B%5Csqrt%7Bn%21%7D%7D%5CBig%29+%7Cn%5Crangle&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle |\\alpha\\rangle = \\mathrm{exp}\\Big( -\\frac{|\\alpha|^2}{2}\\sum_{n=0}^{\\infty}\\frac{\\alpha^2}{\\sqrt{n!}}\\Big) |n\\rangle' title='\\displaystyle |\\alpha\\rangle = \\mathrm{exp}\\Big( -\\frac{|\\alpha|^2}{2}\\sum_{n=0}^{\\infty}\\frac{\\alpha^2}{\\sqrt{n!}}\\Big) |n\\rangle' class='latex' \/><\/p>\n<p><strong>Continuous-Variable gates<\/strong><\/p>\n<p>Unitary operations can always be associated with a generating Hamiltonian <img src='https:\/\/s0.wp.com\/latex.php?latex=H&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='H' title='H' class='latex' \/><\/p>\n<p style=\"text-align: center;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+U+%3D+%5Cmathrm%7Bexp%7D%28-itH%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle U = \\mathrm{exp}(-itH)' title='\\displaystyle U = \\mathrm{exp}(-itH)' class='latex' \/><\/p>\n<p>A CV quantum computer is said to be universal if it can implement with a finite number of steps any unitary which is polynomial in the annihilation \/ creation operator. Two kinds of universal gates can be distinguished:<\/p>\n<ul>\n<li><strong>Gaussian gates<\/strong>: these are gates that are quadratic, such as displacement, rotation, squeezing, and beam splitter&nbsp;gates. These are equivalent to the <a href=\"https:\/\/en.wikipedia.org\/wiki\/Gottesman%E2%80%93Knill_theorem\" target=\"_blank\" rel=\"noopener\">Clifford quantum gates<\/a> from the qubit model (<a class=\"mw-redirect\" title=\"Hadamard gate\" href=\"https:\/\/en.wikipedia.org\/wiki\/Hadamard_gate\" target=\"_blank\" rel=\"noopener\">Hadamard gates<\/a>,&nbsp;<a title=\"Controlled NOT gate\" href=\"https:\/\/en.wikipedia.org\/wiki\/Controlled_NOT_gate\" target=\"_blank\" rel=\"noopener\">controlled NOT gates<\/a>, Phase Gate);<\/li>\n<li><strong>Non-Gaussian gates<\/strong>: these are gates which are of degree 3 or higher, e.g., the&nbsp;cubic phase gate.&nbsp;<\/li>\n<\/ul>\n<p>The following table lists a few fundamental CV gates :<\/p>\n<table width=\"450\" cellspacing=\"0\" cellpadding=\"0\" border=\"1\">\n<tbody>\n<tr style=\"border-top: none !important;\">\n<td style=\"text-align: left; border-top: none !important;\"><b>Gate<\/b><\/td>\n<td style=\"text-align: left; border-top: none !important;\"><b>Unitary<\/b><\/td>\n<td style=\"text-align: middle; border-top: none !important;\"><b>Symbol<\/b><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left;\"><span style=\"color: #999999;\"><strong>Displacement<\/strong><\/span><\/td>\n<td style=\"text-align: left;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+D_i%28%5Calpha%29%3D%5Cmathrm%7Bexp%7D%28%5Calpha+%5Chat%7Ba%7D%5E%5Cdagger_i+-+%5Calpha%5E%2A+%5Chat%7Ba%7D_i%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle D_i(\\alpha)=\\mathrm{exp}(\\alpha \\hat{a}^\\dagger_i - \\alpha^* \\hat{a}_i)' title='\\displaystyle D_i(\\alpha)=\\mathrm{exp}(\\alpha \\hat{a}^\\dagger_i - \\alpha^* \\hat{a}_i)' class='latex' \/><\/td>\n<td style=\"text-align: middle; vertical-align: middle;\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-2602\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/dgate.png\" alt=\"\" width=\"50\" height=\"26\"><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left;\"><span style=\"color: #999999;\"><strong>Rotation<\/strong><\/span><\/td>\n<td style=\"text-align: left;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+R_i%28%5Cphi%29%3D%5Cmathrm%7Bexp%7D%28i%5Cphi+N_i%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle R_i(\\phi)=\\mathrm{exp}(i\\phi N_i)' title='\\displaystyle R_i(\\phi)=\\mathrm{exp}(i\\phi N_i)' class='latex' \/><\/td>\n<td style=\"text-align: middle; vertical-align: middle;\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-2607\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/rgate.png\" alt=\"\" width=\"50\" height=\"27\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/rgate.png 155w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/rgate-150x85.png 150w\" sizes=\"(max-width: 50px) 100vw, 50px\" \/><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left;\"><span style=\"color: #999999;\"><strong>Squeezing<\/strong><\/span><\/td>\n<td style=\"text-align: left;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+S_i%28z%29%3D+%5Cmathrm%7Bexp%7D+%5CBig%28+%5Cfrac%7B1%7D%7B2%7D+%28z%5E%2A+%5Chat%7Ba%7D_i%5E2+-+z%5Chat%7Ba%7D%5E%7B%5Cdagger+2%7D_i%29+%5CBig%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle S_i(z)= \\mathrm{exp} \\Big( \\frac{1}{2} (z^* \\hat{a}_i^2 - z\\hat{a}^{\\dagger 2}_i) \\Big)' title='\\displaystyle S_i(z)= \\mathrm{exp} \\Big( \\frac{1}{2} (z^* \\hat{a}_i^2 - z\\hat{a}^{\\dagger 2}_i) \\Big)' class='latex' \/><\/td>\n<td style=\"text-align: middle; vertical-align: middle;\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-2610\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sgate.png\" alt=\"\" width=\"50\" height=\"29\"><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left;\"><span style=\"color: #999999;\"><strong>Beam Splitter<\/strong><\/span><\/td>\n<td style=\"text-align: left;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+B_%7Bi%2Cj%7D%28%5Ctheta%2C+%5Cphi%29%3D+%5Cmathrm%7Bexp%7D+%5CBig%28+%5Ctheta+%28e%5E%7Bi%5Cphi%7D%5Chat%7Ba%7D_i%5E%7B%5Cdagger%7D+%5Chat%7Ba%7D_j+-+e%5E%7B-i%5Cphi%7D+%5Chat%7Ba%7D_i+%5Chat%7Ba%7D_j%5E%7B%5Cdagger%7D+%29+%5CBig%29+&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle B_{i,j}(\\theta, \\phi)= \\mathrm{exp} \\Big( \\theta (e^{i\\phi}\\hat{a}_i^{\\dagger} \\hat{a}_j - e^{-i\\phi} \\hat{a}_i \\hat{a}_j^{\\dagger} ) \\Big) ' title='\\displaystyle B_{i,j}(\\theta, \\phi)= \\mathrm{exp} \\Big( \\theta (e^{i\\phi}\\hat{a}_i^{\\dagger} \\hat{a}_j - e^{-i\\phi} \\hat{a}_i \\hat{a}_j^{\\dagger} ) \\Big) ' class='latex' \/><\/td>\n<td style=\"text-align: middle; vertical-align: middle;\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-2612\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/bsgate.png\" alt=\"\" width=\"50\" height=\"36\"><\/td>\n<\/tr>\n<tr style=\"vertical-align: middle;\">\n<td style=\"text-align: left; vertical-align: middle;\"><span style=\"color: #999999;\"><strong>Cubic Phase<\/strong><\/span><\/td>\n<td style=\"text-align: left; vertical-align: middle;\"><img src='https:\/\/s0.wp.com\/latex.php?latex=%5Cdisplaystyle+V_i%28%5Cgamma%29%3D+%5Cmathrm%7Bexp%7D%5Cbig%28+i+%5Cfrac%7B%5Cgamma%7D%7B6%7D%5Chat%7Bx%7D%5E3_i+%5Cbig%29&#038;bg=ffffff&#038;fg=000000&#038;s=0' alt='\\displaystyle V_i(\\gamma)= \\mathrm{exp}\\big( i \\frac{\\gamma}{6}\\hat{x}^3_i \\big)' title='\\displaystyle V_i(\\gamma)= \\mathrm{exp}\\big( i \\frac{\\gamma}{6}\\hat{x}^3_i \\big)' class='latex' \/><\/td>\n<td style=\"text-align: middle; vertical-align: middle;\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone wp-image-2613\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/vgate.png\" alt=\"\" width=\"50\" height=\"27\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/vgate.png 158w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/vgate-150x85.png 150w\" sizes=\"(max-width: 50px) 100vw, 50px\" \/><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><strong>Software architecture for&nbsp;photonic quantum computing<\/strong><\/p>\n<p>The <a href=\"https:\/\/strawberryfields.readthedocs.io\/en\/latest\/index.html\" target=\"_blank\" rel=\"noopener\">Strawberry Fields<\/a> software platform is a full-stack Python library for designing, simulating, and optimizing Continuous Variable quantum optical circuits.<\/p>\n<p>The following sketch outlines its key elements :<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2624\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/strawberryfields.png\" alt=\"\" width=\"700\" height=\"293\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/strawberryfields.png 1049w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/strawberryfields-300x126.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/strawberryfields-768x321.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/strawberryfields-1024x429.png 1024w\" sizes=\"(max-width: 700px) 100vw, 700px\" \/><\/p>\n<p>The platform consists of three main components:<\/p>\n<ul>\n<li>an API for quantum programming based a programming language named Blackbird;<\/li>\n<li>a suite of three virtual quantum computer backends, built in <a href=\"https:\/\/en.wikipedia.org\/wiki\/NumPy\" target=\"_blank\" rel=\"noopener\">NumPy<\/a> and <a href=\"https:\/\/en.wikipedia.org\/wiki\/TensorFlow\" target=\"_blank\" rel=\"noopener\">Tensorflow<\/a>, each targeting specialized uses (Optimization, Quantum Machine Learning, &#8230;);<\/li>\n<li>an engine which can compile Blackbird programs on various backends, including the three built-in simulators, and \u2013 in the near future \u2013 photonic quantum information processors.<\/li>\n<\/ul>\n<p>The frontend encompasses the Strawberry Fields Python API and the Blackbird quantum programming language. These elements provide access points for users to design quantum circuits.<\/p>\n<p>Blackbird is a quantum assembly language, capable of representing the basic continuous-variable (CV) states, gates, and measurements. There are four main types of operations:<\/p>\n<ul>\n<li>State preparation;<\/li>\n<li>Gate application;<\/li>\n<li>Measurements;<\/li>\n<li>Adding and removing subsystems that these operations act on.<\/li>\n<\/ul>\n<p>The following example shows an implementation of <a href=\"https:\/\/www.quantum-bits.org\/?p=1857\" target=\"_blank\" rel=\"noopener\">Quantum Teleportation<\/a> with the Blackbird quantum programming language:<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2637\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/qtwithblackbird-414x1024.png\" alt=\"\" width=\"375\" height=\"927\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/qtwithblackbird-414x1024.png 414w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/qtwithblackbird-121x300.png 121w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/qtwithblackbird.png 599w\" sizes=\"(max-width: 375px) 100vw, 375px\" \/><\/p>\n<p>For a backend, the engine can target one of the included quantum computer simulators:<\/p>\n<ul>\n<li>Gaussian Representation (using Numpy)<\/li>\n<li>Fock Representation (using Numpy and Tensorflow)<\/li>\n<\/ul>\n<p>When CV quantum processors will become available, the engine will also build and run circuits on these devices.<\/p>\n<p>Applications, such as the<a href=\"https:\/\/strawberryfields.ai\/\" target=\"_blank\" rel=\"noopener\"> Strawberry Fields Interactive<\/a> website, can be built by leveraging the frontend API. The following screenshot demonstrates how to perform quantum teleportation using Strawberry Fields Interactive :<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter wp-image-2628\" src=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sfinteractive.png\" alt=\"\" width=\"850\" height=\"330\" srcset=\"https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sfinteractive.png 1215w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sfinteractive-300x116.png 300w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sfinteractive-768x298.png 768w, https:\/\/www.quantum-bits.org\/wp-content\/uploads\/2018\/07\/sfinteractive-1024x397.png 1024w\" sizes=\"(max-width: 850px) 100vw, 850px\" \/><\/p>\n<p>More information on the platform is available here: <a href=\"https:\/\/www.xanadu.ai\/software\/\" target=\"_blank\" rel=\"noopener\">https:\/\/www.xanadu.ai\/software\/<\/a><\/p>\n<p>Yes, the authors are obviously fans of The Beatles \ud83d\ude42<\/p>\n<p><span style=\"font-size: 8pt;\">Note: to speedup the writing process, a few paragraphs and illustrations of this post are based selected papers and Wikipedia articles on the subject of LQOC as well as a few excerpts from the Strawberry Fields Open-source software for photonic quantum computing documentation.<\/span><\/p>\n","protected":false},"excerpt":{"rendered":"<p>The worldwide quest to build practical quantum computers is undergoing a critical period. Fault-tolerant quantum computers will soon provide significant computational speedups for problems like factoring, search, or linear algebra, &#8230;<\/p>\n","protected":false},"author":1,"featured_media":3846,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"ngg_post_thumbnail":0},"categories":[6],"tags":[],"_links":{"self":[{"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/posts\/2513"}],"collection":[{"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=2513"}],"version-history":[{"count":0,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/posts\/2513\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=\/wp\/v2\/media\/3846"}],"wp:attachment":[{"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=2513"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=2513"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.quantum-bits.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=2513"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}