DC Current Generation in the Driven Damped Haldane Model

Basic Information

• Paper ID: 2510.14157
• Title: DC Current Generation in the Driven Damped Haldane Model
• Authors: Konrad Koenigsmann¹, Peter Schauss¹'², Gia-Wei Chern¹
• Classification: cond-mat.dis-nn (Condensed Matter Physics - Disordered Systems and Neural Networks)
• Publication Date: October 17, 2025
• Paper Link: https://arxiv.org/abs/2510.14157

¹Department of Physics, University of Virginia; ²Institute for Quantum Physics, University of Hamburg

Abstract

This paper investigates the nonequilibrium steady states of the Haldane model under continuous-wave laser driving and coupled to a thermal reservoir. Dissipative processes are modeled using the Lindblad formalism applicable to quadratic fermionic systems, enabling the study of relaxation dynamics and the emergence of quasi-stationary states. Although conventional topological invariants and bulk-boundary correspondence cannot be directly applied to nonequilibrium settings, the authors introduce an occupation-weighted Chern number to capture the residual topological features of quasi-stationary states. Furthermore, charge transport in the system under simultaneous driving and damping is investigated, demonstrating that finite DC currents are generated through staggered sublattice potential that breaks inversion symmetry.

Research Background and Motivation

Research Problem

The core problem addressed in this paper is: how do topological phenomena evolve in open quantum systems with both external driving and dissipation, and how can topological properties be characterized and understood in nonequilibrium states?

Significance

1. Emerging Research Frontier: The interplay between topological phenomena and nonequilibrium dynamics in open quantum systems represents a rapidly developing frontier in condensed matter physics
2. Technical Application Value: Deeper understanding of open quantum system dynamics will advance technological progress in quantum information science
3. Theoretical Challenge: The applicability of conventional topological theory in nonequilibrium states requires re-examination

Limitations of Existing Approaches

• Previous studies primarily focused on purely driven or purely dissipative systems
• Lack of systematic exploration of system behavior when driving and dissipation coexist
• Traditional topological invariants and bulk-boundary correspondence no longer apply in nonequilibrium settings

Research Motivation

Through numerical investigation of the driven-damped Haldane model, explore the complex interplay between topology, symmetry, and dissipation, providing a new theoretical framework for topological phenomena in open quantum systems.

Core Contributions

1. Introduction of Occupation-Weighted Chern Number: Proposes a new topological invariant to characterize residual topological features of nonequilibrium quasi-stationary states
2. Revelation of DC Current Generation Mechanism: Demonstrates that finite DC currents can be generated through staggered sublattice potential breaking inversion symmetry
3. Establishment of Driving-Dissipation-Topology Relationship: Systematically investigates the effects of driving strength, frequency, and damping strength on topological properties
4. Discovery of Current Reversal Phenomenon: Observes current direction reversal with varying driving strength and provides theoretical explanation

Methodology Details

Task Definition

Investigate the Haldane model under continuous-wave laser driving and thermal reservoir coupling, analyzing how the system relaxes from an initial topological phase to a nonequilibrium quasi-stationary state, and study the topological properties and transport characteristics of this quasi-stationary state.

Model Architecture

Haldane Model Hamiltonian

The basic Haldane model is defined on a honeycomb lattice:

$\hat{H} = -t_1 \sum_{\langle i,j \rangle} (\hat{c}_i^\dagger \hat{c}_j + \text{h.c.}) + t_2 \sum_{\langle\langle i,j \rangle\rangle} (e^{i\phi_{ij}} \hat{c}_i^\dagger \hat{c}_j + \text{h.c.}) + M \left(\sum_{i \in A} \hat{n}_i - \sum_{i \in B} \hat{n}_i\right)$

where:

• $t_1$: nearest-neighbor hopping strength
• $t_2$: next-nearest-neighbor hopping strength
• $\phi_{ij}$: complex phase breaking time-reversal symmetry
• $M$: staggered sublattice potential (breaking inversion symmetry)

Driving Mechanism: Peierls Substitution

