We consider the quantum magic in systems of dense neutrinos undergoing coherent flavor transformations, relevant for supernova and neutron-star binary mergers. Mapping the three-flavor-neutrino system to qutrits, the evolution of quantum magic is explored in the single scattering angle limit for a selection of initial tensor-product pure states for $N_ν\le 8$ neutrinos. For $|ν_e\rangle^{\otimes N_ν}$ initial states, the magic, as measured by the $α=2$ stabilizer Renyi entropy $M_2$ , is found to decrease with radial distance from the neutrino sphere, reaching a value that lies below the maximum for tensor-product qutrit states. Further, the asymptotic magic per neutrino, $M_2/N_ν$ , decreases with increasing $N_ν$ . In contrast, the magic evolving from states containing all three flavors reaches values only possible with entanglement, with the asymptotic $M_2/N_ν$ increasing with $N_ν$ . These results highlight the connection between the complexity in simulating quantum physical systems and the parameters of the Standard Model.We consider the quantum magic in systems of dense neutrinos undergoing coherent flavor transformations, relevant for supernova and neutron-star binary mergers. Mapping the three-flavor-neutrino system to qutrits, the evolution of quantum magic is explored in the single scattering angle limit for a selection of initial tensor-product pure states for $N_ν\le 8$ neutrinos. For $|ν_e\rangle^{\otimes N_ν}$ initial states, the magic, as measured by the $α=2$ stabilizer Renyi entropy $M_2$, is found to decrease with radial distance from the neutrino sphere, reaching a value that lies below the maximum for tensor-product qutrit states. Further, the asymptotic magic per neutrino, $M_2/N_ν$, decreases with increasing $N_ν$. In contrast, the magic evolving from states containing all three flavors reaches values only possible with entanglement, with the asymptotic $M_2/N_ν$ increasing with $N_ν$. These results highlight the connection between the complexity in simulating quantum physical systems and the parameters of the Standard Model.

arXiv:2411.04203Nov 6, 20240 citations

Classical and quantum evolution in a simple coherent neutrino problem

Joshua D. Martin, Alessandro Roggero, Huaiyu Duan +2 more

The extraordinary neutrino flux produced in extreme astrophysical environments like the early Universe, core-collapse supernovae and neutron star mergers may produce coherent quantum neutrino oscillations on macroscopic length scales. The Hamiltonian describing this evolution can be mapped into quantum spin models with all-to-all couplings arising from neutrino-neutrino forward scattering. To date many studies of these oscillations have been performed in a mean-field limit where the neutrinos time evolve in a product state. In this paper we examine a simple two-beam model evolving from an initial product state and compare the mean-field and many-body evolution. The symmetries in this model allow us to solve the real-time evolution for the quantum many-body system for hundreds or thousands of spins, far beyond what would be possible in a more general case with an exponential number ( ${2}^{N}$ ) of quantum states. We compare mean-field and many-body solutions for different initial product states and ratios of one- and two-body couplings, and find that in all cases in the limit of infinite spins the mean-field (product state) and many-body solutions coincide for simple observables. This agreement can be understood as a consequence of the fact that the typical initial condition represents a very local but dense distribution about a mean energy in the spectrum of the Hamiltonian. We explore quantum information measures like entanglement entropy and purity of the many-body solutions, finding intriguing relationships between the quantum information measures and the dynamical behavior of simple physical observables.The extraordinary neutrino flux produced in extreme astrophysical environments like the early Universe, core-collapse supernovae and neutron star mergers may produce coherent quantum neutrino oscillations on macroscopic length scales. The Hamiltonian describing this evolution can be mapped into quantum spin models with all-to-all couplings arising from neutrino-neutrino forward scattering. To date many studies of these oscillations have been performed in a mean-field limit where the neutrinos time evolve in a product state. In this paper we examine a simple two-beam model evolving from an initial product state and compare the mean-field and many-body evolution. The symmetries in this model allow us to solve the real-time evolution for the quantum many-body system for hundreds or thousands of spins, far beyond what would be possible in a more general case with an exponential number (${2}^{N}$) of quantum states. We compare mean-field and many-body solutions for different initial product states and ratios of one- and two-body couplings, and find that in all cases in the limit of infinite spins the mean-field (product state) and many-body solutions coincide for simple observables. This agreement can be understood as a consequence of the fact that the typical initial condition represents a very local but dense distribution about a mean energy in the spectrum of the Hamiltonian. We explore quantum information measures like entanglement entropy and purity of the many-body solutions, finding intriguing relationships between the quantum information measures and the dynamical behavior of simple physical observables.

