Molecular junctions in the Coulomb blockade regime: rectification and nesting

Bo Song,Dmitry A. Ryndyk,Gianaurelio Cuniberti

DOI: https://doi.org/10.1103/PhysRevB.76.045408

2007-07-18

Abstract:Quantum transport through single molecules is very sensitive to the strength of the molecule-electrode contact. Here, we investigate the behavior of a model molecular junction weakly coupled to external electrodes in the case where charging effects do play an important role (Coulomb blockade regime). As a minimal model we consider a molecular junction with two spatially separated donor and acceptor sites. Depending on their mutual coupling to the electrodes, the resulting transport observables show well defined features such as rectification effects in the I-V characteristics and nesting of the stability diagrams. To be able to accomplish these results, we have developed a theory which allows to explore the charging regime via the nonequilibrium Green function formalism parallel to the widely used master equation technique. Our results, beyond their experimental relevance, offer a transparent framework for the systematic and modular inclusion of a richer physical phenomenology.

Mesoscale and Nanoscale Physics,Materials Science

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **quantum transport properties in molecular junctions due to charging effects in the Coulomb Blockade (CB) regime**. Specifically, the research focuses on:

1. **Behavior of molecular junctions in the case of weak coupling**: The paper explores the influence of charging effects on transport properties when there is weak coupling between molecules and electrodes. This weak - coupling situation poses a challenge to existing experimental techniques.

2. **Two - site molecular junction model**: The author uses a molecular junction model containing two spatially separated donor and acceptor sites to study its transport behavior under different coupling conditions. This includes phenomena such as rectification effects and nesting of the stability diagrams.

3. **Application of the Non - Equilibrium Green's Function (NEGF) method**: In order to describe the charging effects, the author develops a method based on the non - equilibrium Green's function. This method can handle single - electron charging effects and can systematically introduce more abundant physical phenomena, such as electron - vibron coupling.

4. **Improved description of transport properties**: By introducing a new Ansatz (assumption), the author is able to more accurately describe the transport properties of molecular junctions in the case of weak coupling, especially the stability diagrams at small bias voltages and gate voltages.

### Main research content

- **Green's functions of single - site junctions (SSJ) and two - site junctions (DSJ)**: The author derives the non - equilibrium Green's functions of single - site junctions and two - site junctions in the Coulomb blockade regime and calculates the corresponding transport observables.

- **Two - site junctions in serial and parallel configurations**: For two - site junctions in serial configuration (such as donor/acceptor - type rectifiers), the author shows the diode characteristics of the I - V curves and predicts the conductance peak splitting and double splitting of the current maximum due to charging effects. For two - site junctions in parallel configuration, the author observes the nesting feature of the stability diagrams.

### Conclusion

The method proposed in the paper not only theoretically provides a powerful framework for describing the transport properties of single molecules under strong Coulomb interactions, but also its results are consistent with those of experiments and the master equation method. In addition, this method provides a transparent basis for systematically introducing more physical phenomena (such as electron - vibron coupling).

### Formula summary

The main formulas involved in the paper include:

- **Green's function equation for single - site junctions**:
\[
(\omega-\epsilon_\sigma-\Sigma^{r / a}_\sigma)G^{r / a}_{\sigma,\tau}=\delta_{\sigma,\tau}+U G^{(2)r / a}_{\sigma,\tau}
\]
\[
(\omega-\epsilon_\sigma - U-\Sigma^{r / a}_\sigma)G^{(2)r / a}_{\sigma,\tau}=\langle n_{\bar{\sigma}}\rangle\delta_{\sigma,\tau}
\]

- **Improved lesser Green function**:
\[
G^<(\omega)=G^<_0+G^r_0\Sigma^r G^<+G^<_0\Sigma^a G^a
\]

- **Current expression**:
\[
J=\frac{ie}{2\hbar}\int d\epsilon\frac{1}{2\pi}\text{Tr}\{(\Gamma_L-\Gamma_R)G^{(U,t)<}+[f_L(\omega)\Gamma_L - f_R(\omega)\Gamma_R](G^{(U,t)r}