AC electric field driving is introduced via Peierls substitution: $t_i \rightarrow t_i \exp\left(i \int_{\mathbf{r}_a}^{\mathbf{r}_b} \mathbf{A}(t) \cdot d\boldsymbol{\ell}\right)$

where $\mathbf{A}(t) = [A \cos(\omega t)] \hat{\mathbf{r}}$ is the periodic vector potential.

Dissipation Modeling: Lindblad Master Equation

System-reservoir coupling is described through the Lindblad master equation: $\frac{d\hat{\rho}}{dt} = -i[\hat{H}, \hat{\rho}] + \mathcal{D}(\hat{\rho})$

Using the relaxation time approximation: $\mathcal{D}(\hat{\rho}) = 2\gamma[\hat{\rho} - \hat{\rho}_{\text{eq}}]$

where $\gamma = 1/(2\tau)$ is the damping strength.

Technical Innovations

Occupation-Weighted Chern Number

Since the quasi-stationary density matrix is no longer a projection operator ($\hat{\rho}^2 \neq \hat{\rho}$), the conventional TKNN formula fails. The authors propose an occupation-weighted density matrix:

$\hat{\rho}_o = \sum_{\mathbf{k}} \sum_n f_n(\mathbf{k}) \rho_n(\mathbf{k})$

where $f_n(\mathbf{k})$ is the band occupation number calculated from the time-averaged density matrix, and the occupation-weighted Chern number $\nu_o$ is then computed using the modified TKNN formula.

Current Operator

The current operator for an individual bond is defined as: $\hat{J}_{ij} = i\boldsymbol{\delta}_{ij}(t_{ij} \hat{c}_i^\dagger \hat{c}_j - \text{h.c.}) = -2\boldsymbol{\delta}_{ij} \text{Im}(t_{ij} \hat{c}_i^\dagger \hat{c}_j)$

Experimental Setup

Parameter Settings

• $t_1 = 1$ (energy unit)
• $t_2 = 0.2$
• $|\phi_{ij}| = \pi/4$
• $M = 0.1$ (set to 0 in some calculations to preserve inversion symmetry)
• Driving direction: positive x-direction

Numerical Methods

• Time Evolution: Fourth-order Runge-Kutta method for solving the Lindblad master equation
• Real-Space Simulation: $13 \times 13$ lattice
• Momentum-Space Simulation: up to $200 \times 200$ lattice
• Initial Conditions: filled valence band, empty conduction band (half-filling)

Evaluation Metrics

1. Occupation-Weighted Chern Number $\nu_o$: characterizes topological properties of quasi-stationary states
2. Unit-Cell Averaged Net Current: characterizes DC current generation
3. Occupation Number Changes: analyzes occupation changes relative to half-filling

Experimental Results

Main Results

Behavior of Occupation-Weighted Chern Number

Figure 2 shows $\nu_o$ as a function of driving parameters:

• Driving Frequency $\omega$: $\nu_o$ decreases with increasing $\omega$, approaching quantized values as $\omega \rightarrow \infty$
• Driving Strength $A$: $\nu_o$ increases monotonically with $A$
• Damping Strength $\gamma$: $\nu_o$ decreases with increasing $\gamma$, approaching quantized values as $\gamma \rightarrow \infty$

Occupation Number Analysis

Figure 3 shows occupation number changes relative to half-filling:

• After driving is turned on, some electrons in the valence band are excited to the conduction band
• Two partially filled bands form, rather than one full and one empty band
• Occupation number changes vanish in the limits $\omega \rightarrow \infty$, $A \rightarrow 0$, and $\gamma \rightarrow \infty$

DC Current Generation

Figure 4 shows the net current at each lattice site in the quasi-stationary state:

• $M = 0$ (preserving inversion symmetry): although site currents are nonzero, the unit-cell averaged net current is $7.12 \times 10^{-9}$ (essentially zero)
• $M = 0.1$ (breaking inversion symmetry): produces significant unit-cell averaged net current of $1.32 \times 10^{-2}$

Current Reversal Phenomenon

Figure 5 shows the x and y components of unit-cell averaged net current as functions of driving strength $A$:

