Смещение и затухание

Подсказка

Индуцированные на диске заряды можно считать точечными на большом расстоянии; для малых отклонений линеаризуйте силу \(F_e(z)\) около положения равновесия, чтобы найти смещение и изменение частоты, а затухание связано с токами в диске при движении шарика.

Ответ

2.1) \(\nu_0=\dfrac1{2\pi}\sqrt{k/m}\);  2.2) \(F_e=\dfrac{q^2V}{8\pi^2\varepsilon_0z^5}\), \(V=\pi r^2h\);  2.3) \(q_{\max}=\sqrt{\dfrac{5^5}{6^6}\dfrac{8\pi^2\varepsilon_0z_0^5(kz_0)}{V}}\) (\(\gamma^*=5^5/6^6\approx0{,}067\));  2.4) \(\xi\approx\gamma=1{,}0\cdot10^{-2}\);  2.5) \(\dfrac{\Delta\nu}{\nu_0}\approx-\dfrac52\gamma=-2{,}5\%\);  2.6) \(n\approx2{,}9\cdot10^7\) периодов.

Решение

2.1. Выразите частоту \(\nu_0\) колебаний шарика через массу шарика и жесткость пружины. Уравнение движения шарика имеет вид

\[ma = -kx,\tag{1}\]

из которого непосредственно следует формула для частоты

\[\nu _ { 0 } = \frac { 1 } { 2 \pi } \sqrt { \frac { k } { m } } \, .\tag{2}\]

2.2. Получите выражение для силы \(F_e(z)\), действующей на неподвижный шарик со стороны диска, как функцию его расстояния до центра диска \(z\) (для \(z \gg r\)).

Электрическое поле внутри проводящего диска должно отсутствовать. Под действием поля \(\vec{E}_0\) точечного заряда \(q\) на поверхности диска будут индуцироваться заряды с поверхностной плотностью \(\pm\sigma\), которые будут создавать поле \(\vec{E}'\), равное по величине и противоположное по направлению поля точечного заряда. Учитывая, что размеры диска малы по сравнению с расстоянием до заряда, можно пренебречь изменением вектора напряженности \(\vec{E}_0\) в пределах диска. Поэтому и поле \(\vec{E}'\) следует считать однородным, а поверхностную плотность индуцированных зарядов постоянной.

Напряженность поля точеного заряда определяется формулой \(E_0 = \frac{q}{4\pi\varepsilon_0 z^2}\), а напряженность поля индуцированных зарядов связана с их поверхностной плотностью соотношением \(E' = \frac{\sigma}{\varepsilon_0}\). Приравнивая эти величины, определим поверхностную плотность индуцированных зарядов \(\sigma' = \frac{q}{4\pi z^2}\), и величину индуцированного заряда на одной стороне диска \(q' = \sigma' S = \frac{qS}{4\pi z^2}\), (где \(S = \pi r^2\) - площадь диска).

Интересно отметить, что индуцированный дипольный момент диска равен \(p' = q'h = \frac{qSh}{4\pi z^2} = \frac{q}{4\pi z^2}V\). Из сравнения с традиционным выражением для этой величины \(p = \alpha\varepsilon_0 E_0\) следует, что поляризуемость диска в точности равна его объему.

Теперь с помощью закона Кулона легко найти силу, действующую на точечный заряд со стороны двух индуцированных зарядов (который на большом расстоянии также можно считать точечными)

\[F_e = \frac{qq'}{4\pi\varepsilon_0} \left( \frac{1}{z^2} - \frac{1}{(z+h)^2} \right) = \frac{qq'}{4\pi\varepsilon_0} \cdot \frac{(z+h)^2 - z^2}{z^2(z+h)^2}.\]

Учитывая, что \(h \ll z\), окончательно получим

\[F_e = \frac{qq'}{4\pi\varepsilon_0} \cdot \frac{(z+h)^2 - z_0^2}{z^2(z+h)^2} \approx \frac{qq'h}{2\pi\varepsilon_0 z^3} = \frac{q}{2\pi\varepsilon_0 z^3} \cdot \frac{q}{4\pi z^2} V = \frac{q^2 V}{8\pi^2 \varepsilon_0 z^5}.\tag{3}\]

где \(V = \pi r^2 h\) - объем диска.

2.3. Получите формулу для определения максимального заряда шарика, при котором он может совершать колебательное движение.

