Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2112.09005 · Physical Review A

Universe as a nonlinear quantum simulation: Large- n limit of the central-spin model

Michael R. Geller

Published Oct 16, 2023• 1 citation
arXivHTMLPDF

Abstract

We investigate models of nonlinear qubit evolution based on mappings to an $n$-qubit central-spin model (CSM) in the large-$n$ limit, where mean-field theory is exact. Extending a theorem of Erd\"os and Schlein [J. Stat. Phys. 134, 859 (2009)], we establish that the CSM is rigorously dual to a nonlinear qubit when $n\ensuremath{\rightarrow}\ensuremath{\infty}$. The duality supports a type of nonlinear quantum computation in systems, such as a condensate, where a large number of ancillas couple symmetrically to a central qubit. It also enables a gate-model implementation of nonlinear quantum simulation with a rigorous error bound. Two variants of the model, with and without coupling between ancillas, map to effective models with different nonlinearity and symmetry. Without coupling the CSM simulates initial-condition nonlinearity, where the Hamiltonian is a linear combination of $\mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{x}){\ensuremath{\sigma}}^{x}, \mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{y}){\ensuremath{\sigma}}^{y}$, and $\mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{z}){\ensuremath{\sigma}}^{z}$, where ${\ensuremath{\sigma}}^{x}, {\ensuremath{\sigma}}^{y}$, and ${\ensuremath{\sigma}}^{z}$ are Pauli matrices and ${\ensuremath{\rho}}_{0}$ is the initial density matrix. With symmetric ancillas coupling it simulates linear combinations of $\mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{x}){\ensuremath{\sigma}}^{x}, \mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{y}){\ensuremath{\sigma}}^{y}$, and $\mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{z}){\ensuremath{\sigma}}^{z}$, where $\ensuremath{\rho}$ is the current state. This case can simulate qubit torsion, which has been shown by Abrams and Lloyd [Phys. Rev. Lett. 81, 3992 (1998)] to enable an exponential speedup for state discrimination in an idealized setting. The duality discussed here might also be interesting from a quantum foundation perspective. There has long been interest in whether quantum mechanics might possess some type of small unobserved nonlinearity. If not, what is the principle prohibiting it? The duality implies that there is not a sharp distinction between universes evolving according to linear and nonlinear quantum mechanics: A one-qubit universe prepared in a pure state $|\ensuremath{\varphi}\ensuremath{\rangle}$ at the time of the big bang and symmetrically coupled to ancillas prepared in the same state would appear to evolve nonlinearly for any finite time $t>0$ as long as there are exponentially many ancillas $n\ensuremath{\gg}\mathrm{exp}[O(t)]$.

DOI
10.1103/physreva.108.042210

Related papers

Nonlinear Hamiltonians and Boolean satisfiability

Michael R. Geller, Victoria S. Ordonez, Yohannes Abate

We consider an extended model of quantum computation where a scalable fault-tolerant quantum computer is coupled to one or more ancilla qubits that evolve according to a nonlinear Schrödinger equation. Following the approach of Abrams and Lloyd, an efficient quantum circuit evaluating an $n$-bit Boolean function in conjunctive normal form is used to prepare an ancilla encoding its number $s$ of satisfying assignments ($0 \le s \le 2^n$). This is followed by a nonlinear quantum state discrimination gate on the ancilla qubit that is used to learn properties of $s$. Here we consider three types of state discriminators generated by different nonlinear Hamiltonians. First, given a restricted Boolean satisfiability problem with the promise of at most one satisfying assignment ($ 0 \le s \le 1$), we show that a qubit with $\langle σ^z \rangle σ^z$ nonlinearity can be used to efficiently determine whether $s = 0$ or $s = 1$, solving the UNIQUE SAT problem. Here $\langle A \rangle := \langle ψ| A |ψ\rangle $ denotes expectation in the current state. UNIQUE SAT is NP-hard under a randomized polynomial-time reduction (of course any discussion of complexity assumes a scalable, fault-tolerant implementation). Second, for unrestricted satisfiability problems with $ 0 \le s \le 2^n$, a Hamiltonian with $ \langle σ^x \rangle σ^y - \langle σ^y \rangle σ^x$ nonlinearity can be used to efficiently determine whether $s=0$ or $s>0$, thereby solving 3SAT, which is NP-complete. Finally, we show that $ \langle σ^y \rangle \langle σ^z \rangle σ^x - \langle σ^x \rangle \langle σ^z \rangle σ^y $ nonlinearity can be used to efficiently measure $s$ and solve #SAT, which is #P-complete. The nonlinear models are of mean field type and might be simulated with ultracold atoms.

