Skip to content

Stability and electronic structure of NV centers at dislocation cores in diamond

MetadataDetails
Publication Date2022-11-28
JournalPhysical review. B./Physical review. B
AuthorsReyhaneh Ghassemizadeh, Wolfgang Körner, Daniel F. Urban, Christian ElsÀsser
InstitutionsUniversity of Freiburg, Fraunhofer Institute for Mechanics of Materials
Citations16
AnalysisFull AI Review Included

Technical Documentation & Analysis: NV Center Self-Assembly in Dislocated Diamond

Section titled “Technical Documentation & Analysis: NV Center Self-Assembly in Dislocated Diamond”

Reference Paper: Stability and electronic structure of NV centers at dislocation cores in diamond (Ghassemizadeh et al., 2022)


This research provides critical theoretical validation for engineering linear, aligned Nitrogen-Vacancy (NV-) center arrays in diamond, a key requirement for next-generation quantum sensing and computing.

  • Energetic Trapping: DFT analysis confirms that dislocation cores act as powerful traps, lowering the NV- defect formation energy by up to 3 eV compared to bulk diamond, promoting self-assembly.
  • Optimal Structure Identified: The core of the 30° partial glide dislocation is the most energetically favorable site for NV- formation while preserving essential quantum properties.
  • Spin Triplet Preservation: The most stable configuration (NV axis parallel to the dislocation line) successfully maintains the crucial S=1 spin triplet ground state.
  • High Quantum Fidelity: Key magnetic parameters (Axial Zero-Field Splitting, D) deviate by only ~3% from bulk experimental values, confirming the viability of these structurally distorted NV centers for high-coherence applications.
  • Linear Array Potential: This mechanism offers a pathway for controlled patterning and assembly of linear NV center chains, enabling enhanced magnetic sensitivity (scaling with √X) and solid-state qubit registers.
  • 6CCVD Relevance: The findings directly support the need for high-purity, defect-engineered Single Crystal Diamond (SCD) materials, a core specialization of 6CCVD.

ParameterValueUnitContext
Maximum Formation Energy Gain (NV-)Up to 3.5eVEnergy gain relative to bulk NV- formation.
Lowest Formation Energy (30° Core, Triplet)-2.79eVPosition 1/1i, relative to bulk NV-.
Triplet/Singlet Energy Difference (ΔETS)-0.28eVTriplet state is lower in energy (more stable) than the singlet state at the 30° core (Pos 1).
Axial ZFS (D) - 30° Core (Pos 1)3.08GHzZero-Field Splitting component, essential for spin control.
Transversal ZFS (E) - 30° Core (Pos 1)120MHzNon-zero value indicates structural asymmetry/strain field effects.
13C HFS Deviation (Max)< 3%Deviation of hyperfine constants from bulk values at the 30° core (Pos 1).
Minimum Supercell Dimension15ÅRequired distance to avoid electrostatic self-interaction of periodic NV images.
DFT Cutoff Energy420eVPlane-wave basis cutoff energy used in VASP calculations.
Structural Relaxation Tolerance< 0.002eV/ÅResidual force criterion for ionic relaxation steps.

The theoretical investigation relied on advanced Density Functional Theory (DFT) simulations using large atomistic supercell models.

  1. Simulation Package: Calculations performed using the Vienna Ab Initio Simulation Package (VASP).
  2. Dislocation Modeling: Atomistic supercell models were constructed containing two oppositely oriented partial glide dislocations (30°, 90° SP, or 90° DP) to form an approximate quadrupole arrangement, compensating for long-range elastic strain fields.
  3. Supercell Size: Large supercells (1176 Carbon atoms; dimensions 30.5 Å x 17.6 Å x 15.1 Å) were utilized to ensure bulk-like conditions and minimize periodic image interactions.
  4. XC Functional: The Generalized Gradient Approximation (GGA) (Perdew, Burke, and Ernzerhof - PBE) was used for exchange-correlation energy.
  5. Band Gap Correction: The DFT-1/2 method was applied to correct the systematic underestimation of the diamond band gap by PBE, achieving a calculated gap of 5.75 eV (close to the experimental 5.47 eV).
  6. NV Position Analysis: NV defects were placed at “core positions” (structurally disturbed) and “quasi-bulk positions” (minimal strain field) identified via Differential Displacement Map (DDM) analysis.
  7. Quantum Parameter Calculation: Zero-Field Splitting (ZFS) and Hyperfine Structure (HFS) tensor components were computed using VASP subroutines, focusing on the triplet (S=1) and singlet (S=0) electronic configurations.