Необходимым и достаточным условием возможности колебаний является наличие положения устойчивого равновесия.

Условию равновесия соответствует равенство сил упругости и электростатического взаимодействия (3), то есть для координаты точки равновесия справедливо уравнение

\[kx = \frac{q^2}{8\pi^2\varepsilon_0(z_0 - x)^5}.\tag{4}\]

Это уравнение является алгебраически уравнением шестой степени, поэтому не может быть решено элементарными методами. Преобразуем уравнение к виду

\[kz_0 \frac{x}{z_0} = \frac{q^2V}{8\pi^2\varepsilon_0 z_0^5 \left(1 - \frac{x}{z_0}\right)^5},\]

обозначим \(\xi = \frac{x}{z_0}\) и введем безразмерный параметр рассматриваемой системы

\[\gamma = \frac{q^2V}{8\pi^2\varepsilon_0 z_0^5 (kz_0)}.\tag{5}\]

В этих обозначениях уравнение (4) имеет простой вид

\[\xi = \frac{\gamma}{(1 - \xi)^5}.\tag{6}\]

Колебания шарика возможны до тех пор, пока существует положение устойчивого равновесия. Иными словами, пока уравнение (6) имеет решения на интервале \(0 \le \xi \le 1\). Функция \(F(\xi) = \frac{\gamma}{(1 - \xi)^5}\) резко возрастает при приближении аргумента к единице. При малых значениях параметра \(\gamma\) эта функция имеет две точки пересечения с прямой \(F = \xi\) (то есть имеется два корня уравнения), при больших значениях \(\gamma\) корней нет – на рисунке \(\gamma_1 > \gamma_0\).

Из анализа зависимостей сил упругости и электрического притяжения следует, что меньший корень соответствует положению устойчивого положения равновесия, а больший неустойчивого. «Критическим» значением является такое значение \(\gamma^*\), при котором прямая \(F = \xi\) касается графика функции \(F(\xi)\) в одной точке \(\xi^*\). Таким

образом, в этой точке производная от этой функции должна равняться \(F'(\xi)=1\). Итак, для определения этого значения имеем два уравнения

\[\xi^* = \frac{\gamma^*}{(1-\xi^*)^5}; \quad 1 = \frac{5\gamma^*}{(1-\xi^*)^6},\]

Разделив первое из них на второе, получим \(\xi^* = \frac{1-\xi^*}{5}\), откуда следует \(\xi^* = \frac{1}{6}\),

а искомое значение параметра выражается «очень красиво»

\[\gamma^* = \xi^* \cdot (1-\xi^*)^5 = \frac{5^5}{6^6} \approx 0{,}067.\]

Используя выражение для параметра (5), находим максимальное значение заряда

\[q_{\text{max}} = \sqrt{\frac{5^5}{6^6} \frac{8\pi^2 \varepsilon_0 z_0^5 (kz_0)}{V}}.\tag{7}\]

Возможен еще один способ определения условия существования корней уравнения (6). Преобразуем это уравнение к виду

\[\xi(1-\xi)^5 = \gamma.\]

Ясно, что это уравнение будет иметь корни, если максимум функции \(f(\xi) = \xi(1-\xi)^5\) превышает значение параметра \(\gamma\). Этот метод приводит к тому конечному результату.

2.4. Чему равно относительное смещение положения равновесия шарика \(\xi = \frac{x}{z_0}\)?

Заметим, что приведенное в условии численное значение параметра \(\gamma\) меньше критического \(\gamma^* = \frac{5^5}{6^6} \approx 0{,}067\), поэтому при таком заряде существует положение устойчивого равновесия шарика. Его координата является меньшим корня уравнения (6).

С точностью до малых величин первого порядка, его решение элементарно

\[\xi = \frac{x}{z_0} \approx \gamma = 1{,}0 \cdot 10^{-2}.\tag{8}\]

2.5. Чему равно относительно изменение частоты колебаний шарика \(\frac{\Delta\nu}{\nu_0}\)?

Обозначим расстояние от диска до положения равновесия заряженного шарика \(z_1\). Переместим начало отсчета в эту точку. Уравнение движения шарика при малом отклонении \(x\) от положения равновесия имеет вид (при преобразовании функций опять используется малость отклонений \(x \ll z_1\)):

\[m a = - k ( z _ { 1 } + x ) + \frac { q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } ( z _ { 1 } - x ) ^ { 5 } } = - k z _ { 1 } - k x + \frac { q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 5 } \left( 1 - \displaystyle \frac { x } { z _ { 1 } } \right) ^ { 5 } } \approx\]

