[Paper Review] Neutrino transport with Monte Carlo method: II. Quantum Kinetic Equations

C. Kato, Hiroki Nagakura+1 more|arXiv (Cornell University)|Aug 13, 2021

Neutrino Physics ResearchPhysics and Astronomy91 references32 citations

TL;DR

This paper presents a novel Monte Carlo solver for quantum kinetic equations (QKE-MC) that self-consistently models neutrino flavor conversion, transport, and matter collisions in core-collapse supernovae and binary neutron star mergers. By embedding flavor degrees of freedom into Monte Carlo particles and introducing an effective mean free path method to suppress statistical noise, the solver accurately captures collective neutrino oscillations, including fast and slow modes, with validation against analytical solutions and prior simulations.

ABSTRACT

Neutrinos have an unique quantum feature as flavor conversions. Recent studies suggested that collective neutrino oscillations play important roles in high-energy astrophysical phenomena. Quantum kinetic equation (QKE) is capable of describing the neutrino flavor conversion, transport and matter collision self-consistently. However, we have experienced many technical difficulties in their numerical implementation. In this paper, we present a new QKE solver based on Monte Carlo (MC) approach. This is an upgraded version of our classical MC neutrino transport solver; in essence, a flavor degree of freedom including mixing state is added into each MC particle. This extension requires updating numerical treatments of collision terms, in particular for scattering processes. We deal with the technical problem by generating a new MC particle at each scattering event. To reduce statistical noise inherent in MC methods, we develop the effective mean free path method. This suppresses a sudden change of flavor state due to collisions without increasing the number of MC particles. We present a suite of code tests to validate these new modules with comparing to the results reported in previous studies. Our QKE-MC solver is developed with fundamentally different philosophy and design from other deterministic- and mesh methods, suggesting that it will be complementary to others, and potentially provide new insights into physical processes of neutrino dynamics.

Motivation & Objective

• To develop a self-consistent numerical framework for quantum kinetic equations (QKE) that simultaneously models neutrino flavor conversion, transport, and matter collisions.
• To address the technical challenges of implementing QKE in Monte Carlo methods, particularly in handling scattering processes and statistical noise.
• To validate the new QKE-MC solver against known analytical and numerical results for vacuum, matter, and collective neutrino oscillations.
• To investigate the impact of scattering processes on fast flavor conversion, especially in the context of unresolved discrepancies in the literature.
• To provide a complementary, particle-based approach to deterministic and mesh-based QKE solvers, offering new insights into nonlinear neutrino dynamics.

Proposed method

• Extends a classical Monte Carlo neutrino transport solver by adding a flavor degree of freedom to each particle, enabling direct tracking of flavor state evolution.
• Introduces a new collision treatment where a new Monte Carlo particle is generated at each scattering event to properly handle momentum-exchanged scattering and flavor transitions.
• Develops an effective mean free path (EMFP) method to suppress statistical noise from sudden flavor changes due to collisions without increasing particle count.
• Employs a cross-validation strategy using benchmark problems from the literature, including vacuum and matter oscillations, fast and slow mode collective oscillations, and scattering-influenced instabilities.
• Uses a finite-difference code as a consistency check for fast flavor conversion without scatterings, confirming agreement with the QKE-MC results.
• Performs convergence tests on angular resolution (Nθ), particle count per grid (N), and EMFP parameter (a), demonstrating robustness and accuracy.

Experimental results

Research questions

• RQ1 Can a Monte Carlo method be effectively extended to solve quantum kinetic equations for neutrino flavor dynamics, including collective oscillations and collisions?
• RQ2 How can statistical noise from collision-induced flavor transitions be suppressed in Monte Carlo simulations without increasing computational cost?
• RQ3 What is the impact of scattering processes on the onset and evolution of fast flavor conversion in dense neutrino media?
• RQ4 How does the QKE-MC solver compare quantitatively with existing deterministic and mesh-based solvers for benchmark problems?
• RQ5 Can the QKE-MC solver resolve discrepancies in the literature regarding the role of collision terms in inducing flavor instabilities?

Key findings

• The QKE-MC solver successfully reproduces vacuum and matter-induced neutrino oscillations with high accuracy, validating the core flavor evolution module.
• The solver accurately captures fast flavor conversion in the absence of scatterings, with results consistent with a finite-difference code and analytical expectations.
• For cases including scatterings, the QKE-MC solver produces results that agree with an empirical analytic formula, confirming its physical accuracy.
• Convergence tests show that Nθ = 64, N = 500, and a = 10⁻⁵ provide sufficient resolution and noise control, with results stable across different particle counts and EMFP parameters.
• The effective mean free path (EMFP) method significantly reduces statistical fluctuations from collision-induced flavor jumps, improving simulation fidelity without increasing particle count.
• The QKE-MC solver resolves a discrepancy with Shalgar & Tamborra (2021a) by confirming their results through cross-validation, establishing consensus on the physical accuracy of the method.