Research article
UDC 519.6
DOI: 10.22363/2658-4670-2019-27-3-217-230
Charge diffusion in homogeneous molecular chains based on the analysis of generalized frequency spectra in the framework of the Holstein model
Dmitry A. Tikhonov1,2, Egor V. Sobolev1,3, Victor D. Lakhno1
1 Institute of Mathematical Problems of Biology Branch of Keldysh Institute of Applied Mathematics of RAS 1 Professor Vitkevich St., Pushchino, Moscow Region 142290, Russian Federation 2 Institute of Theoretical and Experimental Biophysics of RAS 3 Institutskaya St., Pushchino, Moscow Region 142290, Russian Federation 3 European Molecular Biology Laboratory, Hamburg Unit c/o DESY, Building 25A, Notkestrafie 85, 22607 Hamburg, Germany
(received: November 12, 2019; accepted: December 23, 2019)
We analyzed numerically computed velocity autocorrelation functions and generalized frequency spectra of charge distribution in homogeneous DNA sequences at finite temperature. The autocorrelation function and generalized frequency spectrum (frequency-dependent diffusion coefficient) are phenomenologically introduced based on the functional of mean-square displacement of the charge in DNA. The charge transfer in DNA was modeled in the framework of the semi-classical Holstein model. In this model, DNA is represented by a chain of oscillators placed into thermostat at a given temperature that is provided by the additional Langevin term. Correspondence to the real DNA is provided by choice of the force parameters, which are calculated with quantum-chemical methods. We computed the diffusion coefficient for all homogenous DNA chains with respect to the temperature and found a special scaling of independent variables that the temperature dependence of the diffusion coefficient for different homogenous DNA is almost similar. Our calculations suggest that for all the sequences, only one parameter of the system is mainly responsible for the charge kinetics. The character of individual motions contributing to the charge mobility and temperature-dependent regimes of charge distribution is determined.
Key words and phrases: charge transfer, velocity autocorrelation function, generalized frequency spectrum, DNA, Holstein model
1. Introduction
Kinetics of charge transfer in polarized one-dimensional chains at finite temperature is an attractive and pressing theoretical problem by itself [1]—[4].
© Tikhonov D.A., Sobolev E.V., Lakhno V.D., 2019
This work is licensed under a Creative Commons Attribution 4.0 International License http://creativecommons.org/licenses/by/4.0/
Its particular importance is related to the discovery of conducting properties in DNA [5]-[7]. On the first hand, this transfer may determine the biological functions of DNA. On the second hand, this opens up fresh opportunities in nanobioelectronics [8]-[10].
A basic qualitative simulation of charge transfer in DNA is a pioneered Holstein model for describing the charge transfer process in one-dimensional chains [11]. In the semi-classical approximation, DNA nucleotide pairs are modeled by a sequence of unconnected oscillators arranged in a chain [12]. A charge travels along a chain in the strong coupling approximation. Motions of the sites influence the charge propagation, and conversely, the charge probability density affects the charge motion.
The Holstein model is rather simple and can hardly describe the charge propagation in DNA realistically. Nevertheless, it is nonlinear and demonstrates different regimes. The study of this model can give qualitative insight into the nature of charge transfer in one-dimensional chains, including biological polymers.
The Holstein model is thoroughly investigated. Numerous papers are devoted to the analytical and numerical analysis of the original quantum model and its semi-classical approximation. However, the charge transfer kinetics at finite temperature is still to be understood. Of importance in this respect is the diffusion coefficient of a charge. Earlier, we studied its temperature dependence inhomogeneous chains of different nucleotide structures and found the scales which mainly determine differences in the charge diffusion in different chains [13].
This work continues the earlier started analysis. In order to better appreciate the kinetics of an added charge propagation along an infinite homogeneous chain, find all the motions which contribute into its mobility and determine the diffusion coefficient, we calculate and analyze autocorrelation functions of the charge propagation velocities and a relevant generalized frequency spectrum which is a frequency-dependent diffusion coefficient. We show that the earlier suggested scale is natural for the model, and the diffusion coefficients differ only in the low-frequency range for different sequences. In the limit of high and low temperatures, the nucleotide structure is irrelevant. Analysis of autocorrelation functions suggests that the charge propagation is contributed by dissipation and reflection at long distances as well as reflections from neighboring sites. As the temperature grows, these motions arise, coexist, and alternate. In the limit of finite temperature, the charge diffusion ceases.
A detailed straightforward analysis of the Holstein model kinetics at finite temperature requires direct numerical simulation concerned with a vast amount of computations. Determinate time dependencies should be identified from chaotic trajectories by averaging a large number of such trajectories. An essential tool to accomplish these ends is our unique method [14], [15], which enables speeding up calculations by three orders of magnitude as compared to difference schemes. The approach implies combining Magnus expansion methods used to reveal a charge evolution with a stochastic difference scheme applied to calculate motion trajectories of classical sites. This computational method enabled us to model reasonably sized samples so that to calculate a velocity autocorrelation function smooth enough to be subsequently analyzed.
2. Theory
2.1. Discrete Holstein model in semi-classical approximation.
Modelling a thermostat
The charge transfer in a one-dimensional double-strand chain will be considered in the framework of the Holstein model. An electron travels along a chain whose sites are modeled by unconnected classical harmonical oscillators. The model Hamiltonian reads:
H = ^»nm ln)(ml+a^qn + + 2, (1)
n,m n n n
where unm are matrix elements or energies of electron transitions between the sites n and m, a is a constant of electron coupling with displacement of the n-th site qn, M and K are mass and elasticity coefficient of the oscillators, respectively.
We will deal with homogeneous chains in the nearest neighbor approximation:
_ c
If we seek a solution corresponding to Hamiltonian (1), in the form
№)=^bn (t)ln),
n
then the relevant Schroedinger equation for a charge and the equation for the sites motion at finite temperature is written as:
db
ih~dt =u(bn-i + b^+i ) + ^bn, (2a)
Md^f = —Kqn - a lbnI2 + Vd-^ + An (t). (2b)
To model the temperature of the surrounding medium we place the chain in a Langevin thermostat. For this purpose equation for classical sites (2b) will include a term with friction (where r is a friction coefficient) and a random force An (t) with a normal distribution and the autocorrelation function:
(t) Am (f)) = 2TkBT5ntmS (t-f), where T is temperature and kB is the Bolzmann constant.
2.2. Change-over to dimensionless parameters. Scaling of autocorrelation functions
If we neglect the influence of quantum equation (2a) on classical equation (2b), the latter becomes a Langevin-type equation, and the solutions do not depend on the charge distribution over the chain. In paper [16], we called this variant of the system an "adiabatic approximation". In the complete system, the quantum subsystem affects classical displacements but to a limited extent.
Therefore we explicitly separate a temperature multiplier in front of them, proceeding from the expression for dispersion of the oscillator coordinate qn with the elasticity coefficient K:
(il ) = f.
This expression follows from virial theorem and determines the scale of variation of a classical oscillator amplitude as a function of temperature:
k^T
Qn = \l
The time scale specifies the matrix element t = rh/v. Let us rewrite equations (2) with regard to the chosen scales in dimensionless form:
= bn_i + bn+1 + dun bn, (3a)
^T = -U0Un - ^ K|2 + + Zn (r), (3b)
(Zn (r)Zn (T')) = 21ulS(T-T'),
where
h r 2 h2 na T
Uo=VVM> 1 = -VM> X = a^3M, -
The quantum equation contains only one dimensionless parameter 0, which is the amplitude multiplier of classical displacements. It is mainly responsible for the influence of the classical subsystem on the quantum one, depending on temperature. This fact suggests that the electron distribution kinetics may mainly depend on 6 for all chain sequences. In particular, the influence of the classical system becomes negligible for 0^0, and the dependence on the nucleotide sequence completely disappears. Our results demonstrate that this choice of the scale is appropriate.
Apart from frequency u0 and friction 7, responsible for oscillator motion characteristics, the classical equation involves a parameter x/6, which determines a reverse influence of the quantum subsystem on the classical one. Analysis of x/0 suggests that this influence is neglectable for high temperatures. Besides, we also can neglect it in the case of nucleotide chains with a large matrix element v. If we set a constant x = 0 for these cases, we obtain an adiabatic approximation considered earlier in [16].
2.3. Charge mean-square displacement and generalized frequency
spectrum
The kinetic properties of stochastic system (3) are described by the velocity autocorrelation function and the generalized frequency spectrum. These
functions can be calculated from the density of probability of finding an electron at site n in a chain of length 2N + 1:
P(n,r) = bn (r) b*a (r).
In the semi-classical approximation the root-mean-square displacements in dimensionless form are calculated by the formula:
N
(x2)=$(r) = Y, n2p(n,r)- (4)
n=-N
Then the diffusion coefficient can be found as a slope of the curve £ (t) on its linear segment or as an asymptotic value of the derivative of the displacement with respect to time v(t) = £ (t)' , or via the integral of the second derivative of the displacement with respect to time ^(t) = £(t)" . We will call this coefficient a static one and denote it as D (0):
D(0) = ±J%(r)dT.
The natural generalization of the static diffusion coefficient is a frequency-dependent diffusion coefficient D(fi), which can be found as a generalized frequency spectrum of the second derivative of the displacement function with respect to time:
1 f™
D(4>) = - J ^ (t) cos (2n0r) dr, (5)
2 Jf0
where 0 is dimensionless frequency which is related to frequency in hertz
(/) as:
Function ^ (t) in semi-classical approximation is a velocity autocorrelation function.
3. Results and discussion
At the initial moment a charge was inserted in the center of a homogeneous chain of length 2N+1 with temperature distribution of the degrees of freedom of harmonical oscillators:
dn
K(0) = 6nfi, (0) = N(0,1), Un (0) = N(0,1), (6)
where is Kronecker delta, N (^, a) is a normal distribution with the mean ^ and dispersion a2.
The trajectory of the system motion was calculated from initial states (6), each time with new values of un and un. The calculations were carried out by the method, which combines Magnus expansion with a stochastic difference scheme [14], [15]. For the trajectories to correspond to charge propagation
along an infinite chain, the chain length N was chosen such that during the time of modeling the front of propagation of the charge probability density would not come up to its end.
Then we calculated the second derivative of the root-mean-square displacement with respect to (3), (4) by averaging over samples. After that, we found the frequency spectrum by formula (5) using the fast cosine transform.
We chose the parameters listed in Table 1 as in our previous work [13]. These parameters are typical for DNA.
Table 1
Dimensional parameters of the Holstein model for dsDNA
M=10-21 g K = 0.062 eV/A2 r = 6 • 10-10 g/s a = 0.13 eV/A
uA = 0.030 eV uc = 0.041 eV uG = 0.084 eV uT = 0.158 eV
We repeated the calculations for a set of temperatures and homogeneous dsDNA chains composed of different nucleotides. Table 2 lists dimensionless parameters.
Table 2
Dimensionless parameters for dsDNA as a function of a nucleotide sequence
parameter poly A poly C poly G poly T
0 0.02193 0.01605 0.007837 0.004167
0.01316 0.009631 0.004702 0.0025
0.004219 0.001654 1.925 x 10-4 2.894 x 10-5
A common property of the velocity autocorrelation functions at any temperature is that they are equal to 4 in the small area after an initial time instant. This value corresponds to the velocity of the ballistic motion at zero temperature, which is equal to 2. A comparison of autocorrelation functions for chains composed of different nucleotides suggests that for the same value of 9, the functions differ only in the attenuation asymptotics. The spectra demonstrate differences at low frequencies for finite values of d and, accordingly, different diffusion coefficients D (0) (see Figure 1, dashed lines). Here we use the scale were these differences are insignificant. In particular, they lack in the limit of low and high temperatures.
Earlier in work [13], we carried out a more detailed analysis of the temperature dependencies of the static diffusion coefficient for homogeneous dsDNA chains with different nucleotide sequences and suggested the scales in which the temperature dependencies of D (0) are close for chains with different sequences. Analyzing the spectra in this scale here, we see that differences in the frequency-dependent diffusion coefficient blend as the frequency is increasing, and virtually disappear at the frequency, where an absolute maximum of the spectrum is observed (see Figure 1, solid lines).
3
Figure 1. Temperature dependencies of the frequency-dependent coefficient; A, Poly A dsDNA; B, Poly C dsDNA; C, Poly G dsDNA; D, Poly T dsDNA. Dashed lines for the dependencies of 2D (0) on d; solid lines for the dependencies of (max 2D (^))d on 9
As the temperature changes, the form of autocorrelation functions and spectra change too. Their analysis enables us to identify four different regimes. Since we have shown that any differences between the chains with different nucleotide structures are insignificant in our scale, and their character is determined, we will deal only with Poly A dsDNA chains in what follows. In the limit of infinitely small temperatures, the autocorrelation function tends to attenuate exponent:
iP(T) = ip(0) e-T/T°.
Its generalized frequency spectrum has an analytical form similar to that of the real part of the frequency-dependent conductivity, according to Drude formula:
D( <f>) =
i>(0) Tp 1 + 4 n 2(f)2 Tp
In accordance with fluctuation-dissipative theorem, the mobility fi (oj) is related to D(uS) as [17]:
Ii(2tt(f)) = ^D (2tt<f>).
The dependence of fi (0) on temperature is given in [18].
Notice also that the maximum of the frequency-dependent diffusion coefficient and, accordingly, the mobility at the relevant frequency is greater than the static diffusion coefficient (see Figure 1). It suggests that DNA can serve as a conductor of alternating current.
In the range of extra-low temperatures for near-zero 0, the autocorrelation functions and the spectra keep monotonously attenuating (see Figure 2) until
the dimensionless parameter 0 reaches the value 0 = 0.14 in Poly T dsDNA or 6 = 0.3 in Poly A dsDNA. This boundary is slightly different for chains with different nucleotide sequences (see Figure 1). For the first regime, scattering at long distances exceeding one site on the average is typical.
(a) 4 3.5 3 2.5 2 1.5 1
0.5 0 0.5
(b)
10
,2
Q
10
i
10
,0
=
?
Figure 2. (a) Velocity autocorrelation functions and (b) generalized frequency spectrum of the charge propagation along Poly A dsDNA in the range of low and extra-low temperatures; A, d = 0.1; B, d = 0.2; C, 9 = 0.4; D, d = 0.8
10° 101 102 103
10-4 10-3 10-2 10-1 100
In the range of low temperatures, until the dimensionless parameter becomes 9=1, reflection at long distances arises. In this regime, the velocity autocorrelation function (see Figure 2(a)) is not monotone and has a single minimum in the negative range. The absolute value of the minimum increases with growing 0, while its coordinate decreases. Now the spectra demonstrate a low-frequency maximum whose frequency increases with growing d, while the value decreases (see Figure 2(b)).
Then, up to the value of d = 3.5, a range of moderate temperatures is observed. In this temperature range a high-frequency maximum caused by oscillations of the autocorrelation function ^ (r) occurs after the first minimum (see Figure 3(a)). At the same time, the absolute value of the first minimum of the autocorrelation function keeps growing. High-frequency oscillations arise while the function attenuates after the first minimum. Thus, in the range of moderate temperatures, reflection at long distances, and reflection from neighboring sites coexist.
Looking at the spectrum 2D(u), we see that the low-frequency maximum keeps decreasing, while its frequency decreases and approaches the frequency of the second maximum (Figure 3(b)). The position of the high-frequency peak on the spectrum is independent of temperature. The first and the second maxima of the spectrum are well seen and have nearly the same value for d = 1.9 (see the inset in Figure 3(b)).
(a) 4 3.5 3 2.5 2
^ 1.5
1
0.5 0
-0.5 -1
(b)
10-
10°
10'
Q
CM
10
10
10
10
102
10
10
10-
10u
=
?
Figure 3. (a) Velocity autocorrelation functions and (b) generalized frequency spectrum of the charge propagation along Poly A dsDNA in the range moderate temperatures; A, 9 = 1.2; B, 9 = 1.6; C, 9 = 1.9; D, 9 = 2.4; E, 9 = 3.5
Notice that the value of the parameter 0 = 1.9 for which the high-frequency peak is approximately equal in value to the low-frequency one is the same for different types of chains, which is confirmed by the curves of the maxima shown in Figure 1. The value of the parameter 6 for which the low-frequency maximum of the spectrum disappears while the high-frequency one emerges can reasonably be called a critical value. Dimensional critical temperatures corresponding to the critical parameter are listed in Table 3 for different types of chains.
Table 3
Critical temperature for 9 = 1.9 and position of high-frequency peak for different nucleotide
sequences
units poly A poly C poly G poly T
■fp THz 14.5 19.8 40.6 76.4
Tc ° C -134.76 -14.68 811.80 3565.37
Regularity and homogeneity of the chain cause the high-frequency peak. It occurs at a frequency of (f) = 1/tx. The dimension frequency depends only on matrix element:
The values of this peak frequency for the different nucleotide sequences are presented in Table 3. Obviously, in the case of inhomogeneous regular chains,
these frequencies will be more than one, and the spectra will contain several relevant peaks.
In the range of high temperatures, motion mainly fulfills as hopping between neighboring sites. As distinct from the previous regimes, the absolute value of the first minimum of the autocorrelation function decreases here with growing d (see Figure 4(a)). A decrease of its coordinate is considerably retarded and tends to 1. At the same time, the initial attenuation occurs more sharply. Initially, the function sharply decreases to zero and then slowly goes to the first minimum. Oscillations of the autocorrelation function in this regime occur about the abscissa axis rather than about the low-frequency envelope.
(a) 4 3.5 3 2.5 2 1.5 1
0.5 0 0.5
(b)
Q
CM
10"3 10";
10" =
100 10
10
10
10
10
10
10
10
101
102
?
Figure 4. (a) Velocity autocorrelation functions and (b) generalized frequency spectrum of the charge propagation along Poly A dsDNA in the range of high and ultrahigh temperatures; A, d = 4.4; B, d= 10; C, d = 40; D, 9 = 160; E, d = 640
The spectrum demonstrates only a high-frequency peak, which decreases with growing 0 (see Figure 4(b)). The spectrum smoothly widens.
Obviously, in the limit of high temperatures, the spectrum widens and becomes monotone. The charge propagation ceases.
As is seen from Figures 2-4, for rather high u, at any temperatures, except for extra-low and superhigh ones, D (u) increases with growing temperature. This effect is known for static disorder [17].
4. Conclusion
We considered the frequency spectra of the diffusion coefficient (mobility) of an excess charge in a Holstein molecular chain for different temperatures.
We revealed the character of motions contributing to the charge diffusion in the Holstein model. The presence of temperature causes a dynamical disorder even in an ideal chain. Let us call the main differences between the frequency dependencies of chains with static disorder and those with
the dynamic disorder caused by temperature fluctuations. In the static case, the diffusion coefficient increases with growing frequency [17], while in the dynamic case, more complicated behavior is observed (see Figures 2-4). We demonstrated the existence of asymptotics at low and high temperatures. For finite temperatures, we identified four temperature ranges in which autocorrelation functions and their spectra have similar peculiarities/features. Accordingly, the charge diffusion inside these ranges demonstrates the same character. The character is changing in passing to a different temperature range. Low temperatures cause weak scattering without reflections, as a result of which the ballistic motion changes for the diffusion one. As the temperature increases, at first, a reflection at long distances emerges, then it is added with reflections from neighboring sites. Both these motions slow down the charge diffusion. In the limit of high temperatures, the motion degenerates into hopping between neighboring sites, and the charge diffusion ceases.
We revealed the character of motions contributing to the charge diffusion in the Holstein model.
All the results are presented in dimensionless values, which enabled us to identify the main parameters and recognize the different regimes.
The autocorrelation functions in themselves, as well as their spectra, are also a significant result. Their calculation was made possible only with the use of special unique methods of numerical integration of differential equations.
Presently, DNA is considered to be a promising material for constructing molecular wires. However, even homogeneous synthetic nucleotide chains at room temperatures demonstrate deficient mobility in the case of direct current [19], [20]. A possible solution to this problem, according to the results obtained, is the use of not direct but alternating current of rather high frequency.
Acknowledgments
This article is published on the basis of preprint [21].
References
[1] P. Maniadis, G. Kalosakas, K. 0. Rasmussen, and A. R. Bishop, "AC conductivity in a DNA charge transport model," Physical Review E, vol. 72, p. 021912, 2 Aug. 2005. DOI: 10.1103/PhysRevE.72.021912.
[2] G. L. Goodvin, A. S. Mishchenko, and M. Berciu, "Optical conductivity of the Holstein polaron," Physical Review Letters, vol. 107, p. 076 403, 7 Aug. 2011. DOI: 10.1103/PhysRevLett.107.076403.
[3] L. D. Siebbeles and Y. A. Berlin, "Quantum motion of particles along one-dimensional pathways with static and dynamic energy disorder," Chemical Physics, vol. 238, no. 1, pp. 97-107, 1998. DOI: 10.1016/S0301-0104(98)00311-5.
[4] P. Prins, F. C. Grozema, J. M. Schins, and L. D. A. Siebbeles, "Frequency dependent mobility of charge carriers along polymer chains with finite length," Physica Status Solidi B, vol. 243, no. 2, pp. 382-386, 2006. DOI: 10.1002/pssb.200562719.
[5] C. J. Murphy, M. R. Arkin, Y. Jenkins, N. D. Ghatlia, S. H. Bossmann, N. J. Turro, and J. K. Barton, "Long-range photoinduced electron transfer through a DNA helix," Science, vol. 262, no. 5136, pp. 10251029, 1993. DOI: 10.1126/science.7802858.
[6] P. O'Neill, A. W. Parker, M. A. Plumb, and L. D. A. Siebbeles, "Guanine modifications following ionization of DNA occurs predominantly via intra- and not interstrand charge migration: an experimental and theoretical study," Journal of Physical Chemistry B, vol. 105, no. 22, pp. 5283-5290, 2001. DOI: 10.1021/jp003514t.
[7] G. I. Livshits et al., "Long-range charge transport in single G-quadruplex DNA molecules," Nature Nanotechnology, vol. 9, no. 12, pp. 1040-1046, 2014. DOI: 10.1038/nnano.2014.246.
[8] V. D. Lakhno, "DNA nanobioelectronics," International Journal of Quantum Chemistry, vol. 108, no. 11, pp. 1970-1981, 2008. DOI: 10.1002/ qua.21717.
[9] T. Chakraborty, Charge migration in DNA: perspectives from Physics, Chemistry, and Biology, ser. NanoScience and Technology. Springer Berlin Heidelberg, 2007.
[10] A. Offenhausser and R. Rinaldi, Eds., Nanobioelectronics - for Electronics, Biology, and Medicine. Springer, New York, 2009, 337 pp.
[11] T. Holstein, "Studies of polaron motion," Annals of Physics, vol. 8, no. 3, pp. 325-342, 1959. DOI: 10.1016/0003-4916(59)90002-8.
[12] N. S. Fialko and V. D. Lakhno, "Nonlinear dynamics of excitations in DNA," Physics Letters A, vol. 278, pp. 108-111, 1-2 2000. DOI: 10.1016/S0375-9601(00)00755-6.
[13] D. A. Tikhonov, N. S. Fialko, E. V. Sobolev, and V. D. Lakhno, "Scaling of temperature dependence of charge mobility in molecular Holstein chains," Physical Review E, vol. 89, p. 032 124, 3 Mar. 2014. DOI: 10. 1103/PhysRevE.89.032124.
[14] E. V. Sobolev, D. Tikhonov, and N. S. Fialko, "About Numerical solution of the Holstein's discrete model [O chislennom reshenii uravneniy diskretnoy modeli Kholsteyna]," in Proceedings of the XIX All-Russian Conference "Theoretical bases and generation of numerical algorithms of solving mathematical physics problems", devoted to K. I. Babenko, Durso, Russia, 2012, in Russian, Moscow: Keldysh Institute of Applied Mathematics, 2012, p. 90.
[15] E. V. Sobolev, D. Tikhonov, and N. S. Fialko, "Numerical solution of the Holstein's discrete model in the problem of charge transfer in DNA [Chislennoye resheniye uravneniy diskretnoy modeli Kholsteyna v zadache o modelirovanii perenosa zaryada v DNK]," in Proceedings of the 4th International Conference on Mathematical Biology and Bioin-formatics, Pushchino, Russia, 2012, V. D. Lakhno, Ed., in Russian, Moscow: MAKS Press, 2012, p. 18.
[16] D. A. Tikhonov, E. V. Sobolev, V. D. Lakhno, and N. S. Fialko, "Adia-batic approximation for the calculation of the charge mobility in the DNA Holstein model [Adiabaticheskoye priblizheniye pri raschetakh podvizhnosti zaryada v kholsteynovskoy modeli DNK]," Matematich-eskaya biologiya i bioinformatika, vol. 6, no. 2, pp. 264-272, 2011, in Russian. DOI: 10.17537/2011.6.264.
[17] J. C. Dyre and T. B. Schr0der, "Universality of AC conduction in disordered solids," Reviews of Modern Physics, vol. 72, pp. 873-892, 3 Jul. 2000. DOI: 10.1103/RevModPhys.72.873.
[18] V. D. Lakhno and N. S. Fialko, "Bloch oscillations in a homogeneous nucleotide chain," English, JETP Letters, vol. 79, no. 10, pp. 464-467, 2004. DOI: 10.1134/1.1780553.
[19] D. M. Basko and E. M. Conwell, "Effect of solvation on hole motion in DNA," Physical Review Letters, vol. 88, p. 098 102, 9 Feb. 2002. DOI: 10.1103/PhysRevLett.88.098102.
[20] V. D. Lakhno and N. S. Fialko, "Solvation effects on hole mobility in the poly G/Poly C duplex," Russian Journal of Physical Chemistry A, vol. 86, no. 5, pp. 832-836, 2012. DOI: 10.1134/S0036024412050196.
[21] D. A. Tikhonov, E. V. Sobolev, and V. D. Lakhno, "Charge diffusion in homogeneous molecular chains based on the analysis of generalized frequency spectra in the framework of the Holstein model," Tech. Rep. 70-e, 2018, pp. 1-16. DOI: 10.20948/prepr-2018-70-e.
For citation:
D. A. Tikhonov, E. V. Sobolev, V. D. Lakhno, Charge diffusion in homogeneous molecular chains based on the analysis of generalized frequency spectra in the framework of the Holstein model, Discrete and Continuous Models and Applied Computational Science 27 (3) (2019) 217-230. DOI: 10.22363/26584670-2019-27-3-217-230.
Information about the authors:
Dmitry A. Tikhonov (Russian Federation) — Candidate of Physical and Mathematical Sciences, Senior researcher, Institute of Mathematical Problems of Biology Branch of Keldysh Institute of Applied Mathematics Russian Academy of Sciences (e-mail: dmitry.tikhonov@gmail.com, phone: +7(915)4013703, ORCID: https://orcid.org/0000-0002-1779-464X, ResearcherID: O-2214-2013, Scopus Author ID: 7006659212)
Egor V. Sobolev (Germany) — Candidate of Physical and Mathematical Sciences, Postdoctoral fellow, European Molecular Biology Laboratory, Hamburg Unit (e-mail: egor@embl-hamburg.de, phone: +49(4089)902184, ORCID: https://orcid.org/0000-0003-2478-5685, ResearcherID: O-2216-2013, Scopus Author ID: 55971569900)
Victor D. Lakhno (Russian Federation) — Doctor of Physical and Mathematical Sciences, Scientific Director, Institute of Mathematical Problems of Biology Branch of Keldysh Institute of Applied Mathematics Russian Academy of Sciences (e-mail: lak@impb.ru, phone: +7(4967)318504, ORCID: https://orcid.org/0000-0001-9224-769X, ResearcherID: N-7346-2013, Scopus Author ID: 7003392437)
УДК 519.6
DOI: 10.22363/2658-4670-2019-27-3-217-230
О диффузии заряда в однородных молекулярных цепочках на основе анализа обобщенных частотных спектров в рамках модели Холстейна
Д. А. Тихонов12, Е. В. Соболев13, В. Д. Лахно1
1 Институт математических проблем биологии (ИМПБ РАН) ул. проф. Виткевича, д. 1, г.Пущино, Московская область, 142290, Россия
2 Институт теоретической и экспериментальной биофизики РАН ул. Институтская, д. 3, г. Пущино, Московская область, 142290, Россия 3 Европейская лаборатория молекулярной биологии, отделение в Гамбурге с/о DESY, д. 25А, Ноткештрассе 85, 22607 Гамбург, Германия
В статье проведён анализ автокорреляционных функций скорости и обобщённых частотных спектров распространения заряда в однородных последовательностях ДНК при конечной температуре. Функции рассчитаны численно в рамках квазиклассической модели Холстейна. Показано, что в системе только один параметр главным образом определяет кинетику заряда для всех последовательностей. Анализ позволил определить характер отдельных движений, вносящих вклад в подвижность заряда, и выделить различные режимы распространения заряда в зависимости от температуры.
Ключевые слова: перенос заряда, автокорреляционная функция скорости, обобщённый частотный спектр, ДНК, модель Холстейна