Source author record

Huaiyu Duan

Huaiyu Duan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2024arXiv

Equilibration of quantum many-body fast neutrino flavor oscillations

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 $ν_e {\barν}_e \leftrightarrows ν_μ{\barν_μ} $ yield a rapid equilibrium satisfying $n( ν_e) n({\bar ν}_e) = n( ν_μ) n({\bar ν}_μ) = n( ν_τ) n({\bar ν}_τ)$ in addition to the standard lepton number conservation in regimes where off-diagonal vacuum oscillations are small compared to $ν-ν$ interactions.

preprint2023arXiv

Many-body neutrino flavor entanglement in a simple dynamic model

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.

preprint2022arXiv

Classical and Quantum Evolution in a Simple Coherent Neutrino Problem

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.

preprint2022arXiv

Code Comparison for Fast Flavor Instability Simulation

The fast flavor instability (FFI) is expected to be ubiquitous in core-collapse supernovae and neutron star mergers. It rapidly shuffles neutrino flavor in a way that could impact the explosion mechanism, neutrino signals, mass outflows, and nucleosynthesis. The variety of initial conditions and simulation methods employed in simulations of the FFI prevent an apples-to-apples comparison of the results. We simulate a standardized test problem using five independent codes and verify that they are all faithfully simulating the underlying quantum kinetic equations under the assumptions of axial symmetry and homogeneity in two directions. We quantify the amount of numerical error in each method and demonstrate that each method is superior in at least one metric of this error. We make the results publicly available to serve as a benchmark.

preprint2021arXiv

Fast flavor oscillations in dense neutrino media with collisions

We investigate the impact of the nonzero neutrino splitting and elastic neutrino-nucleon collisions on fast neutrino oscillations. Our calculations confirm that a small neutrino mass splitting and the neutrino mass hierarchy have very little effect on fast oscillation waves. We also demonstrate explicitly that fast oscillations remain largely unaffected for the time/distance scales that are much smaller than the neutrino mean free path but are damped on larger scales. This damping originates from both the direct modification of the dispersion relation of the oscillation waves in the neutrino medium and the flattening of the neutrino angular distributions over time. Our work suggests that fast neutrino oscillation waves produced near the neutrino sphere can propagate essentially unimpeded which may have ramifications in various aspects of the supernova physics.

preprint2020arXiv

Fast Neutrino Flavor Conversion Modes in Multidimensional Core-collapse Supernova Models: the Role of the Asymmetric Neutrino Distributions

A dense neutrino gas, such as the one anticipated in the supernova environment, can experience fast neutrino flavor conversions on scales much shorter than those expected in vacuum probably provided that the angular distributions of $ν_e$ and $\barν_e$ cross each other. We perform a detailed investigation of the neutrino angular distributions obtained by solving the Boltzmann equations for fixed matter profiles of some representative snapshots during the post-bounce phase of core-collapse supernovae in multidimensional calculations of an $11.2\mathrm{M}_{\odot}$ and a $27\mathrm{M}_{\odot}$ progenitor models. Although the $11.2\mathrm{M}_{\odot}$ model features $ν_e - \barν_e$ angular crossings and the associated fast modes at different time snapshots, the $27\mathrm{M}_{\odot}$ model does not show any crossings within the decoupling region. We show that this can be understood by studying the multipole components of the neutrino distributions. In fact, there is a higher chance for the occurrence of $ν_e - \barν_e$ angular crossings for the zones where the multipole components of the neutrino distributions are strong enough. We also show that there can exist more than one crossings between the angular distributions of $ν_e$ and $\barν_e$. In addition, apart from the crossings within the neutrino decoupling region, there is a class of $ν_e-\barν_e$ angular crossings which appears very deep inside the proto-neutron star.

preprint2020arXiv

On the occurrence of fast neutrino flavor conversions in multidimensional supernova models