\[\approx - k z _ { 1 } - k x + \frac { q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 5 } } \left( 1 + 5 \frac { x } { z _ { 1 } } \right) = - k z _ { 1 } + \frac { q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 5 } } - \left( k - \frac { 5 q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 6 } } \right) x\]

Первые два слагаемых взаимно уничтожаются, так как точка \(z _ { 1 }\) - положение равновесия. Таким образом, получаем уравнение гармонических колебаний

\[a = - \frac { k } { m } \left( 1 - \frac { 5 q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } k z _ { 1 } ^ { 5 } } \right) x\]

с частотой \(\omega = { \sqrt { { \frac { k } { m } } \left( 1 - { \frac { 5 q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 6 } k } } \right) } } = { \sqrt { \omega _ { 0 } ^ { 2 } \left( 1 - { \frac { 5 q ^ { 2 } V } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 6 } k } } \right) } } \approx \omega _ { 0 } \left( 1 - { \frac { 5 q ^ { 2 } V } { 1 6 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 6 } k } } \right) .\)

Следовательно, относительное изменение частоты колебаний равно

\[\frac { \Delta \nu } { \nu _ { 0 } } = \frac { \Delta \omega } { \omega _ { 0 } } = - \frac { 5 q ^ { 2 } V } { 1 6 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 1 } ^ { 6 } k } \approx - \frac { 5 } { 2 } \gamma = - 2 { , } 5 \% \, .\tag{8}\]

2.6. Оцените, за сколько периодов колебаний заряженного шарика их амплитуда уменьшится \(\eta = 1 { , } 0 \%\) .

Так как индуцированные на поверхности диска заряды зависят от положения заряженного шарика, то при его движении величины этих зарядов будут изменяться, то есть по диску будут протекать электрические токи. Причиной затухания колебаний является потеря энергии на выделение джоулевой теплоты при протекании электрического тока при «перезарядке» диска. При малом сопротивлении диска можно пользоваться «квазистационарным» приближением, то есть считать, что заряд диска соответствует равновесному значению при заданном положении шарика. Величина равновесного индуцированного заряда была подсчитана ранее и равна \(q ^ { \prime } = \frac { q S } { 4 \pi z ^ { 2 } }\). Сила тока протекающего поперек диска легко вычисляется

\[I = \frac { d q ^ { \prime } } { d t } = \frac { d q ^ { \prime } } { d z } \frac { d z } { d t } = \frac { q S } { 2 \pi z ^ { 3 } } v \, ,\tag{9}\]

где \(v = - \frac { d z } { d t }\) - скорость шарика. Так как колебания шарика малы, то можно пренебречь изменением расстояния \(z\) в знаменателе формулы (9) и положить \(z ( t ) \approx z _ { 1 }\). В этом случае мощность теплоты, выделяющейся в диске по закону Джоуля-Ленца равна

\[P = I ^ { 2 } R = \left( \frac { q S } { 2 \pi z ^ { 3 } } \right) ^ { 2 } \rho \frac { h } { S } v ^ { 2 } \approx \frac { q ^ { 2 } S h \rho } { 4 \pi ^ { 2 } z _ { 1 } ^ { 6 } } v ^ { 2 } \, .\]

Если скорость шарика изменяется по закону \(v ( t ) = v _ { 0 } \cos ( \omega t + \varphi )\), то среднее за период значение квадрата скорости равно \(\left\langle v ^ { 2 } ( t ) \right\rangle = \frac { 1 } { 2 } v _ { 0 } ^ { 2 }\). Потери энергии за один период колебаний равны

\[\Delta W = \langle P \rangle T = \frac { q ^ { 2 } S h \rho } { 4 \pi ^ { 2 } z _ { 1 } ^ { 6 } } \Bigl < v ^ { 2 } \Bigr > T = \frac { q ^ { 2 } S h \rho } { 4 \pi ^ { 2 } z _ { 1 } ^ { 6 } } \cdot \frac 1 2 v _ { 0 } ^ { 2 } \cdot \frac { 2 \pi } { \omega } \, .\tag{10}\]

