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:2506.05757 · Nature Communications

Pathfinding quantum simulations of neutrinoless double-β decay

Ivan A. Chernyshev, Roland C. Farrell, Marc Illa, Martin J. Savage A, Andrii Maksymov, Felix Tripier, Miguel Angel Lopez-Ruiz, Andrew Arrasmith, Yvette de Sereville, Aharon Brodutch, Claudio Girotto, Ananth P. Kaushik, Martin Roetteler

Published Jan 23, 2026• 2 citations
arXivHTMLPDF

Abstract

We present results from co-designed quantum simulations of the neutrinoless double-β decay of a simple nucleus in 1+1D quantum chromodynamics using IonQ's Forte-generation trapped-ion quantum computers. Electrons, neutrinos, and up and down quarks are distributed across two lattice sites and mapped to 32 qubits, with an additional 4 qubits used for flag-based error mitigation. A four-fermion interaction is used to implement weak interactions, and lepton-number violation is induced by a neutrino Majorana mass. Quantum circuits that prepare the initial nucleus and time evolve with the Hamiltonian containing the strong and weak interactions are executed on IonQ Forte Enterprise. Enabled by tuned model parameters, lepton-number violation is observed in real time, providing a clear signal of neutrinoless double-β decay. This was made possible by co-designing the simulation to maximally utilize the all-to-all connectivity and native gate-set available on IonQ's quantum computers. Quantum circuit compilation techniques and co-designed error-mitigation methods, informed from executing benchmarking circuits with up to 2,356 two-qubit gates, enabled observables to be extracted with high precision. We discuss the potential of future quantum simulations to provide yocto-second resolution of the reaction pathways in these, and other, nuclear processes.

DOI
10.1038/s41467-026-68536-8

Related papers

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. II. Single-baryon β -decay in real time

Roland C. Farrell, Ivan A. Chernyshev, Sarah J. M. Powell +3 more

A framework for quantum simulations of real-time weak decays of hadrons and nuclei in a two-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in the dynamics, the latter through four-Fermi effective operators. Quantum circuits that implement time evolution in this lattice theory are developed and run on Quantinuum's H1-1 20-qubit trapped ion system to simulate the $\ensuremath{\beta}$-decay of a single baryon on one lattice site. These simulations include the initial state preparation and are performed for both one and two Trotter time steps. The potential intrinsic error-correction properties of this type of lattice theory are discussed, and the leading lattice Hamiltonian required to simulate $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$-decay of nuclei induced by a neutrino Majorana mass term is provided.

arXiv:2209.10781•Physical review. D/Physical review. D.•Mar 30, 2023•76 citations

Trapped-ion quantum simulation of the Fermi-Hubbard model as a lattice gauge theory using hardware-aware native gates

Dhruv Srinivasan, Beyer, Alex, Daiwei Zhu +8 more

The Fermi-Hubbard model (FHM) is a simple yet rich model of strongly interacting electrons with complex dynamics and a variety of emerging quantum phases. These properties make it a compelling target for digital quantum simulation. Trotterization-based quantum simulations have shown promise, but implementations on current hardware are limited by noise, necessitating error mitigation techniques like circuit optimization and post-selection. A mapping of the FHM to a Z2 LGT was recently proposed that restricts the dynamics to a subspace protected by additional symmetries, and its ability for post-selection error mitigation was verified through noisy classical simulations. In this work, we propose and demonstrate a suite of algorithm-hardware co-design strategies on a trapped-ion quantum computer, targeting two key aspects of NISQ-era quantum simulation: circuit compilation and error mitigation. In particular, a novel combination of iteratively preconditioned gradient descent (IPG) and subsystem von Neumann Entropy compression reduces the 2-qubit gate count of FHM quantum simulation by 35%, consequently doubling the number of simulatable Trotter steps when used in tandem with error mitigation based on conserved symmetries, debiasing and sharpening techniques. Our work demonstrates the value of algorithm-hardware co-design to operate digital quantum simulators at the threshold of maximum circuit depths allowed by current hardware, and is broadly generalizable to strongly correlated systems in quantum chemistry and materials science.

arXiv:2411.07778•Nov 12, 2024•0 citations

Large circuit execution for NMR spectroscopy simulation on NISQ quantum hardware

Artemiy Burov, Julien Baglio, Clément Javerzac-Galy

