This article is part of the series International Conference on Superlattices, Nanostructures and Nanodevices (ICSNN 2010).

Open Access Nano Express

Effect of phonons on the ac conductance of molecular junctions

Akiko Ueda1*, Ora Entin-Wohlman12 and Amnon Aharony12

Author affiliations

1 Department of Physics, Ben Gurion University, Beer Sheva 84105, Israel

2 Tel Aviv University, Tel Aviv 69978, Israel

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2011, 6:204  doi:10.1186/1556-276X-6-204

The electronic version of this article is the complete one and can be found online at:

Received:16 August 2010
Accepted:9 March 2011
Published:9 March 2011

© 2011 Ueda et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We theoretically examine the effect of a single phonon mode on the structure of the frequency dependence of the ac conductance of molecular junctions, in the linear response regime. The conductance is enhanced (suppressed) by the electron-phonon interaction when the chemical potential is below (above) the energy of the electronic state on the molecule.

PACS numbers: 71.38.-k, 73.21.La, 73.23.-b


Molecular junctions, made of a single molecule (or a few molecules) attached to metal electrodes, seem rather well established experimentally. An interesting property that one can investigate in such systems is the interplay between the electrical and the vibrational degrees of freedom as is manifested in the I-V characteristics [1,2].

To a certain extent, this system can be modeled by a quantum dot with a single effective level ε0, connected to two leads. When electrons pass through the quantum dot, they are coupled to a single phonon mode of frequency ω0. The dc conductance of the system has been investigated theoretically before, leading to some distinct hallmarks of the electron- phonon (e-ph) interaction [3-6]. For example, the Breit-Wigner resonance of the dc linear conductance (as a function of the chemical potential μ, and at very low temperatures) is narrowed down by the e-ph interaction due to the renormalization of the tunnel coupling between the dot and the leads (the Frank-Condon blockade) [4,5]. On the other hand, the e-ph interaction does not lead to subphonon peaks in the linear response conductance when plotted as a function of the chemical potential. In the nonlinear response regime, in particular for voltages exceeding the frequency ω0 of the vibrational mode, the opening of the inelastic channels gives rise to a sharp structure in the I-V characteristics. In this article, we consider the ac linear conductance to examine phonon-induced structures on transport properties when the ac field is present.

Model and calculation method

We consider two reservoirs (L and R), connected via a single level quantum dot. The reservoirs have different chemical potentials, μL = μ+Re[δμLeiωt] and μR = μ+Re[δμReiωt]. When electrons pass through the quantum dot, they are coupled to a single phonon mode of frequency ω0. In its simplest formulation, the Hamiltonian of the electron-phonon (e-ph) interaction can be written as , where b (c0) and b() are the annihilation and the creation operators of phonons (electrons in the dot), and γ is the coupling strength of the e-ph interaction. The broadening of the resonant level on the molecule is given by Γ = ΓL + ΓR, with , where ν is the density of states of the electrons in the leads and tL(R) is the tunneling matrix element coupling the dot to the left (right) lead.

The ac conductance of the system is derived by the Kubo formula. In the linear response regime, the current is given by I = (IL - IR)/2, where


Here, is the Fourier transform of the two particle Green function,


where , with and ck(p) denoting the creation and annihilation operators of an electron of momentum k(p) in the left (right) lead. The ac conductance is then given by


In this article we consider the case of the symmetric tunnel coupling, ΓL = ΓR. We also assume δμL = - δμR = δμ/2. The e-ph interaction is treated by the perturbation expansion, to order γ2. The resulting conductance includes the self-energies stemming from the Hartree and from the exchange terms of the e-ph interaction, while the vertex corrections of the e-ph interaction vanish when the tunnel coupling is symmetric. We also take into account the RPA type dressing of the phonon, resulting from its coupling with electrons in the leads [3].


The total conductance is given by G = G0 + Gint, where G0 is the ac conductance without the e-ph interaction, while Gint GH + Gex contains the Hartree contribution GH and the exchange term Gex. Figure 1 shows the conductance G as a function of ε0 - μ, for a fixed ac frequency ω = 0.5Γ. The solid line indicates G0. The dotted line shows the full conductance G, with γ = 0.3Γ. The peak becomes somewhat narrower, and it is shifted to higher energy, which implies a lower (higher) conductance for ε0 < μ (ε0 > μ). However, no additional peak structure appears.

thumbnailFigure 1. The ac conductance as a function of (ε0 - μ). The ac frequency ω = Γ. ΓL = ΓR and δμL = -δμR. Solid line: without e-ph interaction. Dotted line: γ = 0.3Γ and ω0 = Γ.