arXiv:2112.12686Physical review. D/Physical review. D.Apr 27, 202235 citations

Qutrit and qubit circuits for three-flavor collective neutrino oscillations

Francesco Turro, Ivan A. Chernyshev, Ramya Bhaskar +1 more

We explore the utility of qutrits and qubits for simulating the flavor dynamics of dense neutrino systems. The evolution of such systems impacts some important astrophysical processes, such as core-collapse supernovae and the nucleosynthesis of heavy nuclei. Many-body simulations require classical resources beyond current computing capabilities for physically relevant system sizes. Quantum computers are therefore a promising candidate to efficiently simulate the many-body dynamics of collective neutrino oscillations. Previous quantum simulation efforts have primarily focused on properties of the two-flavor approximation due to their direct mapping to qubits. Here, we present new quantum circuits for simulating three-flavor neutrino systems on qutrit- and qubit-based platforms, and demonstrate their feasibility by simulating systems of two, four, and eight neutrinos on IBM and Quantinuum quantum computers.We explore the utility of qutrits and qubits for simulating the flavor dynamics of dense neutrino systems. The evolution of such systems impacts some important astrophysical processes, such as core-collapse supernovae and the nucleosynthesis of heavy nuclei. Many-body simulations require classical resources beyond current computing capabilities for physically relevant system sizes. Quantum computers are therefore a promising candidate to efficiently simulate the many-body dynamics of collective neutrino oscillations. Previous quantum simulation efforts have primarily focused on properties of the two-flavor approximation due to their direct mapping to qubits. Here, we present new quantum circuits for simulating three-flavor neutrino systems on qutrit- and qubit-based platforms, and demonstrate their feasibility by simulating systems of two, four, and eight neutrinos on IBM and Quantinuum quantum computers.

arXiv:2407.13914Physical review. D/Physical review. D.Feb 18, 202516 citations

Many-body neutrino flavor entanglement in a simple dynamic model

Martin, Joshua D., A. Roggero, Huaiyu Duan +1 more

Dense neutrino gases form in extreme astrophysical sites, and the flavor content of the neutrinos likely has an important impact on the subsequent dynamical evolution of their environment. Through coherent forward scattering among neutrinos, the flavor content of the gas evolves under a time-dependent potential which can be modeled in a quantum many-body formalism as an all-to-all coupled spin-spin interaction. This two-body potential generically introduces entanglement and greatly complicates the study of these systems. In this work we study the evolution of the quantum many-body problem as well as the typically employed mean-field approximation to it for a small number of neutrinos ( $N = 16$ ). We consider randomly chosen one- and two-body couplings in the Hamiltonian, and the resulting evolution of several initial product states. We subsequently compare many-body and mean-field predictions for one-body observables, and we consider one- and two-body entanglement to assess under what conditions the many-body and mean-field predictions are likely to disagree. Except for a special category of prototypical initial conditions, we find that the typically employed mean-field approximation is insufficient to capture the evolution of one-body operators in the systems we consider. We also observe a loss of coherence in one- and two-body trace-reduced subsystems which suggests that the evolution may be well approximated as a classical mixture of separable states.Dense neutrino gases form in extreme astrophysical sites, and the flavor content of the neutrinos likely has an important impact on the subsequent dynamical evolution of their environment. Through coherent forward scattering among neutrinos, the flavor content of the gas evolves under a time-dependent potential which can be modeled in a quantum many-body formalism as an all-to-all coupled spin-spin interaction. This two-body potential generically introduces entanglement and greatly complicates the study of these systems. In this work we study the evolution of the quantum many-body problem as well as the typically employed mean-field approximation to it for a small number of neutrinos ($N = 16$). We consider randomly chosen one- and two-body couplings in the Hamiltonian, and the resulting evolution of several initial product states. We subsequently compare many-body and mean-field predictions for one-body observables, and we consider one- and two-body entanglement to assess under what conditions the many-body and mean-field predictions are likely to disagree. Except for a special category of prototypical initial conditions, we find that the typically employed mean-field approximation is insufficient to capture the evolution of one-body operators in the systems we consider. We also observe a loss of coherence in one- and two-body trace-reduced subsystems which suggests that the evolution may be well approximated as a classical mixture of separable states.

