Приборы и техника эксперимента, 2019, № 3, стр. 150-156

ХОЛЛОВСКИЙ МИКРОСКОП ДЛЯ ИССЛЕДОВАНИЯ ВЫСОКОТЕМПЕРАТУРНЫХ СВЕРХПРОВОДНИКОВ

Х. Р. Ростами *

Фрязинский филиал Института радиотехники и электроники им. В.А. Котельникова РАН
141190 Фрязино, Московской обл., пл. Введенского, 1, Россия

* E-mail: rostami@ms.ire.rssi.ru

Поступила в редакцию 09.04.2018
После доработки 26.11.2018
Принята к публикации 27.12.2018

Полный текст (PDF)

Аннотация

На основе преобразователей Холла разработан трехмерный растровый микроскоп с магнитной чувствительностью ~2.5 ⋅ 10–3 Гс/Гц1/2. Динамический диапазон микроскопа по магнитному полю составляет от 10–3 до ±3 ⋅ 103 Гс. Стабильность величины заданного магнитного поля лучше чем 10–5 Гс. Минимальная величина магнитного поля при его ступенчатом задании составляет 10–3 Гс. Максимальный размер площади обзора исследуемых объектов 5 × 5 мм2. Минимальные размеры шага сканирования по осям Х, Y составляют 10 и 1 мкм соответственно при грубом и плавном перемещении столиков. Максимальный размер перемещения по оси Z равен 25 мм с минимальным шагом сканирования 10 и 1 мкм соответственно при грубом и плавном перемещении. Микроскоп позволяет создавать аппаратную функцию произвольной формы за счет комбинации заданных величин: температуры, транспортного тока, постоянного магнитного поля и осциллирующего, затухающего во времени, переменного магнитного поля. Приведены примеры применения микроскопа для исследования высокотемпературных сверхпроводников BSCCO и YBCO.

В диагностических целях широко применяются холловские микроскопы [14], которые обеспечивают непосредственное измерение индукции магнитного поля, обладают широкими функциональными возможностями, позволяют визуализировать магнитное микросостояние протяженных объектов, просты в изготовлении и имеют низкую стоимость. Однако у них сравнительно низкие пространственное разрешение и чувствительность к регистрируемому сигналу, что не позволяет отчетливо выделить спектры сигнала отклика объектов на слабое магнитное поле.

Как показано в [5, 6], при исключении влияния на объект спада импульса магнитного поля и использовании в качестве задающего поля только фронта переменного затухающего магнитного поля H(t) удается существенно повысить пространственное разрешение микроскопа при исследовании динамики магнитного потока в высокотемпературных сверхпроводниках (ВТСП). При приложении к образцу затухающего магнитного поля, амплитуда которого с высокой стабильностью ступенчато изменяется, были обнаружены и селективно выделены незначительные пространственные выбросы локальных характеристик образца. Таким образом, в отличие от традиционных методик, возможно одновременное исследование кристаллических и магнитных микросостояний объектов [5].

Разработанные к настоящему времени микроскопы измеряют двумерную картину распределения магнитного поля на поверхности образца. Однако при исследовании пространственных распределений неоднородных полей, полей размагничивания и т.д. требуется получить трехмерную картину распределения магнитного поля вблизи поверхности и вокруг образцов [57]. Исходя из вышеизложенного, разработка трехмерного растрового микроскопа с повышенным пространственным разрешением, позволяющего локальным способом проводить исследование объектов в однородных и локальных осциллирующих затухающих магнитных полях, представляется актуальной задачей.

Для комплексного исследования пространственного распределения магнитного поля вблизи поверхности образца в микроскопе предусмотрено два режима сканирования.

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

Во втором режиме фиксируется поле и производится сканирование по координатам. Разрешение в данном случае определяется размерами измерительного зонда и точностью перемещения координатных столиков.

Образец, закрепленный на подвижном координатном столике, в ходе измерений можно перемещать в трех взаимно перпендикулярных направлениях по отношению к датчику. При этом при каждом значении Bz для пространственных изменений магнитного поля объекта по осям X, Y, превышающих $S_{{eff}}^{{1/2}}$ (Seff – эффективная площадь рабочей поверхности датчика), может быть установлена истинная картина распределения магнитного поля Вz(X, Y, Z) над образцом.

