CC BY
70
UDC 537.21, 539.217, 620.3+004.942 Doi: 10.31772/2587-6066-2018-19-4-677-682
For citation: Zalizniak V. E., Zolotov O. A., Zolotova O. P. [Mathematical model of conducting nanopore for molecular dynamics simulations]. Siberian Journal of Science and Technology. 2018, Vol. 19, No. 4, P. 677-682. Doi: 10.31772/2587-6066-2018-19-4-677-682
Для цитирования: Зализняк В. Е., Золотов О. А., Золотова О. П. Математическая модель проводящей нано-поры для расчетов методом молекулярной динамики // Сибирский журнал науки и технологий. 2018. Т. 19, № 4. С. 677-682. Doi: 10.31772/2587-6066-2018-19-4-677-682
MATHEMATICAL MODEL OF CONDUCTING NANOPORE FOR MOLECULAR
DYNAMICS SIMULATIONS
V. E. Zalizniak1, O. A. Zolotov1, O. P. Zolotova2*
1 Siberian Federal University 79, Svobodny Av., Krasnoyarsk, 660041, Russian Federation 2 Reshetnev Siberian State University of Science and Technology 31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation
* E-mail: zolotova@sibsau.ru
An electrostatic model of conducting nanopore is presented in the paper. The model does not require solution of the Poisson equation for the potential. This model is intended for use in simulation of transport phenomena of charged particles in conducting nanopores by the method of molecular dynamics. This method is based on Newton's equations of motion and it allows one to determine the variation of position, velocity and acceleration of particles with time. The electric field from the charge distributed over the nanopore surface is approximated by the field from fictitious point charges on the same surface. To verify the proposed model of fictitious charges system capacitance is calculated. The obtained values of capacitance are compared with classical results for conducting tubule and with the results obtained by the other similar method. The comparison shows that relative discrepancy between results is less than 10 %. There is a need to further develop the proposed model both in case of a large number of fictitious charges and in case when charged particles are in close proximity to the nanopore surface. The proposed method can be easily applied to an arbitrary shape nanopore. The model can be used in the development of various nanodevices, among them the devices used in life support systems of manned space vehicles.
Keywords: conducting nanopore, ion transport, molecular dynamics.
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОВОДЯЩЕЙ НАНОПОРЫ ДЛЯ РАСЧЕТОВ МЕТОДОМ МОЛЕКУЛЯРНОЙ ДИНАМИКИ
В. Е. Зализняк1, О. А. Золотов1, О. П. Золотова2*
1 Сибирский федеральный университет Российская Федерация, 660041, г. Красноярск, просп. Свободный, 79
2Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева Российская Федерация, 660037, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31
* E-mail: zolotova@sibsau.ru
Представлена электростатическая математическая модель проводящей нанопоры. Расчет потенциала проводящей нанопоры не требует решения уравнения Пуассона. Модель предназначена для компьютерного моделирования транспорта заряженных частиц в проводящих нанопорах методом молекулярной динамики. Этот вычислительный метод основан на втором законе Ньютона и позволяет получить траекторию, которая описывает положения, скорости и ускорения частиц со временем. Электрическое поле распределенного по поверхности нанопоры заряда аппроксимировано полем фиктивных точечных зарядов, расположенных на той же поверхности. Выполнено сравнение электроемкости системы фиктивных зарядов, рассчитанной предложенным методом, с классическими решениями для проводящего тубулена и с решениями, полученными другим подобным методом. Предложенный метод показал достаточную точность (около 10 %) и предсказуемость поведения. Необходимо его дальнейшее исследование и развитие как для случая достаточно большого числа фиктивных зарядов, так и для случая достаточно сильного приближения ионов к поверхности нанопоры. Рассмотренный метод введения фиктивных зарядов может быть легко применен к нанопорам произвольной
формы. Результаты выполненного моделирования могут быть использованы при разработке различных нано-технологических устройств, в том числе используемых в качестве компонентов систем жизнеобеспечения пилотируемых космических аппаратов.
Ключевые слова: проводящая нанопора, транспорт ионов, молекулярная динамика.
Introduction. Nanopores represent a type o pores with a diameter of about several nanometers. Channels of such size can be both natural and artificially made objects. Currently researchers take great interest in the issues of studying and modeling of nanopore properties due to their broad application in science and technology [1; 2]. Artificial nanopores can be used, for example, for imitation of biological nanopore functions. They also have attractive potential applications in nanofluid fields [3]. Among these applications are desalination, ionic sieves [4], sensors for biological agents and sequencing of DNA or RNA [5].
Various forms of artificial nanopores using different production technologies and materials, including polymers, inorganic substances, biotic and composite materials have been developed [6; 7].
Conducting nanopores attract considerable attention due to their important potential applications. In a number of works production methods of similar structures are given. For example, such nanopores can be made via integration of carbon structures into the corresponding nanoporous material. A multilayered carbon nanotube was made inside a pore and on the surface of an anode aluminum oxidic membrane, applying the method of chemical sedimentation from a steam phase [8]. Strengthened single-layer metal carbon nanotubes were grown on a silicon plate via chemical sedimentation from a steam phase [9]. Synthesis of conducting carbon tubules inside the pores of the anode aluminum and oxidic membrane was investigated [10]. Similar structures can be of interest in terms of production of selective membranes switched by electric field.
In recent years mathematical and electrostatic models of ion transport in various types of nanopores and nano-channels have been vigorously developed, computer modeling of processes have been carried out as well [11-13].
At computer modeling of ion transport, applying methods of molecular dynamics, the problem of electric field defining, created in the system, has to be solved at each step in time. In presence of electro-conductive elements in the system, solving such a task with classical methods becomes rather resource-intensive. Therefore, problem solving methods of such electrostatic tasks should be fairly simple, although they are not as accurate as classical ways. In this regard there is a problem of such methods accuracy assessment in comparison with classical ones. A rather typical assessment might be a comparison of integrated electric parameters of the system, obtained by the offered methods, with classical decisions. In this work electric capacitance parameter will be compared with the data from the reference book [14].
The purpose of this work is to simulate the conducting nanopore electric field in presence of external charges (ions). This model is intended for use in fundamental research of charged particles transfer phenomenon (ions and charged proteinaceous residue) in conducting nanopores.
In the work [15] a method of the conducting nanopore electric field in the presence of external charges (ions) modeling in relation to problems of molecular dynamics was offered. In the present work an alternative for further development of this method is being investigated.
The results of simulation can be used when developing various devices, including, for instance, a component of life support systems in piloted spacecrafts.
Problem definition. The membrane containing a nanopore is simulated in the form of L wide uniform dielectric plate. The nanopore is a cylindrical hole with radius a and length L (fig. 1).
Fig. 1. Model of nanopore. Dots denote fictitious electric charges
Рис. 1. Модель нанопоры. Точками обозначены фиктивные электрические заряды
We shall consider the problem of ion transport in the context of continual electrostatic model. We shall assume that the ion is located in random position r* . Electrostatic potential of ф(г) in and outside the nanopore satisfies the Poisson's equation
1 N
Аф = - 4ЛТ £ 9"5(r- r*)' (1)
^'^О n=1
where А - Laplasian operator; qn - ion charge; N - ion number; 5(r - rn*) - delta-Dirac function. The equation
(1) has to be solved across all the space, except the area of membrane material with the potential set of boundary conditions on the tubulen surface. It is also necessary to demand the continuity of the membrane, normal to the surface, components of electric induction vector and, whenever possible, decrease of potential while being removed from the system.
Under consideration of charged particles transport phenomenon, we are interested in the force affecting the ion from the conducting nanopore which is expressed
as f = -q -grad(ф (r*)).
The solution of the equation (1) with the set boundary conditions is a difficult task. The problem becomes even
more complicated if to consider not only cylindrical, but also other forms of nanopores. Further we will try to present simple and rather universal model of conducting nanopores which does not require solving the Poisson equation for the potential.
Electrostatic model of conducting nanopore. We shall approximate electric field distributed on the body (tubulen) of the charge by the field of the dot charges located on the same surface. We shall place these fictitious charges in the chosen in advance surface dots. We will determine their size from the condition of maximum approach to mechanical balance. Assuming that fictitious charges could move on a tubulen surface, then equality to zero projection on the turbulen surface of total force affecting each charge from other charges would be condition of balance. This case is equivalent to the requirement of equality to zero projections on the tubulen surface of total electric field strengths created by all other charges in the point of each fictitious charge. In our case as positions of charges are fixed, the above-stated requirement is generally impracticable. Therefore, we shall demand the minimal sum of these projections squares.
Thus, we shall place the fictitious charges Qk, k = 1, ..., K, at the point with coordinates rk. Then we obtain electric-field vector Projection EIX at point r¿:
E (r ) _ 1 V Q (ri - rk) - ni (ni * (ri - rk)) Eft (ri) _ V Qk-^-
G k*i
r - r,.
(2)
where n,- - normal vector to the surface at the point r,-; * is dot product of vectors.
By summing squares (2) where all i = 1, ..., K, we obtain
VE2(ri)
: VVVQkQl
16n2- G
(ri - rk) * (ri - rl) -- (ni * (ri - rk ))(ni * (ri - r¡))
I ! - г, I31 Г, - Г, I3 .
(3)
q _V Qi.
(4)
V E,- (ri ) _ -A-2 Q'X A X Q,
i_1 16n -„
(5)
A * _V
(ri- rk) * (ri- r,) - (ni * (ri- rk))(ni * (ri- r,))
I r - r, I I r - r,
i*k,l I li lk
Under condition (4):
q = 1'x Q.
(6)
notes transposing operation of a vector column in a vector line; 1 = (1,1,1,...1) - made of units dimension vector K.
Solution to variational problem (5), (6), on condition that nondegenerate matrix A, will be vector
Q_q
A-1 x l
l'x A-1 x l.
(7)
In case of nondegenerate matrics A, under condition that a0 - evector A, corresponding to zero evalue, nonorthogonal 1, we obtain
Q_q
l x aG
(8)
We shall look for the minimum (3) on condition of specified total charge q of conducting tubulen:
Considering that a set of unknown charges is vector Q = (Qk), than we have a standard problem of square form minimization of (3):
Here x denotes dot product of vectors of K dimension and vectors and dimensional matrix K x K; stroke de-
The result (8) can be easily summarized and in case of greater 1 dimension of the evector subspace, corresponding to zero evalue A: in this case a0 - any non orthogonal to 1 vector from that subspace. However, in case a0 is orthogonal to 1, the task should be reduced, while decreasing A size and leaving behind only zero evalues. It should be noted that from the method A it is clear that A is symmetrical non-negative matrix.
We shall consider the alternative where A is a degenerate matrix. Generally this is true: we shall place nonzero dot charges of one sign on the tubulen surface and allow them to move on a tubulen surface. It is obvious that such system has a stationary position of charges. If to arrange fictitious charges similarly, then the matrix A will degenerate as there are charge values at which vector projection of electric field strength affecting each charge is equal to zero. Moreover, another stationary position, mostly continuously depending on charges, corresponds to a different set of value sizes which are freely moving on the tubulen surface. Therefore, we have a rather big area of coordinates of fictitious charges at which the matrix degenerate. At numerical calculations it will mean that if we chose the "wrong" positions of fictitious charges, and the matrix is not degenerated, then, while increasing in number of fictitious charges, generally, con-ditionality ratio of this matrix will grow.
The present work is devoted to preliminary assessment of the offered method practical applicability under molecular and dynamic calculations, therefore further cases of small number of fictitious charges at which in the reviewed examples matrix A is not degenerated are investigated.
The offered method is a further development of [15]. In the work [15] comparative analysis of the presented method of electrostatic task solution with other methods of electrostatics is carried out. As a comparative criterion, the results of electric capacitance calculation given in the reference book [14] were used. In this work for the above method applicability assessment we will use the results of electric capacitance of the conducting tubulen calculations and compare them with the results [14] and [15].
We shall continue the above suggested calculations in order to obtain estimates of conducting tubulen electric capacitance. We shall receive these estimates, comparing
a
0
i_1
1=1
the energy of the system considered with the expression for the energy of the conducting body:
w =—у
4пе, ^
QQj
0 i* j I ri rj
. q
2C,
(9)
where C - conducting tubulen electric capacitance. Adding (7) to (9), finally obtain the expression for electric capacitance:
C = 2ne,
IB „
(1'x A-1 X 1)2
0 1'x A"1 X B X A"1 X 1' 1
(10)
ri" rj 0
' * j; i = /
It shall be noted that in [15] nondegenerate matrics B was proved
Model validation. In order to validate the accuracy of the model suggested, we shall calculate the self-capacitance of the cylindrical hole of radius a and length L and compare it with the data given in the table [14].
We shall place the fictitious charges in points (zn, a, 0m) of the cylindrical system of coordinates, where
^n = (n " 1)
em = (m - 1)
L
N - Г 2n
N -1'
, n = 1,...,N,
(11)
m = 1,..., N.
Then K = N. Calculations are given in table for 5 < LI a < 20 h N = 10, 20, 30. Electrical capacitance value form the reference book [14] are marked C0, discrepancy between calculated and table-based [14] capacitance value d = (C - C0)IC0.
In fig. 2 the obtained results are presented graphically, the best results from the work [15] obtained under the
same algorithm of fictitious charges arrangement and N = 80 are given.
Results presented in table and fig. 2 show that when Li a > 8, discrepancy amounts to less than 15 % already for N = 20, and smoothly decreases under increase both LIa , and N, which is also an advantage of the suggested method. For comparison we must note that given in fig. 2 works [15] do not possess this property. From the practical point of view it might be expected that at Li a > 10, N = 30 parameters of charged particles interaction with the conducting nanopore are defined to a precision of less than 10 %.
Conclusion. One electrostatic model of conducting nanopore is presented in the work. The model is intended for use when modeling a nanopore, applying the method of molecular dynamics. This computing method is based on the second law of Newton. In case force intensity operating on each of atoms is known, integration of the movement equations allows to obtain a trajectory which describes positions, velocities and accelerations of particles over the time.
If fictitious charges are not entered, then to find the force operating on charged particle from the conducting nanopore, at first it is necessary to solve the equation (1) in each timepoint. In [15] the method allowing a single invert of the matrix of equation system by fictitious charges introduction is offered. The method presented in this work is similar, but demonstrates a better accuracy as it is shown in table and fig. 2. It is obvious that the offered model can be easily applied to nanopores of any form.
Thus, the offered method showed the sufficient accuracy and predictability. Therefore, its further research and development is necessary for both large number of fictitious charges leading to significant growth of matrix con-ditionality and for the case of rather close approximation of ions to a nanopore surface.
The results of calculations of the electrical capacitance of a conducting nanopore
L/a, relative units. N = 10 N = 20 N = 30 C0/2^80a, [14], relative units.
C/2^e0a, relative units. d, % C/2^e0a, relative units. d, % C/2^e0a, relative units. d, %
5 2.294 35.58 2.129 25.83 2.073 22.52 1.692
6 2.457 30.90 2.283 21.63 2.223 18.43 1.877
7 2.604 28.34 2.418 19.17 2.354 16.02 2.029
8 2.740 25.75 2.543 16.70 2.475 13.58 2.179
9 2.869 23.45 2.660 14.46 2.588 11.36 2.324
10 2.993 21.47 2.772 12.50 2.696 9.42 2.464
11 3.113 19.68 2.879 10.69 2.799 7.61 2.601
12 3.229 18.06 2.983 9.07 2.899 6.00 2.735
13 3.342 16.65 3.083 7.61 2.995 4.54 2.865
14 3.451 15.30 3.180 6.25 3.088 3.17 2.993
15 3.557 14.08 3.276 5.07 3.179 1.96 3.118
16 3.659 12.90 3.369 3.95 3.268 0.83 3.241
17 3.759 11.84 3.461 2.98 3.355 -0.18 3.361
18 3.855 10.78 3.551 2.04 3.441 -1.12 3.480
19 3.948 9.73 3.639 1.14 3.524 -2.06 3.598
20 4.038 8.75 3.726 0.35 3.607 -2.85 3.713
С, relative units
[15]
5 10 15 20
L/a, relative units
Fig. 2. Relationship between capacitance of a nanopore and the value of Lia for various numbers of fictitious charges. Solid line - tabulated data on capacitance [14]
Рис. 2. Зависимость емкости нанопоры от величины Lia при различном количестве фиктивных зарядов. Сплошная жирная линия - табличные значения емкости [14]
References
1. Dekker C. Solid-state nanopores. Nature Nanotech-nology. 2007, Vol. 2, No. 4, P. 209-215. DOI: 10.1038/nnano.2007.27.
2. Thomas M., Corry B., Hilder T. A. What have we learnt about the mechanisms of rapid water transport, ion rejection and selectivity in nanopores from molecular simulation? Small. 2014, Vol. 10, No. 8, P. 1453-1465. DOI: 10.1002/smll.201302968.
3. Daiguji H. Ion transport in nanofluidic channels. Chemical Society Reviews. 2010, Vol. 39, No. 3, P. 901911. DOI: 10.1039/B820556F.
4. Majd S., Yusko E. C., Billeh Y. N. et. al. Applications of biological pores in nanomedicine, sensing, and nanoelectronics. Current Opinion in Biotechnology. 2010, Vol. 21, No. 4, P. 439-476. DOI: 10.1016/j.copbio.2010.05.002.
5. Singer A., Rapireddy S., Ly D. H., Meller A. Electronic barcoding of a viral gene at the single-molecule level. Nano Letters. 2012, Vol. 12, No. 3, P. 1722-1728. DOI: 10.1021/nl300372a.
6. Hou X., Zhang H., Jiang L. Building bio-inspired artificial functional nanochannels: from symmetric to asymmetric modification. Angewandte Chemie International Edition. 2012, Vol. 51, No. 22, P. 5296-5307. DOI: 10.1002/anie.201104904.
7. Zhang Y., Kong X. Y., Gao L. et. al. Fabrication of nanochannels. Materials. 2015, Vol. 8, No. 9, P. 6277-6308. DOI: 10.3390/ma8095304.
8. Kasi J. K., Kasi A. K., Wongwiriyapan W. et. al. Synthesis of carbon nanotube and carbon nanofiber in nanopore of anodic aluminum oxide template by chemical vapor deposition at atmospheric pressure.
Advanced Materials Research - Trans Tech Publications. 2012, Vol. 557-559, P. 544-549. DOI: 10.4028/www.scientific.net/AMR.557-559.544.
9. Choi W., Ulissi Z. W., Shimizu S. F. et. al. Diameter-dependent ion transport through the interior of isolated single-walled carbon nanotubes. Nature Communications. 2013, Vol. 4, P. 2397. DOI: 10.1038/ncomms3397.
10. Simunin M. M., Khartov S. V., Shiverskii A. V. et. al. Structures based on graphitized nanotubulenes with a common electrode in a matrix of porous anodic alumina for the purpose of forming electrically switchable membranes. Technical Physics Letters. 2015, Vol. 41, No. 11, P. 1047-1050. DOI: 10.1134/S1063785015110103.
11. Levin Y. Electrostatics of ions inside the nanopores and trans-membrane channels. Europhysics Letters. 2006, Vol. 76, No. 1, P. 163-169. DOI: 10.1209/epl/i2006-10240-4.
12. Modi N., Winterhalter M., Kleinekathoefer U. Computational modeling of ion transport through nanopores. Nanoscale. 2012, Vol. 4, No. 20, P. 61666180. DOI: 10.1039/b000000x.
13. Hilder T. A., Gordon D., Chung S. H. Computational modeling of transport in synthetic nanotubes. Nanomedicine: Nanotechnology, Biology and Medicine. 2011, Vol. 7, No. 6, P. 702-709. DOI: 10.1016/j.nano.2011.02.011.
14. Iossel Yu. Ya., Kochanov E. S., Strunsky M. G. Raschet elektricheskoy emkosti [Calculation of Electric Power]. Leningrad, Energoizdat, 1981, 288 p. (In Russ.).
15. Zolotov O. A, Zalizniak V. E. Simple mathematical model of conducting nanopore. Zhurnal Sibirskogo federal'nogo universiteta. Seriya Matematika i fizika. 2018, Vol. 11, No. 4, P. 505-512. DOI: 10.17516/19971397-2018-11-4-505-512.
Библиографические ссылки
1. Dekker C. Solid-state nanopores // Nature Nanotechnology. 2007. Vol. 2. P. 209-215.
2. Thomas M., Corry B., Hilder T. A. What have we learnt about the mechanisms of rapid water transport, ion rejection and selectivity in nanopores from molecular simulation? // Small. 2014. Vol. 10. P. 1453-1465.
3. Daiguji H. Ion transport in nanofluidic channels // Chemical Society Reviews. 2010. Vol. 39. P. 901-911.
4. Applications of biological pores in nanomedicine, sensing, and nanoelectronics / S. Majd [et al.] // Current Opinion in Biotechnology. 2010. Vol. 21. P. 439-476.
5. Electronic barcoding of a viral gene at the single-molecule level / A. Singer [et al.] // Nano Letters. 2012. Vol. 12. P. 1722-1728.
6. Hou X., Zhang H., Jiang L. Building bio-inspired artificial functional nanochannels: from symmetric to asymmetric modification // Angewandte Chemie International Edition. 2012. Vol. 51. P. 5296-5307.
7. Fabrication of nanochannels / Y. Zhang [et al.] // Materials. 2015. Vol. 8. P. 6277-6308.
8. Synthesis of carbon nanotube and carbon nanofiber in nanopore of anodic aluminum oxide template by chemical vapor deposition at atmospheric pressure / J. K. Kasi [et al.] // Advanced Materials Research -Trans Tech Publications. 2012. Vol. 557-559. P. 544-549.
9. Diameter-dependent ion transport through the interior of isolated single-walled carbon nanotubes / W. Choi [et al.] // Nature Communications. 2013. Vol. 4. P. 2397.
10. Структуры на основе графитизированных на-нотубуленов с общим электродом в матрице пористого анодного оксида алюминия для задачи формирования переключаемых электрическим полем мембран / М. М. Симунин [и др.] // Письма в журнал технической физики. 2015. № 21 (41). С. 52-59.
11. Levin Y. Electrostatics of ions inside the nanopores and trans-membrane channels // Europhysics Letters. 2006. Vol. 76. P. 163-169.
12. Modi N., Winterhalter M., Kleinekathoefer U. Computational modeling of ion transport through nanopores // Nanoscale. 2012. Vol. 4. P. 6166-6180.
13. Hilder T. A., Gordon D., Chung S. H. Computational modeling of transport in synthetic nanotubes // Nanomedicine: Nanotechnology, Biology and Medicine. 2011. Vol. 7. P. 702-709.
14. Иоссель Ю. Я., Кочанов Э. С., Струнский М. Г. Расчет электрической емкости. Л. : Энергоиздат, 1981. 288 с.
15. Золотов О. А., Зализняк В. Е. Простая математическая модель проводящей нанопоры // Журнал Сибирского федерального университета. Серия «Математика и физика». 2018. № 4 (11). С. 505-512.
© Zalizniak В. Е., Zolotov О. А., Zolotova О. P., 2018
