{"id":361053,"date":"2025-11-07T01:15:30","date_gmt":"2025-11-07T01:15:30","guid":{"rendered":"https:\/\/www.europesays.com\/us\/361053\/"},"modified":"2025-11-07T01:15:30","modified_gmt":"2025-11-07T01:15:30","slug":"effective-descriptions-of-bosonic-systems-can-be-considered-complete","status":"publish","type":"post","link":"https:\/\/www.europesays.com\/us\/361053\/","title":{"rendered":"Effective descriptions of bosonic systems can be considered complete"},"content":{"rendered":"<p>Bosonic quantum information theory<\/p>\n<p>Quantum information theory with bosonic systems typically involves modeling unbounded operators over infinite-dimensional Hilbert spaces with continuous spectrum, and is thus often referred to as continuous-variable quantum information theory<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Braunstein, S. L. &amp; Van Loock, P. Quantum information with continuous variables. Rev. Mod. Phys. 77, 513 (2005).\" href=\"#ref-CR26\" id=\"ref-link-section-d59843725e790\">26<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" title=\"Weedbrook, C. et al. Gaussian quantum information. Rev. Mod. Phys. 84, 621 (2012).\" href=\"#ref-CR27\" id=\"ref-link-section-d59843725e790_1\">27<\/a>,<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 28\" title=\"Adesso, G., Ragy, S. &amp; Lee, A. R. Continuous variable quantum information: Gaussian states and beyond. Open Syst. Inf. Dyn. 21, 1440001 (2014).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR28\" id=\"ref-link-section-d59843725e793\" rel=\"nofollow noopener\" target=\"_blank\">28<\/a>. In this framework, systems corresponding to individual quantum harmonic oscillators are called qumodes, or simply modes, and quantum states of m-mode bosonic systems are elements of an infinite-dimensional Hilbert space \\({{{\\mathcal{H}}}}^{\\otimes m}\\), with each single-mode Hilbert space spanned by a countably infinite (Fock) basis \\({\\{\\left\\vert n\\right\\rangle \\}}_{n\\in {\\mathbb{N}}}\\).<\/p>\n<p>Creation and annihilation operators \\({\\hat{a}}^{{\\dagger} }\\) and \\(\\hat{a}\\) for each mode are defined by their action on the Fock basis as \\({\\hat{a}}^{{\\dagger} }\\left\\vert n\\right\\rangle=\\sqrt{n+1}\\left\\vert n+1\\right\\rangle\\), \\(\\hat{a}\\left\\vert n+1\\right\\rangle=\\sqrt{n+1}\\left\\vert n\\right\\rangle\\) and \\(\\hat{a}\\left\\vert 0\\right\\rangle=0\\). The (unbounded) canonical bosonic operators are then given by \\(\\hat{q}=\\frac{1}{\\sqrt{2}}(\\hat{a}+{\\hat{a}}^{{\\dagger} })\\) and \\(\\hat{p}=\\frac{1}{i\\sqrt{2}}(\\hat{a}-{\\hat{a}}^{{\\dagger} })\\) and satisfy the canonical commutation relation \\([\\hat{q},\\hat{p}]=i\\hat{I}\\), where \\(\\hat{I}\\) is the identity operator. Polynomial Hamiltonians over m modes are the Hermitian operators of the form \\(P({\\hat{q}}_{1},{\\hat{p}}_{1},\\ldots,{\\hat{q}}_{m},{\\hat{p}}_{m})\\), where P is a polynomial and where \\(({\\hat{q}}_{1},{\\hat{p}}_{1},\\ldots,{\\hat{q}}_{m},{\\hat{p}}_{m})\\) are the canonical operators of the modes 1,\u00a0\u2026,\u00a0m. The number operator is defined as \\({\\hat{n}}={\\hat{a}}^{{\\dagger} }{\\hat{a}}\\), and its expectation value is the average particle number of a bosonic mode, which we refer to as energy hereafter. For all \\(d\\in {\\mathbb{N}}\\), we denote by \\({{{\\mathcal{H}}}}_{d}: \\!\\!={{\\rm{span}}}{\\{\\left\\vert n\\right\\rangle \\}}_{0\\le n\\le d}\\) the subspace of states with a number of particles at most d.<\/p>\n<p>In some sense, not all states in the Hilbert space are valid physical states, even if they are normalized: for instance, some may have infinite energy, such as \\(\\sqrt{6}{\\sum }_{n\\ge 1}\\frac{1}{\\pi n}\\left\\vert n\\right\\rangle\\). Similarly, there are normalized states with infinite fluctuations in either position or momentum or both. To avoid inconsistencies and ensure that the action of polynomial Hamiltonians leads to physical states, it is common to define a set of physical states as a dense subspace of the Hilbert space known as the Schwartz space<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 20\" title=\"Hall, B. C. Quantum theory for mathematicians, 267 (Springer, 2013).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR20\" id=\"ref-link-section-d59843725e1946\" rel=\"nofollow noopener\" target=\"_blank\">20<\/a>\u2014informally, the set of states with bounded canonical operator moments, with a position wave function and all of its derivatives decaying sufficiently fast at infinity, which we denote by \\({{\\mathcal{S}}}\\subset {{\\mathcal{H}}}\\) (see the Methods section). Similarly, denoting by \\({{\\mathcal{U}}}({{\\mathcal{H}}})\\) the group of unitary operators over \\({{\\mathcal{H}}}\\), we define<\/p>\n<p>$${{\\mathcal{U}}}({{\\mathcal{S}}}) : \\!\\!=\\left\\{{\\hat{U}}\\in {{\\mathcal{U}}({{\\mathcal{H}}})\\,| \\,\\forall \\left\\vert \\psi \\right\\rangle \\in {{\\mathcal{S}}},\\,{\\hat{U}}\\left\\vert \\psi \\right\\rangle \\in {{\\mathcal{S}}}}\\right\\}.$$<\/p>\n<p>\n                    (1)\n                <\/p>\n<p>These form a group of physical unitary operators, i.e., which map physical states to physical states. This set can be readily extended to the multimode setting. Physical unitary operators have the property that they map bosonic states with a finite number of particles to states of finite energy: for all \\(n\\in {\\mathbb{N}}\\) and all \\(\\hat{U}\\in {{\\mathcal{U}}}({{\\mathcal{S}}})\\),<\/p>\n<p>$${E}_{\\hat{U}}(n) : \\!\\!={\\sup }_{\\left\\vert \\psi \\right\\rangle \\in {{{\\mathcal{H}}}}_{n}}\\left\\langle \\psi \\right\\vert {\\hat{U}}^{{\\dagger} }\\hat{n}\\hat{U}\\left\\vert \\psi \\right\\rangle <\/p>\n<p>\n                    (2)\n                <\/p>\n<p>In that case, the supremum is a maximum and the quantity \\({E}_{\\hat{U}}(n)\\) can be thought of as the maximal amount of energy involved when implementing the unitary evolution \\(\\hat{U}\\) of an initial state with a number of particles at most n. This is a slightly weaker property than energy-limited quantum dynamics<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 29\" title=\"van Luijk, L. Energy-limited quantum dynamics. Commun. Math. Phys. 406, 120 (2025).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR29\" id=\"ref-link-section-d59843725e2491\" rel=\"nofollow noopener\" target=\"_blank\">29<\/a> which map state of finite energy to states of finite energy.<\/p>\n<p>We denote by \\(D(\\rho,\\sigma )=\\frac{1}{2}\\parallel \\rho -\\sigma {\\parallel }_{1}\\) the trace norm distance between density operators<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 22\" title=\"Nielsen, M. A. &amp; Chuang, I. L. Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, 2011), 10th edn.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR22\" id=\"ref-link-section-d59843725e2571\" rel=\"nofollow noopener\" target=\"_blank\">22<\/a>, which is used to define the diamond norm for quantum maps. The diamond norm induces a distinguishability measure between maps that is often too stringent for quantum mechanical applications involving unbounded operators, as physical unitary channels may be always maximally separated for this norm, even though they may have similar effects on all states below a certain energy<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 30\" title=\"Winter, A. Energy-constrained diamond norm with applications to the uniform continuity of continuous variable channel capacities. arXiv preprint arXiv:1712.10267 (2017).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR30\" id=\"ref-link-section-d59843725e2575\" rel=\"nofollow noopener\" target=\"_blank\">30<\/a>. Instead, the strong topology of pointwise convergence provides a suitable notion of continuity in infinite dimensions, which is implied by the closeness in energy-constrained diamond norm (see the Methods section), defined as<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 31\" title=\"Shirokov, M. E. On the energy-constrained diamond norm and its application in quantum information theory. Probl. Inf. Transm. 54, 20&#x2013;33 (2018).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR31\" id=\"ref-link-section-d59843725e2579\" rel=\"nofollow noopener\" target=\"_blank\">31<\/a>:<\/p>\n<p>$$\\parallel {{\\mathcal{E}}}{\\parallel }_{\\diamond }^{E} : \\!\\!={\\sup }_{{{\\rm{Tr}}}[\\rho (\\hat{H}\\otimes {\\hat{I}})]\\le E}\\parallel ({{\\mathcal{E}}}\\otimes {{\\rm{id}}})\\rho {\\parallel }_{1}.$$<\/p>\n<p>\n                    (3)\n                <\/p>\n<p>Here, \\({{\\mathcal{E}}}\\) is a Hermitian-preserving map, \u03c1 is a density operator over \\({{\\mathcal{H}}}\\otimes {{{\\mathcal{H}}}}^{{\\prime} }\\), with \\({{{\\mathcal{H}}}}^{{\\prime} }\\) isomorphic to \\({{\\mathcal{H}}}\\), and \\(\\hat{H}\\) is a Hamiltonian over \\({{\\mathcal{H}}}\\) with ground state energy equal to 0 used to specify the energy bound E\u00a0&gt;\u00a00. We employ this topology hereafter, with \\(\\hat{H}=\\hat{n}\\) the number operator.<\/p>\n<p>Effective dimension of physical unitary evolutions<\/p>\n<p>Having introduced the necessary notation, we show in this section that any physical unitary channel can be approximated to arbitrary precision by a finite-dimensional unitary channel, where the quality of the approximation is measured by the energy-constrained diamond norm. This is captured by the notion of approximate effective dimension, which we define as follows:<\/p>\n<p>                  Definition 1<\/p>\n<p>(Approximate effective dimension). Let \\(\\hat{U}\\) be a single-mode unitary operator, let E\u22650 be an energy parameter, and let \u03f5\u00a0&gt;\u00a00 be an approximation parameter. Let also \\(d\\in {\\mathbb{N}}\\). The unitary operator \\(\\hat{U}\\) has (E,\u00a0\u03f5)-approximate effective dimension d\u00a0+\u00a01 if there exists a unitary operator \\({\\hat{V}}_{d}\\) over \\({{{\\mathcal{H}}}}_{d}\\) which approximates \\(\\hat{U}\\) in the following sense: for any unitary operator \\({\\hat{V}}^{{\\prime} }\\) over \\({{\\rm{span}}}{(\\left\\vert n\\right\\rangle )}_{n\\ &gt; \\ d}\\), denoting \\(\\hat{V}={\\hat{V}}_{d}\\oplus {\\hat{V}}^{{\\prime} }\\), \\({{\\mathcal{V}}}=\\hat{V}{{\\boldsymbol{\\cdot }}}{\\hat{V}}^{{\\dagger} }\\) and \\({{\\mathcal{U}}}=\\hat{U}{{\\boldsymbol{\\cdot }}}{\\hat{U}}^{{\\dagger} }\\),<\/p>\n<p>$$\\parallel {{\\mathcal{U}}}-{{\\mathcal{V}}}{\\parallel }_{\\diamond }^{E}\\le \\epsilon .$$<\/p>\n<p>\n                    (4)\n                <\/p>\n<p>When that is the case, we say that \\({\\hat{V}}_{d}\\) is a (d\u00a0+\u00a01)-dimensional (E,\u00a0\u03f5)-approximation of \\(\\hat{U}\\).<\/p>\n<p>When a unitary dynamics has a finite approximate effective dimension, it can effectively be modeled by a finite-dimensional unitary evolution. We now give our first main result:<\/p>\n<p>                  Theorem 1<\/p>\n<p>(Effective dimension of physical unitary channels). Let \\({{\\mathcal{U}}}=\\hat{U}{{\\boldsymbol{\\cdot }}}{\\hat{U}}^{{\\dagger} }\\) be a single-mode unitary channel, with \\(\\hat{U}\\in {{\\mathcal{U}}}({{\\mathcal{S}}})\\) a physical unitary operator. Let E\u22650 be an energy parameter, and let \u03f5\u00a0&gt;\u00a00 be an approximation parameter. Then, there exists \\(d={{\\mathcal{O}}}\\left(\\frac{E}{{\\epsilon }^{4}}{E}_{\\hat{U}}\\left(\\frac{64E}{{\\epsilon }^{2}}\\right)\\right)\\in {\\mathbb{N}}\\) such that \\(\\hat{U}\\) has (E,\u00a0\u03f5)-approximate effective dimension d\u00a0+\u00a01. Moreover, a finite-dimensional (E,\u00a0\u03f5)-approximation of \\(\\hat{U}\\) can be computed in polynomial time in d, E, and \\(\\frac{1}{\\epsilon }\\).<\/p>\n<p>The proof of this result is based on constructing a good unitary operator approximation of a truncated unitary operator, while carefully controlling truncation parameters for the state space, before and after applying the unitary channel, so that the truncated unitary operator is itself a good approximation of the original unitary operator. We refer to the Methods section for the main techniques and to the\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">Supplementary Information<\/a> for the proof.<\/p>\n<p>A direct consequence of Theorem 1 is that the corresponding sequence of (E,\u00a0\u03f5)-approximations strongly converges to \\(\\hat{U}\\) as \u03f5 goes to 0 in trace distance (see the Methods section). From an operational standpoint, this implies that we can effectively truncate the Hilbert space of physical bosonic computations to a certain effective dimension, as long as the cut-off is high enough. As such, Theorem 1 provides rigorous bounds on the effective dimension required to simulate bosonic dynamics with provable accuracy<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 9\" title=\"Tong, Y., Albert, V. V., McClean, J. R., Preskill, J. &amp; Su, Y. Provably accurate simulation of gauge theories and bosonic systems. Quantum 6, 816 (2022).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR9\" id=\"ref-link-section-d59843725e3975\" rel=\"nofollow noopener\" target=\"_blank\">9<\/a>. In other words, there will be no catastrophic loss of information in the measurements statistics when restricting to the cut-off evolution: by the operational property of the trace distance, results of a cut-off bosonic computation will be indistinguishable from results of the original computation. However, this cut-off comes at a price: Theorem 1 identifies a direct relation between the effective dimension sufficient to embed a unitary evolution to a good precision and the energy \\({E}_{\\hat{U}}\\) of the unitary evolution, as defined in Eq. (<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"equation anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#Equ2\" rel=\"nofollow noopener\" target=\"_blank\">2<\/a>). Such an energy bound may grow very fast for seemingly simple dynamics, such as alternating Gaussian and non-Gaussian unitary evolutions<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 25\" title=\"Chabaud, U., Joseph, M., Mehraban, S. &amp; Motamedi, A. Bosonic Quantum Computational Complexity. In Meka, R. (ed.) 16th Innovations in Theoretical Computer Science Conference (ITCS 2025), 325 of Leibniz International Proceedings in Informatics (LIPIcs), 33:1&#x2013;33:19 (Schloss Dagstuhl &#x2013; Leibniz-Zentrum f&#xFC;r Informatik, Dagstuhl, Germany, 2025). arXiv:2410.04274.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR25\" id=\"ref-link-section-d59843725e4022\" rel=\"nofollow noopener\" target=\"_blank\">25<\/a>. We give an explicit construction in the\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">Supplementary Information<\/a> showing that this growth may be arbitrarily fast.<\/p>\n<p>Universality of polynomial Hamiltonians<\/p>\n<p>We have shown that bosonic unitary evolutions may be well approximated by finite-dimensional ones. In this section, we give our second main result, namely that any finite-dimensional Hamiltonian generating such a finite-dimensional unitary evolution can be realized exactly by a polynomial Hamiltonian:<\/p>\n<p>                  Theorem 2<\/p>\n<p>(Finite-dimensional universality of polynomial Hamiltonians). Let \\(d\\in {\\mathbb{N}}\\) and let \\(\\hat{H}\\) be a Hermitian operator over \\({{{\\mathcal{H}}}}_{d}\\). There exists a polynomial Hamiltonian \\({P}_{\\hat{H}}(\\hat{q},\\hat{p})\\) of degree at most 3d over \\({{\\mathcal{H}}}\\) such that<\/p>\n<p>$${P}_{\\hat{H}}(\\hat{q},\\hat{p})=\\hat{H}\\oplus {\\hat{H}}^{{\\prime} },$$<\/p>\n<p>\n                    (5)\n                <\/p>\n<p>where \\({\\hat{H}}^{{\\prime} }\\) is a Hermitian operator over \\({{\\rm{span}}}{(\\left\\vert n\\right\\rangle )}_{n\\ &gt; \\ d}\\). In particular, \\({P}_{\\hat{H}}(\\hat{q},\\hat{p})\\) generates an evolution in \\({{\\mathcal{H}}}\\) given for all \\(\\left\\vert \\psi \\right\\rangle \\in {{{\\mathcal{H}}}}_{d}\\) by<\/p>\n<p>$${e}^{i{P}_{\\hat{H}}(\\hat{q},\\hat{p})}\\left\\vert \\psi \\right\\rangle={e}^{i\\hat{H}}\\left\\vert \\psi \\right\\rangle,$$<\/p>\n<p>\n                    (6)\n                <\/p>\n<p>Moreover, the polynomial \\({P}_{\\hat{H}}\\) can be computed in polynomial time in d.<\/p>\n<p>The proof of this result is based on the use of interpolation polynomials for reproducing exactly a target finite-dimensional operator on a subspace of the Fock basis, taking advantage of the sparsity of the canonical operators in that basis. We refer to the Methods section for the main techniques and to the\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">Supplementary Information<\/a> for the multimode generalisation of Theorem 2 and its proof. In particular, defining the unitary evolution operator \\({e}^{i{P}_{\\hat{H}}(\\hat{q},\\hat{p})}\\) on \\({{\\mathcal{H}}}\\) is highly non-trivial: we restrict its action to states in \\({{{\\mathcal{H}}}}_{d}\\) to avoid such mathematical technicalities.<\/p>\n<p>An important consequence of this theorem is that polynomial Hamiltonians allow us to explore the full Hilbert space, i.e., there are no inaccessible regions of the Hilbert space when starting from a fixed reference state, such as the vacuum state, and restricting to using gates generated by polynomial Hamiltonians:<\/p>\n<p>                  Corollary 1<\/p>\n<p>(Universal quantum controllability of polynomial Hamiltonians in infinite dimension). Let \\(\\left\\vert \\psi \\right\\rangle \\in {{\\mathcal{H}}}\\). For all \u03f5\u00a0&gt;\u00a00, there exists a polynomial Hamiltonian \\({P}_{\\epsilon }(\\hat{q},\\hat{p})\\) such that evolving the vacuum state \\(\\left\\vert 0\\right\\rangle\\) under the Schr\u00f6dinger equation with Hamiltonian \\({P}_{\\epsilon }(\\hat{q},\\hat{p})\\) for constant time yields a state \\(\\left\\vert {\\psi }_{\\epsilon }\\right\\rangle\\) satisfying \\(D(\\left\\vert \\psi \\right\\rangle,\\left\\vert {\\psi }_{\\epsilon }\\right\\rangle )\\le \\epsilon\\).<\/p>\n<p>Corollary 1 proves the universal quantum controllability of polynomial Hamiltonians in infinite dimensions, resolving the open question in<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 21\" title=\"Wu, R.-B., Tarn, T.-J. &amp; Li, C.-W. Smooth controllability of infinite-dimensional quantum-mechanical systems. Phys. Rev. A 73, 012719 (2006).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR21\" id=\"ref-link-section-d59843725e5204\" rel=\"nofollow noopener\" target=\"_blank\">21<\/a> and [ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 25\" title=\"Chabaud, U., Joseph, M., Mehraban, S. &amp; Motamedi, A. Bosonic Quantum Computational Complexity. In Meka, R. (ed.) 16th Innovations in Theoretical Computer Science Conference (ITCS 2025), 325 of Leibniz International Proceedings in Informatics (LIPIcs), 33:1&#x2013;33:19 (Schloss Dagstuhl &#x2013; Leibniz-Zentrum f&#xFC;r Informatik, Dagstuhl, Germany, 2025). arXiv:2410.04274.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR25\" id=\"ref-link-section-d59843725e5208\" rel=\"nofollow noopener\" target=\"_blank\">25<\/a>, Open question 3]. We give a quick proof in the single-mode case hereafter, the multimode case being analogous.<\/p>\n<p>                  Proof<\/p>\n<p>Since the state \\(\\left\\vert \\psi \\right\\rangle\\) is normalized, there exists d\u03f5 such it is \u03f5-close in trace distance to the state truncated at d\u03f5 and renormalized, which we denote by \\(\\left\\vert {\\psi }_{\\epsilon }\\right\\rangle\\). Then, there exists a finite-dimensional unitary operator \\({\\hat{U}}_{\\epsilon }={e}^{i{\\hat{H}}_{\\epsilon }}\\) over \\({{{\\mathcal{H}}}}_{{d}_{\\epsilon }}\\) which maps \\(\\left\\vert 0\\right\\rangle\\) to \\(\\left\\vert {\\psi }_{\\epsilon }\\right\\rangle\\). By Theorem 2, the polynomial Hamiltonial \\({P}_{\\epsilon }({\\hat{q}}, {\\hat{p}}) : \\!\\!={P}_{{\\hat{H}}_{\\epsilon }}({\\hat{q}},{\\hat{p}})\\) satisfies \\({P}_{\\epsilon }({\\hat{q}},{\\hat{p}})={\\hat{H}}_{\\epsilon }\\oplus {\\hat{H}}_{\\epsilon }^{{\\prime} }\\). Evolving the vacuum under the Schr\u00f6dinger equation with Hamiltonian \\({P}_{\\epsilon }(\\hat{q},\\hat{p})\\) thus leads to the same state as if evolving the vacuum under the (d\u03f5\u00a0+\u00a01)-dimensional Schr\u00f6dinger equation with Hamiltonian \\({\\hat{H}}_{\\epsilon }\\), which concludes the proof. \u25a1<\/p>\n<p>Beyond quantum state preparation, Theorem 2 provides a universal method for exactly reproducing any finite-dimensional Hamiltonian using polynomial Hamiltonians. This method can be used to engineer universal bosonic simulators capable of emulating the unitary evolution of any discrete-variable quantum system.<\/p>\n<p>Furthermore, Theorem 2 also implies that bosonic circuits composed of input vacuum state, unitary gates generated by polynomial Hamiltonians, and number measurements, can simulate universal qudit computations efficiently (without requiring feed-forward of measurement outcomes), generalizing a recent result for qubit computations [ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 25\" title=\"Chabaud, U., Joseph, M., Mehraban, S. &amp; Motamedi, A. Bosonic Quantum Computational Complexity. In Meka, R. (ed.) 16th Innovations in Theoretical Computer Science Conference (ITCS 2025), 325 of Leibniz International Proceedings in Informatics (LIPIcs), 33:1&#x2013;33:19 (Schloss Dagstuhl &#x2013; Leibniz-Zentrum f&#xFC;r Informatik, Dagstuhl, Germany, 2025). arXiv:2410.04274.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR25\" id=\"ref-link-section-d59843725e5863\" rel=\"nofollow noopener\" target=\"_blank\">25<\/a>, Theorem 4.1].<\/p>\n<p>Effective bosonic computations<\/p>\n<p>In this section, we explore the combined consequences of Theorem 1 (effective dimensions) and Theorem 2 (effective Hamiltonians) for the validity of effective descriptions of bosonic computations. At a fundamental level, both results show that effective descriptions are capturing bosonic computations well, i.e., one can always truncate the dimension of physical unitary evolutions up to an arbitrarily small error and use polynomial Hamiltonians to model finite-dimensional evolutions exactly. Together, these two results imply that polynomial Hamiltonians can generate any physical unitary evolution approximately to arbitrary precision:<\/p>\n<p>                  Theorem 3<\/p>\n<p>(Infinite-dimensional universality of polynomial Hamiltonians). Let \\({{\\mathcal{U}}}=\\hat{U}{{\\boldsymbol{\\cdot }}}{\\hat{U}}^{{\\dagger} }\\) be a single-mode unitary channel, with \\(\\hat{U}\\in {{\\mathcal{U}}}({{\\mathcal{S}}})\\) a physical unitary operator. Let E\u22650 be an energy parameter, and let \u03f5\u00a0&gt;\u00a00 be an approximation parameter. Then, there exists \\(d={{\\mathcal{O}}}\\left(\\frac{E}{{\\epsilon }^{4}}{E}_{\\hat{U}}\\left(\\frac{64E}{{\\epsilon }^{2}}\\right)\\right)\\in {\\mathbb{N}}\\) and a polynomial Hamiltonian \\(P(\\hat{q},\\hat{p})\\) of degree at most 3d such that, writing \\({{\\mathcal{V}}}\\) the corresponding unitary evolution induced on \\({{{\\mathcal{H}}}}_{d}\\),<\/p>\n<p>$$\\parallel {{\\mathcal{U}}}-{{\\mathcal{V}}}{\\parallel }_{\\diamond }^{E}\\le \\epsilon .$$<\/p>\n<p>\n                    (7)\n                <\/p>\n<p>Moreover, the polynomial P can be computed from \\(\\hat{U}\\) in polynomial time in d, E, and \\(\\frac{1}{\\epsilon }\\).<\/p>\n<p>As we explain in the Methods section, this result is a direct consequence of Theorem 1 and Theorem 2. Formally, Theorem 3 shows that the group of unitary gates generated by polynomial Hamiltonians is dense in the group of physical unitary operators in the strong operator topology (see the Methods section). This places the definition of universality for continuous-variable quantum gate sets, based on the ability to approximate to arbitrary precision any unitary evolution generated by a polynomial Hamiltonian<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 11\" title=\"Lloyd, S. &amp; Braunstein, S. L. Quantum computation over continuous variables. Phys. Rev. Lett. 82, 1784&#x2013;1787 (1999).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR11\" id=\"ref-link-section-d59843725e6344\" rel=\"nofollow noopener\" target=\"_blank\">11<\/a>, on an equal footing with its discrete-variable counterpart, based on the property to approximate to arbitrary precision any unitary operation<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 22\" title=\"Nielsen, M. A. &amp; Chuang, I. L. Quantum Computation and Quantum Information: 10th Anniversary Edition (Cambridge University Press, 2011), 10th edn.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR22\" id=\"ref-link-section-d59843725e6348\" rel=\"nofollow noopener\" target=\"_blank\">22<\/a>.<\/p>\n<p>As such, this shows that polynomial Hamiltonians generate truly universal gate sets and paves the way for a Solovay\u2013Kitaev theorem<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 23\" title=\"Dawson, C. M. &amp; Nielsen, M. A. The Solovay-Kitaev algorithm. arXiv preprint quant-ph\/0505030 (2005).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR23\" id=\"ref-link-section-d59843725e6355\" rel=\"nofollow noopener\" target=\"_blank\">23<\/a> in the infinite-dimensional setting. We obtain a version of this theorem hereafter, for universal gate sets generated by polynomial Hamiltonians as in Theorem 2:<\/p>\n<p>                  Theorem 4<\/p>\n<p>(Solovay\u2013Kitaev theorem for polynomial Hamiltonians). Let E\u00a0&gt;\u00a00, \u03f5\u00a0&gt;\u00a00 and \\(d\\ge \\frac{64E}{{\\epsilon }^{2}}\\in {\\mathbb{N}}\\). Let \\({{\\mathcal{G}}}\\) be a finite set of unitary operators over \\({{{\\mathcal{H}}}}_{d}\\) generating a dense subset of \\({{\\mathcal{U}}}({{{\\mathcal{H}}}}_{d})\\) and let \\({{\\mathcal{P}}}\\) be its realization with polynomial Hamiltonians from Theorem 2. There is a constant c such that for any physical unitary operator \\(\\hat{U}\\in {{\\mathcal{U}}}({{\\mathcal{S}}})\\) with (E,\u00a0\u03f5)-approximate effective dimension d\u00a0+\u00a01, there exists a finite sequence \\(\\hat{V}\\) of gates from \\({{\\mathcal{P}}}\\) of length \\({{\\mathcal{O}}}({\\log }^{c}(1\/\\epsilon ))\\) and such that \\(\\parallel {{\\mathcal{U}}}-{{\\mathcal{V}}}{\\parallel }_{\\diamond }^{E}\\le 2\\epsilon\\), where \\({{\\mathcal{U}}}=\\hat{U}{{\\boldsymbol{\\cdot }}}{\\hat{U}}^{{\\dagger} }\\) and \\({{\\mathcal{V}}}=\\hat{V}{{\\boldsymbol{\\cdot }}}{\\hat{V}}^{{\\dagger} }\\).<\/p>\n<p>Recall that Theorem 1 ensures the existence of an (E,\u00a0\u03f5)-approximate effective dimension for all physical unitary operators. The proof of Theorem 4 combines Theorems 1 and 2 with the Solovay\u2013Kitaev theorem for qudits<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 23\" title=\"Dawson, C. M. &amp; Nielsen, M. A. The Solovay-Kitaev algorithm. arXiv preprint quant-ph\/0505030 (2005).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR23\" id=\"ref-link-section-d59843725e6883\" rel=\"nofollow noopener\" target=\"_blank\">23<\/a> of dimension d\u00a0+\u00a01, and provides an explicit algorithm for producing the sequence of polynomial Hamiltonians, given the target unitary \\(\\hat{U}\\in {{\\mathcal{U}}}({{\\mathcal{S}}})\\). We refer to the Methods section for the main techniques and to the\u00a0<a data-track=\"click\" data-track-label=\"link\" data-track-action=\"supplementary material anchor\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#MOESM1\" rel=\"nofollow noopener\" target=\"_blank\">Supplementary Information<\/a> for the proof.<\/p>\n<p>Theorem 4 generalises the Gaussian Solovay\u2013Kitaev theorem from ref. <a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 24\" title=\"Becker, S., Datta, N., Lami, L. &amp; Rouz&#xE9;, C. Energy-constrained discrimination of unitaries, quantum speed limits, and a Gaussian Solovay-Kitaev theorem. Phys. Rev. Lett. 126, 190504 (2021).\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR24\" id=\"ref-link-section-d59843725e6944\" rel=\"nofollow noopener\" target=\"_blank\">24<\/a>, thus providing an answer to [<a data-track=\"click\" data-track-action=\"reference anchor\" data-track-label=\"link\" data-test=\"citation-ref\" aria-label=\"Reference 25\" title=\"Chabaud, U., Joseph, M., Mehraban, S. &amp; Motamedi, A. Bosonic Quantum Computational Complexity. In Meka, R. (ed.) 16th Innovations in Theoretical Computer Science Conference (ITCS 2025), 325 of Leibniz International Proceedings in Informatics (LIPIcs), 33:1&#x2013;33:19 (Schloss Dagstuhl &#x2013; Leibniz-Zentrum f&#xFC;r Informatik, Dagstuhl, Germany, 2025). arXiv:2410.04274.\" href=\"http:\/\/www.nature.com\/articles\/s41467-025-64872-3#ref-CR25\" id=\"ref-link-section-d59843725e6948\" rel=\"nofollow noopener\" target=\"_blank\">25<\/a>, Open question 4] for the gate sets based on the polynomial Hamiltonians used in Theorem 2. In particular, it implies the computational equivalence of bosonic quantum circuits based on these universal gate sets.<\/p>\n","protected":false},"excerpt":{"rendered":"Bosonic quantum information theory Quantum information theory with bosonic systems typically involves modeling unbounded operators over infinite-dimensional Hilbert&hellip;\n","protected":false},"author":3,"featured_media":361054,"comment_status":"","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[25],"tags":[26265,10046,10047,492,13632,5649,18678,159,67,132,68],"class_list":{"0":"post-361053","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-physics","8":"tag-computer-science","9":"tag-humanities-and-social-sciences","10":"tag-multidisciplinary","11":"tag-physics","12":"tag-quantum-information","13":"tag-quantum-optics","14":"tag-quantum-simulation","15":"tag-science","16":"tag-united-states","17":"tag-unitedstates","18":"tag-us"},"share_on_mastodon":{"url":"https:\/\/pubeurope.com\/@us\/115505772958272022","error":""},"_links":{"self":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/361053","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/comments?post=361053"}],"version-history":[{"count":0,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/posts\/361053\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media\/361054"}],"wp:attachment":[{"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/media?parent=361053"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/categories?post=361053"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.europesays.com\/us\/wp-json\/wp\/v2\/tags?post=361053"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}