Схема конструкции микроскопа представлена на рис. 1. Перед измерениями держатель 1 (см. рис. 1) столика (2, 3) при помощи дифференциального микровинта поднимается вверх по оси Z, и через окно в несущей 4 в столик ввинчивается медная оправа 5 с образцом 6. Затем столик с образцом подводится на минимальное расстояние ∼1 мкм к медному диску 7 с преобразователями Холла (п.Х.) 8, 9. Фланец 10 жестко связан с штангой 11. Внутри штанги проходят поворотные тяги 12, обеспечивающие посредством червячных передач 13, 14 соответствующие линейные перемещения координатного столика в направлениях X, Y. Фиксаторы 15 служат для фиксации начала и конца работы червячных передач. Вся система перемещения находится внутри проточного микрокриостата 16.

Рис. 1.

Структурная схема микроскопа. 1 – держатель; 2, 3 – (X, Y)-столик; 4 – несущая; 5 – медная оправа; 6 – образец; 7 – медный диск; 8, 9 – преобразователи Холла; 10 – фланец; 11 – штанга; 12 – поворотные тяги; 13, 14 – червячные передачи; 15 – фиксаторы; 16 – проточный микрокриостат; 17 – микрокапилляр; 18 – тонкая капиллярная сетка; 19 – термометр; 20 – соосные медные соленоиды; 21 – гелиевый криостат КГ-100; 22 – экспресс-нагреватель; 23 – биполярный источник тока; 24, 25, 26 и 27 – токовые ключи.

Жидкий гелий или азот, поступающий через спиралеобразный микрокапилляр 17 (бифилярно намотанный поверх нагревателя), проходит через тонкую капиллярную сетку 18 и подается в микрокриостат. Термометр 19, приклеенный клеем БФ2 к образцу, совместно с нагревателем 17 при помощи модернизированной схемы на основе терморегулятора РТП-3М обеспечивают температуру в камере с точностью лучше 0.01% в диапазоне 2.5–300 К. Токовые и холловские контакты п.Х., разработанного на основе датчика ПХЭ 607118А1, утоплены в глубь подложки, а провода подведены с ее боковых сторон. Измерительные провода проходят между соосных тонких нержавеющих трубок центральной штанги 11. Проточный микрокриостат размещен внутри двух соосных резистивных (с медными обмотками) соленоидов 20. Вся система находится внутри гелиевого криостата ГК-100 (21) и окружена цилиндрическим магнитным экраном, изготовленным из алюминия и пермаллоя. Нагреватель 22 служит для быстрого изменения температуры.

Сигналы положения образца по координатам X, Y вырабатываются многооборотными потенциометрами, связанными с поворотными тягами 12 и червячными передачами (в схеме не показаны), синхронно работающими с червячными передачами 13, 14. При заполнении криостата жидким азотом большая секция катушек 20 создает поле в диапазоне от 10–3 до ±1.3 ⋅ 103 Гс. Вторая секция катушек 20 служит для компенсации постоянной составляющей поля Земли по координате Z.

На рис. 2 приведена структурная схема системы позиционирования микроскопа.

Рис. 2.

Структурная схема системы позиционирования микроскопа. 1, 2X-, Y-платформа; 3, 4 – ходовые гайки; 5, 6 – червячные колеса; 7, 8 – червяки; 9, 10 – пластины, регулирующие степень соприкосновения рабочих поверхностей платформы и ходовых гаек; 11 – нагреватель; 12 – оправа; 13 – образец.

Для обеспечения однородной температуры образца и уменьшения влияния на результаты измерения остаточных магнитных полей (X, Y)-платформы 1, 2 и ходовые гайки 3, 4 изготовлены из меди, а длинные винты для перемещения ходовых гаек – из латуни. Для уменьшения трения, меняющегося из-за изменения коэффициентов теплового расширения медных координатных столиков, и снижения нагрузок на винты вращения выбран угол наклона 60° между рабочими поверхностями платформ и ходовых гаек.