arXiv:2301.07049Jan 17, 20231 citation

Spectral splits and entanglement entropy in collective neutrino oscillations

Amol V. Patwardhan, Michael J. Cervia, A. B. Balantekin

In environments such as core-collapse supernovae, neutron star mergers, or the early universe, where the neutrino fluxes can be extremely high, neutrino-neutrino interactions are appreciable and contribute substantially to their flavor evolution. Such a system of interacting neutrinos can be regarded as a quantum many-body system, and prospects for nontrivial quantum correlations, i.e., entanglement, developing in a gas of interacting neutrinos have been investigated previously. In this work, we uncover an intriguing connection between the entropy of entanglement of individual neutrinos with the rest of the ensemble, and the occurrence of spectral splits in the energy spectra of these neutrinos, which develop as a result of collective neutrino oscillations. In particular, for various types of neutrino spectra, we demonstrate that the entanglement entropy is highest for the neutrinos whose locations in the energy spectrum are closest to the spectral split(s). This trend demonstrates that the quantum entanglement is strongest among the neutrinos that are close to these splits, a behavior that seems to persist even as the size of the many-body system is increased.In environments such as core-collapse supernovae, neutron star mergers, or the early universe, where the neutrino fluxes can be extremely high, neutrino-neutrino interactions are appreciable and contribute substantially to their flavor evolution. Such a system of interacting neutrinos can be regarded as a quantum many-body system, and prospects for nontrivial quantum correlations, i.e., entanglement, developing in a gas of interacting neutrinos have been investigated previously. In this work, we uncover an intriguing connection between the entropy of entanglement of individual neutrinos with the rest of the ensemble, and the occurrence of spectral splits in the energy spectra of these neutrinos, which develop as a result of collective neutrino oscillations. In particular, for various types of neutrino spectra, we demonstrate that the entanglement entropy is highest for the neutrinos whose locations in the energy spectrum are closest to the spectral split(s). This trend demonstrates that the quantum entanglement is strongest among the neutrinos that are close to these splits, a behavior that seems to persist even as the size of the many-body system is increased.

arXiv:2109.08995Physical review. D/Physical review. D.Dec 27, 202150 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.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.12334Oct 9, 20260 citations

N_ν\le 8

neutrinos. For

|ν_e\rangle^{\otimes N_ν}

initial states, the magic, as measured by the

α=2

stabilizer Renyi entropy

M_2

, is found to decrease with radial distance from the neutrino sphere, reaching a value that lies below the maximum for tensor-product qutrit states. Further, the asymptotic magic per neutrino,

M_2/N_ν

, decreases with increasing

N_ν

. In contrast, the magic evolving from states containing all three flavors reaches values only possible with entanglement, with the asymptotic

M_2/N_ν

increasing with

N_ν

. These results highlight the connection between the complexity in simulating quantum physical systems and the parameters of the Standard Model.

{2}^{N}

) of quantum states. We compare mean-field and many-body solutions for different initial product states and ratios of one- and two-body couplings, and find that in all cases in the limit of infinite spins the mean-field (product state) and many-body solutions coincide for simple observables. This agreement can be understood as a consequence of the fact that the typical initial condition represents a very local but dense distribution about a mean energy in the spectrum of the Hamiltonian. We explore quantum information measures like entanglement entropy and purity of the many-body solutions, finding intriguing relationships between the quantum information measures and the dynamical behavior of simple physical observables.

N = 16

). We consider randomly chosen one- and two-body couplings in the Hamiltonian, and the resulting evolution of several initial product states. We subsequently compare many-body and mean-field predictions for one-body observables, and we consider one- and two-body entanglement to assess under what conditions the many-body and mean-field predictions are likely to disagree. Except for a special category of prototypical initial conditions, we find that the typically employed mean-field approximation is insufficient to capture the evolution of one-body operators in the systems we consider. We also observe a loss of coherence in one- and two-body trace-reduced subsystems which suggests that the evolution may be well approximated as a classical mixture of separable states.

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.