quant-pharXiv:2605.14822•May 15, 2026•0 citations

Quantum mechanical treatment of large spin baths

Robin Röhrig, Philipp Schering, Lars B. Gravert +2 more

The electronic spin in quantum dots can be described by central spin models (CSMs) with a very large number ${N}_{\text{eff}}\ensuremath{\approx}{10}^{4}$ to ${10}^{6}$ of bath spins posing a tremendous challenge to theoretical simulations. Here, a fully quantum mechanical theory is developed for the limit ${N}_{\text{eff}}\ensuremath{\rightarrow}\ensuremath{\infty}$ by means of iterated equations of motion (iEoM). We find that the CSM can be mapped to a four-dimensional impurity coupled to a noninteracting bosonic bath in this limit. Remarkably, even for infinite bath the CSM does not become completely classical. The data obtained by the proposed iEoM approach are tested successfully against data from other, established approaches. Thus the iEoM mapping extends the set of theoretical tools that can be used to understand the spin dynamics in large CSMs.

arXiv:1711.08919•Physical review. B./Physical review. B•Apr 23, 2018•9 citations

Emergent Non-Markovian Nonlinear Qubit From Collective Spin Interactions

Gregory T. Carroll, Michael R. Geller, Andre Erpenbeck

Open-system descriptions are typically introduced by coupling a quantum system to an external environment. Here we show that a closed interacting many-body system can itself generate a controlled non-Markovian quantum channel acting on a reduced nonlinear qubit through finite-size corrections to a nonlinear mean-field limit. We demonstrate this using the Kitagawa-Ueda one-axis twisting model, $H=χJ_z^2$, a paradigmatic model of collective spin dynamics, spin squeezing, and two-component Bose-Einstein condensates. Although the large-$N$ regime of this model has been extensively studied, the conventional fixed-$χ$ scaling does not yield a nontrivial dynamical large-$N$ limit. In this paper, we investigate a complementary large-$N$ formulation obtained from the double limit $N\rightarrow\infty$ and $χ\rightarrow O(g/N)$, where $g$ is a coupling constant. We derive the leading finite-$N$ corrections to this limit and show that they correspond to an emergent non-Markovian dephasing process, producing a Gaussian decay of the Bloch-vector coherence with characteristic timescale $t_\varphi\geq\sqrt{N}/(2g)$. Exact finite-$N$ calculations demonstrate that this effective open-system description becomes quantitatively accurate for systems containing on the order of one hundred qubits. The resulting framework provides a microscopic realization of non-Markovian dephasing generated intrinsically by a closed many-body system and enables efficient simulation of collective quantum dynamics beyond unitary mean-field theory. These results link the long-studied phenomenon of phase diffusion in atomic ensembles and Bose-Einstein condensates to the growing effort to characterize non-Markovian, beyond-Lindblad noise in quantum computing hardware, providing a rare case in which such a noise channel is derived from microscopic dynamics rather than fit phenomenologically.

quant-pharXiv:2608.07723•Aug 11, 2026•0 citations

Quantum algorithm for Valiant-Vazirani reduction

Patrick Kelly, Victoria S. Ordonez, Michael R. Geller +1 more