The dense neutrino medium in a core-collapse supernova or a neutron-star merger event can experience fast flavor conversions on time/distance scales that are much smaller than those of vacuum oscillations. It is believed that fast neutrino flavor transformation occurs in the region where the angular distributions of $ν_e$ and $\barν_e$ cross each other. We present the first study of this crossing phenomenon and the fast neutrino flavor conversions in multidimensional (multi-D) supernova models. We examine the neutrino distributions obtained by solving the Boltzmann transport equation for several fixed profiles which are representative snapshots taken from separate 2D and 3D supernova simulations with an $11.2 M_\odot$ progenitor model. Our research shows that the spherically asymmetric patterns of the $ν_e$ and $\barν_e$ fluxes in multi-D models can assist the appearance of the crossing between the $ν_e$ and $\barν_e$ angular distributions. In the models that we have studied, there exist unstable neutrino oscillation modes in and beyond the neutrino decoupling region which have amplitude growth rates as large as an $e$-fold per nanosecond. This finding can have important consequences for the explosion mechanism, nucleosynthesis, and neutrino signals of core-collapse supernovae.

preprint2020arXiv

Spectral swaps in a two-dimensional neutrino ring model

Neutrinos emitted deep within a supernova explosion experience a self-induced index of refraction. In the stationary, one-dimensional (1D) supernova "bulb model", this self-induced refraction can lead to a collective flavor transformation which is coherent among different neutrino momentum modes. Such collective oscillations can produce partial swaps of the energy spectra of different neutrino species as the neutrinos stream away from the proto-neutron star. However, it has been demonstrated that the spatial symmetries (such as the spherical symmetry in the bulb model) can be broken spontaneously by collective neutrino oscillations in multi-dimensional models. Using a stationary, 2D neutrino ring model we demonstrate that there exist two limiting scenarios where collective oscillations may occur. In one limit, the collective flavor transformation begins at a radius with relatively high neutrino densities and develops small-scale flavor structures. The loss of the spatial correlation in the neutrino flavor field results in similar (average) energy spectra for the anti-neutrinos of almost all energies and the neutrinos of relatively high energies. In the other limit, the flavor transformation starts at a radius where the neutrino densities are smaller (e.g., due to the suppression of the high matter density near the proto-neutron star). Although the spatial symmetry is broken initially, it is restored as the neutrino densities decrease, and the neutrinos of different flavors partially swap their energy spectra as in the 1D bulb model. This finding may have interesting ramifications in other aspects of supernova physics.

preprint2015arXiv

Collective neutrino oscillations and spontaneous symmetry breaking

Neutrino oscillations in a hot and dense astrophysical environment such as a core-collapse supernova pose a challenging, seven-dimensional flavor transport problem. To make the problem even more difficult (and interesting), neutrinos can experience collective oscillations through nonlinear refraction in the dense neutrino medium in this environment. Significant progress has been made in the last decade towards the understanding of collective neutrino oscillations in various simplified neutrino gas models with imposed symmetries and reduced dimensions. However, a series of recent studies seem to have "reset" this progress by showing that these models may not be compatible with collective neutrino oscillations because the latter can break the symmetries spontaneously if they are not imposed. We review some of the key concepts of collective neutrino oscillations by using a few simple toy models. We also elucidate the breaking of spatial and directional symmetries in these models because of collective oscillations.

preprint2015arXiv

Flavor instabilities in the multi-angle neutrino line model

Neutrino flavor oscillations in the presence of ambient neutrinos is nonlinear in nature which leads to interesting phenomenology that has not been well understood. It was recently shown that, in the two-dimensional, two-beam neutrino Line model, the inhomogeneous neutrino oscillation modes on small distance scales can become unstable at larger neutrino densities than the homogeneous mode does. We develop a numerical code to solve neutrino oscillations in the multi-angle/beam Line model with a continuous neutrino angular distribution. We show that the inhomogeneous oscillation modes can occur at even higher neutrino densities in the multi-angle model than in the two-beam model. We also find that the inhomogeneous modes on sufficiently small scales can be unstable at smaller neutrino densities with ambient matter than without, although a larger matter density does shift the instability region of the homogeneous mode to higher neutrino densities in the Line model as it does in the one-dimensional supernova Bulb model. Our results suggest that the inhomogeneous neutrino oscillation modes can be difficult to treat numerically because the problem of spurious oscillations becomes more severe for oscillations on smaller scales.

