Search

The Phase-Space approach (PSA), which was originally introduced in [Lacroix et al., Phys. Rev. D 106, 123006 (2022)] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beam. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonian. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers.

We propose a general method embedded in the ab initio nuclear framework to reconstruct linear response functions and calculate sum rules. Within our approach, based on the Gaussian integral transform, we consistently treat the groundstate and the excited spectrum. Crucially, the method allows for a robust uncertainty estimation of the spectral reconstruction. We showcase it for the spin response in neutron matter. Our calculations are performed using state-of-the-art many-body coupled-cluster method and Hamiltonians derived in the chiral effective field theory, emphasizing the analysis of finite-size effects. This work serves as a stepping stone towards further studies of neutrino interactions in astrophysical environments from first principles.

The auxiliary field diffusion Monte Carlo method uses imaginary-time projection techniques to accurately solve the ground-state wave function of atomic nuclei and infinite nuclear matter. In this work, we present a novel representation of the imaginary-time propagator based on restricted Boltzmann machines. We test its accuracy against the routinely employed Hubbard-Stratonovich transformations by evaluating ground-state energies and single-particle densities of selected light nuclei. This analysis paves the way for incorporating more realistic nuclear potentials in the auxiliary field diffusion Monte Carlo method, including isospin-dependent spin-orbit terms and cubic spin-isospin operators, which characterize accurate phenomenological and chiral effective field theory Hamiltonians., Comment: 8 pages, 3 figures

Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. This work reports the development of an algorithm that makes it possible to obtain phase shifts for generic non-relativistic elastic scattering processes on a quantum computer. This algorithm is based on extracting phase shifts from the direct implementation of the real-time evolution. The algorithm is improved by a variational procedure, making it more accurate and resistant to the quantum noise. The reliability of the algorithm is first demonstrated by means of classical numerical simulations for different potentials, and later tested on existing quantum hardware, specifically on IBM quantum processors., Comment: 17 pages, 15 figures, 4 tables

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with $N$ degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to $N$, non-trivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large-$N$ partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counter-examples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with $N$. Indeed, we derive the amusing fact that these \emph{thermodynamic} partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the "flavor-momentum-locked" state., Comment: 26 pg., 3 app, multiple figs

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other non symmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the non symmetric GEVP by means of a quantum annealer, i.e to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the proprietary simulated-anneaing package and the D-Wave Advantage 4.1 system. The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features., Comment: 20 pages, 11 figures

We show in this paper that a strong and easy connection exists between quantum error correction and Lattice Gauge Theories (LGT) by using the Gauge symmetry to construct an efficient error-correcting code for Abelian LGTs. We identify the logical operations on this gauge covariant code and show that the corresponding Hamiltonian can be expressed in terms of these logical operations while preserving the locality of the interactions. Furthermore, we demonstrate that these substitutions actually give a new way of writing the LGT as an equivalent hardcore boson model. Finally we demonstrate a method to perform fault-tolerant time evolution of the Hamiltonian within the gauge covariant code using both product formulas and qubitization approaches. This opens up the possibility of inexpensive end to end dynamical simulations that save physical qubits by blurring the lines between simulation algorithms and quantum error correcting codes.

We explore the practicality of early fault-tolerant quantum algorithms, focusing on ground-state energy estimation problems. Specifically, we address the computation of the cumulative distribution function (CDF) of the spectral measure of the Hamiltonian and the identification of its discontinuities. Scaling to bigger system sizes unveils three challenges: the smoothness of the CDF for large supports, the absence of tight lower bounds on the overlap with the actual ground state, and the complexity of preparing high-quality initial states. To tackle these challenges, we introduce a signal processing technique for identifying the inflection point of the CDF. We argue that this change of paradigm significantly simplifies the problem, making it more accessible while still being accurate. Hence, instead of trying to find the exact ground-state energy, we advocate improving on the classical estimate by aiming at the low-energy support of the initial state. Furthermore, we offer quantitative resource estimates for the maximum number of samples required to identify an increase in the CDF of a given size. Finally, we conduct numerical experiments on a 26-qubit fully-connected Heisenberg model using a truncated density-matrix renormalization group (DMRG) initial state of low bond dimension. Results show that the prediction obtained with the quantum algorithm aligns well with the DMRG-converged energy at large bond dimensions and requires several orders of magnitude fewer samples than predicted by the theory. Hence, we argue that CDF-based quantum algorithms are a viable, practical alternative to quantum phase estimation in resource-limited scenarios., Comment: 16 pages, 9 figures

