Редактирование: Разработка электро-механических моделей углеродных нановискеров как детекторов масс

'''МАГИСТЕРСКАЯ РАБОТА'''<br>
''Автор работы'': [[Ванюшкина Валентина]]<br>
''Научный руководитель'': [[Л.В. Штукин]]<br>

==Введение==
Наноструктуры - объекты, созданные человеком и имеющие малые размеры (от
нескольких нанометров до сотен микрометров). Эта область науки открывает много новых
и интересных возможностей для экспериментов и теоретических расчётов. Наноструктуры
могут быть очень полезны для медицины, механики, инженерии и биологии: могут быть ис-
пользованы для разрезания клеток, взвешивания объектов на нано- и микро- уровнях, а так
же для построения более сложных структур.
Одной из технологий для создания наноструктур является напыление аморфного углерода
на металлическую иглу.
Целью данной работы является подробное исследование двух углеродных наноструктур
Одна из них это прямолинейный одиночный вискер. Длина такого вискера порядка несколь-
ких микрометров, а максимальная толщина всего лишь несколько десятков нанометров. Вис-
кер выращен на вершине металлической иглы и закреплен на ней. Вторая исследуемая кон-
струкция, это система из двух вискеров, связанных между собой упругой перемычкой. Пере-
мычка так же создана из углерода путём напыления.
Задачи данной работы заключаются в исследовании и описании одного из экспериментов
с одиночным вискером на игле с помощью применения приближённого метода решения; опи-
сание модели и изучение дифференциального резонатора из нановискеров для возможного
применения его в качестве детектора масс; изучения одного из способов точного измерения
колебаний нановискера и дифференциального резонатора с помощью различных емкостных
датчиков.
<gallery>
File:Ishape1.png|Одиночный вискер
File:Doublewhisker.png|Двойной вискер
</gallery>


==Задача о вискере, колеблющемся на игле==
Первая исследуемая конструкция представляет из себя одиночный прямой вискер, колеблющийся на игле. В ходе одного из экспериментов над такой системой было обнаружено, что существует две близкие собственные частоты этой системы с похожими формами колебаний.
[[File:Sameshapes.png|thumb|400px|left|Колебания одиночного вискера на игле]]

Уравнения каждой из таких балок: <math> \rho F\frac{\partial^2 w}{\partial t^2}=-E J\frac{\partial^4 w}{\partial x^4} </math>. Это уравнения колебаний балки Бернулли-Эйлера. 
Для построения приближённого решения задачи о двух балках используется метод, предложенный В.М.Фридманом, заключающемся в использовании спектральных свойств отдельных элементов системы. В соответствии с работой уравнения движения системы
балок могут быть для удобства записаны в операторном виде.
<math>\xi=\left( \begin{array}{cc}
Q\\
M
\end{array} \right) </math><math>\eta=\left( \begin{array}{cc} w\\ \theta \end{array} \right) </math>
Так же вводятся <math>D</math> и <math>D^*</math> - дифференциальные операторы вида, а <math>R</math> и <math>B</math> - алгебраические операторы инерции и упругости: <math>D=\left( \begin{array}{cc}
-\frac{\partial }{\partial x} & 0\\
1 & -\frac{\partial }{\partial x}
\end{array} \right) \qquad D^*=\left( \begin{array}{cc} \frac{\partial }{\partial x} & 1\\ 0 & \frac{\partial }{\partial x} \end{array} \right)</math>
,<math>R=\left( \begin{array}{cc}
\rho & 0\\
0 & j
\end{array} \right) \qquad B=\left( \begin{array}{cc} b & 0\\ 0 & \beta \end{array} \right)</math>
Важно заметить, что для <math>D</math> и <math>D^*</math> верно следующее соотношение, оно позволяет связать эти два дифференциальных оператора.
<math>\int\limits_0^l D \xi \eta \,dx +\xi\eta\bigg|_{x=0}^{x=l} = \int\limits_0^l \xi D^* \eta \,dx</math>
Тогда уравнения колебаний балок примут следующий вид:<math>\begin{array}{l}
D \xi - \lambda^2 R \eta = 0
\\
D^* \eta - B \xi = 0
\end{array}</math>

Так же необходимо учесть разложение решений на две составляющие:<math> \xi = \widetilde{\xi} + \bar{\xi} \qquad \eta = \widetilde{\eta} + \bar{\eta}</math>.
Здесь составляющая с волной обозначает разложение по всем собственным формам каждой из балок, а слагаемое с чертой отвечает за сопряжение вискера с иглой. Решение с волной это сумма собственных форм колебаний с некоторыми коэффициентами:<math> \widetilde{\xi} = \sum_{i=1}^\infty \alpha_i \widetilde{\xi}_i \widetilde{\eta} = \sum_{i=1}^\infty \beta_i \widetilde{\eta}_i  </math>