preprint2015arXiv

Flavor instabilities in the neutrino line model

A dense neutrino medium can experience collective flavor oscillations through nonlinear neutrino-neutrino refraction. To make this multi-dimensional flavor transport problem more tractable, all existing studies have assumed certain symmetries (e.g., the spatial homogeneity and directional isotropy in the early universe) to reduce the dimensionality of the problem. In this work we show that, if both the directional and spatial symmetries are not enforced in the neutrino line model, collective oscillations can develop in the physical regimes where the symmetry-preserving oscillation modes are stable. Our results suggest that collective neutrino oscillations in real astrophysical environments (such as core-collapse supernovae and black-hole accretion discs) can be qualitatively different from the predictions based on existing models in which spatial and directional symmetries are artificially imposed.

preprint2015arXiv

Neutrino flavor instabilities in a time-dependent supernova model

A dense neutrino medium such as that inside a core-collapse supernova can experience collective flavor conversion or oscillations because of the neutral-current weak interaction among the neutrinos. This phenomenon has been studied in a restricted, stationary supernova model which possesses the (spatial) spherical symmetry about the center of the supernova and the (directional) axial symmetry around the radial direction. Recently it has been shown that these spatial and directional symmetries can be broken spontaneously by collective neutrino oscillations. In this paper we analyze the neutrino flavor instabilities in a time-dependent supernova model. Our results show that collective neutrino oscillations start at approximately the same radius in both the stationary and time-dependent supernova models unless there exist very rapid variations in local physical conditions on timescales of a few microseconds or shorter. Our results also suggest that collective neutrino oscillations can vary rapidly with time in the regimes where they do occur which need to be studied in time-dependent supernova models.

preprint2015arXiv

Performance Analysis of an Astrophysical Simulation Code on the Intel Xeon Phi Architecture

We have developed the astrophysical simulation code XFLAT to study neutrino oscillations in supernovae. XFLAT is designed to utilize multiple levels of parallelism through MPI, OpenMP, and SIMD instructions (vectorization). It can run on both CPU and Xeon Phi co-processors based on the Intel Many Integrated Core Architecture (MIC). We analyze the performance of XFLAT on configurations with CPU only, Xeon Phi only and both CPU and Xeon Phi. We also investigate the impact of I/O and the multi-node performance of XFLAT on the Xeon Phi-equipped Stampede supercomputer at the Texas Advanced Computing Center (TACC).

preprint2015arXiv

Physics of neutrino flavor transformation through matter-neutrino resonances

In astrophysical environments such as core-collapse supernovae and neutron star-neutron star or neutron star-black hole mergers where dense neutrino media are present, matter-neutrino resonances (MNRs) can occur when the neutrino propagation potentials due to neutrino-electron and neutrino-neutrino forward scattering nearly cancel each other. We show that neutrino flavor transformation through MNRs can be explained by multiple adiabatic solutions similar to the Mikheyev-Smirnov-Wolfenstein mechanism. We find that for the normal neutrino mass hierarchy, neutrino flavor evolution through MNRs can be sensitive to the shape of neutrino spectra and the adiabaticity of the system, but such sensitivity is absent for the inverted hierarchy.

preprint2015arXiv

Quantum Monte Carlo calculations of the thermal conductivity of neutron star crusts

