Bipul Pandey,Peter B. Littlewood
Abstract:Previously, we introduced a method for systematically correcting a quasiparticle green's function via a power series expansion. Here we present an ODE based formalisms of power series correction that goes beyond the cumulant approximation and implement it to 1D Holstein chain for a wide range of coupling strengths in a scalable and inexpensive fashion at both zero and finite temperature. We show that this first differential formalism of the power series is both qualitatively and quantitatively in excellent agreement with exact diagonalization results on 1D Holstein chain with dispersive bosons for a large range of electron-boson coupling strength. We investigate carrier mass growth rate and carrier energy displacement across a wide range of coupling strength. Finally, we present a heuristic argument which predicts most of the rich satellite structure without explicit calculation.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of how to correct the single - particle Green's function by power series beyond the cumulant approximation in the 1D Holstein chain model. Specifically, the paper proposes a power - series correction form based on ordinary differential equations (ODEs) to deal with the electron - boson interaction problem under different coupling strengths. This method is applicable not only to the zero - temperature case but also to the finite - temperature case, and can be extended to macroscopic systems with low computational cost.
#### Main objectives of the paper:
1. **Propose a new method**: Introduce a power - series correction method based on ODEs to go beyond the traditional cumulant approximation.
2. **Scope of application**: Apply this method to the 1D Holstein chain model, covering different cases from weak coupling to strong coupling.
3. **Verify accuracy**: Verify the accuracy and effectiveness of this method under different coupling strengths by comparing with the exact diagonalization results.
4. **Explore physical phenomena**: Study the changes of carrier effective mass and energy shift with coupling strength, and explain the formation mechanism of satellite structures.
#### Specific background and problem description:
- **Limitations of traditional methods**: Although the traditional cumulant approximation method can handle weak - coupling cases, it may produce divergent or inaccurate results in strong - coupling cases.
- **Advantages of the new method**: The proposed power - series correction method can not only handle strong - coupling cases but also has an advantage in computational efficiency, especially for large - scale systems.
- **Physical significance**: This method helps to understand the influence of electron - boson coupling on material properties, such as carrier effective mass, energy shift, and lifetime.
#### Explanation of mathematical formulas:
The key formulas involved in the paper include:
- **1D Holstein Hamiltonian**:
\[
H = H_e + H_{\text{bos}} + H_{e-\text{bos}}
\]
where,
\[
H_e = \epsilon_0 \sum_i c_i^\dagger c_i - t_{\text{el}} \sum_i (c_i^\dagger c_{i+1} + c_{i+1}^\dagger c_i)
\]
\[
H_{\text{bos}} = \tilde{\omega} \sum_j b_j^\dagger b_j + t_d \sum_j (b_j^\dagger b_{j+1} + b_{j+1}^\dagger b_j)
\]
\[
H_{e-\text{bos}} = g \sum_i c_i^\dagger c_i (b_i^\dagger + b_i)
\]
- **Green's function and self - energy**:
\[
G(k,t) = G_0(k,t) \sum_{n = 0}^\infty \left(\frac{g^2}{N}\right)^n C_n(k,t)
\]
\[
-i \Sigma(k,t) = \frac{g^2}{N} \sum_q D_0(q,t) G(k - q,t)
\]
- **Differential form of power - series correction**:
\[
\frac{dP_k(t)}{dt} = (-i \frac{g^2}{N}) \left[ \int_0^t d\tau e^{i \epsilon_k \tau} \Sigma_{k,q = 0}(\tau) P_k(t) + \sum_{q \neq 0} \int_0^t d\tau e^{i \epsilon_k \tau} \Sigma_{k,q}(\tau) P_{k - q}(\tau) P_k(t-\tau) \right]
\]
Through these formulas, the paper shows how to use power - series correction.