This research confirms that controlled crystallographic defects, specifically 30° partial glide dislocations, are highly effective for the self-assembly and alignment of high-quality NV centers. 6CCVD is uniquely positioned to supply the foundational materials and engineering support required to translate these theoretical findings into practical quantum devices.

The successful replication and extension of this research depend on ultra-high-purity diamond with controlled defect structures.

  • Optical Grade Single Crystal Diamond (SCD): 6CCVD recommends our high-purity MPCVD SCD. This material offers the necessary low background nitrogen concentration, which is critical for controlling the subsequent NV formation process (typically via irradiation and annealing).
  • Substrate Engineering: While the paper models inherent growth defects, 6CCVD can assist in sourcing or preparing SCD substrates designed to promote specific dislocation types or densities, enabling targeted experimental validation of the 30° core trapping mechanism.

Translating theoretical models into functional quantum devices requires precision material customization, which is a core strength of 6CCVD.

Requirement from Research6CCVD CapabilityTechnical Specification
Large-Scale IntegrationCustom DimensionsPlates/wafers up to 125mm (PCD) and large SCD plates.
Device FabricationUltra-Low Roughness PolishingRa < 1nm for SCD, essential for high-fidelity optical coupling and waveguide integration.
Qubit Control / ReadoutCustom MetalizationInternal capability to deposit Au, Pt, Pd, Ti, W, Cu for microwave delivery lines, electrodes, and contacts.
Defect Depth ControlPrecise Thickness ControlSCD layers available from 0.1”m up to 500”m, allowing for shallow or deep NV layer creation.

The stability of the NV triplet state is highly sensitive to the local strain field (evidenced by the non-zero E component of ZFS). 6CCVD’s in-house PhD team specializes in defect engineering and material optimization for quantum applications.

  • Material Selection: We provide consultation on selecting the optimal SCD growth parameters and post-growth processing (e.g., irradiation and annealing recipes) required to maximize NV yield and alignment along engineered dislocation lines for quantum magnetometry projects.
  • Strain Management: Our team can advise on polishing and etching techniques to minimize surface strain effects that could degrade the coherence time (T2) of the aligned NV centers.

For custom specifications or material consultation, visit 6ccvd.com or contact our engineering team directly.

View Original Abstract

We present a density functional theory analysis of the negatively charged\nnitrogen-vacancy (NV) defect complex located at or close to the core of\n30$^\circ$ and 90$^\circ$ partial glide dislocations in diamond. Formation\nenergies, electronic densities of states, structural deformations, hyperfine\nstructure and zero-field splitting parameters of NV centers in such\nstructurally distorted environments are analyzed. The formation energies of the\nNV centers are up to 3 eV lower at the dislocation cores compared to the bulk\nvalues of crystalline diamond. We found that the lowest energy configuration of\nthe NV center at the core of a 30$^\circ$ partial glide dislocation is realized\nwhen the axis of the NV center is oriented parallel to the dislocation line.\nThis special configuration has a stable triplet ground state. Its hyperfine\nconstants and zero field splitting parameters deviate by only 3% from values of\nthe bulk NV center. Hence, this is an interesting candidate for a self-assembly\nof a linear array of NV centers along the dislocation line.\n