При настройке механической части установки минимальное расстояние между строго параллельными, тщательно отполированными поверхностями оправы образца и медного диска, в котором установлены п.Х., составляло 1 ± 0.1 мкм. Такое расстояние задавалось следующим образом. Сначала с помощью подложки из NdGaO3 толщиной 1 мм, установленной в оправе, внутри которой имелась выемка глубиной 1 мм, при Z‑перемещении обеспечивалось касание винтиком оправы. Винтик, создающий омический контакт с оправой, ввинчивался через изолятор в медном диске. Затем, после фиксации разрыва контакта, что соответствовало минимальному зазору между медным диском и оправой, с помощью дифференциального микровинта Z-перемещения оправа отодвигалась на 2 мкм. После этого ввинчиванием винтика устанавливался омический контакт с оправой. Далее оправа отодвигалась на 10 мм, и установленная подложка заменялась подложкой из NdGaO3 толщиной 1 мм, на поверхность которой лазерным распылением была нанесена эпитаксиальная пленка YBa2Cu3O7 – x (YBCO) толщиной ~1 мкм. Далее осуществлялось Z-перемещение оправы с образцом к винтику до фиксации омического контакта. Затем винтик отвинчивался до разрыва контакта, после чего проводилось измерение.

Размеры и толщина зеркально отполированной с обеих сторон подложки задавались заранее, исходя из технологических требований напыления пленок. Медная оправа с углублением строго по толщине подложки также позволяла выполнить замену подложки на пленку с боковой стороны, не меняя заранее установленный зазор между образцом и п.Х. Кроме того, толщина подложки и образца дополнительно измерялись с помощью микроскопа. Для этого сначала с помощью двух омических контактов по оси Z (сверху и снизу) устанавливалась толщина подложки (один из винтиков фиксировался, а другой регулировался). Точность зазора 2 мкм контролировалась также по зависимости ${{B}_{{tr}}}(Z){\text{/}}B_{{tr}}^{{\max }}(0)$ с помощью экстраполяции кривой к оси ординат, смещенной по оси Z на –2 мкм.

Такую процедуру настройки необходимо проводить перед измерениями и для образцов другой толщины. Это позволяет компенсировать люфты резьбы микровинта при прямом и обратном перемещении образца по оси Z.

Чтобы уменьшить усилие вращающих механизмов, способное привести к деформации винтиков при детектировании омических контактов, сопротивление контактов регистрировалось с помощью высокочувствительного преобразователя сопротивление–напряжение. С помощью омических контактов таким же образом задавались и контролировались начало и конец перемещения (X, Y)-столиков.

Для контроля степени параллельности плоскостей, поверхности образца и (X, Y)-столиков при сканировании были использованы радиальные распределения плотности захваченного магнитного потока эпитаксиальных пленок YBCO и намагниченности однородно намагниченной магнитной пленки. При изменении угла наклона между поверхностями образца и п.Х. изменялось показание п.Х. Для контроля неизменности угла наклона при перемещении ходовых гаек, а также отсутствия гистерезиса были построены зависимости Bz(X) и Bz(Y) в прямом и в обратном направлениях сканирования. Таким образом, по показаниям п.Х. определялась степень параллельности сканирования образца по отношению к поверхности п.Х.

Для исключения гистерезиса, который мог появиться при перемещении по оси Z, были построены нормированные осевые распределения ${{B}_{{tr}}}(Z){\text{/}}B_{{tr}}^{{\max }}(0)$ в прямом и обратном направлениях сканирования.

После установки образца координатный столик выводился в положение, соответствующее краю линии, проходящей примерно посередине длинной стороны плоскости образца. Устанавливалось соответствие между механическими перемещениями образца в направлениях X, Y, Z и соответствующими показаниями многооборотных потенциометров позиционирования. Затем задавалась требуемая температура в измерительной камере и проводилась настройка электронной схемы магнитометра. Быстрое и точное перемещение столиков в трех направлениях обеспечивалось дифференциальными микровинтами с точностью 10 и 1 мкм соответственно в грубом и плавном диапазонах. В грубом интервале после выбора размера поля обзора интересующего участка исследуемого образца и задания значения координат Y, Z образец перемещался с шагом 1 мкм относительно п.Х. вдоль оси Х. Для каждого значения координаты Х измерялась величина поля, и данные автоматически считывались компьютером. После окончания сканирования по оси Х вводилось новое значение координаты Y, и для него измерения повторялись и т.д. В результате была получена двумерная картина распределения поля над образцом [8]. Далее образец перемещали по оси Z вверх на заданный шаг и заново снимали распределение поля над образцом. На основе полученных данных строили трехмерное изображение Bz(X, Y, Z). Экспресс-выбор соответствующей области поверхности образца (для сканирования) также можно осуществлять с помощью винтов, на которых установлены омические контакты.