Из предыдущего рассмотрения следует, что изменения положения равновесия и частоты малы, поэтому в данной формуле можно положить \(z _ { 1 } \approx z _ { 0 }\) , \(\omega \approx \omega _ { 0 }\) . Начальная энергия колебаний \(W _ { 0 } = \frac { m v _ { 0 } ^ { 2 } } { 2 }\) , поэтому формулу (10) можно записать в виде

\[\Delta W = \frac { q ^ { 2 } S h \rho } { 4 \pi ^ { 2 } z _ { 0 } ^ { 6 } } \cdot \frac 1 2 m v _ { 0 } ^ { 2 } \cdot \frac { 2 \pi } { m \omega _ { 0 } } = \frac { q ^ { 2 } S h \rho } { 8 \pi ^ { 2 } \varepsilon _ { 0 } z _ { 0 } ^ { 6 } k } \cdot \frac 1 2 m v _ { 0 } ^ { 2 } \cdot \frac { 4 \pi \rho \varepsilon _ { 0 } } { \omega _ { 0 } } \cdot \frac { k } { m } = 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma W _ { 0 } \, .\tag{10a}\]

Как следует из этого выражения, потери энергии за период пропорциональны начальной энергии, следовательно, энергия как функция числа колебаний будет убывать в геометрической прогрессии, действительно, за один период энергия уменьшится до величины

\[\Delta W = W _ { 0 } - W _ { 1 } = 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma W _ { 0 } \quad \Rightarrow \quad W _ { 1 } = \left( 1 - 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \right) W _ { 0 } \, ,\]

после \(n\) колебаний энергия станет равной

\[W _ { n } = \left( 1 - 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \right) ^ { n } W _ { 0 } \, .\tag{11}\]

Так как энергия колебаний пропорциональна квадрату амплитуды, то из формулы (11) следует, что амплитуда колебаний зависит от числа колебаний по закону

\[A _ { n } = A _ { 0 } \left( 1 - 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \right) ^ { \frac { n } { 2 } } .\]

Учитывая, что рассматриваемое изменение амплитуды колебаний мало, то очередной раз можно воспользоваться приближением

\[A _ { n } = A _ { 0 } \left( 1 - 4 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \right) ^ { \frac { n } { 2 } } \approx A _ { 0 } \left( 1 - 2 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \cdot n \right) ,\]

или относительное уменьшение амплитуды

\[\frac { A _ { 0 } - A _ { n } } { A _ { 0 } } = \eta = 2 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma \cdot n \, .\]

Таким образом, искомое число колебаний равно

\[n = \frac { \eta } { 2 \pi \rho \varepsilon _ { 0 } \omega _ { 0 } \gamma } \approx 2 { , } 9 \cdot 1 0 ^ { 7 } \, ,\tag{12}\]

что соответствует времени более одного месяца.

2.7. Оцените, при каком удельном электрическом сопротивлении материала диска затухание колебаний шарика будет максимальным.

Проведенное выше описание затухание колебаний не полно.

Оно справедливо только для малых удельных сопротивлений. Действительно, из полученных выражений следует, что для непроводящей пластинки \((\rho \to \infty)\), колебания затухают мгновенно, хотя в этом случае отсутствуют механизмы потери энергии, следовательно, колебания должны быть вообще незатухающими.

Для более точного описания затухания необходимо более точно рассмотреть динамику изменения

поверхностной плотности индуцированных зарядов на поверхности диска. Для этого воспользуемся законом Ома\(^1\) \(j = \frac{1}{\rho} E\) и связью плотности тока с изменением поверхностной плотности заряда \(j = \frac{d\sigma}{dt}\). Учитывая, что поле внутри диска создается как полем заряда шарика, так и индуцированными зарядами, запишем уравнение для изменения поверхностной плотности индуцированных зарядов

\[\frac{d\sigma}{dt} = \frac{1}{\rho} \left( E - \frac{\sigma}{\varepsilon_0} \right) = \frac{1}{\rho \varepsilon_0} \left( \frac{q}{4\pi (z_1 - x)^2} - \sigma \right).\tag{13}\]

Обозначим \(\tau = \rho \varepsilon_0\) и упростим это уравнение, используя стандартную методику

\[\tau \frac{d\sigma}{dt} = \frac{q}{4\pi (z_1 - x)^2} - \sigma = \frac{q}{4\pi z_1^2} \left( 1 - \frac{x}{z_1} \right)^{-2} - \sigma = \frac{q}{4\pi z_1^2} \left( 1 + 2 \frac{x}{z_1} \right) - \sigma =\]

