PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nanoreslettspringer open web sitethis articlemanuscript submissionregistrationjournal front pagespringer open web site
 
Nanoscale Res Lett. 2011; 6(1): 204.
Published online Mar 9, 2011. doi:  10.1186/1556-276X-6-204
PMCID: PMC3211260
Effect of phonons on the ac conductance of molecular junctions
Akiko Ueda,corresponding author1 Ora Entin-Wohlman,1,2 and Amnon Aharony1,2
1Department of Physics, Ben Gurion University, Beer Sheva 84105, Israel
2Tel Aviv University, Tel Aviv 69978, Israel
corresponding authorCorresponding author.
Akiko Ueda: akiko/at/bgu.ac.il; Ora Entin-Wohlman: oraentin/at/bgu.ac.il; Amnon Aharony: aaharony/at/bgu.ac.il
Received August 16, 2010; Accepted March 9, 2011.
Abstract
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.
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 An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i1.gif, where b (c0) and b(An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i2.gif) 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 An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i3.gif, 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
A mathematical equation, expression, or formula.
 Object name is 1556-276X-6-204-i4.gif
(1)
Here, An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i5.gif is the Fourier transform of the two particle Green function,
A mathematical equation, expression, or formula.
 Object name is 1556-276X-6-204-i6.gif
(2)
where An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i7.gif, with An external file that holds a picture, illustration, etc.
Object name is 1556-276X-6-204-i8.gif 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
A mathematical equation, expression, or formula.
 Object name is 1556-276X-6-204-i9.gif
(3)
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 [equivalent] GH + Gex contains the Hartree contribution GH and the exchange term Gex. Figure Figure11 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.
Figure 1
Figure 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 (more ...)
Next, Figure Figure2a2a shows the full ac conductance G as a function of the ac frequency ω, when ε0 - μ = Γ. The solid line in Figure Figure2a2a 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 Figure2b2b 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 Figure3.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.
Figure 2
Figure 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 (more ...)
Figure 3
Figure 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 (more ...)
Conclusion
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 < μ.
Abbreviations
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.
Acknowledgements
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).
  • Park H, Park J, Lim AKL, Anderson EH, Alivisatos AP, MacEuen PL. Nanomechanical oscillations in a single-C60 transistor. Nature (London) 2000;407:57. doi: 10.1038/35024031. [PubMed] [Cross Ref]
  • 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. doi: 10.1103/PhysRevLett.100.196804. [PubMed] [Cross Ref]
  • Mitra A, Aleiner I, Millis AJ. Phonon effects in molecular transistors: Quantal and classical treatment. Phys Rev B. 2004;69:245302. doi: 10.1103/PhysRevB.69.245302. [Cross Ref]
  • Koch J, von Oppen F. Franck-Condon Blockade and Giant Fano Factors in Transport through Single Molecules. Phys Rev Lett. 2005;94:206804. doi: 10.1103/PhysRevLett.94.206804. [PubMed] [Cross Ref]
  • Entin-Wohlman O, Imry Y, Aharony A. Voltage-induced singularities in transport through molecular junctions. Phys Rev B. 2009;80:035417. doi: 10.1103/PhysRevB.80.035417. [Cross Ref]
  • Entin-Wohlman O, Imry Y, Aharony A. Transport through molecular junctions with a nonequilibrium phonon population. Phys Rev B. 2010;81:113408. doi: 10.1103/PhysRevB.81.113408. [Cross Ref]
Articles from Nanoscale Research Letters are provided here courtesy of
Springer