Предполагается,что объяснение этого эксперимента заключается в том, что нельзя рассматривать вискер, как консольную балку на неподвижном основании, необходимо учитывать колебания иглы, хотя она много превышает размерами вискер. Более того, можно предположить, что игла колеблется со своей высокой частотой, соответствующей одной из многоузловых форм колебаний. Поэтому во всем решении была учтена высокая частота балки и первая частота вискера.  То есть разложение по собственным формам колебаний проводилось до слагаемого, равного по номеру количеству учтенных собственных частот. 
[[File:Blabaa.png|thumbnail|500px|Форма колебаний иглы и вискера в системе]]
Результатом такого метода решения стал набор собственных чисел всей системы игла-вискер. Колебания системы при полученных частотах имеют вид, изображенный на рисунке.

==Дифференциальный резонатор==
Под дифференциальным резонатором из углеродных вискеров понимается система, состоящая из двух параллельно расположенных вискеров, соединённых между собой упругой перемычкой и закреплённых на жёстком основании. Свободные колебания системы с двумя степенями свободы и близкими собственными частотами имеют вид биений. Период огибающей сильно зависит даже от малых изменений упругих и инерционных свойств системы, что может позволить повысить точность измерения малых добавленных масс.
В качестве модели, для исследования такого резонатора были взяты две балки Бернулли-Эйлера, с заделкой на нижнем конце и свободные с верхнего. При этом упругая перемычка моделируется, как пружина, соединяющая два вискера. При частотах колебаний, близких к резонансным уравнения будут иметь следующий вид: <math>\begin{array}{c}
m \ddot{x_1} + b \dot{x_1} +c x_1 + c \varepsilon (x_1 - x_2) = F(t)
\\
m \ddot{x_2} + b \dot{x_2} +c x_2 + c \varepsilon (x_2 - x_1) = F(t)
\end{array}</math> В правой части стоит вынуждающая сила колебаний, которая выражает колебания подложки под вискерами. 
[[File:2whiskera.png|thumb|300px|left|Колебания одиночного вискера на игле]]
Для измерения масс в приведенных выше уравнениях необходимо взять разные массы вискеров. Одна из них считается равной елинице, а вторая отличается от неё на несколько процентов. Тогда при увеличении разницы между массами вискеров, частота огибающей будет возрастать.
<gallery mode=""traditional"">
File:2-2.png|Колебания одного из вискеров при разнице масс 3%
File:3-2.png|Колебания одного из вискеров при разнице масс 5%
File:6-2.png|Колебания одного из вискеров при разнице масс 15%
</gallery>
[[File:Graph11.png|thumb|300px|right|Зависимость частоты огибающей от приложенной массы]]
Зависимость частоты огибающей от массы, прикрепленной к одному из вискеров является практически линейной, если добавочная масса более 2%. То есть возможно производить градуировку нановесов на основе 2-3 измерений. Это значительно упрощает производство таких приборов.

== Емкостные датчики колебаний==
В работе рассмотрены три вида емкостных датчиков колебаний. 
* Первый самый простой датчик представляет из себя электрическую цепь, включающую в себя конденсатор с переменным сопротивлением. Одной из обкладок конденсатора является вискер. Колебания вискера вызывают колебания напряжения в цепи. Проблемой данного датчика является нелинейность выходного сигнала. 
<gallery mode=""traditional"">
File:Condensator.png|Устройство емкостного датчика колебаний одиночного вискера
File:2condens.png|Устройство емкостного датчика колебаний двойного вискера
File:Lineinii.png|Устройство улучшенного емкостного датчика колебаний одиночного вискера
</gallery>
Размерные уравнения колебаний такого вискера и соответствующего ему напряжения выглядят следующим образом:<math> \begin{array}{l}
U-R i-u_c=0
\\
m \ddot{x}+b \dot{x} + cx - \frac{C u_c ^2}{d}= F(t)
\end{array} </math>. Первое уравнение показывает распределение напряжений в цепи. То есть напряжение от источника частично тратится на сопротивление цепи, а частично на перезарядку конденсатора. Второе уравнение это уравнение механических колебаний вискера с демпфером, но с учетом действующей на него электрической силы. 
В результате численного решения этих уравнений получаются результаты, представленные на рисунке. Видно, что вискер имеет гармоническую форму колебаний, а напряжения нет. Но период и частота совпадают. 
*Второй датчик рассчитан на измерение колебаний двойного вискера и на детектирование масс. Он состоит из двух вискеров, связанных между собой упругой перемычкой, каждый из которых подключен к своему конденсатору.  Уравнение колебаний будут выглядеть следующим образом:  <math> \begin{array}{l}
m \ddot{x_1} + b \ddot{x_1} - \frac{C_1 u_1^2}{d}+\epsilon c (x_1 - x_2) = F(t)
\\
m \ddot{x_2} + b \ddot{x_2} - \frac{C_2 u_2^2}{d}+\epsilon c (x_2 - x_1) = F(t)
\\
U-R i_1-u_1c=0
\\
U-R i_2-u_2c=0
\end{array}</math>
Эти уравнения отличаются от предыдущих тем, что в механической части добавилось слагаемое, отвечающее за упругую связь между вискерами. В результате интегрирования видно, что график напряжения полностью повторяет график механических колебаний, главное, что возможно точно измерить период огибающей биений, то есть померить массу с помощью напряжения. 