\[= \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x - \sigma\]

Плотность тока перезарядки (следовательно, и потери энергии) определяются переменной составляющей плотности заряда, поэтому представим \(\sigma = \frac{q}{4\pi z_1^2} + s(t)\), где \(s(t)\) - изменяющаяся часть поверхностной плотности заряда. Для последней справедливо уравнение

\[\tau \frac{ds}{dt} = \frac{q}{2\pi z_1^3} x - s.\]

Как мы показали выше, затухание колебаний является крайне малым, поэтому координату шарика с высокой точностью можно представить в виде\(^2\)

\[x = A \cos \omega t,\]

тогда для определения \(s(t)\) получим уравнение

\[\tau \frac{ds}{dt} + s = \frac{q}{2\pi z_1^3} A \cos \omega t,\tag{14}\]

которое проще всего решить методом векторных диаграмм. Так, при \(s = B \cos(\omega t + \varphi)\)

\[\tau \frac{ds}{dt} = -B \omega \sin(\omega t + \varphi) = B \omega \cos\left( \omega t + \varphi + \frac{\pi}{2} \right).\]

Тогда

векторная диаграмма, описывающее уравнение (14) имеет вид, показанный на рисунке (где обозначено \(A_1 = \frac{q}{2\pi z_1^3} A\)). Из данной диаграммы следует

— \(^1\) Можно использовать закон Ома и в «обычной» форме – «для участка цепи». \(^2\) Конечно, в качестве частоты колебаний можно взять собственную частоту колебаний \(\omega_0\)

\[B ^ { 2 } \Big ( 1 + \big ( \omega \tau \big ) ^ { 2 } \Big ) = A _ { 1 } ^ { 2 } ,\]

\[t g \varphi = - \omega \tau\]

\[\mathrm { \Pi } \mathrm { \Pi } \mathrm { \Pi } \ s = \frac { A } { \sqrt { 1 + \left( \omega \tau \right) ^ { 2 } } } \frac { q } { 2 \pi z _ { 1 } ^ { 3 } } \cos ( \omega t + \varphi ) .\]

Далее, аналогично выводу соотношения (10), можно утверждать, что потери энергии за период пропорциональны (оставлены только сомножители, зависящие от удельного электрического сопротивления)

\[\Delta W \propto B ^ { 2 } \omega \tau \propto \frac { \omega \tau } { 1 + \left( \omega \tau \right) ^ { 2 } } .\]

Данная функция достигает максимума при \(\omega \tau = 1\), следовательно, при этом условии затухание и будет максимальным. Таким образом, получаем \(2 \pi \nu _ { 0 } \rho \varepsilon _ { 0 } = 1\), или окончательно

\[\rho = \frac { 1 } { 2 \pi \nu _ { 0 } \varepsilon _ { 0 } } \approx 1 { , } 8 \cdot 1 0 ^ { 9 } \ O M \cdot { \scriptscriptstyle M } .\tag{15}\]

2.8. Качественно опишите, какие именно свойства индуцированных зарядов (и их изменения) обеспечивают появление сил, действующих на шарик, приводящих к а) смещению положения равновесия; б) изменению частоты колебаний; в) затуханию колебаний.

На основании предыдущих рассуждений следует, что а) постоянная составляющая поверхностной плотности заряда приводит к смещению положения равновесия; б) изменение поверхностной плотности в «квазистационарном приближении» приводит к изменению частоты колебаний; в) Заметим, что если все силы, действующие на шарик, зависят только от координат, то затухание колебаний отсутствует. Поэтому следует искать причины, приводящие к появлению составляющих сил, зависящих от скорости движения шарика. Причиной появления таких слагаемых является «запаздывание» - распределение индуцированных зарядов отстает от равновесного распределения, соответствующего положению шарика в «данный момент времени». Наличие «запаздывания» распределения зарядов, по отношению к положению заряда приводят к затуханию колебаний.

Дополнение, не входящее в оцениваемое решение.} Эти утверждения можно строго подтвердить на «языке уравнений» движения шарика и изменения поверхностной плотности индуцированных зарядов.

