Diagrammatic quantum Monte Carlo toward the calculation of transport properties in disordered semiconductors

This paper proposes a numerically exact, system-size-independent diagrammatic quantum Monte Carlo method that unifies dynamic and static disorder to efficiently calculate transport properties, such as carrier mobility, in realistic disordered semiconductors within the thermodynamic limit.

The Gist

Imagine you are trying to predict how fast a crowd of people can run through a giant, chaotic maze. This maze represents a semiconductor (the material inside your computer chips or solar panels). The people are electrons (the carriers of electricity).

In a perfect world, the maze would be a straight, empty hallway. The people would run at the speed of light. But in the real world, the maze is messy. There are two main types of messiness:

1. Static Disorder (The Broken Tiles): Imagine some tiles in the hallway are permanently cracked, or some walls are slightly crooked because of bad construction. These are impurities, vacancies, or defects in the material. They don't move; they are just there, permanently blocking or slowing down the runners.
2. Dynamic Disorder (The Dancing Floor): Now imagine the floor itself is vibrating, shaking, and rippling because of heat or sound waves. The people have to dodge these moving obstacles. In physics, these vibrations are called phonons (sound/heat waves), and the interaction between the runners and the shaking floor is electron-phonon interaction.

The Problem

For decades, scientists have struggled to calculate exactly how fast these electrons can move (their mobility) when both types of messiness are happening at the same time.

• If they only looked at the broken tiles, they missed the shaking floor.
• If they only looked at the shaking floor, they missed the broken tiles.
• If they tried to simulate the whole maze on a computer, the math got so huge and complex that the computer crashed, or the results were just guesses.

The Solution: A New "Magic Lens"

The authors of this paper, Yu-Chen Wang and Yi Zhao, have invented a new mathematical tool called Diagrammatic Quantum Monte Carlo (DQMC).

Think of this tool as a super-powered, magical lens that lets you see the future of the electron's journey without having to simulate every single step of every single person in the crowd.

Here is how their "magic lens" works, using simple analogies:

1. The "Recipe Book" (The Diagrams)

Instead of trying to track every electron, the authors break the problem down into a "recipe book" of possible paths (called Feynman diagrams).

• Imagine you want to know the average time it takes to get from point A to point B.
• Instead of timing one person, you write down every possible way they could get there: "Run straight," "Jump over a crack," "Dance around a vibration," "Wait for a tile to settle."
• Their method creates a massive list of these "recipes" (diagrams) and uses a clever sampling technique (Monte Carlo) to pick the most important ones to calculate, ignoring the ones that don't matter.

2. The "Universal Translator" (Reciprocal Space)

The biggest breakthrough in this paper is how they handle the Static Disorder (the broken tiles).

• Usually, broken tiles are described in "Real Space" (e.g., "Tile #42 is cracked"). This is messy to calculate when the maze is infinite.
• The authors invented a way to translate this "Real Space" messiness into "Reciprocal Space" (a fancy math language that describes patterns and frequencies).
• The Analogy: Imagine instead of looking at a messy pile of Lego bricks, you look at the pattern of how they fit together. By translating the broken tiles into this pattern language, they found that the "broken tiles" behave mathematically very similarly to the "shaking floor."
• This allowed them to treat both the permanent defects and the moving vibrations using the same set of rules. It's like realizing that both a pothole and a speed bump can be solved with the same suspension system.

3. The "Infinite Maze" (Thermodynamic Limit)

Most computer simulations have to stop at a certain size because they run out of memory. They simulate a small chunk of the material and hope it represents the whole.

• This new method is special because it doesn't care how big the maze is.
• Whether the semiconductor is the size of a grain of sand or the size of a planet, the calculation cost stays the same. They can calculate the behavior of an infinite material, which is exactly what real materials are (effectively).

Why Does This Matter?

This new tool is like giving scientists a crystal ball for designing better electronics.

• Better Solar Panels: It helps us understand why some materials capture energy better than others, even when they are full of defects.
• Faster Computers: It helps predict how fast electrons can move in new types of organic semiconductors (plastic-like electronics), which could lead to flexible, cheap, and super-fast devices.
• Solving Mysteries: For years, scientists have been confused by materials that act like they are "frozen" (localized) but still conduct electricity like a fluid. This tool helps explain that weird behavior by showing how the "shaking floor" and "broken tiles" work together to create a unique dance for the electrons.

The Bottom Line

The authors have built a universal, exact, and efficient calculator that can handle the messy reality of the real world. They figured out how to mix the math of "permanent damage" (static disorder) with "moving chaos" (dynamic disorder) into one smooth formula.

This means we can finally stop guessing how new materials will behave and start designing them with confidence, leading to faster phones, better batteries, and more efficient energy systems.

Technical Summary

1. Problem Statement

The transport properties of charge carriers (electrons, holes, excitons) in disordered semiconductors and light-harvesting systems are governed by a complex interplay of dynamic disorder (electron-phonon interactions) and static disorder (structural imperfections like impurities, vacancies, and irregular stacking).

• Challenges: Existing theoretical models often fail to capture the full complexity of these systems.
  • Hopping models neglect quantum coherence.
  • Mixed quantum-classical methods often violate detailed balance and fail to describe quantum coherence accurately.
  • Real-time quantum dynamics are computationally prohibitive for large systems and suffer from the "dynamical sign problem."
  • Semiclassical Boltzmann equations fail when electron-phonon interactions are strong (non-perturbative).