Hamiltonian moments in Fourier space - expectation values of the unitary evolution operator under a Hamiltonian at different times - provide a convenient framework to understand quantum systems. They offer insights into the energy distribution, higher-order dynamics, response functions, correlation information and physical properties. This paper focuses on the computation of Fourier moments within the context of a nuclear effective field theory on superconducting quantum hardware. The study integrates echo verification and noise renormalization into Hadamard tests using control reversal gates. These techniques, combined with purification and error suppression methods, effectively address quantum hardware decoherence. The analysis, conducted using noise models, reveals a significant reduction in noise strength by two orders of magnitude. Moreover, quantum circuits involving up to 266 CNOT gates over five qubits demonstrate high accuracy under these methodologies when run on IBM superconducting quantum devices., Comment: 12 pages, 7 pages appendix, comments are welcome

Time scales associated with many-body fast neutrino flavor conversions in core-collapse supernova are explored in the context of an effective two-flavor model with axial symmetry. We present a preliminary study of time scales obtained from a linear stability analysis and from the distributions of Loschmidt echo crossing times (intimately connected to dynamical phase transitions in non-equilibrium systems) determined by time evolution with the exact many-body Hamiltonian. Starting from a tensor-product initial state describing systems of $N$ neutrinos, with $N/2$ electron-type and $N/2$ heavy-type, with uniform angular distributions, the Loschmidt echo crossing times, $t_{\mathcal{L}_{\times}}$, are found to exhibit two distinct time scales that are exponentially separated. The second peak structure at longer times, effectively absent for $N=4$, develops with increasing $N$. When re-scaled in terms of $\log t_{\mathcal{L}_{\times}}$, the distributions are found to become increasingly well described by the sum of two stable distributions. The distribution of Loschmidt echo crossing times differs somewhat from the results of the (numerical) linear stability analysis, which exhibits a peak at finite frequency and a second peak consistent with zero frequency. The exact analysis suggests that the zero-frequency instability manifests itself as a modest flavor-conversion time scale., Comment: 13 pages, 5 figures; Updated references, and acknowledgments section

Neutrino gases are expected to form in high density astrophysical environments, and accurately modeling their flavor evolution is critical to understanding such environments. In this work we study a simplified model of such a dense neutrino gas in the regime for which neutrino-neutrino coherent forward scattering is the dominant mechanism contributing to the flavor evolution. We show evidence that the generic potential induced by this effect is non-integrable and that the statistics of its energy level spaces are in good agreement with the Wigner surmise. We also find that individual neutrinos rapidly entangle with all of the others present which results in an equilibration of the flavor content of individual neutrinos. We show that the average neutrino flavor content can be predicted utilizing a thermodynamic partition function. A random phase approximation to the evolution gives a simple picture of this equilibration. In the case of neutrinos and antineutrinos, processes like $\nu_e {\bar{\nu}}_e \leftrightarrows \nu_\mu {\bar{\nu}_\mu} $ yield a rapid equilibrium satisfying $n( \nu_e) n({\bar \nu}_e) = n( \nu_\mu) n({\bar \nu}_\mu) = n( \nu_\tau) n({\bar \nu}_\tau)$ in addition to the standard lepton number conservation in regimes where off-diagonal vacuum oscillations are small compared to $\nu-\nu$ interactions., Comment: 16 pages, 8 figures, 1 appendix

Quantum computers offer an intriguing path for a paradigmatic change of computing in the natural sciences and beyond, with the potential for achieving a so-called quantum advantage, namely a significant (in some cases exponential) speed-up of numerical simulations. The rapid development of hardware devices with various realizations of qubits enables the execution of small scale but representative applications on quantum computers. In particular, the high-energy physics community plays a pivotal role in accessing the power of quantum computing, since the field is a driving source for challenging computational problems. This concerns, on the theoretical side, the exploration of models which are very hard or even impossible to address with classical techniques and, on the experimental side, the enormous data challenge of newly emerging experiments, such as the upgrade of the Large Hadron Collider. In this roadmap paper, led by CERN, DESY and IBM, we provide the status of high-energy physics quantum computations and give examples for theoretical and experimental target benchmark applications, which can be addressed in the near future. Having the IBM 100 x 100 challenge in mind, where possible, we also provide resource estimates for the examples given using error mitigated quantum computing.

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., Comment: 6 pages, 5 figures

Simulating many-body quantum systems is a promising task for quantum computers. However, the depth of most algorithms, such as product formulas, scales with the number of terms in the Hamiltonian, and can therefore be challenging to implement on near-term, as well as early fault-tolerant quantum devices. An efficient solution is given by the stochastic compilation protocol known as qDrift, which builds random product formulas by sampling from the Hamiltonian according to the coefficients. In this work, we unify the qDrift protocol with importance sampling, allowing us to sample from arbitrary probability distributions, while controlling both the bias, as well as the statistical fluctuations. We show that the simulation cost can be reduced while achieving the same accuracy, by considering the individual simulation cost during the sampling stage. Moreover, we incorporate recent work on composite channel and compute rigorous bounds on the bias and variance, showing how to choose the number of samples, experiments, and time steps for a given target accuracy. These results lead to a more efficient implementation of the qDrift protocol, both with and without the use of composite channels. Theoretical results are confirmed by numerical simulations performed on a lattice nuclear effective field theory., Comment: 17 pages, 9 pages supplemental material

In addition to the need for stable and precisely controllable qubits, quantum computers take advantage of good readout schemes. Superconducting qubit states can be inferred from the readout signal transmitted through a dispersively coupled resonator. This work proposes a novel readout classification method for superconducting qubits based on a neural network pre-trained with an autoencoder approach. A neural network is pre-trained with qubit readout signals as autoencoders in order to extract relevant features from the data set. Afterwards, the pre-trained network inner layer values are used to perform a classification of the inputs in a supervised manner. We demonstrate that this method can enhance classification performance, particularly for short and long time measurements where more traditional methods present lower performance., Comment: 16 pages, 23 figures

Accurate predictions of inclusive scattering cross sections in the linear response regime require efficient and controllable methods to calculate the spectral density in a strongly-correlated many-body system. In this work we reformulate the recently proposed Gaussian Integral Transform technique in terms of Fourier moments of the system Hamiltonian which can be computed efficiently on a quantum computer. One of the main advantages of this framework is that it allows for an important reduction of the computational cost by exploiting previous knowledge about the energy moments of the spectral density. For a simple model of medium mass nucleus like $^{40}$Ca and target energy resolution of $1$ MeV we find an expected speed-up of $\approx 125$ times for the calculation of the giant dipole response and of $\approx 50$ times for the simulation of quasi-elastic electron scattering at typical momentum transfers., Comment: 13 pages, 4 figures

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 non-linear 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 noisy 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 formulae 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-ions qubits., Comment: Updated version accepted for publication

It is for the first time that Quantum Simulation for High Energy Physics (HEP) is studied in the U.S. decadal particle-physics community planning, and in fact until recently, this was not considered a mainstream topic in the community. This fact speaks of a remarkable rate of growth of this subfield over the past few years, stimulated by the impressive advancements in Quantum Information Sciences (QIS) and associated technologies over the past decade, and the significant investment in this area by the government and private sectors in the U.S. and other countries. High-energy physicists have quickly identified problems of importance to our understanding of nature at the most fundamental level, from tiniest distances to cosmological extents, that are intractable with classical computers but may benefit from quantum advantage. They have initiated, and continue to carry out, a vigorous program in theory, algorithm, and hardware co-design for simulations of relevance to the HEP mission. This community whitepaper is an attempt to bring this exciting and yet challenging area of research to the spotlight, and to elaborate on what the promises, requirements, challenges, and potential solutions are over the next decade and beyond., Comment: This is a whitepaper prepared for the topical groups CompF6 (Quantum computing), TF05 (Lattice Gauge Theory), and TF10 (Quantum Information Science) within the Computational Frontier and Theory Frontier of the U.S. Community Study on the Future of Particle Physics (Snowmass 2021). 103 pages and 1 figure