Для измерения магнитного момента образцов были использованы два близких по техническим характеристикам п.Х. марки ПХЭ 607118А1 (НПО “Сенсор”, Санкт-Петербург). Оба п.Х. находились в одной плоскости, при этом один из них был размещен в центре медного диска под образцом, а другой – опорный п.Х. – на расстоянии 9 мм от центра, в невозмущенном образцом магнитном поле. Чувствительности п.Х. выравнивались путем регулировки тока опорного п.Х. Каждый раз перед измерениями после задания необходимой температуры образец сначала перемещали по оси Z на расстояние 25 мм, после чего с помощью электронной схемы компенсировали сигнал разбаланса п.Х., который имел разные значения для установленных температур. Температурный разбаланс преобразователей Холла особенно сильно сказывался при регистрации слабых магнитных полей.

Для изменения параметров неоднородного поля и их измерения три одинаковых по техническим характеристикам п.Х. размещали на трех строго взаимно перпендикулярных гранях кубика размером 5 × 5 × 5 мм, который закрепляли на трехкоординатном столике микроскопа. В результате измеряли три компоненты Bx, By, Bz и величину поля $B = {{(B_{x}^{2} + B_{y}^{2} + B_{z}^{2})}^{{1/2}}}$ произвольного направления. Этим же способом были измерены все тензоры магнитного поля – основные элементы матрицы градиентов магнитного поля:

${{B}_{{x,y,z}}} = \left| {\begin{array}{*{20}{c}} {\frac{{\partial {{B}_{x}}}}{{\partial x}}}&{\frac{{\partial {{B}_{x}}}}{{\partial y}}}&{\frac{{\partial {{B}_{x}}}}{{\partial z}}} \\ {\frac{{\partial {{B}_{y}}}}{{\partial x}}}&{\frac{{\partial {{B}_{y}}}}{{\partial y}}}&{\frac{{\partial {{B}_{y}}}}{{\partial z}}} \\ {\frac{{\partial {{B}_{z}}}}{{\partial x}}}&{\frac{{\partial {{B}_{z}}}}{{\partial y}}}&{\frac{{\partial {{B}_{z}}}}{{\partial z}}} \end{array}} \right|,$
а также параметры созданного поля в отверстии соленоида.

Вышеприведенные принципиальные решения позволили существенно расширить технические возможности предлагаемого микроскопа по сравнению с ранее разработанным нами микроскопом [9]. Принцип действия входящего в состав микроскопа высокочувствительного холловского магнитометра подробно описан в [10, 11].

Предложенный микроскоп позволяет проводить исследования во внешних магнитных полях величиной от 10–3 до ±3 ⋅ 103 Гс в диапазоне температур 2.5–300 К. Нестабильность и нелинейность во всем динамическом диапазоне развертки поля не хуже 0.01%. Температуру также можно сканировать по заданному закону во всем указанном диапазоне с нестабильностью не хуже 0.01% [12, 13].

Испытания показали, что при прецизионном изготовлении деталей микроскопа разрешение системы позиционирования по обеим координатам, благодаря использованию фиксаторов 15 (см. рис. 1) начала и конца перемещения столиков и применению принципа дифференциального перемещения, было лучше 1 мкм. Полный оборот тяги обеспечивал линейное перемещение в соответствующем направлении <1 мкм. Полное линейное перемещение по осям X и Y составляло 5 мм. По оси Z система обеспечивала линейное перемещение на расстояние 25 мм с разрешением не хуже 1 мкм. Это на порядок лучше ранее достигнутого нами в [9] разрешения. В микроскопе, описанном в [9], перемещения по осям Z, X, Y обеспечивались с помощью головок, имеющих минимальный шаг 10 мкм.

