О потенциале взаимодействия открытой углеродной нанотрубки с молекулами-компонентами природного газа
,
Национальный исследовательский Томский государственный университет, 634050, проспект Ленина 36, Томск, Российская Федерация,
E-mail: *****@***ru
Аннотация. В работе изучается взаимодействие молекул гелия (Не) и метана (СH4) с энергетическим барьером, создаваемым однослойной углеродной нанотрубкой. Получены поля потенциала взаимодействия трубки как нанообъекта с одиночными молекулами или атомами. Расчетами была установлена зависимость скорости молекулы на оси симметрии однослойной углеродной нанотрубки от осевой координаты. Исследовано влияние радиуса трубки на характер прохождения молекул через открытую трубку.
В работе [1] было показано, как от общих соотношений для потенциала Леннард-Джонса, которому авторы придали вид классического источника, получены выражения для плоской пластины, сферы и бесконечной тонкой нанонити. В данной статье рассмотрено потенциальное поле, создаваемое однослойной углеродной нанотрубкой. Построенное поле справедливо как для внутренней области трубки, так и для её внешней области.
Рис. 1. Нанотрубки различной хиральности
Углеродные нанотрубки представляют для задач фильтрации и адсорбции газов из смесей. С точки зрения построения модели, однослойная нанотрубка является объектом, один из линейных размеров которого много меньше двух других, и составляет значения порядка размера (диаметра) атома углерода. Это позволяет рассматривать их как двумерные объекты. В работе [1] представлена модификация потенциала Леннарда-Джонса, имеющая интегрируемую особенность в нуле:
(1)
Здесь 
– межатомное или межмолекулярное расстояние; σ, ε – параметры потенциала Леннарда-Джонса. Значения констант ε и σ, входящих в потенциал Леннарда-Джонса, для некоторых пар веществ приведены в следующей таблице.
Таблица 1 Константы потенциала Леннарда-Джонса
Пары взаимодействующих молекул | Относительная глубина потенциальной ямы | Радиус влияния потенциала взаимодействия |
C — C | ε/k = 51,2 [K] | σ = 0,3029 [нм] |
He — He | ε/k = 10,2 [K] | σ = 0,257 [нм] |
CH4 — CH4 | ε/k = 148 [K] | σ = 0,38 [нм] |
Если исследуемая система состоит из разнородных молекул (атомов), то для параметров ε и σ справедливы следующие правила усреднения Лорентца-Бертло (Lorentz-Berthelot mixing rule):
(2)
Мы попробуем найти интегральное силовое воздействие от открытой нанотрубки как отдельной структуры, являющейся кристаллической наночастицей.
В цилиндрических координатах 
, квадрат расстояния между двумя точками (M и M', где M' – точка на поверхности) есть:
(3)
Потенциал (1) будем интегрировать по поверхности открытой нанотрубки, т. е. по боковой поверхности цилиндра, тогда с учетом (3) можно записать:
. (4)
Если (4) умножить на плотность распределения источников силового взаимодействия (плотность распределения атомов углерода на поверхности трубки) 
(здесь S – общая площадь поверхности графеновой трубки, n – число атомов на этой поверхности, q ≈ 28 нм–2) и одновременно положить r' = a – радиус нанотрубки, то получим выражение для потенциала воздействия от всех атомов углерода нанотрубки:
. (5)
Здесь l – длина трубки. Причем, выражение (5) будет иметь размерность потенциала (энергии). Заметим, что если геометрические размеры нанотрубки l и a представлять в нанометрах, тогда с учетом размерности q, потенциал будет определяться размерностью одного параметра 
. Выражение (5) еще модно переписать следующим образом:
, (6)
где 
.
Из (6) видно, что потенциал имеет ту же размерность, что и ε. Выражая ε в градусах Кельвина, получаем значения потенциала в той же величине. Т. к. результат интегрирования для осесимметричного объекта не должен зависеть от угла φ, и этот угол, в сущности, является координатой точки наблюдения, то в (4) целесообразно положить φ=0. Что и сделано в (5) и (6).
Рис. 2 Поле потенциала, выраженного в К°, созданное однослойной углеродной нанотрубкой относительно молекулы He.
Рис. 3 Поле потенциала, выраженного в К°, созданное однослойной углеродной нанотрубкой относительно молекулы CH4.
Имея потенциал нанотрубки (5), с помощью методов классической молекулярной динамики можно изучать проницаемость отдельной нанотрубки, а также проницаемость материалов, полученных компактированием однородных или разнородных трубок. Кроме этого, важными являются задачи по изучению сорбционных свойств нанотрубок. Из рис. 2 и 3 видно, что в центральной части нанотрубки имеется потенциальная яма с отрицательными значениями энергии взаимодействия. Наличие этой ямы означает возможность накопления молекул или атомов какого-либо вещества. Таким образом, приосевую зону трубки можно назвать сорбционной зоной. Так же имеется кольцевая сорбционная зона около внешней поверхности трубки. Все эти зоны сформированы силами притяжения потенциала Леннарда-Джонса, то есть силами дальнего действия. Для описания движения молекулы вещества в области потенциала, созданного нанотрубкой, используем следующую форму закона сохранения.
(7)
На основе закона сохранения (7) построены зависимости скоростей молекулы на оси симметрии однослойной углеродной нанотрубки от осевой координаты при различных значениях радиуса трубки. Здесь индексом «0» отмечены значения величин на таком удалении от трубки, при котором 
с учетом этого из закона сохранения (7) получим:
(8)
Рис 4. – Зависимость скорости молекулы CH4 на оси симметрии однослойной углеродной нанотрубки от осевой координаты при различных значениях радиуса трубки.
Рис 5. – Зависимость скорости молекулы He на оси симметрии однослойной углеродной нанотрубки от осевой координаты при различных значениях радиуса трубки.
Попадая во внутреннюю область, молекулы продолжают свое движение, испытывая ускорение на участке до середины трубки и замедление после значения z=0 (см. рис. 4,5) и выходят из трубки, имея значения скорости равное исходной величине. Существуют режимы (отвечающие малым значениям радиуса трубки) при которых определяющей является часть потенциала Леннарда-Джонса, отвечающая за отталкивание молекул. При таких режимах молекула вообще не может войти в трубку, то есть она является не проходимой для данного вещества.
Работа выполнена при поддержке программы повышения конкурентоспособности Томского государственного университета.
Список литературы:
1. , , Потенциальное поле углеродных тел как основа сорбционных свойств барьерных газовых систем // Известия ВУЗов. Физика. – 2015. – Т. 58, № 4. – С.
2. А. I. Potekaev, А. M. Bubenchikov, M. A. Bubenchikov, New physical ideas and method of description and Calculation of resistance to motion of small Particles in a gaseous medium, Russian physics journal Vol. 55 (12) (2013) 1434-1443.
3. M. A. Bubenchikov, А. I. Potekaev, А. M. Bubenchikov, Three fundamental problems of molecular Statistics, Russian physics journal Vol. 56 (3) (2013) 341-348.
4. M. A. Bubenchikov, А. I. Potekaev, А. M. Bubenchikov, Thermophoresis of ultrafine and nanosized Particles, Russian physics journal Vol. 56 (7) (2013) 785-790.
5. M. A. Bubenchikov, А. I. Potekaev, А. M. Bubenchikov, A scheme for hydrogen evolution from natural gas, Applied Mechanics and Materials Vol. 698 (2015) 525-528
6. M. A. Bubenchikov, А. I. Potekaev, А. M. Bubenchikov, Migration of nanoparticles in operating media of heat exchange systems, Applied Mechanics and Materials Vol. 698 (2015) 627-630