We use the quantum Monte Carlo (QMC) techniques to calculate the static structure function $S(q)$ of a one-component ion lattice and use it to calculate the thermal conductivity $κ$ of high-density solid matter expected in the neutron star crust. By making detailed comparisons with the results for the thermal conductivity obtained using standard techniques based on the one-phonon approximation (OPA) valid at low temperature, and the multi-phonon harmonic approximation expected to be valid over a wide range of temperatures, we asses the temperature regime where $S(q)$ from QMC can be used directly to calculate $κ$. We also compare the QMC results to those obtained using classical Monte Carlo to quantitatively asses the magnitude of the quantum corrections. We find that quantum effects became relevant for the calculation of $κ$ at temperature $T \lesssim 0.3 ~Ω_\mathrm{P}$, where $Ω_\mathrm{P}$ is the ion plasma frequency. At $T \simeq 0.1 ~Ω_\mathrm{P}$ the quantum effects suppress $κ$ by about $30\%$. The comparison with the results of the OPA indicates that dynamical information beyond the static structure is needed when $T \lesssim 0.1~ Ω_\mathrm{P}$. These quantitative comparisons help to establish QMC as a viable technique to calculate $κ$ at moderate temperatures in the range $T=0.1-1~Ω_\mathrm{P}$ of relevance to the study of accreting neutron stars. This finding is especially important because QMC is the only viable technique so far for calculating $κ$ in multi-component systems at low-temperatures.

preprint2014arXiv

Multipole expansion method for supernova neutrino oscillations

We demonstrate a multipole expansion method to calculate collective neutrino oscillations in supernovae using the neutrino bulb model. We show that it is much more efficient to solve multi-angle neutrino oscillations in multipole basis than in angle basis. The multipole expansion method also provides interesting insights into multi-angle calculations that were accomplished previously in angle basis.

preprint2013arXiv

Flavor Oscillation Modes In Dense Neutrino Media

We study two-flavor neutrino oscillations in homogeneous neutrino gases in which neutrinos and anti-neutrinos are in nearly pure weak interaction states initially. We find that the monopole and dipole oscillation modes can trigger flavor instabilities in the opposite neutrino mass hierarchies in a nearly isotropic neutrino gas. For a class of simple neutrino systems we are able to identify the normal modes of neutrino oscillations in the linear regime. Our results provide new insights into the recently discovered multi-azimuthal angle (MAA) instability of neutrino oscillations in supernovae.

preprint2011arXiv

Density Fluctuation Effects on Collective Neutrino Oscillations in O-Ne-Mg Core-Collapse Supernovae

We investigate the effect of matter density fluctuations on supernova collective neutrino flavor oscillations. In particular, we use full multi-angle, 3-flavor, self-consistent simulations of the evolution of the neutrino flavor field in the envelope of an O-Ne-Mg core collapse supernova at shock break-out (neutrino neutronization burst) to study the effect of the matter density "bump" left by the He-burning shell. We find a seemingly counterintuitive increase in the overall electron neutrino survival probability created by this matter density feature. We discuss this behavior in terms of the interplay between the matter density profile and neutrino collective effects. While our results give new insights into this interplay, they also suggest an immediate consequence for supernova neutrino burst detection: it will be difficult to use a burst signal to extract information on fossil burning shells or other fluctuations of this scale in the matter density profile. Consistent with previous studies, our results also show that the interplay of neutrino self-coupling and matter fluctuation could cause a significant increase in the electron neutrino survival probability at very low energy

preprint2011arXiv

Neutrino Luminosity and Matter-Induced Modification of Collective Neutrino Flavor Oscillations in Supernovae

We show that the bump in the electron number density profile at the base of the hydrogen envelope in O-Ne-Mg core-collapse supernovae causes an interesting interplay between neutrino-electron and neutrino-neutrino forward scattering effects in the flavor evolution of low-energy nu_e in the neutronization burst. The bump allows a significant fraction of the low-energy nu_e to survive by rendering their flavor evolution nonadiabatic. Increasing the luminosity of the neutronization burst shifts the bump-affected nu_e to lower energy with reduced survival probability. Similarly, lowering the luminosity shifts the bump-affected neutrinos to higher energies. While these low energy neutrinos lie near the edge of detectability, the population of bump-affected neutrinos has direct influence on the spectral swap formation in the neutrino signal at higher energies.

preprint2011arXiv

Self-induced suppression of collective neutrino oscillations in a supernova