There is growing interest in extensions of the standard model of gate-based quantum computation to include auxiliary degrees of freedom evolving according to a nonlinear Schrödinger equation. By reducing the Boolean satisfiability problem SAT to quantum state discrimination, Abrams and Lloyd argued that the right type of nonlinearity can be used to solve NP and #P problems in polynomial time, at least in an idealized noise-free limit. For practical implementation, however, we are restricted to simulated and emergent nonlinearities, such as that appearing in mean field models for ultracold atoms and similar ensembles. A prominent example is the torsion model, which arises in two-component Bose-Einstein condensates and spin models with all-to-all Ising interaction. But torsion-based state discrimination appears to fall short of solving SAT. Here we close this gap by constructing the filtered oracle of the Valiant-Vazirani theorem, providing a randomized polynomial-time reduction from SAT to UNIQUE SAT, a promise problem where there is at most 1 satisfying assignment. In the noise-free limit, the UNIQUE SAT problem can be solved in polynomial time using torsion nonlinearity. Quantum Valiant-Vazirani reduction is no faster than the efficient classical version, but a fault-tolerant implementation coupled to a nonlinear quantum coprocessor simulating torsion would enable polynomial time solution to NP (but not #P) problems.

quant-pharXiv:2606.18428•Jun 18, 2026•0 citations

Quantum simulation from the bottom up: the case of rebits

Dax Enshan Koh, Murphy Yuezhen Niu, Theodore J. Yoder

Abstract Typically, quantum mechanics is thought of as a linear theory with unitary evolution governed by the Schrödinger equation. While this is technically true and useful for a physicist, with regards to computation it is an unfortunately narrow point of view. Just as a classical computer can simulate highly nonlinear functions of classical states, so too can the more general quantum computer simulate nonlinear evolutions of quantum states. We detail one particular simulation of nonlinearity on a quantum computer, showing how the entire class of -unitary evolutions (on n qubits) can be simulated using a unitary, real-amplitude quantum computer (consisting of n + 1 qubits in total). These operators can be represented as the sum of a linear and antilinear operator, and add an intriguing new set of nonlinear quantum gates to the toolbox of the quantum algorithm designer. Furthermore, a subgroup of these nonlinear evolutions, called the -Cliffords , can be efficiently classically simulated, by making use of the fact that Clifford operators can simulate non-Clifford (in fact, non-linear) operators. This perspective of using the physical operators that we have to simulate non-physical ones that we do not is what we call bottom-up simulation, and we give some examples of its broader implications.

arXiv:1708.09355•Journal of Physics A Mathematical and Theoretical•Mar 27, 2018•8 citations

Dynamics of the homogeneous two-qubit X X Z central spin model with the spin bath prepared in superpositions of symmetric Dicke states

Zejiang Li, Pei Yang, Wen‐Long You +1 more

We obtain exact dynamics of a two-qubit central spin model (CSM) consisting of two interacting qubits homogeneously coupled to a spin bath via the $XXZ$-type coupling, with the bath initially prepared in linear superpositions of the symmetric Dicke states. Using the interaction picture Hamiltonian with respect to the non-spin-flipping part of the model, we derive a sequence of equations of motion within each magnetization sector satisfied by the probability amplitudes of the time-evolved state. These equations of motion admit analytical solutions for the single-qubit CSM in which one of the two central qubits decouples from the rest of the system. Based on this, we provide a quantitative interpretation to the observed collapse-revival phenomena in the single-qubit Rabi oscillations when the bath is prepared in the spin coherent state. We then study the disentanglement and coherence dynamics of two initially entangled noninteracting qubits when the two qubits interact with individual baths or with a common bath. For individual baths, the coherent dynamics is found to be positively correlated to the single-qubit purity dynamics, and entanglement sudden disappearance and revivals are observed in both cases. The entanglement creation of two initially separable qubits coupled to a common bath is also studied and collapse and revival behaviors in the entanglement dynamics are observed. Choosing the equally weighted state and the $W$-class states as the bath initial states, we finally study the dynamics of entanglement between two individual bath spins and demonstrate the entanglement sharing mechanism in such a system.

