Skip to content

Describing the Migdal effect with a bremsstrahlung-like process and many-body effects

MetadataDetails
Publication Date2021-09-09
JournalPhysical review. D/Physical review. D.
AuthorsZheng-Liang Liang, Chongjie Mo, Fawei Zheng, Ping Zhang
InstitutionsBeijing Institute of Technology, Qufu Normal University
Citations33
AnalysisFull AI Review Included

6CCVD Technical Documentation: MPCVD Diamond for Sub-GeV Dark Matter Detection

Section titled “6CCVD Technical Documentation: MPCVD Diamond for Sub-GeV Dark Matter Detection”

Research Paper Analysis: Describing the Migdal effect with a bremsstrahlung-like process and many-body effects (arXiv:2011.13352v3)


This research validates the use of crystalline diamond and silicon as highly sensitive semiconductor targets for the direct detection of sub-GeV Dark Matter (DM) via the Migdal effect.

  • Core Achievement: The study successfully integrates a quantum field theory (QFT) description, incorporating bremsstrahlung-like processes and electronic many-body effects (Random Phase Approximation, RPA), to model the Migdal excitation event rate in bulk semiconductors.
  • Material Validation: High-purity crystalline Diamond (Band Gap Eg = 5.47 eV) and Silicon (Eg = 1.12 eV) are confirmed as optimal targets, offering low thresholds for electron-hole pair production.
  • Enhanced Sensitivity: The proposed methodology yields significantly larger event rates in the low-energy regime compared to previous models, attributed to a strong $\omega^{-4}$ scaling factor.
  • Computational Rigor: Calculations utilized Density Functional Theory (DFT) and RPA (via Quantum Espresso and YAMBO codes) with rigorous convergence testing (e.g., 6x6x6 k-point mesh for Diamond) to ensure quantitative accuracy.
  • Detector Feasibility: The analysis provides critical input for designing next-generation detectors, estimating cross-section sensitivities based on single- and two-electron ionization bins for 1 kg-yr exposure.
  • Key Parameter: The low average energy required to produce an electron-hole pair ($\varepsilon$) in these materials (Diamond: 13 eV; Silicon: 3.6 eV) is crucial for maximizing sensitivity to low-mass DM.

The following hard data points were extracted from the computational and experimental parameters used in the analysis of the Migdal effect in bulk semiconductors.

ParameterValueUnitContext
Target Materials AnalyzedDiamond, SiliconN/ABulk semiconductor targets for DM detection.
Diamond Band Gap (Eg)5.47eVExperimental value matched via scissor correction.
Silicon Band Gap (Eg)1.12eVExperimental value matched via scissor correction.
Diamond Lattice Constant3.577ÅUsed in DFT calculation.
Silicon Lattice Constant5.429ÅUsed in DFT calculation.
Reference DM Mass Range (mχ)10 MeV, 100 MeV, 1 GeVN/AMasses used for differential event rate calculations.
Benchmark Cross Section (σχn)10-38cm2Reference value for event rate plots.
Diamond Energy Cutoff (Ecut)50RyUsed for convergence of the energy loss function $F(\omega)$.
Silicon Energy Cutoff (Ecut)20RyUsed for convergence of the energy loss function $F(\omega)$.
Diamond Average Energy per e-h pair ($\varepsilon$)13eVUsed for estimating ionization charge Q.
Silicon Average Energy per e-h pair ($\varepsilon$)3.6eVUsed for estimating ionization charge Q.
Energy Bin Width ($\Delta\omega$)0.05eVAdopted for calculating the dielectric matrix.

The Migdal excitation event rates were calculated by integrating advanced QFT techniques with condensed matter physics simulations.

  1. QFT Framework Integration: The Migdal effect was modeled by combining the bremsstrahlung-like process (accounting for the drag force on electrons by the recoiling ion) and electronic many-body effects (screening and collective behavior).
  2. Dielectric Function (RPA): The central component, the inverse microscopic dielectric matrix $\epsilon^{-1}(\mathbf{k}, \omega)$, was calculated using the Random Phase Approximation (RPA) via the YAMBO code. This accounts for screening and collective plasma oscillations.
  3. DFT Input Generation: Density Functional Theory (DFT) calculations were performed using the Quantum Espresso package to obtain the necessary Bloch eigenfunctions and eigenvalues for the crystalline structures.
  4. Band Gap Matching: A “scissor correction” was applied to the DFT results to align the calculated band gaps with the known experimental values (Diamond: 5.47 eV; Silicon: 1.12 eV), ensuring accurate modeling of electron excitation thresholds.
  5. Convergence Parameters: Computational parameters were rigorously tested for convergence (within 5% difference in event rates), utilizing a 6x6x6 k-point mesh and 50 Ry cutoff for diamond, and a 5x5x5 mesh and 20 Ry cutoff for silicon.
  6. Event Rate Calculation: The total event rate $R$ was derived from the velocity-averaged energy spectra, modulated by the $\omega^{-4}$ factor, and used to estimate cross-section sensitivities for single- and two-electron ionization bins.