*Последний датчик предназначен для точного измерения колебаний одиночного вискера. Его особенность заключается в том, что один и тот же вискер выступает в роли обкладки для двух конденсаторов. Тогда выходной сигнал каждого из напряжений будет нелинейным, как и в первом случае, однако, если в качестве результата рассматривать разность напряжений на конденсаторах, то выходной сигнал получается гармоническим, точно совпадающим по частоте с колебаниями вискера, это дает возможность измерять частоту колебаний одиночного вискера с большой точностью. Уравнения колебаний такого датчика выглядят следующим образом: <math>m \ddot{x} + b \ddot{x} - \frac{C_1 u_1^2}{d}- \frac{C_2 u_2^2}{d} = F(t)
\\
U-R i_1-u_1c=0
\\
U-R i_2-u_2c=0</math>

<gallery mode=""traditional"">
File:UstCondens.png|Колебания вискера и напряжения в простейшем датчике
File:2n.png|Колебания вискера и напряжения в датчике двойного вискера
File:Lineinii2.png|Колебания вискера и напряжений в двойном датчике
File:Lineinii3.png|Колебания вискера и разности напряжений в двойном датчике
</gallery>

==Выводы==
В работе было рассмотрен и объяснён эксперимент с одиночным вискером, колеблющемся
на игле; исследована и смоделирована система, состоящая из двух вискеров, соединённых
упругим элементом; рассмотрены различные емкостные датчики колебаний.

Решение задачи об одиночном вискере, колеблющемся на игле, объясняет природу явления, зафиксированного в эксперименте, а так же предлагает метод приближённого решения
подобных проблем. Метод решения (использования операторного подхода и разложения Галёркина) не увеличивает размерность задачи при увеличении количества искомых собственных частот, боле того, он учитывает особенности системы и позволяет находить только нужные собственные частоты и формы. Применить использованный подход можно так же и ко многим другим
задачам. 

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

Последним этапом проведённого исследования было изучение емкостных датчиков, с помощью которых можно производить измерения колебаний нановискеров. Изучено три вида датчиков. Все они основаны на использовании вискера в качестве одной из обкладок конденсатора, подключенного в электрическую цепь. Такие устройства позволяют измерять напряжение в сети, вместо измерения амплитуды или частоты колебаний самого вискера. 

Все проведённое исследование основано на экспериментах с нановискерами и направлено на их улучшение и объяснение. Понимание процессов, происходящих при натурном исследовании нанообъектов позволяет упростить и усовершенствовать эксперименты, а также
сократить количество неудачных, что ведет к значительной экономии денежных средств,
материалов и времени учёных-эксперементаторов.


