Skip to content

Dynamical downfolding for localized quantum states

MetadataDetails
Publication Date2023-07-19
Journalnpj Computational Materials
AuthorsMariya Romanova, Guorong Weng, Arsineh Apelian, Vojtěch Vlček
InstitutionsUniversity of California, Santa Barbara
Citations14
AnalysisFull AI Review Included
  • Core Value Proposition: Introduction of Dynamical Downfolding (DD) methodology, an efficient theoretical framework for accurately modeling localized, strongly correlated electronic states (quantum defects) embedded within massive, weakly correlated host materials.
  • Key Achievement: The DD method successfully reproduces the experimentally measured optical excitations (Zero Phonon Lines) of the negatively charged Nitrogen-Vacancy (NV-) center in diamond.
  • Critical Dynamic Effects: The methodology fully incorporates dynamical screening effects from the host environment, which is essential for accurate prediction of the singlet-singlet transition (1.18 eV calculated vs. 1.19 eV experimental). Static approximations yield significant errors (0.7 eV).
  • Computational Efficiency: Utilizes an efficient stochastic constrained Random Phase Approximation (s-cRPA) approach, allowing the host environment (containing 1020 valence states) to be treated at minimal computational cost.
  • Mechanism: The environment’s influence is captured by renormalizing the one-body (t) and two-body (W) quasiparticle interaction terms within a small, explicitly correlated subspace (4 localized orbitals).
  • General Applicability: Provides a practical and reliable starting point for future simulations of electronic excitations and observables in localized quantum states and correlated materials.
ParameterValueUnitContext
NV- Triplet-Triplet Transition (Calculated, Dynamic)1.92eVVertical excitation energy (Exp. ZPL: 1.95 eV)
NV- Singlet-Singlet Transition (Calculated, Dynamic)1.18eVVertical excitation energy (Exp. ZPL: 1.19 eV)
NV- Singlet-Singlet Transition (Calculated, Static)0.7eVStatic screening limit (significant deviation)
Pristine Diamond QP Band Gap (GoW0)5.6eVCalculated using 4096-atom supercell
Defect Supercell Size (Converged)511atomsUsed for NV- defect calculations
Correlated Subspace Size4orbitalsMinimal model (N atom + 3 C atoms)
Host Environment Valence States (Sampled)1020statesTotal number of states stochastically sampled
Stochastic Samples (t and W terms)3200samplesUsed for convergence of interaction parameters
Screening Reduction (On-site Wiiii, Dynamic)50% to 52%%Reduction compared to bare Coulomb interaction
Time Propagation (s-cRPA)50a.u.Maximum time for induced charge density evolution
Kinetic Energy Cutoff (DFT)26EhUsed for real-space DFT implementation
Diamond Lattice Parameter3.543AngstromEquilibrium lattice constant used in simulations

The methodology combines Density Functional Theory (DFT), Many-Body Perturbation Theory (MBPT), stochastic sampling, and Exact Diagonalization (ED) in a multi-step workflow:

  1. Structural Relaxation: Atomic positions of the NV- defect in periodic diamond supercells (up to 999 atoms) were relaxed using QuantumESPRESSO, incorporating Tkatchenko-Scheffler total energy corrections.
  2. Mean-Field Starting Point: Initial electronic structure was generated using real-space, non-spin-polarized DFT calculations with the PBE functional.
  3. Subspace Localization: The correlated subspace was defined using Maximally Localized Functions (Wannier orbitals), centered on the nitrogen atom and the three nearest carbon atoms surrounding the vacancy (4 orbitals total).
  4. Environment Screening (s-cRPA): The dynamical response of the weakly correlated host environment was calculated using the stochastic constrained Random Phase Approximation (s-cRPA). This method efficiently samples the induced charge density fluctuations (polarization) using random vectors.
  5. Dynamical Downfolding:
    • One-Body Terms (t): The single-quasiparticle (QP) terms were renormalized by the environment self-energy (ÎŁenv), evaluated at the QP excitation energies (ωqp).
    • Two-Body Terms (W): The explicit QP-QP interaction terms were dynamically screened by the environment polarization (Wenv), evaluated at the two-QP excitation energies (ω2qp) obtained from auxiliary Bethe-Salpeter Equation (BSE) solutions.
  6. Hamiltonian Construction: The effective Hamiltonian (H) for the 4-orbital subspace was constructed using the dynamically renormalized t and W parameters, ensuring symmetry (Hermiticity) was restored a posteriori.
  7. Excitation Spectrum: The final many-body excited states (optical transitions) of the NV- center were obtained by performing Exact Diagonalization of the downfolded Hamiltonian.

The Dynamical Downfolding methodology, validated by its success in modeling the NV- center, is highly relevant for industries relying on precise control and understanding of localized quantum states:

  • Quantum Sensing and Metrology:
    • Directly applicable to optimizing the performance and coherence properties of solid-state quantum sensors (like NV centers) used for high-resolution magnetic field, temperature, and pressure measurements.
    • Enables predictive modeling of new defect centers in wide-bandgap materials (e.g., SiC, AlN) for next-generation sensors.
  • Solid-State Quantum Computing:
    • Crucial for the design and characterization of robust, localized qubits based on point defects, where accurate knowledge of spin-dependent energy levels and optical initialization/readout transitions is paramount.
  • Optoelectronics and Photonics:
    • Modeling of localized excitons and charge transfer dynamics in semiconductor heterostructures and quantum dots, improving the design of efficient light emitters and detectors.
  • Materials Science and Defect Engineering:
    • Serves as a high-throughput computational tool for screening and engineering novel quantum defects in bulk materials, allowing researchers to predict excited state properties and stability before costly experimental synthesis.
  • Correlated Materials Research:
    • The general DD framework can be extended to study other systems where localized correlation effects dominate, such as impurities in metals or localized magnetic moments in complex oxides.