Photogalvanic currents from first-principles real-time density-matrix dynamics

This paper presents a first-principles real-time density-matrix formalism that predicts photogalvanic currents in both transient and steady regimes by fully incorporating quantum scatterings mediated by phonons and photons, thereby elucidating the roles of shift and injection currents and their connection to quantum-geometric quantities in non-centrosymmetric materials.

The Gist

Imagine you have a special kind of solar panel, but instead of just making electricity when the sun shines on it, it can generate a current simply because the material itself is "lopsided" (non-centrosymmetric). This is called the Photogalvanic Effect. It's like a one-way street for electrons: when light hits the material, the electrons are forced to run in a specific direction, creating a direct current without needing any batteries or external wires.

For a long time, scientists tried to predict how strong this current would be, but their calculations were like looking at a race from a helicopter and only seeing the runners at the starting line. They knew how the light kicked the electrons into motion (photo-excitation), but they ignored what happened after the race started. They forgot about the obstacles, the fatigue, and the interactions between the runners and the track itself.

This paper introduces a new, super-detailed simulation tool called FPDMD (First-Principles Real-Time Density Matrix Dynamics). Think of this tool as a high-speed, 4K slow-motion camera that doesn't just watch the start of the race, but tracks every single electron, every collision, and every step they take in real-time.

Here is a breakdown of what they discovered using this new "camera":

1. The "Bounciness" of the Track (Phonons Matter!)

In the old models, scientists thought the current was generated purely by the initial "kick" from the light. They assumed the electrons just ran straight to the finish line.

The New Discovery: The authors found that the "track" itself is alive. The atoms in the material are constantly vibrating (these vibrations are called phonons). When an electron runs, it bumps into these vibrating atoms.

• The Analogy: Imagine running on a trampoline. If you just jump (the light), you go up. But if the trampoline fabric is also bouncing and pushing you (the phonons), it changes your path and speed significantly.
• The Result: In materials like Barium Titanate (a common piezoelectric used in sensors), these "bumps" with vibrating atoms actually contribute just as much to the electricity as the initial light kick did. Ignoring them was like trying to predict a runner's time without accounting for wind resistance or a bumpy road.

2. The Two Types of Currents

The paper distinguishes between two ways light creates electricity, depending on the "color" of the light (its polarization):

• The "Shift" Current (Linear Light): Imagine a crowd of people being pushed by a wind blowing from the left. They all shuffle to the right. This is the Shift Current. The team found that the "bumps" from the vibrating atoms (phonons) help push this crowd even harder than previously thought.
• The "Injection" Current (Circular Light): Imagine a spinning top. If you spin it clockwise, it moves one way; counter-clockwise, it moves the other. This is the Injection Current. The new theory calculates exactly how fast these "spinning tops" (electrons) slow down or speed up due to collisions with the vibrating atoms, giving a much more accurate prediction of the current.

3. The Mystery of the "Bipolar" Pulse

Scientists have been confused by experiments using ultra-fast laser pulses (THz spectroscopy). Sometimes, the electricity generated flips direction instantly: it goes positive, then immediately negative, like a heartbeat (bipolar). Other times, it just goes one way (unipolar).

The Explanation: The new simulation shows this is a relay race between two different forces:

1. The Light Kick: Happens almost instantly (in a fraction of a femtosecond).
2. The Phonon Push: Happens slightly slower (about 25 femtoseconds later).

If the "Light Kick" pushes electrons one way, and the "Phonon Push" (the track vibrations) pushes them the other way, you get that confusing flip-flop (bipolar) signal. If they push in the same direction, you get a steady signal. The paper explains that the "flip" happens because the two forces take different amounts of time to kick in.

Why Does This Matter?

• Better Solar Cells: By understanding exactly how electrons interact with the vibrating atoms, we can design materials that harvest sunlight more efficiently, even from low-energy light (like infrared).
• Super-Fast Sensors: This helps us build detectors that can sense light polarization (like 3D glasses) or detect tiny changes in magnetic fields.
• Solving the Puzzle: It clears up years of confusion about why some experiments didn't match the old theories. The old theories were missing the "vibrating track" part of the story.

In a nutshell: This paper built a time-machine simulation that shows us that generating electricity from light isn't just about the light hitting the material; it's a complex dance between the light, the electrons, and the vibrating atoms of the material itself. By counting every step of that dance, we can finally predict and design better energy technologies.

Technical Summary

Here is a detailed technical summary of the paper "Photogalvanic currents from first-principles real-time density-matrix dynamics":

1. Problem Statement

The photogalvanic effect (PGE), or bulk photovoltaic effect, generates a direct current (DC) in non-centrosymmetric materials under illumination without external fields. While recent research has linked PGE to quantum geometric properties (Berry curvature, quantum metric) in topological materials, existing ab initio studies face significant limitations:

• Focus on Photo-excitation: Previous theories largely focused only on the initial photo-excitation process, neglecting subsequent scattering events.
• Inadequate Scattering Models: Realistic materials involve complex scattering mediated by bosons (photons and phonons), including intra/interband relaxation and electron-hole recombination. Common approximations, such as a single state-independent relaxation time ($\tau_0$), fail to self-consistently describe these mechanisms.
• Transient Dynamics: There is a lack of rigorous theoretical frameworks to explain transient photogalvanic currents (observed via THz emission), particularly the origin of "bipolar" current responses (currents that reverse direction over time).
• Misconception on Shift Current: It is commonly assumed that the shift current (under linearly polarized light) is purely a photo-excitation effect, ignoring the potential contribution of phonon-mediated scattering.

