Quantum-circuit simulation of three-flavor neutrino oscillations: vacuum, matter, and CP diagnostics

We simulate three-flavor neutrino oscillations on two-qubit quantum circuits. The flavor space is embedded as $|00⟩={\nu }_e$, $|01⟩={\nu }_{\mu }$, $|10⟩={\nu }_{\tau }$ with the $|11⟩$ channel kept inert, and the input mixing parameters (NuFIT~5.3, normal ordering) follow the global fit. We verify three physical settings against analytic or exact references: vacuum propagation, reproduced by a mass-eigenstate phase circuit to $\sim 10^{-16}$; constant-density Mikheev--Smirnov--Wolfenstein (MSW) matter effects, showing resonant flavor conversion when the matter potential $A$ crosses $\Delta m_{31}^2$; and varying-density (supernova shock-shell) propagation by slicing plus second-order Suzuki--Trotter splitting, agreeing with exact evolution to $\sim 10^{-6}$. We report a Suzuki--Trotter precision--resource calibration varying order, step number, and matrix/probability error. It delimits brute-force Trotterization on solar baselines: we quantify both standard alternatives, coherence averaging and adiabatic MSW. We also compute quantum-information diagnostics from the same model parameters, namely mode entanglement and the CP asymmetry. Finally, we extend the framework to open quantum systems via a Lindblad master equation with dephasing and absorptive channels. These results provide a reproducible link between low-energy oscillation observables and quantum-information measures.