• Current direction reversal is observed near $A \approx 3$
• This can be explained through the effective hopping strength $t_{\text{eff}} \sim t_i J_0(\alpha)$ in the high-frequency approximation, where $J_0$ is the zeroth-order Bessel function

Physical Mechanism Analysis

DC Current Generation Mechanism

1. Symmetry Requirements: Simultaneous breaking of time-reversal and inversion symmetry is required
2. Physical Picture: Staggered sublattice potential makes nearest-neighbor hoppings A→B and B→A inequivalent
3. Charge Flow: When the driving field modifies occupation numbers, charge preferentially flows from one sublattice to the other
4. Effective Interpretation: Can be understood as the motion of polarization between two sublattices (shift current)

Related Work

Main Research Directions

1. Open Quantum Systems: nonequilibrium steady states, dissipative state preparation, many-body dissipative physics
2. Topological Driven Systems: Floquet topological phases, optically-driven topological materials
3. Quantum Transport: quantized transport in topological systems, edge state currents

Advantages of This Work

• First systematic study of the Haldane model with simultaneous driving and dissipation
• Proposes topological invariants applicable to nonequilibrium states
• Reveals complex interplay between symmetry, topology, and dissipation

Conclusions and Discussion

Main Conclusions

1. Preservation of Topological Features: Although quasi-stationary states are topologically trivial, residual topological features can still be captured through the occupation-weighted Chern number
2. Conditions for DC Current: Breaking inversion symmetry is a necessary condition for generating net DC current
3. Parameter Sensitivity: Current magnitude and direction are highly sensitive to driving strength
4. Theoretical Consistency: Current reversal phenomenon can be explained through Bessel function zeros

Limitations

1. Parameter Range: Study is limited to specific parameter regions, not covering the complete phase diagram
2. System Size: Real-space calculations are limited to relatively small lattices
3. Boundary Effects: Does not consider the interaction between edge states and bulk DC current at interfaces
4. Approximation Validity: The applicable range of high-frequency approximation requires further verification

Future Directions

1. Interface Studies: Introduce interfaces to study the interplay between edge currents and bulk DC currents
2. Phase Diagram Mapping: Extend parameter space for systematic phase diagram studies
3. Experimental Realization: Explore experimental implementation in cold atom or solid-state systems
4. Strong Correlation Effects: Consider the effects of electron-electron interactions

In-Depth Evaluation

Strengths

1. Methodological Innovation: Occupation-weighted Chern number provides a new tool for characterizing nonequilibrium topological states
2. Physical Insight: Deeply reveals the relationship between symmetry breaking and DC current generation
3. Numerical Rigor: Multiple numerical methods are employed to verify result reliability
4. Theoretical Completeness: Forms a complete closed loop from phenomenological observation to theoretical explanation

Weaknesses

1. Theoretical Depth: The deeper physical meaning of occupation-weighted Chern number awaits further clarification
2. Experimental Feasibility: Lacks discussion of experimental implementation feasibility
3. Universality: The applicability of the method to other topological models has not been verified
4. Quantum Effects: Insufficient consideration of quantum coherence effects in long-time evolution

Impact

1. Theoretical Contribution: Provides a new theoretical framework for nonequilibrium topological physics
2. Methodological Value: Occupation-weighted Chern number may become a standard tool for studying driven-dissipative systems
3. Application Prospects: Has potential applications in quantum device design and topological quantum computing

Applicable Scenarios

1. Cold Atom Systems: Artificial topological materials in optical lattices
2. Solid-State Devices: Optically-driven topological phase transitions in 2D materials
3. Quantum Simulators: Digital quantum simulation of open quantum systems
4. Theoretical Research: Study of topological properties in nonequilibrium many-body systems

References

The paper cites 61 relevant references covering multiple fields including open quantum systems, topological physics, and quantum transport, providing a solid theoretical foundation for the research. Key references include Haldane's original work 16, Lindblad master equation theory 17, TKNN formula 21, and other classical literature.

Overall Assessment: This is a high-quality theoretical physics paper that makes important contributions to the frontier field of nonequilibrium topological physics. The paper features methodological innovation, in-depth analysis, and reliable results, providing new perspectives and tools for understanding topological phenomena in open quantum systems.