The research confirms that high-purity crystalline diamond is a critical material for next-generation sub-GeV DM detectors. 6CCVD is uniquely positioned to supply the necessary MPCVD diamond materials and customization required to replicate and advance this research.

Research Requirement6CCVD Material RecommendationTechnical Rationale
High-Purity Crystalline TargetOptical Grade Single Crystal Diamond (SCD)SCD offers the highest purity (Nitrogen < 1 ppb), minimizing defects that could act as charge traps or generate spurious background signals, crucial for low-energy detection.
Large-Area Detector SubstratesHigh-Purity Polycrystalline Diamond (PCD)For scaling up to kg-scale detectors (e.g., 1 kg-yr exposure), PCD provides large-area wafers up to 125mm diameter, maintaining high thermal and electronic quality.
Semiconductor Device FabricationBoron-Doped Diamond (BDD)BDD can be used for creating highly conductive contacts or p-type layers, essential for fabricating robust semiconductor junctions and electrodes in detector architectures.

6CCVD’s in-house capabilities directly address the specific engineering challenges of building high-performance diamond detectors:

  • Custom Dimensions and Thickness: We provide SCD and PCD plates in custom dimensions up to 125mm (PCD) and thicknesses ranging from 0.1 ”m to 500 ”m for active layers, and substrates up to 10 mm for bulk targets.
  • Precision Polishing: Achieving the necessary surface quality for low-noise operation is critical. We offer ultra-low roughness polishing: Ra < 1 nm for SCD and Ra < 5 nm for inch-size PCD, minimizing surface leakage currents.
  • Integrated Metalization: Detector fabrication requires precise electrode deposition. 6CCVD offers internal metalization services, including Au, Pt, Pd, Ti, W, and Cu, allowing researchers to specify complex electrode patterns (e.g., Ti/Pt/Au stacks) directly on the diamond surface.
  • Global Supply Chain: We ensure reliable, global shipping (DDU default, DDP available) of sensitive materials, supporting international collaborations in fundamental physics.

6CCVD’s in-house team of PhD material scientists specializes in optimizing MPCVD diamond properties for extreme applications, including high-energy physics and quantum sensing. We offer consultation on:

  • Material Selection: Determining the optimal balance between purity (SCD) and size (PCD) for specific detector mass requirements.
  • Doping Profiles: Customizing Boron doping levels for specific conductivity and junction requirements.
  • Surface Preparation: Advising on polishing and termination methods to maximize charge collection efficiency and minimize surface defects.

For custom specifications or material consultation for similar Dark Matter detection projects, visit 6ccvd.com or contact our engineering team directly.

View Original Abstract

Recent theoretical studies have suggested that the suddenly recoiled atom\nstruck by dark matter (DM) particle is much more likely to excite or lose its\nelectrons than expected. Such Migdal effect provides a new avenue for exploring\nthe sub-GeV DM particles. There have been various attempts to describe the\nMigdal effect in liquids and semiconductor targets. In this paper we\nincorporate the treatment of the bremsstrahlung process and the electronic\nmany-body effects to give a full description of the Migdal effect in bulk\nsemiconductor targets diamond and silicon. Compared with the results obtained\nwith the atom-centered localized Wannier functions (WFs) under the framework of\nthe tight-binding (TB) approximation, the method proposed in this study yields\nmuch larger event rates in the low energy regime, due to a $\omega^{-4}$\nscaling. We also find that the effect of the bremsstrahlung photon mediating\nthe Coulomb interaction between recoiled ion and the electron-hole pair is\nequivalent to that of the exchange of a single phonon.\n