Neutron matter, through its connection to neutron stars as well as systems like cold atom gases, is one of the most interesting yet computationally accessible systems in nuclear physics. The Configuration-Interaction Monte Carlo (CIMC) method is a stochastic many-body technique allowing to tackle strongly coupled systems. In contrast to other Quantum Monte Carlo methods employed in nuclear physics, the CIMC method can be formulated directly in momentum space allowing for an efficient use of non-local interactions. In this work we extend CIMC method to include three-nucleon interactions through the normal-ordered two-body approximation. We present results for the equation of state of neutron matter in line with other many-body calculations that employ low resolution chiral interactions, and provide predictions for the momentum distribution and the static structure factor., Comment: 9 pages, 6 figures

Compared with time independent Hamiltonians, the dynamics of generic quantum Hamiltonians $H(t)$ are complicated by the presence of time ordering in the evolution operator. In the context of digital quantum simulation, this difficulty prevents a direct adaptation of time independent simulation algorithms for time dependent simulation. However, there exists a framework within the theory of dynamical systems which eliminates time ordering by adding a "clock" degree of freedom. In this work, we provide a computational framework, based on this reduction, for encoding time dependent dynamics as time independent systems. As a result, we make two advances in digital Hamiltonian simulation. First, we create a time dependent simulation algorithm based on performing qubitization on the augmented clock system, and in doing so, provide the first qubitization-based approach to time dependent Hamiltonians that goes beyond Trotterization of the ordered exponential. Second, we define a natural generalization of multiproduct formulas for time-ordered exponentials, then propose and analyze an algorithm based on these formulas. Unlike other algorithms of similar accuracy, the multiproduct approach achieves commutator scaling, meaning that this method outperforms existing methods for physically-local time dependent Hamiltonians. Our work reduces the disparity between time dependent and time independent simulation and indicates a step towards optimal quantum simulation of time dependent Hamiltonians., Comment: 58 pages, 5 figures

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales $t\propto\mu^{-1}$, with $\mu$ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number $N$ of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales $\tau\approx\mu^{-1}\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales $\tau\approx\mu^{-1}\log(N)$. We present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

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., Comment: 16 pages, 10 figures

Quantum simulations of Lattice Gauge Theories (LGTs) are often formulated on an enlarged Hilbert space containing both physical and unphysical sectors in order to retain a local Hamiltonian. We provide simple fault-tolerant procedures that exploit such redundancy by combining a phase flip error correction code with the Gauss' law constraint to correct one-qubit errors for a $\mathbb{Z}_2$ or truncated U(1) LGT in 1+1 and 2+1 dimensions with a link flux cutoff of $1$. Unlike existing work on detecting violations of Gauss' law, our circuits are fault tolerant and the overall error correction scheme outperforms a na\"{i}ve application of the $[5,1,3]$ code. The constructions outlined can be extended to LGT systems with larger cutoffs and may be of use in understanding how to hybridize error correction and quantum simulation for LGTs in higher space-time dimensions and with different symmetry groups.

We investigate the dynamics of heavy impurities embedded in an ultra-cold Fermi gas by using a Generalized Langevin equation. The latter -- derived by means of influence functional theory -- describes the stochastic classical dynamics of the impurities and the quantum nature of the fermionic bath manifests in the emergent interaction between the impurities and in the viscosity tensor. By focusing on the two-impurity case, we predict the existence of bound states, in different conditions of coupling and temperature, and whose life-time can be analytically estimated. Our predictions should be testable using cold-gases platforms within current technology., Comment: 10 pages, 4 figures