Сила электрического взаимодействия \(F \approx \frac{qq'h}{2\pi\varepsilon_0 z^3} = \frac{qV}{2\pi\varepsilon_0 z^3} \cdot \sigma\)

Уравнение колебаний шарика

\[ma = -k(z_1 + x) + \frac{qV}{2\pi\varepsilon_0 (z_1 - x)^3} \cdot \sigma = -k(z_1 + x) + \frac{qV\sigma}{2\pi\varepsilon_0 z_1^3} \cdot \left(1 - \frac{x}{z_1}\right)^{-3} =\]

\[= -k(z_1 + x) + \frac{qV}{2\pi\varepsilon_0 z_1^3} \sigma + \frac{3qV}{2\pi\varepsilon_0 z_1^4} \sigma x\]

Уравнение для изменения поверхностной плотности индуцированных зарядов

\[\frac{d\sigma}{dt} = \frac{1}{\rho} \left( E - \frac{\sigma}{\varepsilon_0} \right) = \frac{1}{\rho\varepsilon_0} \left( \frac{q}{4\pi(z_1 - x)^2} - \sigma \right),\]

или, с учетом использованных обозначений,

\[\tau \frac{d\sigma}{dt} = \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x - \sigma\]

Из этого уравнения выразим \(\sigma = \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x - \tau \frac{d\sigma}{dt}\), и подставим в уравнение колебаний шарика

\[ma = -k(z_1 + x) + \frac{qV}{2\pi\varepsilon_0 z_1^3} \sigma + \frac{3qV}{2\pi\varepsilon_0 z_1^4} \sigma x =\]

\[= -k(z_1 + x) + \frac{qV}{2\pi\varepsilon_0 z_1^3} \left( \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x - \tau \frac{d\sigma}{dt} \right) + \frac{3qV}{2\pi\varepsilon_0 z_1^4} x \left( \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x - \tau \frac{d\sigma}{dt} \right) =\]

\[= -kz_1 + \frac{q^2 V}{8\pi^2 \varepsilon_0 z_1^5} - \left( k - \frac{q^2 V}{4\pi^2 \varepsilon_0 z_1^6} - \frac{3q^2 V}{8\pi^2 \varepsilon_0 z_1^6} \right) x - \frac{qV}{2\pi\varepsilon_0 z_1^3} \tau \frac{d\sigma}{dt} + \frac{3qV}{2\pi\varepsilon_0 z_1^4} x \left( \frac{q}{2\pi z_1^3} x - \tau \frac{d\sigma}{dt} \right)\]

Первые два слагаемых взаимно уничтожаются \(-kz_1 + \frac{q^2 V}{8\pi^2 \varepsilon_0 z_1^5} = 0\) (уравнение для определения смещенного положения равновесия).

Множитель при третьем слагаемом определяет изменение частоты

\[-\left( k - \frac{q^2 V}{4\pi^2 \varepsilon_0 z_1^6} - \frac{3q^2 V}{8\pi^2 \varepsilon_0 z_1^6} \right) x = -\left( k - \frac{5q^2 V}{8\pi^2 \varepsilon_0 z_1^6} \right) x = -m\omega_1^2 x\]

Следующее, содержащее скорость изменения поверхностной плотности, ответственно за затухание. Остальными следует пренебречь, как имеющие более высокий порядок малости.

Итого имеем уравнение колебаний.

\[ma = -m\omega_1^2 x - \frac{qV}{2\pi\varepsilon_0 z_1^3} \tau \frac{d\sigma}{dt}.\]

В использованном приближении, уравнение для плотности заряда можно решать методом последовательных приближений. Пренебрегая «запаздыванием» получаем

\[\sigma_{(0)} \approx \frac{q}{4\pi z_1^2} + \frac{q}{2\pi z_1^3} x,\]

тогда \(\tau \frac{d\sigma}{dt} \approx \tau \frac{d\sigma_{(0)}}{dt} = \frac{q}{2\pi z_1^3} \tau \frac{dx}{dt}\). Подставляя в уравнение движения шарика, получаем уравнение

\[ma = -m\omega_1^2 x - \frac{qV}{2\pi\varepsilon_0 z_1^3} \tau \frac{d\sigma}{dt} \approx -m\omega_1^2 x - \frac{qV}{2\pi\varepsilon_0 z_1^3} \cdot \frac{q}{2\pi z_1^3} \tau \frac{dx}{dt} = -m\omega_1^2 x - \frac{q^2 V \rho}{4\pi^2 z_1^6} \cdot \frac{dx}{dt},\]

в котором в явном виде присутствует слагаемое, «ответственное» за затухание.

Прочитали решение — спросите себя, как можно было до этого додуматься. А через день-два решите задачу заново, не подглядывая.