Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oscillation mitigation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Heat load measurements for the PIP-II pHB650 cryomodule

Phase-3 testing of the pHB650 cryomodule at the PIP-II Injector Test Facility was conducted to evaluate the effectiveness of heat load mitigations performed after earlier phases of testing and to continue pinpointing any sources of unexpectedly high heat loads.. The programme measured HTTS, LTTS, and 2 K isothermal/non-isothermal loads under "standard", "linac", and "simulated dynamic" operating modes, recording data both inside the cryomodule and across the bayonet can circuits. Thermal-acoustic oscillations were eliminated by replacing the original G10 cooldown-valve stem with a stainless-steel stem fitted with wipers. A newly developed Python script automated acquisition of ACNET data, performed real-time heat-load calculations, and generated plots and tables that were posted to the electronic logbook within minutes, vastly reducing manual effort and accelerating feedback between SRF and cryogenics teams. Analysis showed that JT heat-exchanger effectiveness and temperature stratification in the two-phase and relief piping strongly influence the observed loads and helped isolate sources of excess heat. The campaign demonstrates that rigorous pre-test planning, real-time diagnostics, and automated reporting can improve both accuracy and efficiency, providing a template for future PIP-II cryomodule tests and for implementing targeted heat-load mitigations.

Porwisiak, D. [Fermilab; Wroclaw Tech. U.] (ORCID:↗
DOE OSTI2025

Detection of high-frequency gravitational waves using SRF cavities

Today, apart from some isolated R\&D efforts, there are no GW experiments, yet which explore a large part of the vast frequency range above the LIGO/Virgo band. It is planned to establish an experiment at DESY and Fermilab to search for high-frequency GWs in the frequency range of 10\,kHz to 100\,MHz. The basic idea is to use superconducting radiofrequency (SRF) cavities to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. This paper summarizes the challenging environmental boundary requirements, and the R\&D to operate a cavity using a LLRF system which pushes beyond state-of-the-art accuracy and resolutions and a seismic noise mitigated cryostat at 1.8\,K. The focus of this paper is the warm and cold commissioning of a prototype cavity, built 20 years ago during the MAGO collaboration, and its first measurement in our collaborative research project.

Giaccone, Bianca [Fermilab]↗
DOE OSTI2025

Exploring the Nature of Neutrinos with the Deep Underground Neutrino Experiment (DUNE)

The Deep Underground Neutrino Experiment (DUNE) is a next-generation experimental program designed to study the behaviour of neutrino oscillation. DUNE will utilize a neutrino beam originating at Fermilab, near Chicago, and will leverage a detector at Fermilab (Near Detector) and a detector 1300 km away in South Dakota (Far Detector), south of Saskatchewan. In the first phase of DUNE, its Far Detector will comprise of two 10,000 ton (fiducial) liquid argon (LAr) time-projection chamber (TPC) modules – powerful tracking calorimeter detectors – placed nearly a mile underground. With this large, sensitive, underground detector, DUNE aims to collect a high statistics and pure sample of neutrinos at the Far Detector. This setup also offers the potential to study non-beam physical processes via e.g. neutrinos produced in the atmosphere, supernova neutrino bursts, and/or solar neutrinos, etc. A second phase will aim to add more detector mass and expand the program. The Near Detector will consist of a LAr TPC module as well: critical to constraining systematic uncertainties in the oscillation analysis. However, this LAr TPC will have a novel design using a pixel-based readout instead of the traditional wire-based readout. This and the segmentation of the LAr TPC into multiple units are crucial in mitigating the high multiplicity of neutrino interactions expected in any readout window given its proximity to the beam. The Near Detector will feature additional components and capability beyond the LAr TPC, allowing one to deeply characterize the neutrino flux. Due to the complexity of this experimental program, several smaller-scale prototype detectors have been operating to test, validate, and improve both the technical designs and software for processing and analyzing events. By operating in charged particle test beams or neutrino beams, several of the prototypes are also capable of producing valuable results. Canadian institutions are involved in the realization of the DUNE through efforts with both the Near and Far Detectors and prototypes. DUNE is anticipated to begin operating near the end of this decade/the beginning of the next. This talk will focus on the overall DUNE program, for example its ultimate plans, status, and the efforts with prototypes.

Howard, Bruce [York U., Canada; Fermilab]↗
DOE OSTI2025

Suppressing the sample variance of DESI-like galaxy clustering with fast simulations

Ongoing and upcoming galaxy redshift surveys, such as the Dark Energy Spectroscopic Instrument (DESI) survey, will observe vast regions of sky and a wide range of redshifts. In order to model the observations and address various systematic uncertainties, N-body simulations are routinely adopted, however, the number of large simulations with sufficiently high mass resolution is usually limited by available computing time. Therefore, achieving a simulation volume with the effective statistical errors significantly smaller than those of the observations becomes prohibitively expensive. In this study, we apply the Convergence Acceleration by Regression and Pooling (CARPool) method to mitigate the sample variance of the DESI-like galaxy clustering in the AbacusSummit simulations, with the assistance of the quasi-N-body simulations FastPM. Based on the halo occupation distribution (HOD) models, we construct different FastPM galaxy catalogs, including the luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars, with their number densities and two-point clustering statistics well matched to those of AbacusSummit. We also employ the same initial conditions between AbacusSummit and FastPM to achieve high cross-correlation, as it is useful in effectively suppressing the variance. Our method of reducing noise in clustering is equivalent to performing a simulation with volume larger by a factor of 5 and 4 for LRGs and ELGs, respectively. We also mitigate the standard deviation of the LRG bispectrum with the triangular configurations k 2 = 2k 1 = 0.2 h Mpc -1 by a factor of 1.6. With smaller sample variance on galaxy clustering, we are able to constrain the baryon acoustic oscillations (BAO) scale parameters to higher precision. The CARPool method will be beneficial to better constrain the theoretical systematics of BAO, redshift space distortions (RSD) and primordial non-Gaussianity (NG).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Detection of high-frequency gravitational waves using SRF cavities

Today, apart from some isolated R&D efforts, there are no gravitational wave (GW) experiments, yet which explore a large part of the vast frequency range above the LIGO/Virgo band. It is planned to establish an experiment at Deutsches Elektronen-Synchrotron (DESY) and at the Superconducting Quantum Materials and Systems (SQMS) Center at Fermi National Accelerator Laboratory (Fermilab) to search for high-frequency GWs in the frequency range of 10 kHz to 100 MHz. The basic idea is to use superconducting radiofrequency (SRF) cavities to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. This paper summarizes the challenging environmental boundary requirements, and the R&D to operate a cavity using a low level RF (LLRF) system which pushes beyond state-of-the-art accuracy and resolutions and a seismic noise mitigated cryostat at 1.8 K. The focus of this paper is the warm and cold commissioning of a prototype cavity, built 20 years ago during the MAGO collaboration, and its first measurement in our collaborative research project.

Wenskat, M. [DESY; Hamburg U.] (ORCID:000000016546↗
DOE OSTI2025

The Dark Energy Camera All Data Everywhere cosmic shear project V: Constraints on cosmology and astrophysics from 270 million galaxies across 13,000 deg$^2$ of the sky

We present constraints on models of cosmology and astrophysics using cosmic shear data vectors from three datasets: the northern and southern Galactic cap of the Dark Energy Camera All Data Everywhere (DECADE) project, and the Dark Energy Survey (DES) Year 3. These data vectors combined consist of 270 million galaxies spread across 13,000 ${\rm deg}^2$ of the sky. We first extract constraints for $Λ$CDM cosmology and find $S_8= 0.805^{+0.019}_{-0.019}$ and $Ω_{\rm m} = 0.262^{+0.023}_{-0.036}$, which is consistent within $1.9 σ$ of constraints from the Planck satellite. Extending our analysis to dynamical dark energy models shows that lensing provides some (but still minor) improvements to existing constraints from supernovae and baryon acoustic oscillations. Finally, we study six different models for the impact of baryons on the matter power spectrum. We show the different models provide consistent constraints on baryon suppression, and associated cosmology, once the astrophysical priors are sufficiently wide. Current scale-cut approaches for mitigating baryon contamination result in a residual bias of $\approx 0.3σ$ in the $S_8, Ω_{\rm m}$ posterior. Using all scales with dedicated baryon modeling leads to negligible improvement as the new information is used solely to self-calibrate the baryon model on small scales. Additional non-lensing datasets, and/or calibrations of the baryon model, will be required to access the full statistical power of the lensing measurements. The combined dataset in this work represents the largest lensing dataset to date (most galaxies, largest area) and provides an apt testing ground for analyses of upcoming datasets from Stage IV surveys. The DECADE shear catalogs, data vectors, and likelihoods are made publicly available.

Anbajagane, D. [Chicago U., Astron. Astrophys. Ctr↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2026

Toward the origins of binding energy shifts and “satellites” formation during plasma-XPS measurements

In-plasma x-ray photoelectron spectroscopy (plasma-XPS) emerges as a powerful platform for real-time, in situ chemical analysis under conditions relevant to semiconductor processing and other plasma-enabled technologies. This study investigates the origins of binding energy (BE) shifts and the formation of “satellite” peaks observed during plasma-XPS measurements across conductive, dielectric, and gas-phase systems. Using a standard laboratory-based ambient pressure XPS apparatus coupled with an alternating current (AC)-driven capacitively coupled plasma source, we demonstrate that metastable surface species, such as transient Au oxides, can be detected during plasma exposure, revealing chemical states that are hardly accessible using conventional ultrahigh vacuum (UHV) XPS. In dielectric samples (e.g., undoped diamond, sapphire), we observe pressure- and plasma-type-dependent BE shifts of more than 50 eV, attributed to x-ray-induced and plasma-mediated surface charging. These shifts are mitigated at higher pressures/plasmas or in electronegative plasmas (e.g., O 2 ), the latter due to enhanced charge compensation mechanisms involving slow negative ions. For gas-phase species, AC-plasma excitation leads to spectral broadening and the emergence of “satellite” peaks with energy separations of a few electron volts, linked to oscillating local plasma potentials in the probing volume. Furthermore, these findings highlight the complex and important interplay between plasma parameters, surface charging, and local electric fields in shaping XPS spectra. Overall, plasma-XPS emerges as a critical metrological tool for probing transient surface chemistry, with implications for semiconductor processing, material synthesis, and plasma diagnostics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Resonant Emittance Mixing of Flat Beams in Plasma Accelerators

Linear colliders rely on high-quality flat beams to achieve the desired event rate while avoiding potentially deleterious beamstrahlung effects. Here, we show that flat beams in plasma accelerators can be subject to quality degradation due to emittance mixing. This effect occurs when the beam particles' betatron oscillations in a nonlinearly coupled wakefield become resonant in the horizontal and vertical planes. Emittance mixing can lead to a substantial decrease of the luminosity, the main quantity determining the event rate. In some cases, the use of laser drivers or flat particle beam drivers may decrease the fraction of resonant particles and, hence, mitigate emittance deterioration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Enhancing DESI DR1 full-shape analyses using HOD-informed priors

We present an analysis of DESI Data Release 1 (DR1) that incorporates Halo Occupation Distribution (HOD)-informed priors into Full-Shape (FS) modeling of the power spectrum based on cosmological perturbation theory (PT). By leveraging physical insights from the galaxy-halo connection, these HOD-informed priors on nuisance parameters substantially mitigate projection effects in extended cosmological models that allow for dynamical dark energy. The resulting credible intervals now encompass the posterior maximum from the baseline analysis using gaussian priors, eliminating a significant posterior shift observed in baseline studies. In the ΛCDM framework, a combined DESI DR1 FS information and constraints from the DESI DR1 baryon acoustic oscillations (BAO) — including Big Bang Nucleosynthesis (BBN) constraints and a weak prior on the scalar spectral index — yields Ω m = 0.2994 ± 0.0090 and σ 8 = 0.836$^{+0.024}_{-0.027}$, representing improvements of approximately 4% and 23% over the baseline analysis, respectively. For the w 0 w a CDM model, our results from various data combinations are highly consistent, with all configurations converging to a region with w 0 > -1 and w a < 0. This convergence not only suggests intriguing hints of dynamical dark energy but also underscores the robustness of our HOD-informed prior approach in delivering reliable cosmological constraints.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Overview of physics results from the ADITYA-U tokamak and future experiments

The ADITYA upgrade (ADITYA-U), a medium-sized (R 0 = 75 cm, a = 25 cm) conventional tokamak facility in India, has been consistently producing experiments findings by usingcircular and shaped-plasmas. Recognizing the plasma parameters aligning closely with the design parameters of circular limited plasmas, ADITYA-U shifted its focus toward exploring the operational regime for experimentation on saw-tooth and MHD phenomena. Moreover, ADITYA-U has made consistent advancements toward conducting preliminary plasma shaping experiments through the activation of top and bottom divertor coils utilizing hydrogen as well as deuterium fuels. Confinement is improved by a factor of ~1.5 in D 2 plasmas when compared to H 2 plasmas of ADITYA-U. Further, ADITYA-U operations emphasize preventing disruptions and runaway electrons (REs) to ensure safe operations for future fusion devices. Significant suppression of REs has been achieved in ADITYA-U with the application of pulsed localized vertical magnetic field (LVF) perturbation, thereby establishing the technique's independence from the tokamak device. The successful RE mitigation requires a critical threshold of LVF pulse magnitude, which is approximately 1% of the toroidal magnetic field, and a minimum duration of ~5 ms. Apart from this, several novel findings have been achieved in the ADITYA-U experiments, including the modification of sawtooth duration through gas-puff, the emergence of MHD-induced geodesic acoustic mode-like oscillations, the propagation of fast heat pulses induced by MHD activity, the control of RE dynamics through Gas-puffs, the propagation of pinch-driven cold-pulses, the transport and core accumulations of argon impurities, the mass dependency of plasma toroidal rotation and the detection of 'RICE' scaling, as well as the characterization of edge plasma using wall conditioning methods, such as glow discharge cleaning using a combination of Ar-H 2 mixture, localized wall cleaning by electron cyclotron resonant plasma, and the development of machine learning-based disruption predictions, will be discussed in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Edge treatment for spurious mode suppression in thin-film lithium niobate resonators

Abstract Thin-film lithium niobate is an attractive material for RF acoustic devices because of its high electromechanical coupling. However, due to the large coupling and the high anisotropy, thin-film lithium niobate resonators are prone to accidental resonances called spurious modes. These modes compromise the frequency response of the resonators, limiting their use in filter and oscillator applications. In this work, we present a novel method of spurious mode suppression through a special edge treatment etch process. Two thin-film lithium niobate resonators were fabricated, one with smooth sidewalls and one with the edge treatment. It was found that the edge-treated resonators show a weaker spurious mode response. This is potentially a new way to mitigate spurious resonances, a major issue in lithium niobate Lamb wave devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-speed fiber-coupled terahertz time-domain spectroscopy system using a bias-free plasmonic source

We demonstrate a high-speed fiber-coupled terahertz time-domain spectroscopy (THz-TDS) system using a bias-free plasmonic terahertz source. The system is powered by a compact, low-power femtosecond fiber laser oscillator that drives a dual-branch fiber amplifier to pump and probe a plasmonic terahertz source and detector, respectively. A polarization-maintaining, fiber-coupled optical delay line is employed to precisely control the temporal offset between the pump and probe beams through a compact, robust platform. Operating in the self-similar amplification regime, the fiber amplifiers are engineered to mitigate nonlinear pulse distortion and dispersion while delivering robust, high-energy pulses at a 1560-nm wavelength, with an output power exceeding 326 mW. In reflection mode, the THz-TDS system achieves a peak dynamic range of 81 dB and a bandwidth greater than 3.5 THz at a scan rate of 3.75 Hz. The utility of the system for rapid sensing is further demonstrated by quantifying the porosity and density of pharmaceutical tablets through measurements of their physical thickness and effective refractive index.

Jiang, Xinghe [University of California, Los Angel↗
DOE OSTI2025

DESI 2024 VI: cosmological constraints from the measurements of baryon acoustic oscillations

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-α forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range 0.1 < z < 4.2. To mitigate confirmation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data alone are consistent with the standard flat ΛCDM cosmological model with a matter density Ω m =0.295±0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and the robustly measured acoustic angular scale from the cosmic microwave background (CMB), DESI requires H 0 =(68.52±0.62) km s -1 Mpc -1 . In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find Ω m =0.307± 0.005 and H 0 =(67.97±0.38) km s -1 Mpc -1 . Extending the baseline model with a constant dark energy equation of state parameter w, DESI BAO alone requirew=-0.99 +0.15 -0.13 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI with CMB or with type Ia supernovae (SN Ia) individually prefer w 0 > -1 and w a < 0. This preference is 2.6σ for the DESI+CMB combination, and persists or grows when SN Ia are added in, giving results discrepant with the ΛCDM model at the 2.5σ, 3.5σ or 3.9σ levels for the addition of the Pantheon+, Union3, or DES-SN5YR supernova datasets respectively. For the flat ΛCDM model with the sum of neutrino mass ∑ m ν free, combining the DESI and CMB data yields an upper limit ∑ m ν < 0.072 (0.113) eV at 95% confidence for a ∑ m ν > 0 (∑ m ν > 0.059) eV prior. These neutrino-mass constraints are substantially relaxed if the background dynamics are allowed to deviate from flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Measuring and unbiasing the BAO shift in the Ly α forest with AbacusSummit

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) places sub- per cent constraints on measurements of the Baryon Acoustic Oscillation (BAO) scaling parameters from the Ly $\alpha$ forest. However, no systematic error budget stemming from non-linearities in the three-dimensional clustering of the Ly $\alpha$ forest is included in the DESI-Ly $\alpha$ analysis. In this work, we measure the size of the shift of the BAO peak using large Ly $\alpha$ forest mocks produced on the N-body simulation suite AbacusSummit, which adopt the Fluctuating–Gunn–Peterson Approximation (FGPA). Specifically, we measure the Ly $\alpha$ autocorrelation and the Ly $\alpha$-quasar cross-correlation functions. To mitigate the noise, we adopt a linear control variates technique, reducing the error bars by a factor of up to $\sim \sqrt{50}$ on large scales. From the autocorrelation, we detect a small positive shift in radial direction of $\Delta \alpha _{\parallel }= 0.35~{{\ \rm per\ cent}}$ at the 3$\sigma$ level and virtually no shift in the transverse direction, $\alpha _\perp$. From the cross-correlation, we see a similar shift to $\Delta \alpha _\parallel$, albeit with larger error bars, and a small negative shift, $\Delta \alpha _{\perp }=\sim$0.25 per cent, at the 2$\sigma$ level. We also make a connection with the Ly $\alpha$ forest effective field theory (EFT) framework and find that the one-loop EFT power spectrum yields unbiased measurements of the BAO shift parameters in radial and transverse direction for Ly $\alpha$ auto- and the Ly $\alpha$-quasar cross-correlation measurements. When using the one-loop EFT framework, we find that we can recover the BAO parameters without a shift, which has important implications for future Ly $\alpha$ forest analyses based on EFT. This work paves the way for novel full-shape analyses of the currently observing DESI and future surveys such as the PFS, WEAVE-QSO, and 4MOST.

Hadzhiyska, Boryana↗
DOE OSTI2025

Collective neutrino oscillations in three flavors on qubit and qutrit processors

Collective neutrino flavor oscillations are of primary importance in understanding the dynamic evolution of core-collapse supernovae and subsequent terrestrial detection, but also among the most challenging aspects of numerical simulations. This situation is complicated by the quantum many-body nature of the problem due to neutrino-neutrino interactions, which demands a quantum treatment. An additional complication is the presence of three flavors, which often is approximated by the electron flavor and a heavy lepton flavor. In this work, we provide both qubit and qutrit encodings for all three flavors, and develop optimized quantum circuits for the time evolution and analyze the Trotter error. We conclude our study with a hardware experiment of a system of two neutrinos with superconducting hardware: the IBM Torino device for qubits and Advanced Quantum Testbed device at the Lawrence Berkeley National Laboratory for qutrits. We find that error mitigation greatly helps in obtaining a signal consistent with simulations. Finally, while hardware results are comparable at this stage, we expect the qutrit setup to be more convenient for large-scale simulations since it does not suffer from probability leakage into nonphysical qubit space, unlike the qubit setup.

Neutrino oscillations↗