With the latest advances in quantum computing technology, we are gradually moving from the noisy intermediate-scale quantum (NISQ) era characterized by hardware limited in the number of qubits and plagued with quantum noise, to the age of quantum utility where both the newest hardware and software methods allow for tackling problems which have been deemed difficult or intractable with conventional classical methods. One of these difficult problems is the simulation of one-dimensional (1D) nuclear magnetic resonance (NMR) spectra, a major tool to learn about the structure of molecules, helping the design of new materials or drugs. Using advanced error mitigation and error suppression techniques from Q-CTRL together with the latest commercially available superconducting-qubit quantum computer from IBM and trapped-ion quantum computer from IonQ, we present the quantum Hamiltonian simulation of liquid-state 1D NMR spectra in the high-field regime for spin systems up to 34 spins. Our pipeline has a major impact on the ability to execute deep quantum circuits with the reduction of quantum noise, improving mean square error by a factor of 22. It allows for the execution of deep quantum circuits and obtaining salient features of the 1D NMR spectra for both 16-spin and 22-spin systems, as well as a 34-spin system, which lies beyond the regime where unrestricted full Liouvillespace simulations are practical (32 spins, the Liouville limit). Our work is a step toward near-term quantum utility in NMR spectroscopy.

quant-pharXiv:2512.14513•Dec 17, 2025•0 citations

Developing techniques for Simulation of SU(3) Quantum Field Theories on State-of-the-Art Quantum Devices

Chernyshev, Ivan

Quantum computing has long been an experimental technology with the potential to simulate, at scale, phenomena which on classical devices would be too expensive to simulate at any but the smallest scales. Over the last several years, however, it has entered the NISQ era, where the number of qubits are sufficient for quantum advantage but substantial noise on hardware stands in the way of this achievement. This thesis details NISQ device-centered improvements to techniques of quantum simulation of the out-of-equilbrium real-time dynamics of lattice quantum chromodynamics (LQCD) and of dense 3-flavor neutrino systems on digital quantum devices. The first project concerning LQCD is a comparison of methods for implementing the variational quantum eigensolver (VQE) that initializes the ground state of an SU(3) plaquette-chain. The thesis then pivots to a 1+1D lattice of quarks interacting with an SU(3) gauge-field. A VQE-based state-preparation for the vacua and a Trotterized time-evolution circuit is designed and applied to the problems of simulating beta and neutrinoless double beta decay. Finally, these circuits are adapted to a version useable on quantum devices with nearest-neighbor connectivity with minimal overhead, with an eye towards utilizing the higher qubit count of such devices for hadron dynamics and scattering. This thesis covers two projects that concern dense 3-flavor neutrino systems. The first details design and testing of Trotterized time-evolution circuits on state-of-the-art quantum devices. The second, motivated by the Gottesman-Knill theorem's result that deviation from stabilizer states ("magic") is necessary for a problem to exhibit quantum advantage, details results with implications for the Standard Model in general that the 3 flavor ultradense neutrino systems with the highest, most-persistent magic are those that start with neutrinos in all 3 flavors.

arXiv:2502.02502•Feb 4, 2025•0 citations

Trapped-ion quantum simulation of collective neutrino oscillations

Valentina Amitrano, Alessandro Roggero, Piero Luchi +3 more

It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of $N$ neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. We finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.

arXiv:2207.03189•Physical review. D/Physical review. D.•Jan 13, 2023•41 citations

Quantum Co-Design of Inhomogeneous Many-Body Neutrino Fast Flavor Transformation

Zoha Laraib, Sherwood Richers, Alessandro Baroni +2 more

Dense-neutrino flavor evolution is a quantum many-body problem whose fully correlated treatment becomes rapidly more expensive with increasing particle number and spatial structure. We develop a \texttt{QCNO} quantum simulation code that maps a inhomogeneous forward-scattering neutrino Hamiltonian with advection and finite-range interactions to TEBD2 product-formula quantum circuits, and connects ideal many-body simulation, backend-aware execution, and fault-tolerant resource estimation within the same physical model. We reproduce the exact many-body evolution of the suppressed mean-field-like transverse fast flavor instability through $N=30$ and show that even a single open-boundary interaction generates Rënyi entanglement and non-stabilizer magic, although noisy backends still produce polarization RMSEs of order $0.1$--$0.3$. The restricted active interaction graph of this problem allows us to estimate the circuit depth and $T$ gate cost through $N=50$ for both NISQ and fault-tolerant approaches. The measured TEBD2 error in our fiducial ideal simulation of order $10^{-3}$ motivates a synthesis tolerance $\varepsilon_{\rm syn}\sim10^{-7}$, corresponding to about $70$ $T$ gates per $R_Z$ rotation. The generated circuit contains $1.4\times10^4$ $T$ gates per qubit for open boundaries (twice that for closed boundaries), requiring an application-level logical-$T$ error target of order $10^{-7}$ for early fault-tolerant architectures.

