Does fermionic entanglement always outperform bosonic entanglement in dilaton black hole?

It has traditionally been believed that fermionic entanglement generally outperforms bosonic entanglement in relativistic frameworks, and that bosonic entanglement experiences sudden death in extreme gravitational environments. In this study, we analyze the genuine N-partite entanglement, measured by negativity, of bosonic and fermionic GHZ states, focusing on scenarios where a subset of $m$ ($m<N$) constituents interacts with Hawking radiation generated by a Garfinkle-Horowitz-Strominger (GHS) dilaton black hole. Surprisingly, we find that quantum entanglement between the non-gravitational and gravitational modes for the bosonic field is stronger than that in the same modes for the fermionic field within dilaton spacetime. This study challenges the traditional belief that ``fermionic entanglement always outperforms bosonic entanglement" in the relativistic framework. However, quantum entanglement between the gravitational modes and the combined gravitational and non-gravitational modes is weaker for the bosonic field than for the fermionic field in the presence of a dilaton black hole. Finally, the connection between the global N-partite entanglement in the bosonic field and that in the fermionic field is influenced by the gravitational field’s intensity. Our study reveals the intrinsic relationship between quantum entanglement of bosonic and fermionic fields in curved spacetime from a new perspective, and provides theoretical guidance for selecting appropriate field-based quantum resources for relativistic quantum information tasks under extreme gravitational conditions.

💡 Research Summary

The paper investigates whether the widely held belief that fermionic entanglement is always more robust than bosonic entanglement holds in the extreme gravitational environment of a Garfinkle‑Horowitz‑Strominger (GHS) dilaton black hole. The authors consider N‑partite Greenberger‑Horne‑Zeilinger (GHZ) states for both a scalar (bosonic) field and a spin‑½ (fermionic) field. A subset of m observers (with m < N) is placed near the black‑hole event horizon while the remaining N − m observers stay in asymptotically flat spacetime.

First, the paper quantizes the bosonic and fermionic fields on the GHS background. Using the Klein‑Gordon equation for the scalar field and the Dirac equation for the fermion, the authors obtain mode solutions in Kruskal coordinates and derive Bogoliubov transformations that relate the creation and annihilation operators in the Kruskal vacuum to those in the exterior (“out”) and interior (“in”) regions of the horizon. For the bosonic field the Bogoliubov coefficients contain the factor (e^{-8\pi\omega(M-\epsilon)}) (Bose‑Einstein statistics), whereas for the fermionic field the coefficients involve ((e^{-8\pi\omega(M-\epsilon)}+1)^{-1}) (Fermi‑Dirac statistics). Tracing over the inaccessible interior modes yields thermal spectra for both fields, with the expected Bose‑Einstein and Fermi‑Dirac distributions.

Next, the N‑partite GHZ state (|\text{GHZ}_N\rangle = \frac{1}{\sqrt{2}}(|0\rangle^{\otimes N}+|1\rangle^{\otimes N})) is expressed in terms of the transformed modes. After tracing out the interior degrees of freedom, the authors obtain mixed density matrices (\rho_B^{(N)}) and (\rho_F^{(N)}) for the bosonic and fermionic cases, respectively.

To quantify multipartite entanglement they employ negativity. They define the one‑tangle (N_{A(BCD)}) (negativity of subsystem A versus the rest), the two‑tangle (N_{AB}) (pairwise negativity), and the residual entanglement (\pi_A = N_{A(BCD)}^2 - N_{AB}^2 - N_{AC}^2 - N_{AD}^2). For four‑partite systems the global measure is the average residual entanglement, the (\pi_4)-tangle: (\pi_4 = \frac{1}{4}\sum_{X=A}^{D}\pi_X). These quantities are computed analytically (or numerically where necessary) as functions of the Hawking temperature (T_H = 1/