We investigate collective flavor oscillations of supernova neutrinos at late stages of the explosion. We first show that the frequently used single-angle (averaged coupling) approximation predicts oscillations close to, or perhaps even inside, the neutrinosphere, potentially invalidating the basic neutrino transport paradigm. Fortunately, we also find that the single-angle approximation breaks down in this regime; in the full multiangle calculation, the oscillations start safely outside the transport region. The new suppression effect is traced to the interplay between the dispersion in the neutrino-neutrino interactions and the vacuum oscillation term.

preprint2010arXiv

Collective Neutrino Oscillations

We review the rich phenomena associated with neutrino flavor transformation in the presence of neutrino self-coupling. Our exposition centers on three collective neutrino oscillation scenarios: a simple bipolar neutrino system that initially consists of mono-energetic electron neutrinos and antineutrinos; a homogeneous and isotropic neutrino gas with multiple neutrino/antineutrino species and continuous energy spectra; and a generic neutrino gas in an anisotropic environment. We use each of these scenarios to illustrate key facets of collective neutrino oscillations. We discuss the implications of collective neutrino flavor oscillations for core collapse supernova physics and for the prospects of obtaining fundamental neutrino properties, e.g., the neutrino mass hierarchy and $θ_{13}$ from a future observed supernova neutrino signal.

preprint2010arXiv

Multi-Angle Simulation of Flavor Evolution in the Neutrino Neutronization Burst From an O-Ne-Mg Core-Collapse Supernova

We report results of the first 3-by-3 "multi-angle" simulation of the evolution of neutrino flavor in the core collapse supernova environment. In particular, we follow neutrino flavor transformation in the neutronization neutrino burst of an O-Ne-Mg core collapse event. Though in a qualitative sense our results are consistent with those obtained in 3-by-3 single-angle simulations, at least in terms of neutrino mass hierarchy dependence, performing multi-angle calculations is found to reduce the adiabaticity of flavor evolution in the normal neutrino mass hierarchy, resulting in lower swap energies. Differences between single-angle and multi-angle results are largest for the normal neutrino mass hierarchy. Our simulations also show that current uncertainties in the measured mass-squared and mixing angle parameters translate into uncertainties in neutrino swap energies. Our results show that at low theta-13 it may be difficult to resolve the neutrino mass hierarchy using the O-Ne-Mg neutronization neutrino burst.

preprint2010arXiv

The influence of collective neutrino oscillations on a supernova r-process

Recently, it has been demonstrated that neutrinos in a supernova oscillate collectively. This process occurs much deeper than the conventional matter-induced MSW effect and hence may have an impact on nucleosynthesis. In this paper we explore the effects of collective neutrino oscillations on the r-process, using representative late-time neutrino spectra and outflow models. We find that accurate modeling of the collective oscillations is essential for this analysis. As an illustration, the often-used "single-angle" approximation makes grossly inaccurate predictions for the yields in our setup. With the proper multiangle treatment, the effect of the oscillations is found to be less dramatic, but still significant. Since the oscillation patterns are sensitive to the details of the emitted fluxes and the sign of the neutrino mass hierarchy, so are the r-process yields. The magnitude of the effect also depends sensitively on the astrophysical conditions - in particular on the interplay between the time when nuclei begin to exist in significant numbers and the time when the collective oscillation begins. A more definitive understanding of the astrophysical conditions, and accurate modeling of the collective oscillations for those conditions, is necessary.

preprint2009arXiv

Neutrino flavour transformation in supernovae

Rapid progress has been made during recent years in the understanding of the flavour oscillations that occur as neutrinos traverse through supernova. The previous paradigm has given way and it is now clear that the neutrino signals we shall receive from future Galactic supernovae will allow us both to peer inside these extraordinary cosmic events and to probe some of the fundamental properties of these elusive particles. In this review we aim to distill the progress that has been made focusing upon the effects of the dynamic density profile and the emergence of collective flavour oscillations due to neutrino self-interactions.