==Список литературы==
*Левичев В.В. Электронные и фотонные устройства: принцип работы, технологии изготов-
ления. //Учебник, учебное пособие, Год: 2015]
* Я.С.Гринберг, Ю.А.Пашкин, Е.В.Ильичев. Наномеханические резонаторы.//Успехи фи-
зических наук. Т.182. №4. 2012.
* Левичев В.В., Жуков М.В., Мухин И.С., Денисюк А.И., Голубок А.О. Об устойчивости
работы сканирующего силового микроскопа с нановискером на вершине зонда.// Журнал
технической физики, 2013, том 83, вып. 7
* Голубок А.О., Ковров А.В., Левичев В.В., Мухин И.С., Приходько О.А. Формирова-
ние одиночных нановискеров на вершинах зондов сканирующего зондового микроскопа.
Научно-технический вестник Санкт-Петербургского государственного университета ин-
формационных технологий, механики и оптики. 2009. № 4(62). С. 82-87. [Тип: Статья,
Год: 2009]
* Жилин П. А. Прикладная механика. Теория тонких упругих стержней: Учеб. пособие.
СПб.: Изд-во Политехн. ун-та, 2007. 101 с.
* Фридман В. М. Теория упругих колебаний: уравнения и методы / В. М. Фридман ; Рос.
акад. наук, Ин-т проблем машиноведения. - Санкт-Петербург : Наука, 2014. - 253, [1] с.;
16 усл. печ. л. - Библиогр.: с. 248-250 (61 назв.). - 500 экз. - ISBN 978-5-02-038375-3 (в пер.)
*Галёркин Б. Г. Стержни и пластинки. Ряды в некоторых вопросах упругого равновесия
стержней и пластинок. // Вестник инженеров. — 1915. — Т. 1. — С. 897—908.
* Штукин Л.В., Беринский И.Е., Индейцев Д.А., Морозов Н.Ф., Скубов Д.Ю. Электроме-
ханические модели нанорезонаторов. //Физическая мезомеханика 19 1 2016г
* И.Е. Бринский, Д.И. Индейцев, Н.Ф. Морозов, Д.Ю. Скубов, Л.В. Штукин. Дифферен-
циальный резонатор как детектор массы. // Механика твердого тела. №2 2015.
* Д.И. Индейцев, О.С. Лобода, Н.Ф. Морозов, Д.Ю. Скубов, Л.В. Штукин. Автоколеба-
тельный режим нанорезонатора. //Физическая мезомеханика 19 5 2016г с. 23-28
*Н. Ф. Морозов, И.Е. Бринский, Д.И. Индейцев, О.В. Привалова, Д.Ю. Скубов, Л.В. Шту-
кин. Срыв колебаний графенового резонатора, как способ определения его спектральных
характеристик. ДОКЛАДЫ АКАДЕМИИ НАУК, 2014, том 456, № 5, с. 1–5
41
@@ Строка 58: / Строка 58: @@
 Так же необходимо учесть разложение решений на две составляющие:<math> \xi = \widetilde{\xi} + \bar{\xi} \qquad \eta = \widetilde{\eta} + \bar{\eta}</math>.
-Здесь составляющая с волной обозначает разложение по всем собственным формам каждой из балок, а слагаемое с чертой отвечает за сопряжение вискера с иглой. Решение с волной это сумма собственных форм колебаний с некоторыми коэффициентами:<math> \widetilde{\xi} = \sum_{i=1}^\infty \alpha_i \widetilde{\xi}_i \qquad \widetilde{\eta} = \sum_{i=1}^\infty \beta_i \widetilde{\eta}_i  </math>
+Здесь составляющая с волной обозначает разложение по всем собственным формам каждой из балок, а слагаемое с чертой отвечает за сопряжение вискера с иглой. Решение с волной это сумма собственных форм колебаний с некоторыми коэффициентами:<math> \widetilde{\xi} = \sum_{i=1}^\infty \alpha_i \widetilde{\xi}_i \widetilde{\eta} = \sum_{i=1}^\infty \beta_i \widetilde{\eta}_i  </math>
-А решение с чертой представимо в виде:<math>\bar{\eta}=b_1 \left( \begin{array}{cc}
-\\
-\end{array} \right) + b_2 \left( \begin{array}{cc}
-x\\
--1
-\end{array} \right) \qquad \bar{\xi}=a_1 \left( \begin{array}{cc}
-\\
-\end{array} \right) + a_2 \left( \begin{array}{cc}
-x\\
-\end{array} \right)</math>
-В ходе дальнейшего решения, находится зависимость коэффициентов разложения по собственным формам от коэффициентов в решении с чертой. Далее с помощью уравнений колебаний балок в операторном виде составляется система уравнений относительно коэффициентов части решения с чертой. В этой зависимости будет присутствовать <math> \lambda</math> - собственная частота всей системы. В результате получается уравнение относительно ее. Весь ход решения здесь не приведен, в силу его громоздкости, но его можно посмотреть в тексте дипломной работы.
 Предполагается,что объяснение этого эксперимента заключается в том, что нельзя рассматривать вискер, как консольную балку на неподвижном основании, необходимо учитывать колебания иглы, хотя она много превышает размерами вискер. Более того, можно предположить, что игла колеблется со своей высокой частотой, соответствующей одной из многоузловых форм колебаний. Поэтому во всем решении была учтена высокая частота балки и первая частота вискера.  То есть разложение по собственным формам колебаний проводилось до слагаемого, равного по номеру количеству учтенных собственных частот.