{"id":86900,"date":"2025-08-21T19:07:10","date_gmt":"2025-08-21T19:07:10","guid":{"rendered":"https:\/\/www.newsbeep.com\/ca\/86900\/"},"modified":"2025-08-21T19:07:10","modified_gmt":"2025-08-21T19:07:10","slug":"universal-quantum-gate-set-for-gottesman-kitaev-preskill-logical-qubits","status":"publish","type":"post","link":"https:\/\/www.newsbeep.com\/ca\/86900\/","title":{"rendered":"Universal quantum gate set for Gottesman\u2013Kitaev\u2013Preskill logical qubits"},"content":{"rendered":"<p>GKP code<\/p>\n<p>Here we give some background on GKP codes, highlighting some of the more salient features in our experiment and simulations. For a more comprehensive summary, see the original proposal<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 6\" title=\"Gottesman, D., Kitaev, A. &amp; Preskill, J. Encoding a qubit in an oscillator. Phys. Rev. A 64, 012310 (2001).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR6\" id=\"ref-link-section-d7536052e5878\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a> and more recent reviews<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 48\" title=\"Grimsmo, A. L. &amp; Puri, S. Quantum error correction with the Gottesman-Kitaev-Preskill code. PRX Quantum 2, 020101 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR48\" id=\"ref-link-section-d7536052e5882\" rel=\"nofollow noopener\" target=\"_blank\">48<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 49\" title=\"Brady, A. J., Eickbusch, A., Singh, S., Wu, J. &amp; Zhuang, Q. Advances in bosonic quantum error correction with Gottesman-Kitaev-Preskill codes: theory, engineering and applications. Prog. Quantum Electron. 93, 100496 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR49\" id=\"ref-link-section-d7536052e5885\" rel=\"nofollow noopener\" target=\"_blank\">49<\/a>. The GKP code is a stabilizer code defined in a bosonic Hilbert space. Its stabilizer generators define a set of lattice points in phase space, which define the code. A family of GKP codes exists, each corresponding to a distinct lattice in phase space<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 46\" title=\"Royer, B., Singh, S. &amp; Girvin, S. Encoding qubits in multimode grid states. PRX Quantum 3, 010335 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR46\" id=\"ref-link-section-d7536052e5889\" rel=\"nofollow noopener\" target=\"_blank\">46<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 50\" title=\"Conrad, J., Eisert, J. &amp; Arzani, F. Gottesman-Kitaev-Preskill codes: a lattice perspective. Quantum 6, 648 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR50\" id=\"ref-link-section-d7536052e5892\" rel=\"nofollow noopener\" target=\"_blank\">50<\/a>. Here we focus on the square GKP code, which is the most fundamental qubit code in this family. The stabilizer generators for this code are given by the two commuting displacement operators<\/p>\n<p>$${\\hat{S}}_{X}=\\hat{D}(2\\alpha ),\\quad {\\hat{S}}_{Z}=\\hat{D}(2\\beta ),$$<\/p>\n<p>\n                    (2)\n                <\/p>\n<p>where \\(\\beta ={\\rm{i}}\\alpha ={\\rm{i}}\\sqrt{\\pi \/2}\\). These admit the following logical Pauli operators<\/p>\n<p>$${\\hat{X}}_{{\\rm{L}}}=\\hat{D}(\\alpha ),\\quad {\\hat{Y}}_{{\\rm{L}}}=\\hat{D}(\\alpha +\\beta ),\\quad {\\hat{Z}}_{{\\rm{L}}}=\\hat{D}(\\;\\beta ),$$<\/p>\n<p>\n                    (3)\n                <\/p>\n<p>which commute with the stabilizer generators. Furthermore, a logical CZL gate takes the simple form<\/p>\n<p>$${\\hat{{\\rm{C}}}}Z_{{\\rm{L}}}={{\\rm{e}}}^{{\\rm{i}}{\\hat{q}}_{1}{\\hat{q}}_{2}},$$<\/p>\n<p>\n                    (4)\n                <\/p>\n<p>where \\({\\hat{q}}_{j}=({\\hat{a}}_{j}+{\\hat{a}}_{j}^{\\dagger })\/\\sqrt{2}\\) denotes the position operator for mode j.<\/p>\n<p>The code words, which are defined to be the \u00b11 eigenstates of the stabilizer generators in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>), may be expressed in phase space as an infinite sum of coherent states localized at the lattice points defined by \u03b1 and \u03b2:<\/p>\n<p>$$\\begin{array}{rcl}\\left\\vert +{Z}_{{\\rm{L}}}^{\\;{\\rm{ideal}}}\\right\\rangle &amp;=&amp;\\mathop{\\sum }\\limits_{k,l=-\\infty }^{\\infty }{{\\rm{e}}}^{-{\\rm{i}}\\pi kl}{\\left\\vert 2k\\alpha +l\\beta \\right\\rangle }_{{\\rm{c}}},\\\\ \\left\\vert -{Z}_{{\\rm{L}}}^{\\;{\\rm{ideal}}}\\right\\rangle &amp;=&amp;\\mathop{\\sum }\\limits_{k,l=-\\infty }^{\\infty }{{\\rm{e}}}^{-{\\rm{i}}\\pi (kl+l\/2)}{\\left\\vert (2k+1)\\alpha +l\\beta \\right\\rangle }_{{\\rm{c}}},\\end{array}$$<\/p>\n<p>\n                    (5)\n                <\/p>\n<p>where \\({\\left\\vert \\gamma \\right\\rangle }_{{\\rm{c}}}=\\hat{D}(\\gamma )\\left\\vert 0\\right\\rangle\\) is a coherent state and |0\u3009 is the vacuum state. These states are the \u00b11 eigenstates of ZL, and similar expressions for the \u00b11 eigenstates of the other logical Pauli operators, \\(\\left\\vert \\pm\\;{X}_{{\\rm{L}}}^{\\;{\\rm{ideal}}}\\right\\rangle \\,\\left\\vert \\pm\\;{Y}_{{\\rm{L}}}^{\\;{\\rm{ideal}}}\\right\\rangle\\), may be obtained by taking superpositions of these states.<\/p>\n<p>Owing to their infinite extent in phase space, these code words have infinite energy. Finite-energy GKP states may be defined by applying the envelope operator onto the infinite energy code words:<\/p>\n<p>$$\\left\\vert \\pm {P}_{{\\rm{L}}}\\right\\rangle =\\frac{{{\\rm{e}}}^{-{\\varDelta }^{2}{a}^{\\dagger }\\hat{a}}}{{\\mathcal{N}}}\\left\\vert \\pm {P}_{{\\rm{L}}}^{\\;{\\rm{ideal}}}\\right\\rangle ,$$<\/p>\n<p>\n                    (6)\n                <\/p>\n<p>where P\u2009=\u2009X, Y, Z, and \\({\\mathcal{N}}=\\sqrt{\\left\\langle \\pm {P}_{{\\rm{L}}}^{{\\rm{ideal}}}\\right\\vert {e}^{-2{\\Delta }^{2}{a}^{\\dagger }\\hat{a}}\\left\\vert \\pm {P}_{{\\rm{L}}}^{{\\rm{ideal}}}\\right\\rangle }\\). These are approximate \u00b11 eigenstates of the stabilizer generators and \u0394\u2009\u2208\u2009[0, 1] parameterizes this approximation, with \u0394\u21920 recovering the infinite energy code words. Crucially, the logical operators in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>) are not exact logical operators for the finite-energy code words. Methods to construct logical operators that are exact for finite-energy code words via the envelope operator, for example, \\({{\\rm{e}}}^{-{\\varDelta }^{2}{\\hat{a}}^{\\dagger }\\hat{a}}{\\hat{Z}}_{{\\rm{L}}}{{\\rm{e}}}^{{\\varDelta }^{2}{\\hat{a}}^{\\dagger }\\hat{a}}\\), have been explored for SQ gates in ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 17\" title=\"Royer, B., Singh, S. &amp; Girvin, S. Stabilization of finite-energy Gottesman-Kitaev-Preskill states. Phys. Rev. Lett. 125, 260509 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR17\" id=\"ref-link-section-d7536052e7498\" rel=\"nofollow noopener\" target=\"_blank\">17<\/a> and for the CZL gate in ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 19\" title=\"Rojkov, I. et al. Two-qubit operations for finite-energy Gottesman&#x2013;Kitaev&#x2013;Preskill encodings. Phys. Rev. Lett. 133, 100601 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR19\" id=\"ref-link-section-d7536052e7507\" rel=\"nofollow noopener\" target=\"_blank\">19<\/a>.<\/p>\n<p>We stress that the construction for finite-energy GKP states in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ6\" rel=\"nofollow noopener\" target=\"_blank\">6<\/a>) is not unique, and that the exact form of the approximation will depend on the details of the experimental platform that the states are prepared in. This motivates our use of an optimizer to find a logical operation that is tailored to the exact form of the approximate GKP states. However, for numerical purposes, we find an alternative definition for the finite-energy GKP code, which was used in refs. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 32\" title=\"Kolesnikow, X. C., Bomantara, R. W., Doherty, A. C. &amp; Grimsmo, A. L. Gottesman-Kitaev-Preskill state preparation using periodic driving. Phys. Rev. Lett. 132, 130605 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR32\" id=\"ref-link-section-d7536052e7517\" rel=\"nofollow noopener\" target=\"_blank\">32<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 51\" title=\"Rymarz, M., Bosco, S., Ciani, A. &amp; DiVincenzo, D. P. Hardware-encoding grid states in a nonreciprocal superconducting circuit. Phys. Rev. X 11, 011032 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR51\" id=\"ref-link-section-d7536052e7520\" rel=\"nofollow noopener\" target=\"_blank\">51<\/a>, particularly useful. Here the code states are defined to be the quasidegenerate ground states of the Hamiltonian<\/p>\n<p>$${\\hat{H}}_{{\\rm{GKP}}}={\\omega }_{0}{\\hat{a}}^{\\dagger }\\hat{a}-J(\\cos (2\\sqrt{\\pi }\\hat{q})+\\cos (2\\sqrt{\\pi }\\hat{p})),$$<\/p>\n<p>\n                    (7)\n                <\/p>\n<p>where \\(\\hat{q}=(\\hat{a}+{\\hat{a}}^{\\dagger })\/\\sqrt{2}\\) and \\(\\hat{p}=-{\\rm{i}}(\\hat{a}-{\\hat{a}}^{\\dagger })\/\\sqrt{2}\\) are the position and momentum operators for a single bosonic mode, respectively; \u03c90 represents the characteristic energy scale of the harmonic oscillator; and J quantifies the stabilizer potential enforcing the periodic structure of the GKP code. In the limit of large J\/\u03c90, the ground states are approximately logical Hadamard eigenstates in the |+ZL\u3009 basis,<\/p>\n<p>$$\\left\\vert +{H}_{{\\rm{L}}}\\right\\rangle =\\cos (\\pi \/8)\\left\\vert +{Z}_{{\\rm{L}}}\\right\\rangle +\\sin (\\pi \/8)\\left\\vert -{Z}_{{\\rm{L}}}\\right\\rangle ,$$<\/p>\n<p>\n                    (8a)\n                <\/p>\n<p>$$\\left\\vert -{H}_{{\\rm{L}}}\\right\\rangle =-\\sin (\\pi \/8)\\left\\vert +{Z}_{{\\rm{L}}}\\right\\rangle +\\cos (\\pi \/8)\\left\\vert -{Z}_{{\\rm{L}}}\\right\\rangle ,$$<\/p>\n<p>\n                    (8b)\n                <\/p>\n<p>and the squeezing parameter \u0394 of these states is related to the energy scales through \u0394\u2009=\u2009(\u03c90\/(4\u03c0J))1\/4. We find numerical diagonalization of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ7\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>) to be an efficient procedure to obtain the finite-energy GKP code words numerically. The relations in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ9\" rel=\"nofollow noopener\" target=\"_blank\">8<\/a>) may be inverted to obtain the |+ZL\u3009 code words, and the other logical Pauli eigenstates obtained from the appropriate superpositions of these states.<\/p>\n<p>Furthermore, the squeezing parameter of these states in each quadrature \u0394X\/Z may be independently calculated from their stabilizer expectation values through the relations<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 52\" title=\"Duivenvoorden, K., Terhal, B. M. &amp; Weigand, D. Single-mode displacement sensor. Phys. Rev. A 95, 012305 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR52\" id=\"ref-link-section-d7536052e8248\" rel=\"nofollow noopener\" target=\"_blank\">52<\/a><\/p>\n<p>$${\\Delta }_{X}=\\sqrt{-\\frac{1}{2\\pi }\\log \\left[| \\langle {\\hat{S}}_{X}\\rangle {| }^{2}\\right]},$$<\/p>\n<p>\n                    (9a)\n                <\/p>\n<p>$${\\Delta }_{Z}=\\sqrt{-\\frac{1}{2\\pi }\\log \\left[| \\langle {\\hat{S}}_{Z}\\rangle {| }^{2}\\right]},$$<\/p>\n<p>\n                    (9b)\n                <\/p>\n<p>which can also be expressed in units of decibels, with \\({\\Delta }_{X\/Z}\\,({\\rm{dB}})=-10{\\log }_{10}[{\\Delta }_{X\/Z}^{2}]\\). In this work, we set J\/\u03c90\u2009=\u20095.95 in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ7\" rel=\"nofollow noopener\" target=\"_blank\">7<\/a>) to target logical states {|+ZL\u3009, |\u2212ZL\u3009, |+XL\u3009, |+YL\u3009} with squeezing parameters [\u0394X, \u0394Z] of {[8.39, 7.90], [8.36, 8.88], [7.90, 8.39], [8.38, 8.38]}\u2009dB, respectively.<\/p>\n<p>Logical measurements on the GKP code using SSSD<\/p>\n<p>SSSD provides a way to divide the infinite-dimensional Hilbert space into a logical subsystem, containing the logical information encoding in a GKP state, and a gauge subsystem, referred to as the stabilizer subsystem<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"Shaw, M. H., Doherty, A. C. &amp; Grimsmo, A. L. Stabilizer subsystem decompositions for single- and multimode Gottesman-Kitaev-Preskill codes. PRX Quantum 5, 010331 (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR29\" id=\"ref-link-section-d7536052e8626\" rel=\"nofollow noopener\" target=\"_blank\">29<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 53\" title=\"Pantaleoni, G., Baragiola, B. Q. &amp; Menicucci, N. C. Modular bosonic subsystem codes. Phys. Rev. Lett. 125, 040501 (2020).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR53\" id=\"ref-link-section-d7536052e8629\" rel=\"nofollow noopener\" target=\"_blank\">53<\/a>. Here we make use of this formalism to accurately and efficiently readout the logical information encoded in our states.<\/p>\n<p>Logical readout for the ideal square GKP code corresponds to measuring the expectation values of the logical Pauli operators, \\(\\{\\langle {\\hat{X}}_{{\\rm{L}}}\\rangle ,\\langle {\\hat{Y}}_{{\\rm{L}}}\\rangle ,\\langle {\\hat{Z}}_{{\\rm{L}}}\\rangle \\}\\). The expectation values of these logical Pauli operators are obtained by applying an SDF interaction, which transfers information from the bosonic mode to the ancilla, followed by an ancilla measurement. Setting the phases \u03d5s and \u03d5m to be time independent with \u03d5s\u2009=\u20090, the application of the SDF interaction from equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) for a time t realizes the operation \\(\\hat{D}(\\gamma {\\hat{\\sigma }}_{x}\/2)\\), where \\(\\gamma =-{\\rm{i}}\\Omega t{{\\rm{e}}}^{-{\\rm{i}}{\\phi }_{{\\rm{m}}}}\\). The magnitude of displacement, |\u03b3|, is controlled by the duration of the SDF pulse, t, and the phase of the displacement, arg(\u03b3), is controlled by the SDF phase, \u03d5m. Following the SDF pulse, the state measurement of the ancilla gives \\(\\langle {\\hat{\\sigma }}_{z}\\rangle ={\\rm{Re}}[\\langle \\hat{D}(\\gamma (t,{\\phi }_{{\\rm{m}}}))\\rangle ]\\), whereas the imaginary part can be obtained by applying an ancilla qubit rotation before the SDF. Setting \u03b3\u2009=\u2009{\u03b1, \u03b1\u2009+\u2009\u03b2, \u03b2} leads to the measurements of \\(\\{\\langle {\\hat{X}}_{{\\rm{L}}}\\rangle ,\\langle {\\hat{Y}}_{{\\rm{L}}}\\rangle ,\\langle {\\hat{Z}}_{{\\rm{L}}}\\rangle \\}\\).<\/p>\n<p>The expectation value of joint Pauli logical operators on two bosonic modes is obtained by applying two sequential SDF operations followed by an ancilla measurement, resulting in \\(\\mathrm{Re}\\,[\\langle \\hat{D}(\\gamma )\\otimes \\hat{D}(\\delta )\\rangle ]\\).<\/p>\n<p>To decode a GKP state, measuring the Pauli expectation values of finite-energy GKP states with the operators of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>) is only valid in the limit \u0394\u21920. As explained in the main text, one can instead measure the expectation value of the Pauli measurement operators, which was done for the SQ gates and Bell state characterization. These operators are defined as a summation over displacement operators on the logical GKP lattice<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 30\" title=\"Shaw, M. H., Doherty, A. C. &amp; Grimsmo, A. L. Logical gates and read-out of superconducting Gottesman-Kitaev-Preskill qubits. Preprint at &#010;                https:\/\/arxiv.org\/abs\/2403.02396&#010;                &#010;               (2024).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR30\" id=\"ref-link-section-d7536052e9255\" rel=\"nofollow noopener\" target=\"_blank\">30<\/a>:<\/p>\n<p>$${\\hat{X}}_{{\\rm{m}}}=\\frac{1}{\\pi }\\mathop{\\sum }\\limits_{n=-N}^{N-1}\\frac{{(-1)}^{n}}{n+\\frac{1}{2}}\\hat{D}\\left(\\alpha \\left(2n+1\\right)\\right),$$<\/p>\n<p>\n                    (10a)\n                <\/p>\n<p>$${\\hat{Y}}_{{\\rm{m}}}=\\frac{1}{{\\pi }^{2}}\\mathop{\\sum }\\limits_{m,n=-N}^{N-1}\\frac{\\hat{D}\\left(\\alpha (2n+1)+\\beta (2m+1)\\right)}{\\left(n+\\frac{1}{2}\\right)\\left(m+\\frac{1}{2}\\right)},$$<\/p>\n<p>\n                    (10b)\n                <\/p>\n<p>$${\\hat{Z}}_{{\\rm{m}}}=\\frac{1}{\\pi }\\mathop{\\sum }\\limits_{m=-N}^{N-1}\\frac{{(-1)}^{m}}{m+\\frac{1}{2}}\\,\\hat{D}\\left(\\beta \\left(2m+1\\right)\\right),$$<\/p>\n<p>\n                    (10c)\n                <\/p>\n<p>where N sets a truncation, and the identity Pauli measurement operator \\({\\hat{I}}_{{\\rm{m}}}\\) coincides with the identity operator on Fock space. The expectation values \\(\\langle {\\hat{X}}_{{\\rm{m}}}\\rangle\\), \\(\\langle {\\hat{Y}}_{{\\rm{m}}}\\rangle\\) and \\(\\langle {\\hat{Z}}_{{\\rm{m}}}\\rangle\\) are obtained by summing over the expectation values of the displacement operator, \\(\\langle \\hat{D}(\\gamma )\\rangle\\). Here \u03b3 is varied to sample the various points determined by equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>). Expectations of joint Pauli measurement operators for two bosonic modes, such as \\(\\langle {\\hat{X}}_{{\\rm{m}}}\\otimes {\\hat{Z}}_{{\\rm{m}}}\\rangle\\), can be obtained by taking the tensor product of any two single-mode Pauli measurement operator of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>). Their joint expectation value will also be a summation over the expectation values of two-mode displacement operators.<\/p>\n<p>The number of required measurements is reduced by using the fact that \\(\\langle \\hat{D}(\\gamma )\\rangle ={\\langle \\hat{D}(-\\gamma )\\rangle }^{* }\\). With this, the total number of measurements to calculate the three single-mode expectation values of the Pauli measurement operators of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>) is Nmeas\u2009=\u20092N\u2009+\u20092N2. For the expectation values of the 15 non-trivial two-mode Pauli measurement operators, the number of measurements is Nmeas\u2009=\u20094N\u2009+\u20098N2\u2009+\u20098N3\u2009+\u20094N4. The finite-energy envelope of the GKP states implies that \\(\\langle \\hat{D}(\\gamma )\\rangle \\approx 0\\) for large enough |\u03b3|, meaning that the sums in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>) converge with N. In both SQ gate and Bell state experiments, we choose a truncation of N\u2009=\u20092 to minimize the number of measurements required, but still achieve a good approximation (in Supplementary Section <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">C<\/a>, we quantify the error arising from this truncation). This results in Nmeas\u2009=\u200912 and Nmeas\u2009=\u2009168 measurements for logical SQ and TQ states, respectively. For the CZL gate experiment, we choose not to measure the Pauli measurement operators because the required number of measurements for all 16 states is 2,688, leading to a much longer experimental runtime. Instead, we choose to measure the expectation values of the logical Pauli operators of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>), which requires a total of 240 measurements, and incorporate the error from using these operators into our error budget (Supplementary Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>).<\/p>\n<p>State tomography<\/p>\n<p>In the SQ experiment, the logical density matrices \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{in}}}\\) and \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{out}}}\\) are reconstructed from the Pauli measurement expectation values using the relation<\/p>\n<p>$${\\hat{\\rho }}_{{\\rm{L}}}=\\frac{1}{2}{\\hat{E}}^{0}+\\frac{1}{2}\\mathop{\\sum }\\limits_{i=1}^{3}\\left\\langle {\\hat{E}}_{{\\rm{m}}}^{i}\\right\\rangle {\\hat{E}}^{i},$$<\/p>\n<p>\n                    (11)\n                <\/p>\n<p>where \\({\\hat{E}}^{i}\\in \\{\\hat{I},{\\hat{\\sigma }}_{x},{\\hat{\\sigma }}_{y},{\\hat{\\sigma }}_{z}\\}\\) are the usual one-qubit Pauli operators and \\({\\hat{E}}_{{\\rm{m}}}^{i}\\in \\{{\\hat{I}}_{{\\rm{m}}},{\\hat{X}}_{{\\rm{m}}},{\\hat{Y}}_{{\\rm{m}}},{\\hat{Z}}_{{\\rm{m}}}\\}\\) are the Pauli measurement operators, which can be obtained from the experiment using equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>). We construct a tomographically complete set of \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{in}}}\\) by preparing one of the four states {|+ZL\u3009, |\u2212ZL\u3009, |+XL\u3009, |+YL\u3009} and taking measurements for each of the expectation values in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ15\" rel=\"nofollow noopener\" target=\"_blank\">11<\/a>), leading to a total of 4\u2009\u00d7\u200912\u2009=\u200948 measurements (Supplementary Section <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">B<\/a>). Reconstructing each \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{out}}}\\) is performed in a similar manner, but with the gate applied to each initial state.<\/p>\n<p>In the TQ experiment, the input and output density matrices are reconstructed from Pauli expectation values using the relation<\/p>\n<p>$${\\hat{\\rho }}_{{\\rm{L}}}=\\frac{1}{4}{\\hat{E}}^{0}+\\frac{1}{4}\\mathop{\\sum }\\limits_{i=1}^{15}\\langle {\\hat{E}}_{{\\rm{L}}}^{i}\\rangle {\\hat{E}}^{i},$$<\/p>\n<p>\n                    (12)\n                <\/p>\n<p>where \\({\\hat{E}}^{i}\\in {\\{\\hat{I},{\\hat{\\sigma }}_{x},{\\hat{\\sigma }}_{y},{\\hat{\\sigma }}_{z}\\}}^{\\otimes 2}\\) are the TQ Pauli operators and \\({\\hat{E}}_{{\\rm{L}}}^{i}\\in {\\{{\\hat{I}}_{{\\rm{L}}},{\\hat{X}}_{{\\rm{L}}},{\\hat{Y}}_{{\\rm{L}}},{\\hat{Z}}_{{\\rm{L}}}\\}}^{\\otimes 2}\\) are the logical Pauli operators of equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ3\" rel=\"nofollow noopener\" target=\"_blank\">3<\/a>). We first reconstruct each \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{in}}}\\) by preparing one of the 16 possible input states from the set {|+ZL\u3009, |\u2212ZL\u3009, |+XL\u3009, |+YL\u3009}\u22972 and measuring the 15 different non-trivial logical Pauli operators, resulting in 16\u2009\u00d7\u200915\u2009=\u2009240 measurements. Similarly, each \\({\\hat{\\rho }}_{{\\rm{L}}}^{{\\rm{out}}}\\) is retrieved by applying the CZL gate after preparing one of the input states.<\/p>\n<p>The Bell state is characterized using logical QST, which aims to reconstruct the logical density matrix of the experimentally prepared state. In the Bell state experiment, we reconstruct the logical density matrix from<\/p>\n<p>$${\\hat{\\rho }}_{{\\rm{L}}}=\\frac{1}{4}{\\hat{E}}^{0}+\\frac{1}{4}\\mathop{\\sum }\\limits_{i=1}^{15}\\langle {\\hat{E}}_{{\\rm{m}}}^{i}\\rangle {\\hat{E}}^{i},$$<\/p>\n<p>\n                    (13)\n                <\/p>\n<p>where \\({\\hat{E}}^{i}\\in {\\{\\hat{I},{\\hat{\\sigma }}_{x},{\\hat{\\sigma }}_{y},{\\hat{\\sigma }}_{z}\\}}^{\\otimes 2}\\) are the TQ Pauli operators and \\({\\hat{E}}_{{\\rm{m}}}^{i}\\in {\\{{\\hat{I}}_{{\\rm{m}}},{\\hat{X}}_{{\\rm{m}}},{\\hat{Y}}_{{\\rm{m}}},{\\hat{Z}}_{{\\rm{m}}}\\}}^{\\otimes 2}\\) are the Pauli measurement operators. The expectation values are calculated from equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ12\" rel=\"nofollow noopener\" target=\"_blank\">10<\/a>), where a truncation of N\u2009=\u20092 results in 168 unique measurements. We apply a post-processing step in the logical Bell state analysis, ensuring that the reconstructed state is physical. After obtaining \\({\\hat{\\rho }}_{{\\rm{L}}}\\) from equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ17\" rel=\"nofollow noopener\" target=\"_blank\">13<\/a>), we use a convex optimizer to obtain a physical density matrix, \\({\\hat{\\rho }}_{{\\rm{L}},{\\rm{Bell}}}\\), which minimizes \\(| | {\\hat{\\rho }}_{{\\rm{L}},{\\rm{Bell}}}-{\\hat{\\rho }}_{{\\rm{L}}}| {| }_{2}\\) (ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 54\" title=\"Riofr&#xED;o, C. A. et al. Experimental quantum compressed sensing for a seven-qubit system. Nat. Commun. 8, 15305 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR54\" id=\"ref-link-section-d7536052e12519\" rel=\"nofollow noopener\" target=\"_blank\">54<\/a>). \\({\\hat{\\rho }}_{{\\rm{L}},{\\rm{Bell}}}\\) is constrained to be non-negative definite and \\(\\mathrm{Tr}\\,({\\hat{\\rho }}_{{\\rm{L}},{\\rm{Bell}}})=1\\). Expectation values of the physical density matrix \\({\\hat{\\rho }}_{{\\rm{L}},{\\rm{Bell}}}\\) are plotted in Fig. <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"figure anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Fig4\" rel=\"nofollow noopener\" target=\"_blank\">4<\/a>.<\/p>\n<p>Uncertainty analysis<\/p>\n<p>Uncertainties on the logical fidelities presented in this work are determined using a non-parametric bootstrap procedure. For each dataset, we resample the raw spin readout counts with replacement to create 1,000 pseudo-datasets of identical size. Each resampled dataset is processed through the full analysis: first extracting the relevant displacement operator expectation values, then reconstructing the logical density matrix or process matrix and finally calculating the fidelity metric. The standard deviation of the resulting distribution of 1,000 fidelity values is taken as the one standard deviation (1\u03c3) statistical uncertainty reported in the main text.<\/p>\n<p>Calibration and experimental drift<\/p>\n<p>To mitigate errors from slow experimental drifts, we implement an interleaved calibration routine throughout the experiment. The motional mode frequencies \u03c9x and \u03c9y, corresponding to the radial-x and radial-y modes, are periodically recalibrated approximately every 30\u2009s using a calibration protocol described in ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 55\" title=\"MacDonell, R. J. et al. Predicting molecular vibronic spectra using time-domain analog quantum simulation. Chem. Sci. 14, 9439&#x2013;9451 (2023).\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#ref-CR55\" id=\"ref-link-section-d7536052e12752\" rel=\"nofollow noopener\" target=\"_blank\">55<\/a>. Independent measurements reveal that the frequency drifts by \u03c3\u03b4\u2009\u2248\u20092\u03c0\u2009\u00d7\u20093.2\u2009Hz between calibrations, and that drifts in \u03c9x and \u03c9y are highly correlated. We estimate errors due to frequency drifts from numerical simulations by repeating a similar analysis, as carried out in Supplementary Section <a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">C<\/a>. The SDF Hamiltonian in equation (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41567-025-03002-8#Equ1\" rel=\"nofollow noopener\" target=\"_blank\">1<\/a>) is modified by adding a noisy Hamiltonian term of the form \\({\\delta }_{x}{\\hat{a}}_{1}^{\\dagger }{\\hat{a}}_{1}+{\\delta }_{y}{\\hat{a}}_{2}^{\\dagger }{\\hat{a}}_{2}\\). We set \u03b4\u2009=\u2009\u03b4x\u2009=\u2009\u03b4y to model-correlated frequency drifts, and \\(\\delta \\approx {\\mathcal{N}}(0,{\\sigma }_{\\delta }^{2})\\) is sampled from a normal distribution. The resulting infidelity is below 5\u2009\u00d7\u200910\u22123 across all experiments.<\/p>\n<p>By contrast, drifts in the qubit frequency and SDF Rabi rate occur on much slower timescales and are not included in the scheduled calibration routine. These parameters do not appreciably drift over the duration of the experiments and have a negligible impact on the fidelities relative to drifts in the motional frequency.<\/p>\n","protected":false},"excerpt":{"rendered":"GKP code Here we give some background on GKP codes, highlighting some of the more salient features in&hellip;\n","protected":false},"author":2,"featured_media":86901,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[24],"tags":[8878,49,48,8877,8882,8881,3673,8876,8879,8880,314,3230,8269,66,8875],"class_list":{"0":"post-86900","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-physics","8":"tag-atomic","9":"tag-ca","10":"tag-canada","11":"tag-classical-and-continuum-physics","12":"tag-complex-systems","13":"tag-condensed-matter-physics","14":"tag-general","15":"tag-mathematical-and-computational-physics","16":"tag-molecular","17":"tag-optical-and-plasma-physics","18":"tag-physics","19":"tag-quantum-information","20":"tag-qubits","21":"tag-science","22":"tag-theoretical"},"_links":{"self":[{"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/posts\/86900","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/comments?post=86900"}],"version-history":[{"count":0,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/posts\/86900\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/media\/86901"}],"wp:attachment":[{"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/media?parent=86900"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/categories?post=86900"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.newsbeep.com\/ca\/wp-json\/wp\/v2\/tags?post=86900"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}