После завершения реконструкции ускорителя ЛЯП ОИЯИ и наладки тракта пучков частиц, начиная с 1987 года на muSR установках комплекса был выполнен большой цикл работ, в основном связанный с проблемой высокотемпературной сверхпроводимости (ВТСП), изучением магнитных фазовых переходов в чистых редкоземельных элементах и поведением мюона и мюония в криокристаллах и криожидкостях.
Энергия протонов (659 +/- 6) МэВ Интенсивность выведенного пучка до 2,0 мкА Длительность макроимпульсов ~ 3,2 мс Частота повторения макроиммпульсов 250 Гц Длительность микроимпульсов ~ 10 нс Частота повторения микроиммпульсов ~ 14,6 МГц
Комплекс включает набор криостатов, рассчитанных на температуры 4,2-300 К, магнитные системы перпендикулярного (к начальной поляризации пучка поля на 5 КГс и продольного - 7 КГс, автоматизированный сцинтилляционный временной спектрометр (ширина канала 1 нс ), систему обеспечения хладоагентами и сбора гелия.
На рисунке 2.1 показано размещение спектрометра в I корпусе Лаборатории ядерных проблем ОИЯИ. Протонный пучок 1 фазотрона ОИЯИ падает на узел пионообразующих мишеней 2. Часть пионов, вылетающих из мишени, после отклонения магнитом (на рисунке не показан) на угол 420 направляются в мюонный канал 3, состоящий из фокусирующих квадрупольных магнитных линз. Проходя мюонный канал и распадаясь на лету, пионы создают на выходе канала в месте расположения muSR-установок пучок мюонов, который поворачивается магнитами 4 и 5, фокусируется квадрупольной линзой 6 на образец, расположенный в центре магнита 7. На рисунке показан магнит с поперечным направлением магнитного поля. Образцы охлаждаются жидким гелием (азотом) из дьюра 8. Вместо магнита 7 может быть установлен магнит с продольным направлением поля 9. Оставшиеся после пионообразующей мишени протоны пучка попадают в ловушку 10.
В данных экспериментах используется пучок мюонов с характеристиками, приведенными в Таблице 2.2:
Импульс | 125 МэВ/c (Emu = 58 МэВ) |
Импульсный разброс Dp/p | 10% |
Поляризация | 80% |
Сечение пучков на выходе | 10 x 6 см2 |
Интенсивность на 1 мкА протонного пучка |
3x105 mu+/с (примесь 3% e+ 105 mu-/с (6% e-; 2,4% pi-) |
Схема спектрометра с магнитом, создающим поперечное спину мюонов направление поля, показана на рисунке 2.2.
Регистрация событий mu -> e распада осуществляется с помощью сцинтилляционных счетчиков 1 - 5. Мюоны через отверстие в панцире магнита 6 проходят через замедлитель, где теряют часть своей энергии. Толщина замедлителя подбирается исходя из плотности образца, для получения максимального числа остановок мюонов в образце. Типичный вид кривой остановок приведен на рисунке 2.3. На этом же рисунке показаны зависимости количества позитронов, вылетающих "вперед" (4*5-2-3) и "назад" (2*3-1-4) по отношению к направлению импульса мюонов. Между счетчиками 1 и 2 установлен свинцовый коллиматор 7, диаметр отверстия которого может меняться в зависимости от размера образца с помощью вставных колец 8. Между счетчиками 2 и 5 устанавливается криостат 9 с образцом 10.
Сцинтилляторы счетчиков 3 и 4 находятся внутри холодной зоны криостата и свет от них выводится на ФЭУ по световодам 11 и 12, расположенным друг над другом. С помощью соединений 13 и 14 криостат пристыковывается к дюарам с жидким гелием (азотом) и к системе сбора гелия. Все элементы спектрометра находятся внутри панциря магнита 15 на платформах 16, перемещающихся по направляющим 17.
Регистрация остановок мюонов в образце и событий mu - > e распада осуществляется с помощью сцинтилляционных счетчиков на основе полистирола с добавкой 2% Р-терфинила и 0.02% РОРОР[32] и фотоэлектронных умножителей XP-2020, XP-2982 фирмы Philips[33]. Сцинтилляторы и световоды из оргстекла типа ПММА были соединены клеем на основе эпоксидной смолы ЭД-20. Для улучшения условий светосбора сцинтилляторы и световоды были обернуты белой бумагой, а поверх нее для светоизоляции - черной. Счетчики 3 и 4 расположены в холодной зоне криостата и описаны в п.2.3.
Характерной особенностью работы сцинтилляционных счетчиков в muSR-экспериментах является наличие сильных магнитных полей. Как известно, магнитное поле уменьшает амплитуду сигнала с ФЭУ, что приводит к изменению эффективности регистрации частиц. Это обстоятельство особенно важно учитывать при регистрации позитронов распада, поскольку экспериментальный коэффициент асимметрии углового распределения позитронов зависит от энергетического порога регистрирующей системы. Для уменьшения влияния магнитного поля на работу ФЭУ были использованы длинные световоды. Для защиты от рассеянных полей, ФЭУ помещались в кожухи из магнитомягкого материала (пермалоя).
Фотоэлектронные умножители запитывались от блоков ВВИП-2, представляющих собой модули КАМАК двойной ширины с диапазоном регулирования выходного напряжения 0 - 2560 В с дискретностью в 10 разрядов[34]. Блоки ВВИП-2 были установленны в отдельном крейте, т.к. их суммарная потребляемая мощность близка к предельно допустимой для блока питания крейта. Для управления этим крейтом использовался последовательный контроллер КК011 [35].
В кожухах ФЭУ кроме делителей, собранных по стандартной схеме [33], были размещены формирователи, аналогичные описанным в [36]. Пороги регистрации импульсов задавались от специально разработанного для этой цели 6-ти канального цифро-аналогового преобразователя (ЦАП). Он представляет собой модуль КАМАК шириной 1М, состоящий из восьми 6-разрядных ЦАП. Напряжение на выходах ЦАП (-3 - 0 В) задается подаваемым на них кодом. Управление блоком производится по магистрали КАМАК следующими командами:
NA(0 - 7)F(0) - чтение в ЭВМ данных о величине установленного порога по шинам W1-W8; NA(0 - 7)F(9) - сброс установленного порога во всех каналах и запись в них кода 0; NA(0 - 7)F(16) - запись в канал с номером 0 - 7 кода по шинам W1-W6; NA(0 - 7)F(25) - добавление единицы к коду, записанному во все каналы.
Сигналы с формирователей поступали на систему отбора событий, структурная схема которой приведена на рис.2.5.
Рис.2.5 Схема отбора событий.
Основу системы отбора событий составляют 5 блоков временной селекции БВС1 - БВС5 [37]. Эти блоки служат для отбора событий в режимах 2-, 3-, 4- кратных совпадений с использованием одного или двух каналов запрета. Сигналы на них поступают после формирователей-разветвителей 1 - 5.
Блок БВС1 фиксирует все частицы, проходящие через счетчики 1+2. В этом потоке присутствуют как все мюоны (останавливающиеся в образце и проходящие все счетчики - т.н. "пролетные"), так и пионы и позитроны. Блок БВС2 выделяет "пролетные" частицы по критерию (2+3)*(4+5) (знак "+" обозначает логическое ИЛИ).
Сигналы с этих блоков,собранные по "ИЛИ", запускают одновибратор продлевающего типа ОВП, который вырабатывает сигнал "Охрана до" - B ("Busy"), предотвращающий запись в спектры события mu->e распада, если ему предшествовала регистрация мюона, распад которого не был зафиксирован, или "пролетной" частицы. Одновибраторы продлевающего типа серийно не выпускаются и этот блок был разработан и изготовлен в механическом и электрическом конструктиве КАМАК на основе микросхемы К155АГ3. Входные и выходные сигналы блока имеют уровни NIM.
Временной селектор БВС3 выделяет мюоны, остановившиеся в мишени с учетом сигнала "Охрана до" по схеме (1+2+3)-B. Сигналы с этого блока поступают на одновибратор ОВС (одновибратор стартовых событий), который вырабатывает ворота G ("Gate") длительностью 10 мкс. Эти ворота запускают время-цифровой преобразователь (ВЦП) типа КА304[38]. Еще один ВЦП используется совместно аналоговым преобразователем временных интервалов (АПИ)[39], позволяющим работать с меньшей шириной канала. Блок БВС4 (4*5) регистрирует позитроны, вылетающие из мишени "вперед" ef , а БВС5 (2*3-1) - позитроны, вылетающие "назад" eb . Сигналы с этих блоков смешиваются и подаются на вход "Стоп" ВЦП. Кроме этого сигнал eb поступает на вход "Признак" ВЦП, так что спектры от позитронов распада "вперед" записывается в каналы с 0 по 2047, а спектры от позитронов "назад" - в каналы 2048 - 4095. Если во время стартовых ворот G приходит более одного мюона или регистрируется более одного позитрона, то формируется сигнал "Сброс", прекращающий работу ВЦП. С выхода ВЦП 12-разрядные коды передаются через блок управления инкрементной записью КЛ018 в оперативную память КЛ020[40], где и накапливаются измеряемые спектры. Синхронизация работы блоков временной селекции осуществляется с помощью 14 линий задержек, не показанных на рисунке.
При работе с многоканальными системами сцинтилляционных счетчиков (годоскопами) значительная часть времени уходит на измерение характеристик детекторов и настройку спектрометра. Благодаря применению многоканального ЦАП установки порогов и многофункциональных блоков БСВ-4A эти операции были автоматизированы. Управление системой регистрации событий осуществлялось с помощью ПЭВМ IBM PC/AT 386, сопряженной с крейтами КАМАК через контроллеры КК 012 (KK 009)[35,41]. Процедура набора спектров контролировалась на экранах мониторов, на которые выводились набираемые спектры в графическом виде или результаты их Фурье-преобразования. Также выводилась информация об интенсивности мюонного пучка, количестве остановок мюонов в образце, скорости счета позитронов, вылетающих вперед и назад относительно начального импульса мюонов и количестве событий, занесенных в спектр.
Набранные спектры передавались на ЭВМ mVAX-3400, которая использовалась в качестве файл-сервера. Обработка спектров производилась на ЭВМ различного типа, подключенных к институтской сети Ethernet. Обработка экспериментальных спектров производилась методом минимизации стандартного \chi2 - функционала временных гистограмм N(t) с различными функциями поляризации спина мюона P(t)[42]
.
Криостаты, используемые во время проведения экспериментов, позволяют поддерживать температуру изучаемых образцов в диапазоне 4,2 - 300 К[43]. Основные требования, предъявляемые к криостатам для mSR-экспериментов, можно сформулировать следующим образом:
- поддержание температуры образца с требуемой точностью; - обеспечение максимально возможного телесного угла для телескопов счетчиков, регистрирующих позитроны, по отношению к образцу;
- снижение фона в mSR-спектрах от мюонов, остановившихся в стенках криостата;
- использование в конструкции только немагнитных материалов;
- малая величина расхода хладоагентов.
Форма криостатов спектрометра "МЮСПИН" обусловлена конструкцией магнитов спектрометра. Криостаты должны помещаться между полюсами магнитов и оставлять место для сцинтилляционных счетчиков. Криостаты должны быть достаточно длинными, чтобы пристыковываться к емкостям с жидким азотом и гелиеми, а ФЭУ внутренних счетчиков криостата должны быть вынесены за панцири магнитов.
Схема криостата приведена на рис.2.6. Наружный разъемный вакуумный корпус 1 сечением 168x75 мм имеет стенки толщиной 1 мм. Внутри корпуса расположен азотный экран 2, который снижает теплоприток к контейнеру 3 для размещения образца 4. Контейнер представляет собой цилиндрическую латунную коробку, имеющую в направлении пучка мюонов окна 12 из бериллиевой бронзы толщиной 120 мкм, припаянных припоем ПСР-45. Герметизация крышки контейнера производится сплавом Вуда. Охлаждение образца осуществляется хладоагентами (гелием или азотом), поступающими по трубке внутри которой находится нагреватель 5. Вакуум внутри корпуса контролируется стандартной лампой 6 типа ПМТ-4. Криостат подсоединяется к системе откачки с помощью вакуумного вентиля 7. К корпусу прикреплены две трубы из нержавеющей стали, в которых закреплены световоды 9 длиной около 500 мм, вакуумно уплотненные в отверстиях бокового фланца. На концы труб навинчиваются кожухи с фотоумножителями 10. Сцинтилляторы 11 и 13 с короткими световодами находятся внутри азотного экрана и имеют с ним тепловой контакт. Для снижения теплопритока в холодную зону криостата между короткими и длинными световодами имеются вакуумные зазоры размером около 1 мм. При сборке криостата контейнер располагается между сцинтилляторами и центры сцинтилляторов совпадают с центром контейнера. Расстояние от центра образца до счетчиков около 15 мм (телесный угол 4pi/10). Общее количество вещества между счетчиками на пути пучка мюонов составляет 0.18 г/см . Расход жидкого гелия при температуре > 10 К и температурном градиенте на мишени около 0,1 К составляет около 1,5 л/час.
После помещения изучаемого образца в контейнер и сборки криостата, производится его откачка форвакуумным и диффузионными насосами. Затем криостат пристыковывается к системе охлаждения азотного экрана. Скорость охлаждения азотного экрана не должна превышать 2 К/мин, во избежания разрушения клеевого соединения сцинтилляторов со световодами. Температура азотного экрана и скорость его охлаждения контролируются прибором, описанном в п.2.4. Время охлаждения азотного экрана составляло около 2-х часов, и оно обычно использовалось для методических работ и для получения mSR спектров при температурах, близких к комнатной. Примерно столько же времени занимает и отогрев криостата для замены образца. Было изготовлено несколько криостатов с совместимыми счетчиковыми и мишенными частями. Это позволяло оперативно менять во время сеанса изучаемые образцы, вмонтированные в различные криостаты.
Одной из основных задач при проведении mSR экспериментов является изучение зависимости свойств вещества от температуры.
Наибольшее распространение при измерении температуры образцов при изучении свойств вещества с помощью mSR метода получили термометры сопротивления: угольные, германиевые, металлические[44]. Они характеризуются хорошей температурной чувствительностью, малой инерционностью, возможностью проведения непрерывных измерений. Фактором, усложняющим проведение точных измерений температур электрическими методами в mSR-экспериментах, является наличие магнитных полей.
Для измерения температуры образцов нами использовались датчики на основе монокристаллов германия с многокомпонентным легированием производства Института полупроводников АН Украины (ТПК-110, ТПК-510), [45]. Отличительной особенностью этих датчиков является слабое изменение сопротивления под влиянием магнитных полей, т.к. легирующие добавки у этих термометров подобраны так, что магнитосопротивление у них отрицательное и имеет минимум вблизи температуры 4,2 K. Они имеют коэффициент магнитосопротивления DR/R примерно в два раза ниже, чем другие полупроводниковые коммерческие термометры.
Требуемые температуры в диапазоне 4,2 - 300 K выдерживались с точностью не хуже 0,1 K в течение времени набора статистики в одной температурной точке. Это обеспечивалось применением автоматизированой системы управления температурой (АСУТ) на основе ПЭВМ типа IBM/PC.
Структурная схема АСУТ показана на рис. 2.9. В каждом криостате датчики устанавливались на образце (D1) и на входе в контейнер (D2) в потоке хладоагентов. С помощью коммутаторов типа КА OO4[46] выбранный датчик подключался к источнику тока БИТ и цифровому вольтметру типа Щ1516. Интерфейсом вольтметра служил 16-разрядный регистр ввода-вывода КИ 015 [46]. Источник измерительных токов состоял из 6 генераторов постоянного тока, выполненных по схеме "токового зеркала" и опорных сопротивлений для измерения величин токов. Источники обеспечивали стабильность измерительных токов не хуже 0.1% в диапазоне от 10 мкА до 10 мА. Величина измерительных токов контролировалась с помощью того же цифрового вольтметра через коммутатор.
Падение напряжения на датчиках температуры не превышало 550 мВ. По величине измерительного тока и падению напряжения на датчике вычислялось его сопротивление, а затем, по занесенной в ПЭВМ градуировочной кривой - температура.
Регулирование температуры осуществлялось по ПИД-закону[47], т.е. при выработке управляющего кода, подаваемого на цифро-аналоговый преобразователь (ЦАП) 6 учитывалось не только отличие измеряемой температуры от требуемой, но и скорость ее изменения. ЦАП на основе широтно-импульсной модуляции (ШИМ) состоял из двух блоков: генератора импульсов частотой 100 Гц и управляемого по магистрали КАМАК одновибратора с регулируемой длительностью импульсов. Максимальный код, который может быть записан в одновибратор 4095, соответствует длительности импульсов 10 мс, т.е. полному перекрытию импульсов. Импульсы с ЦАП поступают на усилитель мощности УМ, служащий для питания нагревателя криостата. Амплитуда импульсов, подаваемых на нагреватель, равна 15 В, сопротивление нагревателя - 20 Ом.
Требуемая температура поддерживалась с точностью не хуже 0,1 K.
Для измерения температуры азотного экрана использовались диодные датчики. Их действие основанно на температурной зависимости падения напряжения Ud на p-n переходе, через который протекает постоянный ток в прямом направлении[48]. Функция Ud при температурах 50 - 350 K имеет вид
Ud =U0 -g*T, (2.1)
где U0 - ширина запрещенной зоны при 0 K, g - величина определяющая чувствительность диодного датчика (2 - 3 мВ/K).
Мы используем датчики температуры этого типа на основе серийных диодов Д9К, прошедших предварительный отбор, термоциклирование и калибровку. Для измерения температуры с помощью диодных датчиков был разработан и изготовлен прибор, схема которого показана на рисунке 2.10. Питание диодных датчиков осуществляется током 100 мкA от источника тока I .
Напряжение, снимаемое с датчиков U , подается на усилитель-преобразователь, который выполняет преобразование, обратное (2.1), т.е. выходное напряжение этой схемы пропорционально температуре датчика Ub =к1 -к2Ud , где к1=U0 /g, к2 =1/g. Резисторы Rk1 - Rk2 , задающие коэффициенты к1 и к2 подбираются для каждого датчика индивидуально при его калибровке. Они расположены в отдельном корпусе и подключаются к прибору во время работы с теми или иными датчиками. Выходное напряжение с усилителя поступает на стрелочный измерительный прибор, шкала которого проградуирована в градусах Кельвина (0 .. 300 K) и на схему дифференциатора, которая служит для измерения скорости изменения температуры, т.е. выходное напряжение этой схемы пропорционально dT/dt. Скорость изменения температуры показывается другим стрелочным прибором в диапазоне от -5 до +5 K/мин. Переключателем на передней панели к прибору может быть подключен один из двух датчиков или опорное сопротивление R , служащее для контроля работоспособности. Измерительный прибор выполнен в конструктиве "Вишня" со встроенными блоками питания.
Вторым, после температуры, важным параметром при проведении mSR экспериментов является внешнее магнитное поле.
Магнитное поле на изучаемых образцах создавалось одним из двух панцирных электромагнитов с водяным охлаждением. Магнит для создания полей, перпендикулярных начальной поляризации мюонов, позволяет проводить исследования в диапазоне 0 - 6000 Э. Другой магнит использовался для проведения экспериментов в продольном магнитном поле напряженностью до 7000 Э. В этом магните дополнительно были установлены кольца Гельмгольца для измерения поляризации мюонного пучка. Они позволяют получать магнитные поля напряженностью до 50 Э, перпендикулярные начальной поляризации мюонов. Питание магнитов током до 20 A осуществлялось блоком ВСТН-20 производства РНЦ "Курчатовский институт". Относительная стабильность тока была не хуже 5x10 , что позволяло с такой же стабильностью получать поля до 800 Э в "перпендикулярном" магните и до 850 Э - в "продольном". Во время экспериментов с полями большей величины использовались мотор-генераторы системы питания трактов пучков фазотрона. Они обеспечивали стабильность тока не хуже 10-3 . Неоднородность поля в объеме, используемом для размещения образцов (диаметром 60 мм и толщиной 10 мм), составляла 10-4 и была определена muSR-методом с помощью веществ с малой скоростью деполяризации спина мюона.
Для измерения магнитных полей использовался метод ядерного магнитного резонанса в проточной жидкости[49]. Достоинствами такого способа измерения полей являются высокая точность измерений и широкий диапазон исследуемых магнитных полей (в нашем случае 0,1 - 7000 Э). Это обеспечивается тем, что сигнал ЯМР наблюдается в поле постоянного магнита, и измерительный генератор магнитометра не требует перестройки по частоте. Действие такого магнитометра основано на том, что в магнитных полях порядка нескольких килогаусс вода поляризуется. Возникшая разность заселенностей спиновых уровней протонов воды сохраняется в течение 1 - 5 сек. В измеряемом магнитном поле B при совпадении частот ларморовской прецессии спина протонов и наложенного переменного поля, создаваемым перестраиваемым генератором в катушке нутации происходит деполяризация воды (fNMR =gH*B, gH =2,457 КГц/Гс). Степень поляризации воды можно измерять с помощью колебательного контура порогового генератора, помещенного в анализирующий магнит, по уменьшению его добротности. При протекании неполяризованной воды через катушку индуктивности добротность контура падает и генерация срывается, что приводит к исчезновению сигнала ЯМР.
Структурная схема магнитометра представлена на рис.2.11.
Вода пропускается через поляризующий магнит П, создающий магнитное поле около 2 кЭ. Затем по пластиковой трубке диаметром 3 мм жидкость поступает в область измеряемого магнитного поля М, где находится контур нутации L. Катушка нутации (50 вит. ПЭЛ-0,2) намотана на трубку с водой, укрепленную на позитронном счетчике вблизи криостата. Катушка соединена с синтезатором частоты типа Ч6-31, позволяющим получать синусоидальные колебания частотой до 50 МГц с шагом от 0.01 Гц. Синтезатор частоты управлялся релейным регистром КВ 006 в стандарте КАМАК [50]. Частота, подаваемая на контур нутации измеряется частотомером Ч3-34. В качестве интерфейса частотомера используется входной регистр КР 005[51].
Из контура нутации вода поступает в катушку колебательного контура измерительного генератора Г [52]. Генератор, работающий на частоте 5 МГц, помещен в анализирующий магнит А, создающий магнитное поле 1170 Э что соответствует условию ЯМР-резонанса для протонов воды. Поле анализирующего магнита модулируется с помощью дополнительных катушек (на рисунке не показаны) с частотой 50 Гц. Сигнал с генератора Г измеряется цифровым вольтметром В7-34 с интерфейсом КР 005.
При измерении с помощью ЭВМ устанавливается код, подающийся на синтезатор частоты. Нагрузкой синтезатора частоты служит катушка нутации. Амплитуда ЯМР сигнала от анализирующего генератора Г измеряется через 0,2 сек. Это время необходимо для протекания воды от катушки нутации до контура измерительного генератора. Когда частота приближается к резонансной, амплитуда выходного сигнала анализирующего генератора уменьшается. После прохождения резонансной частоты, измеряемое вольтметром напряжение снова увеличивается. Таким образом процедура измерения сводится к нахождению частоты, соответствующей минимуму сигнала анализируещего генератора. Точность измерения магнитного поля не хуже 10-3 при B > 10 Гс, с уменьшением величины поля точность падает, но и в полях до 1 Гс она остается лучше чем 1% .
Круг задач, решаемых muSR-спектрометром, определяется его техническими возможностями. Характеристики пучков поляризованных мюонов, используемых комплексом "МЮСПИН" приведены в Таблице 2.2. Кроме этого mSR-спектрометр имеет ряд параметров, характеризующих его, как физический прибор. Основные параметры спектрометра "МЮСПИН" и основных его узлов приведены в Таблице 2.3.
Таблица 2.3 Основные параметры спектрометра "МЮСПИН"
Диапазон температур 4,2 - 300 K Точность поддержания температуры 0,1 K Диапазон магнитных полей 0 - 0,6 T Неоднородность магнитных полей 10^-4 Количество позитронных телескопов 2 Разрешающее время спектрометра ~ 1 нс Мертвое время спектрометра 20 нс Ширина канала спектрометра 5,5 нс (1,1 нс АПИ) Длинна гистограммы 10 мкс Отношение сигнал/фон 1000
Спектрометр "МЮСПИН" не уступает по своим параметрам mSR-спектрометрам, действующим на других ускорителях и приведенным в Таблице 1.2, а по отношению сигнал/фон превосходит их. Учитывая, что скорость набора информации в muSR спектроскопии временных интервалов ограничена временем жизни мюона из-за требования "1mu-1e", спектрометры на пучках слаботочных ускорителей (сравнительно с мезонными фабриками) также успешно используются для проведения mSR исследований. После реконструкции фазотрона ОИЯИ на установке "МЮСПИН" был проведен большой объем работ по изучению конденсированных сред мюонным методом. По результатам этих экспериментов за период 1988 - 1994 годы сделано более 30 публикаций в ведущих отечественных и международных журналах.