2. Methodology: First-Principles Real-Time Density-Matrix Dynamics (FPDMD)

The authors developed a comprehensive theoretical framework based on Quantum Master Equations (QME) to simulate electron-boson dynamics in real time.

• Hamiltonian: The system is modeled using a total Hamiltonian $H_T = H_e^0 + H_b^0 + H'$, where $H_e^0$ is the unperturbed electronic Hamiltonian (derived from Kohn-Sham Density Functional Theory), $H_b^0$ is the boson (phonon/photon) Hamiltonian, and $H'$ represents the electron-boson interaction.
• Interaction Terms:
  • Electron-Photon: Treated in the velocity gauge with the dipole approximation.
  • Electron-Phonon: Calculated using first-order matrix elements derived from Density-Functional Perturbation Theory (DFPT).
• Dynamics: The time evolution of the reduced density matrix ($\rho$) is governed by the QME:
  $$\frac{d\rho}{dt} = -\frac{i}{\hbar}[H_e^0, \rho] + \left.\frac{d\rho}{dt}\right|_{exc} + \left.\frac{d\rho}{dt}\right|_{ph} + \left.\frac{d\rho}{dt}\right|_{rec}$$
  This explicitly includes terms for photo-excitation, electron-phonon scattering, and electron-hole recombination.
• Current Calculation: The charge current density $J(t)$ is calculated as the trace of the density matrix multiplied by the velocity operator. The current is decomposed into band-diagonal (injection current) and band-off-diagonal (shift current) components.
• Self-Consistency: Unlike perturbative approaches, this method solves for non-equilibrium steady-state electron distributions and state-dependent relaxation times ($\tau_{ks}$) self-consistently.

3. Key Contributions

• Unified Framework: The FPDMD formalism allows for the calculation of photogalvanic currents in all time regimes (transient and steady-state) within a single first-principles framework.
• Explicit Phonon Scattering: The theory explicitly encodes all quantum scatterings mediated by bosons, moving beyond the "photo-excitation only" paradigm.
• Quantum Geometric Connection: The formulation elucidates the connection between macroscopic currents and fundamental quantum-geometric quantities (Berry curvature and quantum metric) in the presence of scattering.
• Transient Resolution: It provides femtosecond-resolution insights into transient dynamics, explaining the crossover between different current mechanisms.

4. Key Results (Case Study: BaTiO$_3$)

The authors applied the FPDMD formalism to tetragonal Barium Titanate (BaTiO$_3$), a prototypical piezoelectric.

• Linear Photogalvanic Effect (Shift Current):

  • Phonon Contribution: Contrary to the belief that shift current is purely an excitation effect, the study found that phonon-mediated scattering contributes significantly to the shift current, particularly at higher light frequencies.
  • Mechanism: The phonon shift current arises from the scattering of electrons from the excitation surface to the band edge. The magnitude of this contribution is comparable to, and at high frequencies dominates, the excitation shift current.
  • Quantum Geometry: The phonon shift vector is linked to the Berry curvature ($\Omega_k$). The circulation of the Berry curvature in BaTiO$_3$ is tied to its ferroelectricity; inverting the circulation inverts the current direction.
  • Agreement with Experiment: Including the phonon contribution brings the theoretical total photocurrent into much better agreement with experimental data across a wider energy range above the band edge.

• Circular Photogalvanic Effect (Injection Current):

  • Developed a self-consistent theory for steady injection current using state-dependent relaxation times ($\tau_{ks}$) derived from the FPDMD electron distribution.
  • The injection current conductivity was shown to be a synthesis of two quantum geometric characteristics: the Berry curvature (on the photo-excitation surface) and the quantum metric (related to intraband electron-phonon scattering).

• Transient Photogalvanic Current:

  • Bipolar vs. Unipolar: The simulation explains the experimentally observed bipolar transient currents ($J(t)$) as a dynamical crossover between two mechanisms with different time lags:
    1. Excitation Shift Current ($J_{exc}$): Responds rapidly ($\tau_{eph} \approx 0.75$ fs).
    2. Phonon Shift Current ($J_{ph}$): Responds with a delay corresponding to energy relaxation ($\tau_{er} \approx 25$ fs).
  • Crossover: If $J_{exc}$ and $J_{ph}$ have opposite signs, the total current is bipolar. If they have the same sign, it is unipolar. This explains the frequency-dependent behavior observed in THz emission spectroscopy.

5. Significance and Impact

• Resolving Controversies: The work resolves confusion regarding the scattering mechanisms contributing to PGE, proving that phonon scattering is not a minor correction but a dominant factor in realistic materials like piezoelectrics.
• Predictive Power: The FPDMD approach is predictive for realistic materials where electron-phonon coupling cannot be ignored, offering a path to design better photodetectors and solar cells (especially for low-energy photons).
• Beyond Charge Currents: The framework is generalizable to other observables, such as spin photocurrents and orbital photocurrents, which are crucial for the emerging fields of spin-optotronics and orbitronics.
• Non-linear Optics: The theory provides a straightforward route to calculate higher-order non-linear optical responses (e.g., Second Harmonic Generation) and characterize topological properties in novel quantum materials.

In summary, this paper establishes a rigorous, first-principles real-time dynamics method that bridges the gap between idealized quantum geometric theories and the complex, scattering-dominated reality of photogalvanic effects in functional materials.