arXiv:2006.00550•Physical Review A•Sep 16, 2020•11 citations
We investigate models of nonlinear qubit evolution based on mappings to an nnn-qubit central-spin model (CSM) in the large-nnn limit, where mean-field theory is exact. Extending a theorem of Erd"os and Schlein [J. Stat. Phys. 134, 859 (2009)], we establish that the CSM is rigorously dual to a nonlinear qubit when n\ensuremath→\ensuremath∞n\ensuremath{\rightarrow}\ensuremath{\infty}n\ensuremath→\ensuremath∞. The duality supports a type of nonlinear quantum computation in systems, such as a condensate, where a large number of ancillas couple symmetrically to a central qubit. It also enables a gate-model implementation of nonlinear quantum simulation with a rigorous error bound. Two variants of the model, with and without coupling between ancillas, map to effective models with different nonlinearity and symmetry. Without coupling the CSM simulates initial-condition nonlinearity, where the Hamiltonian is a linear combination of tr(\ensuremathρ0\ensuremathσx)\ensuremathσx,tr(\ensuremathρ0\ensuremathσy)\ensuremathσy\mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{x}){\ensuremath{\sigma}}^{x}, \mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{y}){\ensuremath{\sigma}}^{y}tr(\ensuremathρ0​\ensuremathσx)\ensuremathσx,tr(\ensuremathρ0​\ensuremathσy)\ensuremathσy, and tr(\ensuremathρ0\ensuremathσz)\ensuremathσz\mathrm{tr}({\ensuremath{\rho}}_{0}{\ensuremath{\sigma}}^{z}){\ensuremath{\sigma}}^{z}tr(\ensuremathρ0​\ensuremathσz)\ensuremathσz, where \ensuremathσx,\ensuremathσy{\ensuremath{\sigma}}^{x}, {\ensuremath{\sigma}}^{y}\ensuremathσx,\ensuremathσy, and \ensuremathσz{\ensuremath{\sigma}}^{z}\ensuremathσz are Pauli matrices and \ensuremathρ0{\ensuremath{\rho}}_{0}\ensuremathρ0​ is the initial density matrix. With symmetric ancillas coupling it simulates linear combinations of tr(\ensuremathρ\ensuremathσx)\ensuremathσx,tr(\ensuremathρ\ensuremathσy)\ensuremathσy\mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{x}){\ensuremath{\sigma}}^{x}, \mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{y}){\ensuremath{\sigma}}^{y}tr(\ensuremathρ\ensuremathσx)\ensuremathσx,tr(\ensuremathρ\ensuremathσy)\ensuremathσy, and tr(\ensuremathρ\ensuremathσz)\ensuremathσz\mathrm{tr}(\ensuremath{\rho}{\ensuremath{\sigma}}^{z}){\ensuremath{\sigma}}^{z}tr(\ensuremathρ\ensuremathσz)\ensuremathσz, where \ensuremathρ\ensuremath{\rho}\ensuremathρ is the current state. This case can simulate qubit torsion, which has been shown by Abrams and Lloyd [Phys. Rev. Lett. 81, 3992 (1998)] to enable an exponential speedup for state discrimination in an idealized setting. The duality discussed here might also be interesting from a quantum foundation perspective. There has long been interest in whether quantum mechanics might possess some type of small unobserved nonlinearity. If not, what is the principle prohibiting it? The duality implies that there is not a sharp distinction between universes evolving according to linear and nonlinear quantum mechanics: A one-qubit universe prepared in a pure state ∣\ensuremathφ\ensuremath⟩|\ensuremath{\varphi}\ensuremath{\rangle}∣\ensuremathφ\ensuremath⟩ at the time of the big bang and symmetrically coupled to ancillas prepared in the same state would appear to evolve nonlinearly for any finite time t>0t>0t>0 as long as there are exponentially many ancillas n\ensuremath≫exp[O(t)]n\ensuremath{\gg}\mathrm{exp}[O(t)]n\ensuremath≫exp[O(t)].
We consider an extended model of quantum computation where a scalable fault-tolerant quantum computer is coupled to one or more ancilla qubits that evolve according to a nonlinear Schrödinger equation. Following the approach of Abrams and Lloyd, an efficient quantum circuit evaluating an nnn-bit Boolean function in conjunctive normal form is used to prepare an ancilla encoding its number sss of satisfying assignments (0≤s≤2n0 \le s \le 2^n0≤s≤2n). This is followed by a nonlinear quantum state discrimination gate on the ancilla qubit that is used to learn properties of sss. Here we consider three types of state discriminators generated by different nonlinear Hamiltonians. First, given a restricted Boolean satisfiability problem with the promise of at most one satisfying assignment ($ 0 \le s \le 1$), we show that a qubit with ⟨σz⟩σz\langle σ^z \rangle σ^z⟨σz⟩σz nonlinearity can be used to efficiently determine whether s=0s = 0s=0 or s=1s = 1s=1, solving the UNIQUE SAT problem. Here $\langle A \rangle := \langle ψ| A |ψ\rangle $ denotes expectation in the current state. UNIQUE SAT is NP-hard under a randomized polynomial-time reduction (of course any discussion of complexity assumes a scalable, fault-tolerant implementation). Second, for unrestricted satisfiability problems with $ 0 \le s \le 2^n$, a Hamiltonian with $ \langle σ^x \rangle σ^y - \langle σ^y \rangle σ^x$ nonlinearity can be used to efficiently determine whether s=0s=0s=0 or s>0s>0s>0, thereby solving 3SAT, which is NP-complete. Finally, we show that $ \langle σ^y \rangle \langle σ^z \rangle σ^x - \langle σ^x \rangle \langle σ^z \rangle σ^y $ nonlinearity can be used to efficiently measure sss and solve #SAT, which is #P-complete. The nonlinear models are of mean field type and might be simulated with ultracold atoms.
The electronic spin in quantum dots can be described by central spin models (CSMs) with a very large number ${N}_{\text{eff}}\ensuremath{\approx}{10}^{4}$ to 106{10}^{6}106 of bath spins posing a tremendous challenge to theoretical simulations. Here, a fully quantum mechanical theory is developed for the limit ${N}_{\text{eff}}\ensuremath{\rightarrow}\ensuremath{\infty}$ by means of iterated equations of motion (iEoM). We find that the CSM can be mapped to a four-dimensional impurity coupled to a noninteracting bosonic bath in this limit. Remarkably, even for infinite bath the CSM does not become completely classical. The data obtained by the proposed iEoM approach are tested successfully against data from other, established approaches. Thus the iEoM mapping extends the set of theoretical tools that can be used to understand the spin dynamics in large CSMs.
Open-system descriptions are typically introduced by coupling a quantum system to an external environment. Here we show that a closed interacting many-body system can itself generate a controlled non-Markovian quantum channel acting on a reduced nonlinear qubit through finite-size corrections to a nonlinear mean-field limit. We demonstrate this using the Kitagawa-Ueda one-axis twisting model, H=χJz2H=χJ_z^2H=χJz2​, a paradigmatic model of collective spin dynamics, spin squeezing, and two-component Bose-Einstein condensates. Although the large-NNN regime of this model has been extensively studied, the conventional fixed-χχχ scaling does not yield a nontrivial dynamical large-NNN limit. In this paper, we investigate a complementary large-NNN formulation obtained from the double limit N→∞N\rightarrow\inftyN→∞ and χ→O(g/N)χ\rightarrow O(g/N)χ→O(g/N), where ggg is a coupling constant. We derive the leading finite-NNN corrections to this limit and show that they correspond to an emergent non-Markovian dephasing process, producing a Gaussian decay of the Bloch-vector coherence with characteristic timescale tφ≥N/(2g)t_\varphi\geq\sqrt{N}/(2g)tφ​≥N​/(2g). Exact finite-NNN calculations demonstrate that this effective open-system description becomes quantitatively accurate for systems containing on the order of one hundred qubits. The resulting framework provides a microscopic realization of non-Markovian dephasing generated intrinsically by a closed many-body system and enables efficient simulation of collective quantum dynamics beyond unitary mean-field theory. These results link the long-studied phenomenon of phase diffusion in atomic ensembles and Bose-Einstein condensates to the growing effort to characterize non-Markovian, beyond-Lindblad noise in quantum computing hardware, providing a rare case in which such a noise channel is derived from microscopic dynamics rather than fit phenomenologically.
There is growing interest in extensions of the standard model of gate-based quantum computation to include auxiliary degrees of freedom evolving according to a nonlinear Schrödinger equation. By reducing the Boolean satisfiability problem SAT to quantum state discrimination, Abrams and Lloyd argued that the right type of nonlinearity can be used to solve NP and #P problems in polynomial time, at least in an idealized noise-free limit. For practical implementation, however, we are restricted to simulated and emergent nonlinearities, such as that appearing in mean field models for ultracold atoms and similar ensembles. A prominent example is the torsion model, which arises in two-component Bose-Einstein condensates and spin models with all-to-all Ising interaction. But torsion-based state discrimination appears to fall short of solving SAT. Here we close this gap by constructing the filtered oracle of the Valiant-Vazirani theorem, providing a randomized polynomial-time reduction from SAT to UNIQUE SAT, a promise problem where there is at most 1 satisfying assignment. In the noise-free limit, the UNIQUE SAT problem can be solved in polynomial time using torsion nonlinearity. Quantum Valiant-Vazirani reduction is no faster than the efficient classical version, but a fault-tolerant implementation coupled to a nonlinear quantum coprocessor simulating torsion would enable polynomial time solution to NP (but not #P) problems.
Abstract Typically, quantum mechanics is thought of as a linear theory with unitary evolution governed by the Schrödinger equation. While this is technically true and useful for a physicist, with regards to computation it is an unfortunately narrow point of view. Just as a classical computer can simulate highly nonlinear functions of classical states, so too can the more general quantum computer simulate nonlinear evolutions of quantum states. We detail one particular simulation of nonlinearity on a quantum computer, showing how the entire class of -unitary evolutions (on n qubits) can be simulated using a unitary, real-amplitude quantum computer (consisting of n + 1 qubits in total). These operators can be represented as the sum of a linear and antilinear operator, and add an intriguing new set of nonlinear quantum gates to the toolbox of the quantum algorithm designer. Furthermore, a subgroup of these nonlinear evolutions, called the -Cliffords , can be efficiently classically simulated, by making use of the fact that Clifford operators can simulate non-Clifford (in fact, non-linear) operators. This perspective of using the physical operators that we have to simulate non-physical ones that we do not is what we call bottom-up simulation, and we give some examples of its broader implications.
We obtain exact dynamics of a two-qubit central spin model (CSM) consisting of two interacting qubits homogeneously coupled to a spin bath via the XXZXXZXXZ-type coupling, with the bath initially prepared in linear superpositions of the symmetric Dicke states. Using the interaction picture Hamiltonian with respect to the non-spin-flipping part of the model, we derive a sequence of equations of motion within each magnetization sector satisfied by the probability amplitudes of the time-evolved state. These equations of motion admit analytical solutions for the single-qubit CSM in which one of the two central qubits decouples from the rest of the system. Based on this, we provide a quantitative interpretation to the observed collapse-revival phenomena in the single-qubit Rabi oscillations when the bath is prepared in the spin coherent state. We then study the disentanglement and coherence dynamics of two initially entangled noninteracting qubits when the two qubits interact with individual baths or with a common bath. For individual baths, the coherent dynamics is found to be positively correlated to the single-qubit purity dynamics, and entanglement sudden disappearance and revivals are observed in both cases. The entanglement creation of two initially separable qubits coupled to a common bath is also studied and collapse and revival behaviors in the entanglement dynamics are observed. Choosing the equally weighted state and the WWW-class states as the bath initial states, we finally study the dynamics of entanglement between two individual bath spins and demonstrate the entanglement sharing mechanism in such a system.