CC BY
64
УДК 536.75
DOI: 10.18698/1812-3368-2019-2-94-103
БРОУНОВСКОЕ ДВИЖЕНИЕ КАК НЕОБРАТИМЫЙ НЕМАРКОВСКИЙ ПРОЦЕСС
А.Н. Морозов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
amor59@mail.ru amor@bmstu.ru
Аннотация
Предложен метод описания броуновского движения в неравновесной среде при протекании необратимых процессов. Рассчитана спектральная плотность флуктуаций скорости броуновской частицы в неравновесной среде и установлено, что в низкочастотной части спектра она представляет собой фликкер-шум. Разработанный метод описания броуновского движения в неравновесной среде применен для расчета флуктуаций тока в малом объеме электролита. Получены оценки константы Хоуге и постоянной времени хаотизации ионов в электролите, совпадающие с экспериментально полученными данными
Ключевые слова
Немарковский процесс, броуновское движение, неравновесное состояние, производство энтропии, фликкер-шум
Поступила 26.02.2018 © Автор(ы), 2019
Экспериментальные исследования флуктуаций электрического тока в различных средах показывают наличие в их спектрах фликкер-шума [1, 2]. В частности, такой шум наблюдается при протекании электрического тока в малых объемах электролита [3-6]. Фликкер-шум может быть объяснен в рамках модели, предполагающей существование флуктуаций кинетических коэффициентов, в частности коэффициента вязкости при броуновском движении [7].
Математическое описание фликкер-шума можно провести с помощью теории немарковских процессов, описывающей случайные процессы с памятью [8]. Так, эта теория описывает броуновское движение [9-11], диффузию [12, 13], теплопроводность [14, 15], флуктуации молекулярного и фотонного газов [16] и другие процессы.
Проведем описание броуновского движения в неравновесной среде [17, 18]. Предположим, что уравнение для скорости движения V броуновской частицы в одномерном случае можно представить в виде модифицированного уравнения Ланжевена
^= р о )+$! V)+х (О, (1)
где у — коэффициент релаксации броуновской частицы; Б (í) — детерминированное внешнее воздействие; £1 ) — случайный процесс, действующий со стороны частиц среды на броуновскую частицу в равновесном случае; X ) — немарковский процесс, описывающий действие частиц среды с учетом ее неравновесного состояния. Уравнение (1) без учета последнего слагаемого в правой части представляет собой традиционное уравнение Ланжевена, которое описывает броуновское движение в равновесной среде.
В уравнении Ланжевена, как правило, случайный процесс ) описывается моделью 5-коррелированного гауссова процесса с корреляционной функцией
Н1 (ь)) = ^5(^2 -^). (2)
т
Здесь к — постоянная Больцмана; Т — температура среды; т — масса броуновской частицы.
Для описания воздействия неравновесной среды на броуновскую частицу в работе [17] предложено ввести в правую часть уравнения Ланжевена дополнительное слагаемое X ), которое может быть представлено в виде интегрального преобразования
X^)= / О(*,т)^2 СО-^^йт, (3)
где О (I, т) — ядро интегрального преобразования,
С {и т)= , / . ; (4)
-т)
5-коррелированный гауссов процесс Ъ,2 ^) описывается корреляционной функцией
<^2 (*2 (*1 )) = —8(*2 " Н ). (5)
т
Здесь ут = 1/ То — интенсивность флуктуаций; То — постоянная времени хаотизации броуновской частицы; Os — производство энтропии при детерминированном движении броуновской частицы,
тБ2 /¿ч
= ——. (6) УТ
Добавление в правую часть уравнения (1) слагаемого X (г) в виде (3) позволяет сохранить для этого уравнения требования флуктуационно-диссипационной теоремы [15]. С учетом формул (3) и (4) уравнение (1) приобретает вид интегродифференциального уравнения
йУ лт \ у йУ (т), -+ уУ + I , -—йт =
йг -X) ¿1
= б (г (г)+ } ] £2 (т)йт. (7)
-т)
Найдем спектральную плотность флуктуаций скорости У (г). Преобразование Лапласа позволяет записать уравнение (7) в изображениях
(5) + УУ (5К^уУ (5) = Б(5) + |1 (5к/"^ («), (8) V VI V
или
у (5 ( б (5 ) +11 (5 + ^\2 (5 ) (9)
Выражение (9) позволяет найти спектральную плотность Су (о) флуктуаций скорости У (г):
Gy(ш) = \V(/ш)|2 ^^- . (Ю)
(vT i ш2 + у2 j + V 2лутш (ш + у) + л:у2ш I
При получении выражения (10) предполагалось, что F (t) = F = const, и были учтены формулы (2) и (5).
Для низкочастотного случая при ш ^ у и ш ^ vT формула (10) принимает вид
( ч 2kT 2-kgsT
GV (ю) =-+-• (11)
my mvTw
Первое слагаемое в (11) описывает флуктуации скорости броуновской частицы в равновесной среде, а второе — в локально-неравновесной с учетом диссипации энергии при движении броуновской частицы.
Из формулы (11) следует, что при движении броуновской частицы в неравновесной среде в низкочастотной части спектра спектральная плотность флуктуаций ее скорости описывается фликкер-шумом.
Применим полученные выражения для описания токовых флуктуаций в экспериментах по измерению фликкер-шума в малых объемах
электролита. В этих экспериментах электролитическая ячейка представляла собой два сосуда с электролитом, в которых размещались электроды, разделенные тонкой пленкой с небольшими отверстиями [19]. Плотность тока в электролите можно вычислить по формуле
1 = ц+п+У ++ -, (12)
где , q~ — заряд положительных и отрицательных ионов; п+, п ~ —
концентрация положительных и отрицательных ионов; V+, V" — скорости движения положительных и отрицательных ионов.
Скорости движения ионов связаны с напряжением и, приложенным к электролитической ячейке, формулой
V+/- = |Д+/- и ■ (13)
Здесь — подвижности положительных и отрицательных ионов; I — толщина пленки. Тогда ток, проходящий через отверстие в пленке, можно представить в виде
I = ^+пц~п (14)
где 5 — площади сечения отверстия в пленке.
Определим относительные флуктуации тока, протекающего через электролитическую ячейку:
§1 ^+п+У + + q ~п~У ")/
Io (q+nV++ q n д )u"o
(15)
где ио — электрическое напряжение, приложенное к электролитической ячейке.
Второе слагаемое в формуле (11) позволяет без учета равновесных флуктуаций найти двустороннюю спектральную плотность флуктуаций тока в малом объеме электролита для низкочастотной части спектра:
~ / ч 2Kl2T
Gbi (®) = -
(q+n+^,+ + q n д ) U(
2
2
f [q+n+)2 üj ^(q n )2 as ^
m+v|N + m
v J
(16)
Ш
Здесь знаки «+» или «-» соответствуют величинам положительного и отрицательного ионов; N+, N" — число положительных и отрицательных ионов в отверстиях пленки,
N+/"= п.
При выводе формулы (16) учтено то, что флуктуации скорости ионов не коррелируют друг с другом.
Формула (6) применительно к рассматриваемому случаю движения ионов в электролите может быть представлена в виде
ll+/-q +/-ТТ 2
j+/_ _ Ц q U 0 S рт
= ^ 0, (17)
где учтено
а +/
у+/- = а
Выражение (16) с учетом формулы (17) можно записать в стандартной форме
(* )=# •
где / = ш/(2тс) — частота флуктуаций тока; N = я") К — число ионов в малом объеме электролита; а — константа Хоуге,
({ , \3 , , / \з Л
n + n
а = ■
(а+n+^++а n р )
(а+) n+^+ | (а ) n д
m v% m vT
(18)
Если в первом приближении для оценочных расчетов принять, что параметры положительных и отрицательных ионов одинаковы, то формула (18) приобретает простой вид
Ч У а = —í— = —.
т|дут ут
Верхняя оценка константы Хоуге прохождения ионов через молекулярные каналы приведена в работе [20]: а<10"5. Тогда величину ут можно оценить как ут > 1018 с-1.
Если осуществить полосовую фильтрацию флуктуаций тока на электролитической ячейке в диапазоне от Ш1 до Ш2, то можно найти дисперсии флуктуаций тока 11 (I) и 12 ), соответствующие равновесному и неравновесному случаям и описываемые первым и вторым слагаемым в формуле (11). Используя формулу для скорости броуновской частицы
Ю2
1 ш2
D = (У2) = -J Gv (ra)d
ш
и выражения (12), (14) и (16), для низкочастотного случая (ш ^ у и ш ^ vT) найдем
D =( I i)
Dl =(12) = 2kT {q+V+(n+f + qv (« )2J(ö)i-®i);
>("+)2(q+)3 -)2(q-)3]ln
Ql
2S 2U 0
l2
rn v^ m vT
ln
Тогда отношение дисперсии флуктуаций тока в неравновесном и равновесном случаях
г2 ||д+(п+)2 (q+)3 m vT[n )2 (q )3 m+v| 1 ln —
§d; 7 ^—' ; v " 7—^(i9)
m m+v|vT
_^ ^ ■ -И2
kTl2 (qV+(n+)2 + q"|д"(п")2\m~m+v+v7 (®2 )
Если параметры положительных и отрицательных ионов в первом приближении полагать одинаковыми, то формула (19) приобретает вид
5d=^q2uLln (ю2/ ю1). (20)
mkTl2vT (ш2 -ш1)
Подстановка в формулу (19) реальных данных для эксперимента, описанного в работе [21], дает оценку интенсивности флуктуаций vT «1018...1019 с-1, что соответствует оценке постоянной времени хаоти-зации ионов в электролите То ~ 10_19...10~18 с.
Заключение. Предложенное описание броуновского движения в неравновесной среде как необратимого немарковского процесса позволило получить непротиворечивые оценки константы Хоуге и постоянной времени хаотизации ионов в электролите.
ЛИТЕРАТУРА
[1] Бочков Г.Н., Кузовлев Ю.Е. Новое в исследованиях 1//-шума. УФН, 1983, т. 141, № 1, с. 151-176.
[2] Кузовлев Ю.Е. Почему природе нужен 1//-шум. УФН, 2015, т. 185, № 7, с. 773-783.
[3] Hooge F.N., Gaal J.L. Fluctuations with a 1// spectrum in the conductance of ionic solutions and in the voltage of concentration cells. Philips Res. Rep., 1971, vol. 26, no. 2, pp. 77-90.
[4] Bezrukov S.M., Pustovoit M.A., Sibilev A.I., et al. Large-scale conductance fluctuations in solutions of strong electrolytes. Physica B: Condens. Matter, 1989, vol. 159, iss. 3, pp. 388-398. DOI: 10.1016/0921-4526(89)90016-1
[5] van den Berg R.J., de Vos A., de Goede J. Electrical noise in solutions of hydrochloric acid in ethanol. Phys. Lett. A, 1989, vol. 139, iss. 5-6, pp. 249-252.
DOI: 10.1016/0375-9601(89)90149-7
[6] Hooge F.N. 1/f noise sources. IEEE Trans. Electron. Devices, 1994, vol. 41, iss. 11, pp. 1926-1935. DOI: 10.1109/16.333808
[7] Морозов А.Н. О броуновском движении в среде с флуктуирующими коэффициентами переноса. Известия вузов. Физика, 1986, № 6, с. 90-91.
[8] Морозов А.Н. Применение теории немарковских процессов при описании броуновского движения. ЖЭТФ, 1996, т. 109, № 4, с. 1304-1315.
[9] Lenzi E.K., Evangelista L.R., Lenzi M.K., et al. Solutions for a non-Markovian diffusion equation. Phys. Lett. A, 2010, vol. 374, iss. 41, pp. 4193-4198.
DOI: 10.1016/j.physleta.2010.08.049
[10] Морозов А.Н., Скрипкин А.В. Применение интегральных преобразований для описания броуновского движения как немарковского случайного процесса. Известия вузов. Физика, 2009, т. 52, № 2, с. 66-74.
[11] Morozov A.N., Skripkin A.V. Spherical particle Brownian motion in viscous medium as non-Markovian random process. Phys. Lett. A, 2011, vol. 375, iss. 46, pp. 4113-4115. DOI: 10.1016/j.physleta.2011.10.001
[12] Морозов А.Н., Скрипкин А.В. Описание испарения сферической частицы жидкости как немарковского случайного процесса с использованием интегральных стохастических уравнений. Известия вузов. Физика, 2010, № 11-2, с. 55-64.
[13] Mura A., Taqqu M.S., Mainardi F. Non-Markovian diffusion equations and processes: analysis and simulations. Physica A, 2008, vol. 387, iss. 21, pp. 5033-5064.
DOI: 10.1016/j.physa.2008.04.035
[14] Морозов А.Н., Скрипкин А.В. Распространение тепла в пространстве вокруг цилиндрической поверхности как немарковский случайный процесс. Инженерно-физический журнал, 2011, т. 84, № 6, с. 1121-1127.
[15] Lisy V., Tóthóva J., Glod L. On the correlation properties of thermal noise in fluids. Int. J. Thermophys, 2013, vol. 34, iss. 4, pp. 629-641. DOI: 10.1007/s10765-012-1290-1
[16] Морозов А.Н., Скрипкин А.В. Равновесные флуктуации температуры молекулярного и фотонного газов в сферической микрополости. Известия вузов. Физика, 2012, № 7, c. 9-18.
[17] Морозов А.Н. Метод описания немарковских процессов, задаваемых системой линейных интегральных уравнений. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2017, № 5, с. 57-66. DOI: 10.18698/1812-3368-2017-5-57-66
[18] Морозов А.Н., Скрипкин А.В. Немарковские физические процессы. М., Физматлит, 2018.
[19] Морозов А.Н. Предварительные результаты измерений меры Кульбака флук-туаций напряжения на электролитической ячейке. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2011, № 2, с. 16-24.
[20] Безруков С.М., Ирхин А.И., Сибилев А.И. Верхняя оценка для интенсивности 1//-шума электролитов: эксперименты с молекулярными каналами. Препринт ЛИЯФ-1190. Л., ЛИЯФ, 1986.
[21] Morozov A.N. Nonlocal influences of natural dissipative processes on the Kullback measure of voltage fluctuations on an electrolytic cell. NeuroQuantology, 2016, vol. 14, no. 3, pp. 477-483. DOI: 10.14704/nq.2016.14.3.920
Морозов Андрей Николаевич — д-р физ.-мат. наук, профессор, заведующий кафедрой «Физика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5, стр. 1).
Просьба ссылаться на эту статью следующим образом:
Морозов А.Н. Броуновское движение как необратимый немарковский процесс. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2019, № 2, с. 94-103. DOI: 10.18698/1812-3368-2019-2-94-103
BROWNIAN MOTION AS AN IRREVERSIBLE NON-MARKOVIAN PROCESS
A.N. Morozov amor59@mail.ru
amor@bmstu.ru
Bauman Moscow State Technical University, Moscow, Russian Federation
Abstract
The paper presents a method of describing Brownian motion in a non-equilibrium medium for the case of irreversible processes. We computed spectral density of velocity fluctuations for a Brownian particle in a non-equilibrium medium and determined that in the low-frequency region it is represented by flicker noise. We employed the method we developed to describe Brownian motion in a non-equilibrium medium to compute fluctuations of current in a small volume of an electrolyte. We derived estimations of the Hooge parameter magnitude and the randomisation time constant for ions in an electrolyte, which match the estimations obtained via experiments
Keywords
Non-Markovian process, Brownian motion, non-equilibrium state, entropy generation, flicker noise
Received 28.02.2018 © Author(s), 2019
REFERENCES
[1] Bochkov G.N., Kuzovlev Yu.E. New aspects in 1/f noise studies. Phys.-Uspekhi, 1983, vol. 26, no. 9, pp. 829-844. DOI: 10.1070/PU1983v026n09ABEH004497
[2] Kuzovlev Yu.E. Why nature needs 1/f noise. Phys.-Uspekhi, 2015, vol. 58, no. 7, pp. 719-729. DOI: 10.3367/UFNe.0185.201507d.0773
[3] Hooge F.N., Gaal J.L. Fluctuations with a 1/f spectrum in the conductance of ionic solutions and in the voltage of concentration cells. Philips Res. Rep., 1971, vol. 2б, no. 2, pp. 77-90.
[4] Bezrukov S.M., Pustovoit M.A., Sibilev A.I., et al. Large-scale conductance fluctuations in solutions of strong electrolytes. Physica B: Condens. Matter, 1989, vol. 159, iss. 3, pp. 388-398. DOI: 10.101б/0921-452б(89)9001б-1
[5] van den Berg R.J., de Vos A., de Goede J. Electrical noise in solutions of hydrochloric acid in ethanol. Phys. Lett. A, 1989, vol. 139, iss. 5-б, pp. 249-252.
DOI: 10.101б/0375-9б01(89)90149-7
[6] Hooge F.N. 1/f noise sources. IEEE Trans. Electron. Devices, 1994, vol. 41, iss. 11, pp. 192б-1935. DOI: 10.1109/1б.333808
[7] Morozov A.N. On Brownian motion in medium with fluctuating transfer coefficient. Izvestiya vuzov. Fizika, 198б, no. б, pp. 90-91 (in Russ.).
[8] Morozov A.N. Use of the theory of non-Markovian processes in the description of Brownian motion. J. Exp. Theor. Phys., 199б, vol. 82, iss. 4, pp. 703-708.
[9] Lenzi E.K., Evangelista L.R., Lenzi M.K., et al. Solutions for a non-Markovian diffusion equation. Phys. Lett. A, 2010, vol. 374, iss. 41, pp. 4193-4198.
DOI: 10.1016/j.physleta.2010.08.049
[10] Morozov A.N., Skripkin A.V. Application of integral transforms to a description of the Brownian motion by a non-Markovian random process. Russ. Phys. J., 2009, vol. 52, iss. 2, pp. 184-195. DOI: 10.1007/s11182-009-9217-4
[11] Morozov A.N., Skripkin A.V. Spherical particle Brownian motion in viscous medium as non-Markovian random process. Phys. Lett. A, 2011, vol. 375, iss. 4б, pp. 4113-4115. DOI: 10.1016/j.physleta.2011.10.001
[12] Morozov A.N., Skripkin A.V. Description of spherical liquid part evaporation as non-Markovian process using stochastic integral equations. Izvestiya vuzov. Fizika, 2010, no. 11-2, pp. 55-б4 (in Russ.).
[13] Mura A., Taqqu M.S., Mainardi F. Non-Markovian diffusion equations and processes: analysis and simulations. Physica A, 2008, vol. 387, iss. 21, pp. 5033-50б4.
DOI: 10.1016/j.physa.2008.04.035
[14] Morozov A.N., Skripkin A.V. Propagation of heat in the space around a cylindrical surface as a non-Markovian random process. J. Eng. Phys. Thermophys., 2011, vol. 84, iss. 6, pp. 1201-1208. DOI: 10.1007/s10891-011-0585-6
[15] Lisy V., Tóthóva J., Glod L. On the correlation properties of thermal noise in fluids. Int. J. Thermophys, 2013, vol. 34, iss. 4, pp. 629-641. DOI: 10.1007/s10765-012-1290-1
[16] Morozov A.N., Skripkin A.V. Equilibrium temperature fluctuations of molecular and photonic gases in a spherical microcavity. Russ. Phys. J., 2012, vol. 55, iss. 7, pp. 736-747. DOI: 10.1007/s11182-012-9875-5
[17] Morozov A.N. Method for describing non-Markovian processes defined by a system of linear integral equations. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2017, no. 5, pp. 57-66 (in Russ.).
DOI: 10.18698/1812-3368-2017-5-57-66
[18] Morozov A.N., Skripkin A.V. Nemarkovskie fizicheskie protsessy [Non-Markovian physical processes]. Moscow, Fizmatlit Publ., 2018.
[19] Morozov A.N. Preliminary results of recording the Kullback measure of voltage fluctuations on electrolytic cell. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2011, no. 2, pp. 16-24 (in Russ.).
[20] Bezrukov S.M., Irkhin A.I., Sibilev A.I. Verkhnyaya otsenka dlya intensivnosti 1/f-shuma elektrolitov: eksperimenty s molekulyarnymi kanalami [Upper estimate for intensity of electrolyte 1/f noise: experiments with molecular channels]. Preprint LIYaF-1190. Leningrad, LIYaF Publ., 1986 (in Russ.).
[21] Morozov A.N. Nonlocal influences of natural dissipative processes on the Kullback measure of voltage fluctuations on an electrolytic cell. NeuroQuantology, 2016, vol. 14, no. 3, pp. 477-483. DOI: 10.14704/nq.2016.14.3.920
Morozov A.N. — Dr. Sc. (Phys.-Math.), Professor, Head of Department of Physics, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Please cite this article in English as:
Morozov A.N. Brownian motion as an irreversible non-Markovian process. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2019, no. 2, pp. 94-103 (in Russ.). DOI: 10.18698/1812-3368-2019-2-94-103