quant-pharXiv:2610.12334•Oct 9, 2026•0 citations
We present results from co-designed quantum simulations of the neutrinoless double-β decay of a simple nucleus in 1+1D quantum chromodynamics using IonQ's Forte-generation trapped-ion quantum computers. Electrons, neutrinos, and up and down quarks are distributed across two lattice sites and mapped to 32 qubits, with an additional 4 qubits used for flag-based error mitigation. A four-fermion interaction is used to implement weak interactions, and lepton-number violation is induced by a neutrino Majorana mass. Quantum circuits that prepare the initial nucleus and time evolve with the Hamiltonian containing the strong and weak interactions are executed on IonQ Forte Enterprise. Enabled by tuned model parameters, lepton-number violation is observed in real time, providing a clear signal of neutrinoless double-β decay. This was made possible by co-designing the simulation to maximally utilize the all-to-all connectivity and native gate-set available on IonQ's quantum computers. Quantum circuit compilation techniques and co-designed error-mitigation methods, informed from executing benchmarking circuits with up to 2,356 two-qubit gates, enabled observables to be extracted with high precision. We discuss the potential of future quantum simulations to provide yocto-second resolution of the reaction pathways in these, and other, nuclear processes.
A framework for quantum simulations of real-time weak decays of hadrons and nuclei in a two-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in the dynamics, the latter through four-Fermi effective operators. Quantum circuits that implement time evolution in this lattice theory are developed and run on Quantinuum's H1-1 20-qubit trapped ion system to simulate the $\ensuremath{\beta}$-decay of a single baryon on one lattice site. These simulations include the initial state preparation and are performed for both one and two Trotter time steps. The potential intrinsic error-correction properties of this type of lattice theory are discussed, and the leading lattice Hamiltonian required to simulate $0\ensuremath{\nu}\ensuremath{\beta}\ensuremath{\beta}$-decay of nuclei induced by a neutrino Majorana mass term is provided.
The Fermi-Hubbard model (FHM) is a simple yet rich model of strongly interacting electrons with complex dynamics and a variety of emerging quantum phases. These properties make it a compelling target for digital quantum simulation. Trotterization-based quantum simulations have shown promise, but implementations on current hardware are limited by noise, necessitating error mitigation techniques like circuit optimization and post-selection. A mapping of the FHM to a Z2 LGT was recently proposed that restricts the dynamics to a subspace protected by additional symmetries, and its ability for post-selection error mitigation was verified through noisy classical simulations. In this work, we propose and demonstrate a suite of algorithm-hardware co-design strategies on a trapped-ion quantum computer, targeting two key aspects of NISQ-era quantum simulation: circuit compilation and error mitigation. In particular, a novel combination of iteratively preconditioned gradient descent (IPG) and subsystem von Neumann Entropy compression reduces the 2-qubit gate count of FHM quantum simulation by 35%, consequently doubling the number of simulatable Trotter steps when used in tandem with error mitigation based on conserved symmetries, debiasing and sharpening techniques. Our work demonstrates the value of algorithm-hardware co-design to operate digital quantum simulators at the threshold of maximum circuit depths allowed by current hardware, and is broadly generalizable to strongly correlated systems in quantum chemistry and materials science.
With the latest advances in quantum computing technology, we are gradually moving from the noisy intermediate-scale quantum (NISQ) era characterized by hardware limited in the number of qubits and plagued with quantum noise, to the age of quantum utility where both the newest hardware and software methods allow for tackling problems which have been deemed difficult or intractable with conventional classical methods. One of these difficult problems is the simulation of one-dimensional (1D) nuclear magnetic resonance (NMR) spectra, a major tool to learn about the structure of molecules, helping the design of new materials or drugs. Using advanced error mitigation and error suppression techniques from Q-CTRL together with the latest commercially available superconducting-qubit quantum computer from IBM and trapped-ion quantum computer from IonQ, we present the quantum Hamiltonian simulation of liquid-state 1D NMR spectra in the high-field regime for spin systems up to 34 spins. Our pipeline has a major impact on the ability to execute deep quantum circuits with the reduction of quantum noise, improving mean square error by a factor of 22. It allows for the execution of deep quantum circuits and obtaining salient features of the 1D NMR spectra for both 16-spin and 22-spin systems, as well as a 34-spin system, which lies beyond the regime where unrestricted full Liouvillespace simulations are practical (32 spins, the Liouville limit). Our work is a step toward near-term quantum utility in NMR spectroscopy.
Quantum computing has long been an experimental technology with the potential to simulate, at scale, phenomena which on classical devices would be too expensive to simulate at any but the smallest scales. Over the last several years, however, it has entered the NISQ era, where the number of qubits are sufficient for quantum advantage but substantial noise on hardware stands in the way of this achievement. This thesis details NISQ device-centered improvements to techniques of quantum simulation of the out-of-equilbrium real-time dynamics of lattice quantum chromodynamics (LQCD) and of dense 3-flavor neutrino systems on digital quantum devices. The first project concerning LQCD is a comparison of methods for implementing the variational quantum eigensolver (VQE) that initializes the ground state of an SU(3) plaquette-chain. The thesis then pivots to a 1+1D lattice of quarks interacting with an SU(3) gauge-field. A VQE-based state-preparation for the vacua and a Trotterized time-evolution circuit is designed and applied to the problems of simulating beta and neutrinoless double beta decay. Finally, these circuits are adapted to a version useable on quantum devices with nearest-neighbor connectivity with minimal overhead, with an eye towards utilizing the higher qubit count of such devices for hadron dynamics and scattering. This thesis covers two projects that concern dense 3-flavor neutrino systems. The first details design and testing of Trotterized time-evolution circuits on state-of-the-art quantum devices. The second, motivated by the Gottesman-Knill theorem's result that deviation from stabilizer states ("magic") is necessary for a problem to exhibit quantum advantage, details results with implications for the Standard Model in general that the 3 flavor ultradense neutrino systems with the highest, most-persistent magic are those that start with neutrinos in all 3 flavors.
It is well known that the neutrino flavor in extreme astrophysical environments changes under the effect of three contributions: the vacuum oscillation, the interaction with the surrounding matter, and the collective oscillations due to interactions between different neutrinos. The latter adds a nonlinear contribution to the equations of motion, making the description of their dynamics complex. In this work we study various strategies to simulate the coherent collective oscillations of a system of NNN neutrinos in the two-flavor approximation using quantum computation. This was achieved by using a pair-neutrino decomposition designed to account for the fact that the flavor Hamiltonian, in the presence of the neutrino-neutrino term, presents an all-to-all interaction that makes the implementation of the evolution dependent on the qubit topology. We analyze the Trotter error caused by the decomposition demonstrating that the complexity of the implementation of time evolution scales polynomially with the number of neutrinos and that the noise from near-term quantum device simulation can be reduced by optimizing the quantum circuit decomposition and exploiting a full-qubit connectivity. We find that the gate complexity using second order Trotter-Suzuki formulas scales better with system size than with other decomposition methods such as quantum signal processing. We finally present the application and the results of our algorithm on a real quantum device based on trapped-ion qubits.
Dense-neutrino flavor evolution is a quantum many-body problem whose fully correlated treatment becomes rapidly more expensive with increasing particle number and spatial structure. We develop a \texttt{QCNO} quantum simulation code that maps a inhomogeneous forward-scattering neutrino Hamiltonian with advection and finite-range interactions to TEBD2 product-formula quantum circuits, and connects ideal many-body simulation, backend-aware execution, and fault-tolerant resource estimation within the same physical model. We reproduce the exact many-body evolution of the suppressed mean-field-like transverse fast flavor instability through N=30N=30N=30 and show that even a single open-boundary interaction generates Rënyi entanglement and non-stabilizer magic, although noisy backends still produce polarization RMSEs of order 0.10.10.1--0.30.30.3. The restricted active interaction graph of this problem allows us to estimate the circuit depth and TTT gate cost through N=50N=50N=50 for both NISQ and fault-tolerant approaches. The measured TEBD2 error in our fiducial ideal simulation of order 10−310^{-3}10−3 motivates a synthesis tolerance εsyn∼10−7\varepsilon_{\rm syn}\sim10^{-7}εsyn​∼10−7, corresponding to about 707070 TTT gates per RZR_ZRZ​ rotation. The generated circuit contains 1.4×1041.4\times10^41.4×104 TTT gates per qubit for open boundaries (twice that for closed boundaries), requiring an application-level logical-TTT error target of order 10−710^{-7}10−7 for early fault-tolerant architectures.