Сравнительно высокое разрешение по (X, Y)-направлениям предусмотрено и в случае использования в микроскопе п.Х. с рабочими поверхностями S ∼ 10 × 10 мкм2 [3] и S ∼ 1 × 1 мкм2 [1] (у использованных нами в установке п.Х. площадь рабочей поверхности S ∼ 50 × 100 мкм2). Исследования плотности захваченного магнитного потока в ВТСП-образцах также показали, что разрешение разработанного холловского микроскопа по индукции магнитного поля составило ≈2.5 ⋅ 10–3 Гс/Гц1/2, динамический диапазон микроскопа по магнитному полю – от 10–3 до ±3 ⋅ 103 Гс, стабильность величины заданного поля не хуже 10–5 Гс. Более широкий динамический диапазон в области слабых полей по сравнению с ранее предложенным в [9] микроскопом необходим для исследования параметров гипервихрей в ВТСП. Такие слабые поля создавались благодаря питанию соленоида высокостабильным биполярным генератором разрядных токов, управляемым двоичным кодом [12, 13]. Для оценки чувствительности микроскопа сначала на сравнительно большом токе с помощью откалиброванного п.Х. устанавливалась постоянная (ток–поле) катушка, а затем при известном токе определялась величина поля и вычислялась магнитная чувствительность микроскопа.

В случае неутопления контактов в глубь подложки п.Х. или при использовании п.Х. промышленного производства из-за токовых и потенциальных контактов, имеющихся на его поверхности, минимальный зазор между рабочей поверхностью п.Х. и поверхностью образца составит около 200 мкм. При этом пространственное разрешение микроскопа не определяется размерами рабочей поверхности п.Х. (в нашем случае 50×100 мкм2) и составит ~200 мкм. В этом случае для установления точного значения величины максимального магнитного поля $B_{{tr}}^{{\max }}(0)$ в центре на поверхности образцов необходимо построить нормированные осевые распределения ${{B}_{{tr}}}(Z){\text{/}}B_{{tr}}^{{\max }}(0)$ от координаты Z. Экстраполяция кривых к оси ординат, смещенной по оси Z на –200 мкм, показала, что потери сигнала на расстоянии 200 мкм от поверхности образцов составляют <5% от $B_{{tr}}^{{\max }}(0)$.

Для создания осциллирующего, затухающего во времени, переменного магнитного поля параллельно медному соленоиду (индуктивностью L ≈ ≈ 4.13 Гн с остаточным активным сопротивлением обмотки R ≈ 92.2 Ом при Т = 77.4 К), создающему внешнее поле, подключалась емкость C = = 0.05 мкФ (см. рис. 1). Параметры соленоида выбирались, исходя из конструкции криостата и требований к созданию максимально возможного однородного магнитного поля в отверстии соленоида. Задающим являлось только значение емкости, которое выбиралось, исходя из обеспечения наиболее высокой разрешающей способности методики, что достигалось воздействием на образец максимального числа колебаний затухающего внешнего магнитного поля.

При одновременном поступлении от компьютера на токовые ключи 24, 25 и 26 (см. рис. 1) импульсов прямой и обратной полярности происходит подсоединение биполярного источника тока 23 к соленоиду L и отключение емкости С и резистора R1. В этом случае магнитная энергия накапливается в соленоиде. При смене полярностей импульсов, поступающих от компьютера в токовые ключи 24, 25 и 26, источник тока 23 отключается, подсоединяется резистор R1 и включается в цепь соленоида L емкость С. В результате в соленоиде возникает переменное затухающее магнитное поле H(t) = H0exp(βt)cosΩt, где β = R/(2L) = 11 c–1 – коэффициент затухания, Ω = (ω2 – β2)1/2 ≈ 352 Гц – собственная частота колебательного контура, ω2 = 1/(LC). LC-контур устойчиво работает в частотном диапазоне до 5 кГц. Увеличение емкости C до 100 мкФ приводит к исчезновению колебаний и к переходу системы в режим работы, аналогичный традиционным методикам, когда параллельно соленоиду вместо емкости подключен резистор. Для создания однородного магнитного поля в токовые ключи 25 и 27 от компьютера подаются разнополярные импульсы, которые отключают емкость С и включают в цепь резистор R2.

На рис. 3 приведено пространственное распределение захваченного магнитного потока Btr(X, Y, Z) вблизи монокристаллического образца Bi2Sr2CaCu2O8 +x (BSCCO), измеренное при температуре 77.4 К с шагом 50 мкм по координатам X, Y для расстояния между образцом и п.Х. в диапазоне Z = 50 мкм–1 мм. Рис. 4 демонстрирует пространственные возможности микроскопа. Радиальное распределение Btr(r) для квазимонокристаллического образца YBCO измерено в режиме ZFC при температуре 77.4 К. Согласно рис. 4а, в образце отчетливо видны границы доменов размером ~3–4 мм, при этом разрешение микроскопа при радиальном сканировании было не хуже 50 мкм.

