Approximate Excited-State Potential Energy Surfaces for Defects in Solids

This paper introduces and benchmarks an efficient approximation technique that quantifies electron-phonon coupling for solid-state defects using only excited-state forces at the ground-state geometry, demonstrating that key optical properties can be accurately predicted with minimal atomic displacements while establishing the one-dimensional accepting-mode Huang-Rhys factor as a strict upper bound.

The Gist

Imagine you are trying to understand how a tiny, broken piece of a crystal (a "defect") behaves when it gets excited by light. This isn't just about the defect itself; it's about how the defect shakes hands with the surrounding atoms in the crystal lattice. When the defect gets excited, it wants to change its shape, and the surrounding atoms have to rearrange themselves to accommodate it. This interaction is called electron-phonon coupling.

Think of the crystal lattice as a trampoline and the defect as a heavy ball sitting on it.

• The Ground State: The ball is sitting still, and the trampoline has a specific dip around it.
• The Excited State: Suddenly, the ball gets a burst of energy and wants to become a different shape (or a different weight). The trampoline has to stretch and warp to fit this new shape.

The problem scientists face is that calculating exactly how the trampoline warps for the new shape is incredibly difficult. Sometimes the math gets stuck (it won't "converge"), and sometimes it's just too expensive to run the simulation.

The Problem: "The Blindfolded Architect"

Usually, to design a building (the excited state), you need to see the blueprints and walk around the construction site to see where the walls need to move. But in this case, the architects (scientists) are blindfolded. They can see the starting point (the ground state), and they can feel the force pushing on the walls when they try to build the new room, but they can't actually walk through the finished room to see where the walls finally settle.

If they can't see the final shape, they can't accurately predict how much energy is lost as heat (vibrations) or how bright the light emitted will be.

The Solution: "The Force-Mode Shortcut"

This paper introduces a clever shortcut. Instead of trying to map the entire, complex 3D warping of the trampoline, the authors say: "Let's just look at the direction the walls are being pushed."

1. The Force Mode: Imagine you push on a wall, and it wants to move in a specific direction. The authors define a "Force Mode" which is just a line pointing in the direction of that push.

   • Analogy: If you push a shopping cart, you know it wants to go forward. You don't need to know the exact friction of every wheel to know it's going to roll forward. You just need to know the direction of your push.
   • Using just this direction, they can estimate the Zero-Phonon Line (ZPL). This is the "pure" color of light the defect emits if no energy is lost to shaking the atoms. It's like guessing the color of a lightbulb just by looking at the voltage, without needing to measure the heat it generates.

2. The Multidimensional Expansion: However, the trampoline doesn't just move in one straight line; it ripples in all directions. To get the full picture of how much energy is lost (the Huang-Rhys factor, which measures how "noisy" the defect is), they add more "ripples" to their model.

   • They start with the main push (the Force Mode).
   • Then, they add small, calculated pushes on the atoms immediately next to the defect (the "First Nearest Neighbors").
   • Then, they add pushes on the next ring of atoms out (the "Second Nearest Neighbors").
   • Analogy: It's like fixing a dent in a car. First, you push the main dent out. Then you smooth out the immediate neighbors. Then you check the next layer. You don't need to fix the entire car to make the dent look good; fixing the local area is usually enough.

The Big Discovery: The "Upper Bound" Secret

The paper also uncovered a fascinating mathematical secret about a method scientists have been using for decades called the "Accepting Mode" approximation.

• The Old Way: Scientists used to assume that all the complex ripples of the trampoline could be squashed into one single, perfect line of movement.
• The New Insight: The authors proved mathematically that this single-line method is actually a strict upper bound.
  • Analogy: Imagine you are trying to guess the weight of a suitcase. The old method says, "It weighs at most 50 lbs." The new method says, "Actually, it weighs 40 lbs, but knowing it's at most 50 lbs is still a very useful safety limit."
  • This means that even though the old method overestimates how much energy is lost to vibrations, it never underestimates it. This is actually great news! If you want to build a quantum computer, you want defects that don't lose energy. Knowing that the "worst-case scenario" (the old method) is still manageable gives you confidence that the real defect might be even better.

Why Does This Matter?

This technique is a game-changer for two reasons:

1. It Saves Time and Money: You don't need to run expensive, complex simulations that might fail. You can get a very good answer by just looking at the forces in the starting state.
2. It Opens New Doors: Scientists can now study "quantum defects" (tiny imperfections used for quantum computing and ultra-secure communication) that were previously too difficult to model. They can predict how bright these defects will shine and how stable they are without needing a supercomputer to do the impossible math.

In a nutshell: The authors found a way to predict how a crystal defect behaves when excited by looking only at the "push" it feels at the start, rather than trying to simulate the entire complex dance of atoms. It's like predicting the path of a falling leaf by looking at the wind direction at the moment it drops, rather than tracking every single twist and turn it makes until it hits the ground.

Technical Summary

1. Problem Statement

Characterizing defects and impurities in solids (e.g., for quantum information or optoelectronics) requires a precise description of electron-phonon coupling. This coupling dictates key properties such as the Zero-Phonon Line (ZPL) energy, the Huang-Rhys factor ($S$), and the efficiency of radiative vs. non-radiative transitions.

• The Bottleneck: Accurate calculation of these properties typically requires determining the excited-state equilibrium geometry and the associated atomic relaxations.
• The Challenge: In many cases, excited-state calculations fail to converge (e.g., in constrained-occupation DFT) or are computationally prohibitive (e.g., in high-throughput screening or using high-level wavefunction methods like TD-DFT or embedding). Consequently, the atomic relaxation vector ($\Delta Q$) is often inaccessible, preventing accurate quantification of electron-phonon coupling.

2. Methodology

The authors propose a novel approximation technique to quantify electron-phonon coupling using only the forces of the excited state evaluated at the ground-state equilibrium geometry. The method relies on the harmonic approximation and proceeds in two stages:

A. The Force-Mode Model (1D Approximation)

• Concept: Inspired by the "accepting-mode" approximation, the authors define a single vibrational mode parallel to the excited-state forces ($F_e$) calculated at the ground-state geometry ($R_g^0$).
• Assumptions:
  • The excited-state vibrational frequency ($\Omega_e$) is approximated by the ground-state frequency ($\Omega_g$) along this mode ("equal-mode approximation").
  • The potential energy surface is harmonic.
• Calculation:
  • The displacement $\Delta Q_F$ is derived directly from the magnitude of the forces: $\Delta Q_F \approx |F_e| / \Omega_g^2$.
  • The relaxation energy ($W_F$) and Huang-Rhys factor ($S_F$) are calculated using this single mode.
  • This allows for the estimation of the ZPL energy ($E_{ZPL} = \Delta E_{vert} - W_F$) without needing the true excited-state geometry.

B. The Multidimensional Model

• Concept: To improve the description of phonon coupling (specifically the Huang-Rhys factor), the authors extend the model to include multiple modes.
• Strategy:
  1. Start with the Force Mode.
  2. Generate additional orthogonal modes using unit Cartesian displacements on atoms near the defect (ordered by distance: 1st nearest neighbors [1NN], 2nd nearest neighbors [2NN], etc.).
  3. Apply Gram-Schmidt orthogonalization to ensure these new modes are independent of the force mode.
• Convergence: The model calculates the total relaxation energy and Huang-Rhys factor by summing contributions from these localized modes.

3. Key Contributions

1. Force-Mode Approximation: A technique to estimate ZPL energy using only ground-state geometry and excited-state forces, bypassing the need for excited-state geometry relaxation.
2. Multidimensional Convergence: Demonstration that electron-phonon coupling converges rapidly when including displacements only up to the second nearest neighbors (2NN) of the defect, rather than requiring the full supercell.
3. Theoretical Insight on Accepting Mode: A rigorous proof that the standard "accepting-mode" Huang-Rhys factor ($S_A$) is a strict upper bound on the true multidimensional Huang-Rhys factor ($S_{tot}$).
4. Anharmonicity Analysis: Identification that while the harmonic/equal-mode approximations introduce errors, they fortuitously cancel with errors from restricting the degrees of freedom to local atoms, leading to accurate results for $E_{ZPL}$ and $S$.

4. Results

The methodology was benchmarked on three well-studied defect systems:

• CN in GaN: A polar material with strong electron-phonon coupling.
• NV Center in Diamond: A non-polar system with moderate coupling.
• Carbon Dimer (C$_B$-C$_N$) in h-BN: A 2D material with weak coupling.

Key Findings:

• ZPL Energy: The single-mode Force-Mode approximation yields a reasonable estimate of the ZPL energy (errors are often < 0.2 eV) because errors in relaxation energy and frequency tend to cancel out.
• Huang-Rhys Factor ($S$): The single-mode approach significantly underestimates $S$. However, including displacements up to 2NN allows the calculated $S$ to converge to within ~10% of the fully relaxed value.
• Luminescence Spectra: The multidimensional model (up to 2NN) successfully reproduces the qualitative features of luminescence spectra, including the phonon sideband structure.
• Upper Bound Proof: For all tested cases, the accepting-mode $S_A$ was found to be larger than the true multidimensional $S_{tot}$ (e.g., overestimating by up to 40% in the carbon dimer case). This validates the accepting mode as a conservative estimate for coupling strength.

5. Significance

• Enabling High-Throughput Screening: This method allows researchers to screen thousands of defects for quantum applications (e.g., single-photon emitters) using only ground-state DFT calculations and a single excited-state force evaluation, avoiding expensive or failed excited-state relaxations.
• Robustness: The technique is effective even when excited-state calculations fail to converge or when using methods where geometry optimization is too costly.
• Theoretical Clarification: By proving $S_A \geq S_{tot}$, the paper provides a rigorous justification for the historical success of the accepting-mode approximation in predicting non-radiative rates (which depend exponentially on $S$).
• Practical Application: The authors provide a clear recipe: use the force mode for ZPL estimation, and include unit displacements up to 2NN for accurate Huang-Rhys factors and luminescence spectra.

In conclusion, this work bridges the gap between computationally expensive excited-state relaxations and the need for quantitative electron-phonon coupling data, offering a practical, accurate, and theoretically grounded approximation for defect physics in solids.