Skip to content

Decoherence of nitrogen-vacancy spin ensembles in a nitrogen electron-nuclear spin bath in diamond

MetadataDetails
Publication Date2022-08-12
Journalnpj Quantum Information
AuthorsHuijin Park, Jung‐Hyun Lee, Sang-Wook Han, Sangwon Oh, Hosung Seo
InstitutionsAjou University, Korea Institute of Science & Technology Information
Citations32
AnalysisFull AI Review Included

As an expert material scientist, I have analyzed the provided research paper on the decoherence of nitrogen-vacancy (NV) spin ensembles in diamond. The following analysis is structured for an engineering audience, adhering strictly to the specified formatting rules.

  • Core Achievement: Performed fully quantum mechanical calculations (Cluster Correlation Expansion, CCE, combined with Density Functional Theory, DFT) to determine the theoretical upper bound for the Hahn-echo spin coherence time (T2) of NV ensembles in a P1 spin bath.
  • Key Finding on Coherence: The computed T2 time exhibits a clear linear dependence on the P1 concentration ([P1]) on a log scale, with a slope of -1.06, which aligns closely with experimental results (-1.07).
  • Decoherence Mechanism: The study confirms that the Jahn-Teller (JT) effect and the resulting anisotropic hyperfine interaction in the P1 center are critical, significantly suppressing the P1 electron spin flip-flop dynamics, thereby enhancing NV coherence.
  • Theoretical vs. Experimental T2: Theoretical T2 values for a pure P1 bath are consistently 2-4 times greater than experimental results (e.g., 98.2 ”s calculated vs. 40 ”s experimental at 4 ppm), indicating the importance of other parasitic electron spins in real diamond samples.
  • Stretched Exponential Parameter (n): The coherence decay curve is found to be close to a single exponential function, with the stretched exponential parameter (n) averaging ~0.9 across the 1-100 ppm range.
  • Practical Utility: The derived T2 vs [P1] relation (T2 = 416.65 × [P1]-1.06 ”s) provides a crucial reference for non-destructive P1 concentration estimation and materials optimization for high-coherence NV devices.
ParameterValueUnitContext
T2 vs [P1] Dependence (Slope)-1.06N/ALog scale linear fit (Theoretical CCE/DFT)
T2 (Theoretical Upper Bound)98.2”sP1 concentration = 4 ppm
T2 (Theoretical Upper Bound)8.4”sP1 concentration = 40 ppm
T2 (Experimental Reference)40”sP1 concentration = 4 ppm
T2 (Experimental Reference)2”sP1 concentration = 40 ppm
Stretched Exponential Parameter (n)~0.9N/AAverage value for [P1] range 1 to 100 ppm
External Magnetic Field (B0)500GApplied parallel to NV [111] symmetry axis
P1 Concentration Range Studied1 to 100ppmRange commonly found in diamond samples
P1 JT Axis Change Time Scale103 to 105sLow temperatures (≀200 K)
DFT Energy Cutoff85RyPlane-wave basis set calculation
DFT Supercell Size1007atomsUsed for simulating isolated P1 center
  1. Central Spin Model Definition: The system was modeled as a central NV electron spin (spin-1) coupled to a bath of substitutional nitrogen impurities (P1 centers). P1 centers include an electron spin (S = 1/2) and a 14N nuclear spin (I = 1).
  2. Coherence Calculation Method: The Cluster Correlation Expansion (CCE) method was employed, specifically using the CCE-2 level of theory combined with the single-sample approach, which treats each bath state independently.
  3. Pulse Sequence: The Hahn-echo pulse sequence (π pulse applied between two free evolution periods, t) was simulated to compute the homogeneous dephasing time T2 of the NV ensemble.
  4. Spin Bath Dynamics: The P1 bath dynamics were determined by the magnetic dipole-dipole interaction between NV and P1 spins, and the internal P1 hyperfine interaction, which includes the Jahn-Teller (JT) effect derived anisotropic hyperfine tensors.
  5. DFT Parameterization: Density Functional Theory (DFT) calculations (using QUANTUM ESPRESSO) were performed to compute the spin Hamiltonian parameters of the P1 center, including the anisotropic hyperfine tensor and quadrupole moment.
  6. Bath Configuration: P1 centers were randomly distributed around the NV center. Four possible JT anisotropy axes were randomly assigned to the P1 centers and fixed during the simulation time scale.
  7. Convergence and Averaging: T2 convergence was achieved by averaging the results over 100 different bath spin configurations.
  • Quantum Sensing and Metrology: NV ensembles are essential for high-sensitivity magnetometry. The derived T2 upper bound guides the selection of diamond materials to maximize sensor coherence and sensitivity (η).
  • Materials Science and Optimization: The T2 vs [P1] relationship allows for non-destructive characterization of diamond quality, enabling the estimation of P1 concentration and the N-to-NV conversion ratio (c), which is crucial for CVD growth optimization.
  • Solid-State Quantum Computing: NV centers are leading solid-state qubit platforms. Maximizing T2 by minimizing P1-induced decoherence is necessary for developing scalable quantum registers and achieving long-lived quantum memory.
  • Quantum Network Development: Highly coherent NV-based devices are required for quantum repeaters and long-distance entanglement distribution, making the control of P1 bath dynamics paramount.
  • Nuclear Spin Hyperpolarization: Applications that require high NV densities or controlled P1 creation (e.g., for polarizing external nuclear spins) directly benefit from understanding the P1-induced coherence limits.
View Original Abstract

Abstract Nitrogen-vacancy (NV) centers in diamond have been developed into essential hardware units for a wide range of solid-state-based quantum technology applications. While such applications require the long spin coherence times of the NV centers, they are often limited due to decoherence. In this study, we theoretically investigate the decoherence of NV-spin ensembles induced by nitrogen impurities (P1 centers), which are one of the most dominant and inevitable magnetic field noise sources in diamond. We combined cluster correlation expansion and density functional theory to compute the Hahn-echo spin-coherence time of the NV centers for a broad range of P1 concentrations. Results indicate a clear linear dependence of T 2 on P1 concentrations on a log scale with a slope of −1.06, which is in excellent agreement with previous experimental results. The interplay between the Jahn-Teller effect and the hyperfine interaction in the P1 center plays a critical role in determining the bath dynamics and the resulting NV decoherence. Our results provide a theoretical upper bound for the NV-spin T 2 over a wide range of P1 densities, serving as a key reference for materials optimization and spin bath characterization to develop highly coherent NV-based devices for quantum information technology.