Рис. 3.

Пространственное распределение плотности захваченного магнитного потока Btr(X, Y, Z) вокруг монокристаллического образца BSCCO, полученное при сканировании с шагом 50 мкм по координатам XY для расстояния между образцом и п. Х. в диапазоне Z = 50 мкм–1 мм. Температура 77.4 К.

Рис. 4.

Радиальные распределения Btr(r) для квазимонокристаллического образца YBCO [5] (а) и эпитаксиальной пленки YBCO (б). Температура 77.4 К.

На рис. 4б приведено радиальное распределение Btr(r) для эпитаксиальной пленки, измеренное в режиме ZFC при температуре 77.4 К. Эта зависимость была использована для юстировки (XY)-столиков. Для демонстрации возможностей микроскопа в режиме затухания поля и сравнения полученных результатов с измеренными с помощью традиционной методики на рис. 5 приведены магнитополевые зависимости плотности захваченного магнитного потока Btr(H0). На рис. 5а приведена зависимость для эпитаксиальной пленки YBCO, измеренная с помощью локального затухающего поля, а на рис. 5б для той же эпитаксиальной пленки YBCO – магнитополевая зависимость, измеренная в обычном режиме для случая, когда параллельно соленоиду подключался резистор [5].

Рис. 5.

Магнитополевая зависимость плотности захваченного магнитного потока Btr(H0) [5]: a – для эпитаксиальной пленки YBCO, измеренная с помощью затухающего магнитного поля; б – для той же эпитаксиальной пленки YBCO, измеренная в обычном режиме, когда параллельно соленоиду подключен резистор. 1 – при росте внешнего магнитного поля, 2 – при его убывании. Температура 77.4 К.

Как показано в [5], в случае применения осциллирующего, затухающего во времени магнитного поля эффективная площадь Seff рабочей поверхности использованного в микроскопе п.Х. становится существенно меньше, чем заданный размер его рабочей поверхности S. Это дополнительно повышает пространственное разрешение микроскопа по сравнению со случаем применения однородного поля. Пространственное разрешение и чувствительность микроскопа можно улучшить также при применении зондов другого типа, например датчиков на основе гетеропереходов с двумерным газом в канале [14].

Список литературы

  1. Oral A., Bending S.J., Henini M. // J. Vac. Sci. Technol. B. 1996. V. 14. P. 1202.

  2. Gregory J.K., Bending S.J., Sandhu A. // Rev. Sci. Instrum. 2002. V. 73. P. 3515.

  3. Cambel V., Fedor J., Gregusova D., Kováč, P., Hušek I. // Supercond. Sci. Technol. 2005. V. 18. P. 417. http:// doi.org/10.1088/0953-2048/18/4/007.

  4. Игнатьев В.К., Орлов А.А., Перченко С.В., Станкевич Д.А. // Письма в ЖТФ. 2017. Т. 43. № 15. С. 3. https://doi.org/doi 10.21883/PJTF.2017.15.44864.16434

  5. Ростами Х.Р. // ФТТ. 2013. Т. 55. № 9. С. 1677.

  6. Ростами Х.Р. // Письма в ЖЭТФ. 2017. Т. 105. № 12. С. 754.

  7. Антонов В.Г., Петров Л.М., Щелкин А.П. // Средства измерений магнитных параметров материалов. Л.: Энергоатомиздат, 1986.

  8. Гудошников С.А., Венгрус И.И., Снигирев О.В. // Сверхпроводимость: физика, химия, техника. 1994. Т. 7. С. 746.

  9. Ростами Х.Р. // Измер. техника. 2014. № 12. С. 51.

  10. Ростами Х.Р. // Измер. техника. 2016. № 12. С. 40.

  11. Ростами Х.Р. // ПТЭ. 2016. № 2. С. 112. https://doi.org/10.7868/ S0032816216010110

  12. Ростами Х.Р. Патент № 2007862 C1 РФ. МПК: H03M1/66 //Бюл. № 3. Опубл. 15.02.94.

  13. Ростами Х.Р. Патент № 94038958 РФ. МПК: H03K005/00 A1 // Бюл. № 20. Опубл. 20.07.96.

  14. Bending S.J. // Physica C. 2010. T. 470. № 19. P. 754.

Дополнительные материалы отсутствуют.