• Gap: There is a lack of a unified, numerically exact method capable of handling multiple electronic bands, high/low-frequency phonons, and both local/nonlocal static and dynamic disorders simultaneously in the thermodynamic limit.

2. Methodology: Diagrammatic Quantum Monte Carlo (DQMC)

The authors propose a new Diagrammatic Quantum Monte Carlo (DQMC) approach formulated in imaginary time to avoid the dynamical sign problem. The method relies on expanding the partition function into Feynman diagrams and performing importance sampling.

Key Theoretical Innovations:

• Unified Hamiltonian in Reciprocal Space:
  The total Hamiltonian is divided into electronic ($\hat{H}_{el}$), phonon ($\hat{H}_{ph}$), electron-phonon interaction ($\hat{H}_{el-ph}$), and static disorder ($\hat{H}_{sd}$) components. Crucially, the static disorder is formulated in reciprocal space using a general expression involving Gaussian random variables ($D_{\mu q}$) and deterministic parameters ($f_{nmkq\mu}$). This allows any Gaussian static disorder with translational symmetric covariance to be represented.
• Generalized Wick's Theorem for Static Disorder:
  A major theoretical contribution is the derivation of a generalized Wick's theorem for static disorder variables. Unlike phonon operators, static disorder variables are complex random numbers. The authors prove that in the thermodynamic limit ($N \to \infty$), the thermal average of products of these variables reduces to a summation over pairwise combinations (similar to standard Wick's theorem), with unpaired terms vanishing as $O(N^{-3/2})$.
• Diagrammatic Expansion:
  The partition function is expanded into an infinite series of interaction vertices. The method utilizes the generalized Wick's theorem to pair disorder variables and standard Wick's theorem for phonons, reducing the problem to a summation over Feynman diagrams.
• Current Operator Formulation:
  The authors derive explicit expressions for the current operator components arising from electron-phonon interactions and static disorder in reciprocal space, enabling the calculation of transport properties like mobility and optical conductivity.
• Numerical Implementation:
  The algorithm performs stochastic sampling over diagram configurations (order $\kappa$, momenta, times, and pairing structures). The numerical cost is independent of system size, allowing calculations in the thermodynamic limit.

3. Key Contributions

1. Unified Framework: Successfully integrates dynamic (electron-phonon) and static (structural) disorders into a single, numerically exact DQMC framework.
2. Generalized Wick's Theorem: Provides the mathematical foundation for treating static disorder as a "coordinate" within the diagrammatic expansion, enabling the handling of correlated disorder.
3. Thermodynamic Limit Capability: By working in reciprocal space and utilizing the generalized Wick's theorem, the method eliminates finite-size effects, a common limitation in real-space simulations.
4. Versatility: The framework is designed to handle multiple electronic bands, arbitrary phonon dispersions (optical and acoustic), and various forms of disorder (local/nonlocal, correlated/uncorrelated).

4. Results and Validation

The method was validated against other numerical techniques (Stochastic Liouville-von Neumann equation - SLN, Polaron theory, Marcus theory, Fermi's Golden Rule) across various parameter regimes:

• Dynamic Disorder:
  • Sign Problem: Local dynamic disorder exhibits no sign problem (average sign $\approx 1$). Nonlocal dynamic disorder shows a mild sign problem at low temperatures but remains manageable (sign $> 0.3$).
  • Coherence: Nonlocal dynamic disorder causes significantly stronger decoherence than local disorder. Interestingly, strong nonlocal disorder can lead to a "resurgence" of next-nearest-neighbor coherence.
• Static Disorder:
  • Anderson Localization: The method successfully captures the exponential localization of ground states in 1D chains with static disorder.
  • Sign Problem: Nonlocal static disorder induces a severe sign problem when the disorder strength $\Delta_{non}/T > 5$, limiting convergence in extreme regimes.
• Transport Properties:
  • Mobility Extraction: By combining DQMC with numerical analytic continuation (Stochastic Optimization Method - SOM), the authors extracted charge carrier mobilities.
  • Regime Validation:
    • Thermally Activated Hopping: DQMC results matched Marcus theory.
    • Phonon-Assisted Transport: Results matched Polaron theory.
    • Nuclear Tunneling: Results matched Fermi's Golden Rule (FGR) calculations, reproducing band-like mobility behavior in organic semiconductors.
  • Current Autocorrelation: The imaginary-time current autocorrelation functions showed expected behaviors (e.g., becoming more concave with increasing disorder strength).

5. Significance and Future Outlook

• Versatile Tool: This approach provides a powerful, first-principles-compatible tool for calculating transport properties in realistic materials where multiple competing mechanisms (localization vs. delocalization, coherent vs. incoherent) are at play.
• Bridging Theory and Experiment: It offers a pathway to explain puzzling experimental observations, such as band-like mobility in organic semiconductors despite strong localization, and long-lived quantum coherence in photosynthetic complexes.
• Scalability: The independence of computational cost from system size makes it suitable for studying large, complex materials (e.g., light-harvesting complexes, organic crystals) that are inaccessible to real-time quantum dynamics.
• Future Application: The authors suggest combining this DQMC framework with Density Functional Theory (DFT) and Wannier interpolation to calculate material-specific parameters, enabling the prediction of mobility in real-world disordered semiconductors.

In conclusion, this paper presents a rigorous, numerically exact method that overcomes the limitations of perturbation theory and finite-size effects, offering a comprehensive solution for understanding charge transport in complex disordered quantum systems.