Next, Figure 2a shows the full ac conductance G as a function of the ac frequency ω, when ε0 - μ = Γ. The solid line in Figure 2a indicates G0. Two broad peaks appear around ω of order ± 1.5(ε0 - μ). The broken lines show G in the presence of the e-ph interaction with ω0 = 2Γ, ω0 = Γ, or ω0 = 0.5Γ. The e-ph interaction increases the conductance in the region between the original peaks, shifting these peaks to lower |ω|, and decreases it slightly outside this region. Figure 2b indicates the additional conductance due to the e-ph interaction, Gint, for the same parameters. Similar results arise for all positive ε0 - μ. Both GH and Gex show two sharp peaks around ω ~ ± (ε0 - μ) (causing the increase in G and the shift in its peaks), and both decay rather fast outside this region. In addition, Gex also exhibits two negative minima, which generate small 'shoulders' in the total G. For ε0 > μ, Gint is dominated by Gex. The exchange term virtually creates a polaron level in the molecule, which enhances the conductance. The amount of increase is more dominant for lower ω0. The situation reverses for ε0 < μ, as seen in Figure 3. Here, G0 remains as before, but the ac conductance is suppressed by the e-ph interaction. Now Gint is always negative, and is dominated by GH. The Hartree term of the e-ph interaction shifts the energy level in the molecule to lower values, resulting in the suppression of G. The amount of decrease is larger for lower ω0.

thumbnailFigure 2. The ac conductance as a function of the ac frequency ω at ε0 - μ = Γ. (a) The total conductance when ΓL = ΓR and δμL = -δμR. The broken lines indicate the conductance in the presence of e-ph interaction with γ = 0.4Γ. ω0 = 2Γ, or 0.5Γ. The solid line is the 'bare' conductance G0, in the absence of e-ph interaction. (b) The additional conductance due to the e-ph interaction, Gint(ω) = GH(ω) + Gex(ω), for the same parameters as in (a).

thumbnailFigure 3. The conductance as a function of the ac frequency ω at ε0 - μ = -Γ. (a) The total conductance when ΓL = ΓR and δμL = -δμR. The broken lines indicate the conductance in the presence of e-ph interaction with γ = 0.3Γ. ω0 = 2Γ, Γ or 0.5Γ. The solid line is the 'bare' conductance G0 in the absence of e-ph interaction. (b) The additional conductance due to the e-ph interaction, Gint(ω) = GH(ω) + Gex(ω), for the same parameters as in (a).


We have studied the additional effect of the e-ph interaction on the ac conductance of a localized level, representing a molecular junction. The e-ph interaction enhances or suppresses the conductance depending on whether ε0 > μ or ε0 < μ.


e-ph: Electron-phonon.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AU carried out the analytical and numerical calculations of the results and drafted the manuscript. OE conceived of the study. AA participated in numerical calculations. All authors discussed the results and commented and approved the manuscript.


This study was partly supported by the German Federal Ministry of Education and Research (BMBF) within the framework of the German-Israeli project cooperation (DIP), and by the US-Israel Binational Science Foundation (BSF).


  1. Park H, Park J, Lim AKL, Anderson EH, Alivisatos AP, MacEuen PL: Nanomechanical oscillations in a single-C60 transistor.

    Nature (London) 2000, 407:57. Publisher Full Text OpenURL

  2. Tal O, Krieger M, Leerink B, van Ruitenbeek JM: Electron- Vibration Interaction in Single-Molecule Junctions: From Contact to Tunneling Regimes.

    Phys Rev Lett 2008, 100:196804. PubMed Abstract | Publisher Full Text OpenURL

  3. Mitra A, Aleiner I, Millis AJ: Phonon effects in molecular transistors: Quantal and classical treatment.

    Phys Rev B 2004, 69:245302. Publisher Full Text OpenURL

  4. Koch J, von Oppen F: Franck-Condon Blockade and Giant Fano Factors in Transport through Single Molecules.

    Phys Rev Lett 2005, 94:206804. PubMed Abstract | Publisher Full Text OpenURL

  5. Entin-Wohlman O, Imry Y, Aharony A: Voltage-induced singularities in transport through molecular junctions.

    Phys Rev B 2009, 80:035417. Publisher Full Text OpenURL

  6. Entin-Wohlman O, Imry Y, Aharony A: Transport through molecular junctions with a nonequilibrium phonon population.

    Phys Rev B 2010, 81:113408. Publisher Full Text OpenURL