Development of non-equilibrium Green's functions for use with full interaction in complex systems

Richard A. Lynn,Robert van Leeuwen

DOI: https://doi.org/10.1088/1742-6596/696/1/012020

2016-02-05

Abstract:We present an ongoing development of an existing code for calculating ground-state, steady-state, and transient properties of many-particle systems. The development involves the addition of the full four-index two electron integrals, which allows for the calculation of transport systems, as well as the extension to multi-level electronic systems, such as atomic and molecular systems and other applications. The necessary derivations are shown, along with some preliminary results and a summary of future plans for the code.

Mesoscale and Nanoscale Physics,Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a non - equilibrium Green’s functions (NEGF) method capable of handling full interactions in complex systems. Specifically, the author aims to improve the existing code to calculate the ground - state, steady - state, and transient properties of multi - particle systems. The following are the main objectives and challenges of this research:

1. **Inclusion of complete four - index two - electron integrals**:
- Traditional methods usually use simplified two - index integrals to approximately describe the interaction terms. Although this improves the calculation speed, it sacrifices accuracy.
- The method proposed in this paper introduces complete four - index two - electron integrals, thereby more accurately describing the interactions in multi - particle systems.

2. **Extension to multi - level electron systems**:
- Existing NEGF methods are mainly applicable to single - level or few - level systems, while the goal of this paper is to extend it to multi - level systems, such as atomic and molecular systems.
- This extension makes it possible to simulate more complex quantum transport phenomena, such as real - time switching between different steady states.

3. **Improve the ability to simulate time - dependent many - body phenomena**:
- With the development of molecular electronics, simulating time - dependent many - body phenomena has become increasingly important.
- The limitation of computing resources makes it difficult to accurately simulate real systems, so approximate techniques need to be adopted. The NEGF method proposed in this paper is in principle exact, but in practical applications, the self - energy term needs to be truncated and the Kadanoff–Baym equations need to be solved by numerical methods.

4. **Study specific physical phenomena**:
- **Shot noise**: Through the new implementation method, the shot noise phenomenon in quantum systems can be more accurately simulated. This involves measuring the change of current over time and analyzing its fluctuation characteristics.
- **Anderson impurity model**: Using the complete interaction terms, study the multi - level behavior at a single site and explore more complex systems. In particular, study the differences between the isolated central site and the contact situation, and possible new physical phenomena, such as the broadening of the Kondo peak.

### Formula summary
- **Hamiltonian**:
\[
\hat{H}(t)=\hat{H}_C(t)+\hat{H}_{\text{leads}}(t)+\hat{H}_T-\mu \hat{N}
\]
where:
\[
\hat{H}_C = \sum_{ij}\sum_\sigma h_{ij}(t)\hat{a}^\dagger_{i\sigma}\hat{a}_{j\sigma}+\frac{1}{2}\sum_{ijkl}\sum_{\sigma\sigma'} w_{ijkl}\hat{a}^\dagger_{i\sigma}\hat{a}^\dagger_{j\sigma'}\hat{a}_{k\sigma'}\hat{a}_{l\sigma}
\]

- **Kadanoff–Baym equations**:
\[
\left[i\partial_{t_1}-\hat{h}(t_1)\right]G^{\gtrless}(t_1,t_2)=I^{\gtrless}_1(t_1,t_2)
\]
\[
G^{\gtrless}(t_1,t_2)\left[-i\leftarrow\partial_{t_2}-\hat{h}(t_2)\right]=I^{\gtrless}_2(t_1,t_2)
\]

- **Self - energy approximation**:
- Hartree - Fock self - energy:
\[
\Sigma^{\text{HF}}_{ij}(t)=-i\sum_{kl}G_{kl}(t,t^+)(2w_{iljk}-w_{ilkj})
\]