Skip to main content
  • Research article
  • Open access
  • Published:

Evidence of phonon-assisted tunnelling in electrical conduction through DNA molecules


We propose a phonon-assisted tunnelling model for explanation of conductivity dependence on temperature and temperature-dependent I-V characteristics in deoxyribonucleic acid (DNA) molecules. The capability of this model for explanation of conductivity peculiarities in DNA is illustrated by comparison of the temperature dependent I-V data extracted from some articles with tunnelling rate dependences on temperature and field strength computed according to the phonon-assisted tunnelling theory.

PACS Codes: 87.15.-v, 71.38.-k, 73.40.Gk


Conductance properties of DNA have recently attracted a lively interest for theoreticians as well as experimenters [1]. Understanding the charge carriers transfer mechanism along DNA double helix is important for possible applications of DNA molecules in nanoelectronic circuit technology [16].

Direct conductivity measurements have shown a very wide range of conducting properties ranging from no conduction [5, 7, 8] to a good linear conductor [2, 9, 10], while in other experiments semiconducting conductivity behaviour emerges [3, 6, 1114].

The wide range of charge transport behaviour seems to arise from different experimental conditions in which the measurements are carried out. These include the nature of the devices used to measure the conductivity, the sequence and length of the DNA, the type of contacts, the environment in different experiments, etc., all can greatly effects the conductivity of the DNA molecules. For instance, Kasumov et al. [4] have shown that strongly deformed DNA molecules deposited on a substrate, whose thickness is less than half the native thickness of the molecule, are insulating, whereas molecules keeping their native thickness are conducting down to very low temperature with a non-ohmic behaviour characteristics of a one dimensional (1D) conductor.

Extensive experimental and theoretical work over the past decade has led to substantial clarification of charge transport mechanisms in DNA. The dominant mechanisms appear to be short-range quantum mechanical tunnelling [1418] or long-range thermally activated hopping [10, 13, 1924]. But these mechanisms are not capable to explain all the field- and temperature-behaviour of experimental data associated with conduction of the DNA molecules. Indeed, the hopping models confronted with difficulties in explaining the observed strong conductivity dependence on the temperature along λ-DNA double helix at high temperatures and a very week dependence at low temperatures [9]. The tunnelling mechanism was excluded in the case of temperature-dependent results [9, 11].

We affirm that in many cases the temperature-dependent conductivity of DNA molecules could be explained by the quantum mechanical tunnelling theories in which the impact of phonons on tunnelling rate is included. In the event, "pure" tunnelling can be observed at low temperatures when the vibrations modes of the molecule are frozen. At moderate temperatures the input of phonon energy to the process of tunnelling must be taken into account and contemporary phonon-assisted tunnelling theories (PhAT) [2527] realise this.

Recently, it has been shown that the PhAT describes well not only the nonlinear I-V curves, but also the temperature-dependent conductivity in conducting polymers [28, 29]. Therefore, we invoke the PhAT theory to describe some the temperature-dependent experimental data on electrical transport through DNA molecules presented by other authors.

Model and comparison with experimental data

We suggest that the thermoactivated current through the DNA molecules is caused by the charge carriers released from localised states located between HOMO and LUMO levels of DNA ones [11]. In the dc case, the said levels are continuously filled from the electrode. Assuming that the electrons are released from these states due to phonon-assisted tunnelling, we will compare the current (the same as the conductance) dependence on the temperature and field strength with the tunnelling rate dependence on these parameters, computed using the PhAT theory. For this purpose we explore the equation (18) in [27] derived by Makram-Ebeid and Lannoo for the phonon-assisted tunnelling of the electrons from the impurity centre. Taking into consideration the fact that this theory has been evaluated using the Condon approximation, it is more suitable for the molecular structures than other ones. For the tunnelling rate dependence on field strength E and temperature T this theory gives:

W ( E , T ) = p = p o + p o R W e ( ε T + p ω ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGxbWvcqGGOaakcqWGgbGrcqGGSaalcqWGubavcqGGPaqkcqGH9aqpdaaeWbqaaiabdkfasjabdEfaxnaaBaaaleaacqWGLbqzaeqaaOGaeiikaGccciGae8xTdu2aaSbaaSqaaiabdsfaubqabaGccqGHRaWkcqWGWbaCcqWIpecAcqWFjpWDcqGGPaqkcqGGSaalaSqaaiabdchaWjabg2da9iabgkHiTiabdchaWnaaBaaameaacqWGVbWBaeqaaaWcbaGaey4kaSIaemiCaa3aaSbaaWqaaiabd+gaVbqabaaaniabggHiLdaaaa@4CBE@


R = exp ( p ω 2 k B T S  cth ω 2 k B T ) I p ( S sh ( ω / 2 k B T ) , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGsbGucqGH9aqpcyGGLbqzcqGG4baEcqGGWbaCdaqadaqaaKqbaoaalaaabaGaemiCaaNaeS4dHGgcciGae8xYdChabaGaeGOmaiJaem4AaS2aaSbaaeaacqWGcbGqaeqaaiabdsfaubaakiabgkHiTiabdofatjabbccaGiabbogaJjabbsha0jabbIgaOLqbaoaalaaabaGaeS4dHGMae8xYdChabaGaeGOmaiJaem4AaS2aaSbaaeaacqWGcbGqaeqaaiabdsfaubaaaOGaayjkaiaawMcaaiabdMeajnaaBaaaleaacqWGWbaCaeqaaOWaaeWaaKqbagaadaWcaaqaaiabdofatbqaaiabbohaZjabbIgaOjabcIcaOiabl+qiOjab=L8a3jabc+caViabikdaYiabdUgaRnaaBaaabaGaemOqaieabeaacqWGubavaaaakiaawIcacaGLPaaacqGGSaalaaa@5E5E@
W e ( ε T ) = ( 2 e F ( 2 m ε T ) 1 / 2 ) exp [ 4 ( 2 m ) 1 / 2 3 e E ] ε T 3 / 2 . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8bkY=MipeYlH8Hipec8Eeeu0xXdbba9frFj0xb9Lqpepeea0xd9q8qiYRWxGi6xij=hbbc9s8aq0=yqpe0xbbG8A8frFve9Fve9Fj0dmeaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGxbWvdaWgaaWcbaGaemyzaugabeaakiabcIcaOGGaciab=v7aLnaaBaaaleaacqWGubavaeqaaOGaeiykaKIaeyypa0ZaaeWaaKqbagaadaWcaaqaaiabikdaYiabdwgaLjabdAeagbqaaiabcIcaOiabikdaYiabd2gaTjabgEHiQiab=v7aLnaaBaaabaGaemivaqfabeaacqGGPaqkdaahaaqabeaacqaIXaqmcqGGVaWlcqaIYaGmaaaaaaGccaGLOaGaayzkaaGagiyzauMaeiiEaGNaeiiCaa3aamWaaeaacqGHsislcqaI0aanjuaGdaWcaaqaaiabcIcaOiabikdaYiabd2gaTjabgEHiQiabcMcaPmaaCaaabeqaaiabigdaXiabc+caViabikdaYaaaaeaacqaIZaWmcqWGLbqzcqWIpecAcqWGgbGraaaakiaawUfacaGLDbaacqWF1oqzdaqhaaWcbaGaemivaqfabaGaeG4mamJaei4la8IaeGOmaidaaOGaeiOla4caaa@5F10@

Here p o = ε T /ħω, ħω is the phonon energy, εT is the centre depth, I P is the modified Bessel function and S is the Huang-Rhys coupling constant.

In the first instance, the comparison of the I-V data measured in the temperature range from 43 K to 294 K for poly(dA)-poly(dT) DNA molecules (from Fig 2a in [11]) with theoretical W(E,T) dependences is presented in Fig. 1. The W(E,T) dependences were calculated using centre depth ε T to be 0.21 eV (the value is slightly higher than the value of the activation energy estimated in Table 1, which is offered in Ref. [11]), and for the electron effective mass the value of 1.5 m e [30] was used. For the phonon energy, the value of 43 meV was selected. This value is similar to the value of 348 cm-1 used for the calculation of the DNA IR active modes in [31]. The coupling constant S was chosen in order to get the best fit of the experimental data with the calculated dependences on the assumption that the field strength for tunnelling is proportional to the applied voltage. As is seen in Fig. 1, the theoretical W (E,T) dependences fit well with the experimental data for entire range of the measured temperatures. The field strength for theoretical curves ranges from 0.16 MV/m to 500 MV/m, which is close to the field strength estimated from the sample thickness (about 20 nm).

Figure 1
figure 1

Comparison of I – V dependences in poly(dA)-poly(dT) DNA molecules extracted from figure 2(a) in [11] (symbols) with theoretical W (E,T) against E dependences (solid curves) calculated for the same T as in the experiment (from top to bottom) using the following parameters: ε T = 0.21 eV,m* = 1.5 me, ħω = 43 meV and S = 12.

Figure 2
figure 2

Experimental σ (T) against 1/T dependence for poly(dA)-poly(dT) DNA molecules extracted from figure 2 (c) in [11] (symbols) fitted to W (E,T) against 1/T dependence, calculated using the same parameters as in Fig. 1.

The judgment on the carriers transfer mechanism is often carried out considering the conductance dependence on the temperature. The conductance measured by Yoo et al. for poly(dA)-poly(dT) was strongly temperature-dependent around room temperature and slightly temperature-dependent at low temperatures [11]. The authors explain such dependence within small polaron hopping model. We note that the W (E,T) versus E characteristics at temperatures below 100 K are weakly dependent on the temperature, and this feature is in excellent agreement with the experimental observation. This circumstance is also evident in Fig. 2 (solid lines) from the plot of ln W (E,T) versus 1/T calculated using the same parameters as in Fig. 1 and for E = 135 MV/m. The symbols in Fig. 2 represent the experimental data extracted from figure 2(c) in [11]. One can see that the theoretical dependences of the phonon-assisted tunnelling rate are in good agreement with the experimental data.

Tran et al. [9] using the resonant cavity technique studied conductivity and its temperature dependence along the λ-DNA double helix at microwave frequencies. They observed similar behaviour of the conductivity on the temperature as in [11], i.e. a strong temperature dependence of conductivity at high temperatures, whereas at low temperatures the conductivity in λ-DNA exhibits a very weak temperature dependence. They explain the temperature-dependent conductivity suggesting two transport mechanisms, i.e. ionic conduction at low temperatures and temperature driven hopping transport processes at high temperatures. The underlying physics of the weak temperature dependence at low temperatures was not understood. The calculation in [32] has also shown that transport through DNA does not have a purely hopping character.

In Fig. 3 the experimental data extracted from Fig. 3 in [9] for λ-DNA (symbols) fitted to the PhAT rate dependences on 1/T for the λ-DNA in buffer (solid curve) and for the dray λ-DNA (dashed curve), are depicted. In this case the theoretical W(T) dependences reflect also the experimental data well.

Figure 3
figure 3

Experimental σ (T) against 1/T dependences for the λ-DNA in buffer and for the dray λ-DNA extracted from figure 3 in [9] (symbols) fitted to theoretical W (E,T) against 1/T dependences calculated for the following parameters: ε T = 0.21 eV, E = 290 MV/m, S = 12 (solid line) and for parameters ε T = 0.17 eV, E = 315 MV/m, S = 8 (dashed line).


In conclusion, the PhAT model is able to explain the peculiarities of field- and temperature-dependent current observed in DNA molecules in a wide region of the electric field strength. A strong temperature dependence of conductivity observed at high temperatures and a very weak temperature dependence at low temperatures of DNA molecules is comprehensible in the framework of the PhAT model. It is worth to note that the W(E,T) dependences at both low and high temperatures are calculated using the same set of parameters, i.e. ε T , m*, T, E, ħω and S. From these ones only the Huang-Rhys coupling constant S is the fitting parameter estimated from the best fitting of the experimental data and theory. Other parameters are known from experiments or from literary sources. An advantage of the PhAT model over the other models used is the possibility to describe the behaviour of the I-V data measured at different temperatures with the same set of parameters characterizing the material.

Thus, the PhAT mechanism, in some cases, could be dominant in the conductance of the DNA molecules.


  1. Endres RG, Cox DL, Singh RRP: Rev Mod Phys. 2004, 76: 195-214. 10.1103/RevModPhys.76.195.

    Article  ADS  Google Scholar 

  2. Fink H-W, Schönenberger Ch: Nature (London). 1999, 398: 407-410. 10.1038/18855.

    Article  ADS  Google Scholar 

  3. Porath D, Bezryadin A, de Vries S, Dekker C: Nature (London). 2000, 403: 635-638. 10.1038/35001029.

    Article  ADS  Google Scholar 

  4. Kasumov AYu, Klinov DV, Roche P-E, Guéron S, Bouchiat H: Appl Phys Lett. 2004, 84: 1007-1009. 10.1063/1.1644909.

    Article  ADS  Google Scholar 

  5. Storm AJ, van Noort J, de Vries S, Dekker C: Appl Phys Lett. 2001, 79: 3881-3883. 10.1063/1.1421086.

    Article  ADS  Google Scholar 

  6. Cai L, Tabata H, Kawai T: Appl Phys Lett. 2000, 77: 3105-3106. 10.1063/1.1323546.

    Article  ADS  Google Scholar 

  7. de Pablo PJ, Moreno-Herrero F, Colchero J, Gómez Herrero J, Herrero P, Baró AM, Ordejón P, Soler JM, Artacho E: Phys Rev Lett. 2000, 85: 4992-4995. 10.1103/PhysRevLett.85.4992.

    Article  ADS  Google Scholar 

  8. Zhang Y, Austin RH, Kraeft J, Cox EC, Ong NP: Phys Rev Lett. 2002, 89: 198102-10.1103/PhysRevLett.89.198102.

    Article  ADS  Google Scholar 

  9. Tran P, Alavi B, Gruner G: Phys Rev Lett. 2000, 85: 1564-1567. 10.1103/PhysRevLett.85.1564.

    Article  ADS  Google Scholar 

  10. Yu ZG, Song X: Phys Rev Leet. 2001, 86: 6018-6021. 10.1103/PhysRevLett.86.6018.

    Article  ADS  Google Scholar 

  11. Yoo K-H, Ha DH, Lee J-O, Park JW, Kim J, Kim JJ, Lee H-Y, Kawai T, Choi HY: Phys Rev Lett. 2001, 87: 198102-10.1103/PhysRevLett.87.198102.

    Article  ADS  Google Scholar 

  12. Rakitin A, Aich P, Papadopoulos C, Kobzar Y, Vedeneev AS, Lee JS, Xu JM: Phys Rev Lett. 2001, 86: 3670-3673. 10.1103/PhysRevLett.86.3670.

    Article  ADS  Google Scholar 

  13. Kutnjak Z, Lahajnar G, Filipiè C: Phys Rev E. 2005, 71: 041901-10.1103/PhysRevE.71.041901.

    Article  ADS  Google Scholar 

  14. Nogues C, Cohen SR, Daube S, Apter N, Naaman R: J Phys Chem B. 2006, 110: 8910-8913. 10.1021/jp060870o.

    Article  Google Scholar 

  15. Zhang W, Ulloa SE: Phys Rev B. 2004, 69: 153203-10.1103/PhysRevB.69.153203.

    Article  ADS  Google Scholar 

  16. Zikic R, Kristiæ PS, Zhang X-G, Fuentes-Cabrera M, Wells J, Zhao X: Phys Rev E. 2006, 74: 011919-10.1103/PhysRevE.74.011919.

    Article  ADS  Google Scholar 

  17. Lagerquist J, Zwolak M, Di Ventra M: Phys Rev E. 2007, 76: 013901-10.1103/PhysRevE.76.013901.

    Article  ADS  Google Scholar 

  18. Roche S, Bicout D, Maciá E, Kats E: Phys Rev Lett. 2003, 91: 109901-10.1103/PhysRevLett.91.109901. erratum: Phys Rev Lett 2004, 92:109901.

    Article  Google Scholar 

  19. Henderson PT, Jones D, Hampikian G, Kan Y, Schuster G: Proc Natl Acad Sci USA. 1999, 96: 8353-8358. 10.1073/pnas.96.15.8353.

    Article  ADS  Google Scholar 

  20. Jortner J, Bixon M, Voityuk AA, Rösch N: J Phys Chem. 2002, 106: 7599-7606.

    Article  Google Scholar 

  21. Bixon M, Jortner J: Chem Phys. 2005, 319: 273-282. 10.1016/j.chemphys.2005.05.013.

    Article  ADS  Google Scholar 

  22. Gutiérrez R, Mandal S, Cuniberti G: Phys Rev B. 2005, 71: 235116-10.1103/PhysRevB.71.235116.

    Article  ADS  Google Scholar 

  23. Klotsa D, Römer RA, Turner MS: Biophys Jour. 2005, 89: 2187-2198. 10.1529/biophysj.105.064014.

    Article  Google Scholar 

  24. Ma S, Xu H, Liu X, Li M: Phycica B. 2007, 391: 98-102. 10.1016/j.physb.2006.09.002.

    Article  ADS  Google Scholar 

  25. Keldysh L: Zh Eksp Teor Fiz [Sov Phys JETP]. 1958, 34: 962-968.

    Google Scholar 

  26. Kiveris A, Kudžmauskas Š, Pipinys P: Phys Status Solidi A. 1976, 37: 321-327. 10.1002/pssa.2210370140.

    Article  ADS  Google Scholar 

  27. Makram-Ebeid S, Lannoo M: Phys Rev B. 1982, 25: 6406-6424. 10.1103/PhysRevB.25.6406.

    Article  ADS  Google Scholar 

  28. Pipinys P, Kiveris A: J Phys: Condens Matter. 2005, 17: 4147-4155. 10.1088/0953-8984/17/26/013.

    ADS  Google Scholar 

  29. Pipinys P, Kiveris A: Physica B. 2005, 370: 1682-171. 10.1016/j.physb.2005.09.007.

    Article  Google Scholar 

  30. Felice RD, Calzolari A, Varsano D, Rubio A: Lect Notes Phys. 2005, 680: 77-116.

    Article  ADS  Google Scholar 

  31. Bykhovskaia M, Gelmont B, Globus T, Woolard DL, Samuels AC, Duong TH, Zakrzewska K: Theor Chem Acc. 2001, 106: 22-27.

    Article  Google Scholar 

  32. Sidowski M, Walczak K: Materials Science-Poland. 2007, 25: 577-582.

    Google Scholar 

Download references

Author information

Authors and Affiliations


Corresponding author

Correspondence to Antanas Kiveris.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License (, which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.

Reprints and permissions

About this article

Cite this article

Pipinys, P., Kiveris, A. Evidence of phonon-assisted tunnelling in electrical conduction through DNA molecules. PMC Phys B 1, 6 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: