Важнейшей задачей современных фундаментальных и прикладных исследований является снижение экологических рисков и минимизация опасности технологических катастроф [1]. Использование инициирующих и бризантных взрывчатых веществ (ВВ) в промышленном производстве (в том числе — в горной промышленности и строительстве) является необходимостью современного времени [1, 2]. Задача оптимизации капсюля оптического детонатора является весьма актуальной. Основными требованиями являются — высокая селективность к оптическому излучению [1, 3] с одновременной стойкостью к другим видам воздействия, минимальная плотность энергии инициирования состава [4]. Инициирующие взрывчатые вещества весьма чувствительны к нагреву, удару, электромагнитному полю [5, 6]. Ряд ученых предпочитают путь сенсибилизации хорошо исследованных штатных вторичных ВВ [7-10]. Применение ВВ понижает стоимость проходческих работ и безусловно будет расширяться в ближайшем будущем. На основе кристаллов азида серебра (инициирующего ВВ) уже созданы оптические детонаторы [11, 12], однако повышенная чувствительность к лазерному импульсу в этих объектах неотделима от опасности несанкционированного инициирования случайными механическими, термическими и электромагнитными наводками [13, 14].
Для создания селективно чувствительных к лазерному излучению материалов в [15] рассмотрено введение в существующие прозрачные ВВ светопоглощающих наночастиц металлов. Минимальная плотность энергии инициирования взрывчатого разложения бризантных ВВ с добавками наночастиц алюминия [10, 15 -17], кобальта [18, 19], никеля [20-22], хрома [23], олова [9, 24], ванадия [24-26], меди [27-31], серебра [32] и золота [32] составляет величину порядка 1 Дж/см2, что на два порядка меньше по сравнению с чистыми прессованными таблетками данного ВВ. Экспериментально проведен цикл исследований оптических и взрывных характеристик перспективных составов, в зависимости от природы металла и матрицы, форм размерных свойств металла для поиска новых материалов для капсюлей оптических детонаторов [15-32]. Модернизирована микроочаговая модель теплового взрыва [33-35], сформулированная изначально для интерпретации закономерностей взрывного разложения инициирующих ВВ [2, 5, 6, 11-14]. Основные направления модернизации модели заключались в учете дополнительных теплофизических процессов [34, 36] и оптических свойств наночастиц металлов [37-41], и образца [42-44]. В работе [24-26] показана перспективность использования наночастиц ванадия в качестве сенсибилизирующих добавок в PETN (штатное бризантное ВВ) для создания оптического детонатора с рекордной чувствительностью к импульсу первой гармоники неодимового лазера. Однако для PETN (вторичное ВВ, селективно чувствительное к лазерному импульсу) и ряда сенсибилизирующих добавок критическая плотность энергии инициирования взрывного разложения второй гармоники меньше, чем для основной [3, 15-26]. При этом важнейший аспект проблемы — рассеивающие свойства наночастиц. В работах [18, 42-44] показана возможность значительного (до 5 раз) увеличения освещенности в прессованных таблетках PETN с наночастицами алюминия, кобальта и никеля в результате многократного рассеяния излучения на наночастицах металлов. Одновременно, эти процессы существенно зависят от коэффициента преломления прозрачной матрицы, в которой находятся наночастицы. Цель настоящей работы: определение зависимостей коэффициента эффективности рассеяния наночастиц ванадия и характеристик индикатрисы рассеяния от радиуса наночастиц в прозрачных средах различной оптической плотности на второй гармонике неодимового лазера.
Любая задача может быть решена теоретически и практически. Методика определения индивидуальных оптических свойств наночастиц металлов в прозрачных матрицах основана на анализе зависимостей коэффициентов отражения и пропускания от толщины и массовой концентрации наночастиц [37] и апробирована для наночастиц алюминия. Задача сложна даже для одного радиуса наночастиц и одной матрицы. Для решения поставленной задачи практическими методами потребуется не один год напряженных исследований. В основном, экспериментальные ресурсы рационально затрачивать после предварительных теоретических исследований, когда выделена перспективная область проявления необходимых для создания исполнительных устройств оптических свойств, как это было ранее сделано для наночастиц золота, серебра, меди, алюминия, кобальта, никеля и ряда других металлов в PETN [15-44] и кристаллов азида серебра [11-14, 45, 46].
В [32] показана возможность теоретического решения поставленной в работе задачи в рамках теории Ми с необходимой точностью даже в спектральном диапазоне проявления плазмонного резонанса. В ряде работ [28-33] сформулирована и апробирована методика расчета оптических свойств наночастиц металлов в вакууме и прозрачных матрицах с показателем преломления ma, созданы пакеты прикладных программ для численного анализа актуальных процессов. Вначале необходима интерполяция комплексного показателя преломления (mi) на длину волны неодимового лазера (532 нм). Следующий этап — расчет оптических свойств наночастиц при выбранном радиусе металла и показателе преломления матрицы (ma). Значение mi = 3.1799 — 3.6083i оценено в работе [47] по методике [48].
Показатели преломления матрицы в настоящих расчетах изменялись от 1 до 2.5 с шагом 0.05. Очевидно, что начальное значение ma соответствует вакууму, ma от 1.3 до 1.8 –практически важным органическим веществам, включая гексоген и PETN [37-44] и кристаллам азида свинца. У более плотного инициирующего ВВ азида серебра показатель преломления около 2 [49]. Для построения зависимостей ряд продолжен до матрицы с показателем преломления 2.5. Отсутствие у ma комплексной составляющей значит пренебрежение возможной примесной составляющей поглощения, что на длине волны 532 нм, соответствующей энергии 2.33 эВ, является достаточно грубым приближением. Однако исследуемые матрицы принимались оптически прозрачными. При длительном хранении ВВ в экстремальных условиях (повышенные температуры, электромагнитное поле, повышенный радиационный фон) возможно существенное разложение матрицы, и разрабатываемая в статье методика неразрушающего контроля должна определять степень деградации матрицы [11 -14, 50]. Разложение кристаллов азида серебра сопровождается появлением твердых продуктов разложения (кластеров серебра), которые в районе второй гармоники неодимового лазера поглощают свет [50]. Радиус наночастицы ванадия выбран равным 50 нм, так как наночастицы такого радиуса целой группы исследованных в нашей лаборатории металлов наилучшим образом поглощают излучение второй гармоники неодимового лазера в матрице PETN [25].
Для решения подобных задач теория Ми была успешно адаптирована в работах [15-23]. В настоящей работе при постоянном значении радиуса наночастицы и длины волны излучения варьируется коэффициент преломления матрицы. Расчеты проводились в лицензионном (№ 824977) математическом пакете MatLab.6.5. Рассчитанные оптические свойства наночастиц ванадия радиуса 50 нм приведены в таблице и на рис. 1 и 2.
ma | Qsca | Scos | Kp | Kb |
1 | 0.4541 | -0.0744 | 0.4460 | 0 |
1.05 | 0.5549 | -0.0687 | 0.4504 | 0.2343 |
1.10 | 0.6672 | -0.0627 | 0.4551 | 0.3289 |
1.15 | 0.7887 | -0.0565 | 0.4599 | 0.3986 |
1.20 | 0.9161 | -0.0502 | 0.4649 | 0.4548 |
1.25 | 1.0452 | -0.0438 | 0.4699 | 0.5019 |
1.30 | 1.1711 | -0.0373 | 0.4750 | 0.5424 |
1.35 | 1.2889 | -0.0308 | 0.4801 | 0.5777 |
1.40 | 1.3941 | -0.0244 | 0.4852 | 0.6087 |
1.45 | 1.4836 | -0.0182 | 0.4902 | 0.6362 |
1.50 | 1.5554 | -0.0122 | 0.4950 | 0.6605 |
1.55 | 1.6091 | -0.0066 | 0.4997 | 0.6826 |
1.60 | 1.6454 | -0.0015 | 0.5040 | 0.7024 |
1.65 | 1.6661 | 0.0028 | 0.5079 | 0.7202 |
1.70 | 1.6733 | 0.0066 | 0.5113 | 0.7366 |
1.75 | 1.6696 | 0.0094 | 0.5142 | 0.7515 |
1.80 | 1.6573 | 0.0112 | 0.5164 | 0.7651 |
1.85 | 1.6387 | 0.0117 | 0.5178 | 0.7775 |
1.90 | 1.6157 | 0.0110 | 0.5183 | 0.7891 |
1.95 | 1.5901 | 0.0088 | 0.5179 | 0.7997 |
2.00 | 1.5632 | 0.0051 | 0.5165 | 0.8095 |
2.05 | 1.5360 | -0.0001 | 0.5141 | 0.8187 |
2.10 | 1.5094 | -0.0070 | 0.5106 | 0.8273 |
2.15 | 1.4839 | -0.0153 | 0.5061 | 0.8354 |
2.20 | 1.4600 | -0.0249 | 0.5007 | 0.8430 |
2.25 | 1.4379 | -0.0357 | 0.4944 | 0.8502 |
2.30 | 1.4176 | -0.0474 | 0.4875 | 0.8570 |
2.35 | 1.3991 | -0.0598 | 0.4801 | 0.8635 |
2.40 | 1.3824 | -0.0724 | 0.4723 | 0.8696 |
2.45 | 1.3673 | -0.0851 | 0.4644 | 0.8754 |
2.50 | 1.3535 | -0.0977 | 0.4565 | 0.8807 |
В первом столбце таблицы приведены значения показателей преломления матриц, при которых проводились расчеты оптических свойств наночастиц ванадия радиуса 50 нм. Во втором — значения коэффициентов эффективности рассеяния (Qsca) наночастиц.

При рассмотрении результатов расчета мы видим нетривиальную зависимость Qsca(ma). В начале рост коэффициента преломления естественно сопровождается увеличением Qsca как и для никеля и кобальта [18-22]. Для наночастиц ванадия в вакууме коэффициент эффективности рассеяния составляет величину 0.4541 и далее быстро увеличивается. До ma = 1.7 Qsca увеличивается практически линейно до значения 1.6733. Последующее увеличение коэффициента преломления до 2.5 приводит к уменьшению Qsca. Это эффект, причины которого предстоит исследовать в дальнейших работах. Уменьшаясь в этом диапазоне ma почто на 30%, Qsca достигает значения 1.3535 на границе рассматриваемой области. Для наглядности зависимость Qsca(ma) построена на рис. 1 полыми кружками и штриховой линией. Звездочками и сплошной линией помечены рассчитанные значения коэффициентов эффективности поглощения (Qabs).
В практически важном диапазоне ma от 1.5 до 2, в котором находятся коэффициенты преломления прозрачных ВВ Qabs остается примерно постоянным в районе 1.75. Это достаточно высокое значение коэффициента эффективности поглощения, которое показывает, что сечение поглощения почти в два раза больше геометрического. Изменение ma в этом диапазоне практически невозможно диагностировать портативными (что важно для внедрения) приборами.

Рассмотрим особенности рассеяния света наночастицами в различных матрицах (столбец 2 таблицы и рис. 1 — штриховая линия). От вакуума до матрицы с ma = 1.702 зависимость Qsca(ma) возрастающая от значения Qsca = 0.4541 до абсолютного максимума 1.6733. Достаточно быстро Qsca увеличивается почти в 5 раз. До коэффициента преломления 1.6 линейное увеличение Qsca сохраняется. На основе этого эффекта возможна разработка методики неразрушающего контроля капсюля оптического детонатора. Деградация или разложение ВВ приводит к изменению коэффициента преломления матрицы, что сопровождается уменьшением коэффициента отражения от образца. Уменьшение Qsca существенно усиливает эффект.
Значительно разнообразнее индикатриса рассеяния наночастиц — нормированное угловое распределение интенсивности рассеяния света. Эта функция определяет единичную плотность энергии, рассеянную на выбранный угол от первоначального направления света. Геометрия задачи определяет сферическую симметрию зависимости индикатрисы, типичную для исследуемого направления [51]. Нормировочный интеграл необходим для проверки корректности расчета. На рисунке 2 приведена индикатриса рассеяния, рассчитанная для наночастиц ванадия радиусов 50 нм в матрице с показателем преломлений 1.55 (как у PETN). Рассчитанная индикатриса рассеяния имеет симметричный вид, почти сферическая форма. Анизотропия однократного рассеяния характеризуется средним значением косинуса угла (столбец 3 таблицы). В анизотропных средах индикатриса рассеяния вытянута вперед и значение Scos ~ 1 [51]. Индикатриса рассеяния зависит от длины волны света, для атмосферы и видимого света (около 2 гармоники неодимового лазера) значение Scos > 0.8 [51]. Ограничение сверху достаточно естественно: среднее значение косинуса угла не может быть больше 1. Знак равенства означает, что вся рассеянная энергия распространяется в первоначальном направлении под углом 0. Этот процесс нельзя трактовать как рассеяние. Поэтому очевидно Scos cos = 0. Преимущественное рассеяние назад определяет отрицательное значение среднего значения косинуса угла, которое для наночастицы ванадия радиусом 50 нм реализуется в начале и конце показателя преломления. На актуальных наночастицах наблюдается эффект смены знака среднего косинуса индикатрисы. До ma 1.65, значение Scos отрицательно, далее до 2.05 — положительное, и — снова отрицательное.

Наглядной характеристикой индикатрисы, определяющей часть рассеянной вперед энергии, является величина Kp, которая определяет часть энергии, рассеянной вперед по отношению к направлению первоначального распространения света. Эта величина представлена в 4 столбце таблице. Две последние характеристики связаны между собой. В некоторых источниках сообщается о достаточности первой их них. В самом деле, при отрицательном среднем косинусе значение Kp меньше 0.5, что достаточно очевидно для индикатрисы, близкой к сферической. Однако реальная индикатриса наночастицы ванадия радиуса 50 нм в различных матрицах уже настолько далека от изотропной, что это очевидное равенство перестает работать. Это наглядно видно из значений последних 2 параметров (средний косинус и доля рассеянной вперед энергии) при показателе преломления матрицы 1.6, когда первый параметр отрицателен (-0.0015), а второй больше 0.5 (0.504). Такая же ситуация с наночастицей ванадия в матрицах с коэффициентами преломления от 2 до 2.2, когда отрицательные значения среднего косинуса «уживаются» с величинами Kp больше 0.5.
Для интерпретации этого эффекта построен рис. 3, на котором представлена рассчитанная индикатриса рассеяния света с длиной волны 532 нм наночастицы ванадия радиуса 50 нм в матрице с показателем преломления 2 (кристаллы азида серебра). Мы видим, что направление максимального рассеяния света в первой четверти соответствует углу примерно 60 о, косинус которого всего 0.5. Для второй четверти (отражение назад) направление максимального рассеяния света соответствует углу примерно 150 о, модуль косинус которого значительно больше 0.5 (0.866). Это обстоятельство делает определение Kp необходимой задачей.
Оптимизация оптических детонаторов идет в направлении снижения необходимой для инициирования взрывного разложения плотности энергии импульса. Повышение коэффициента освещенности в образце (в некоторых случаях в 5 раз) эффективно означает снижение критической плотности энергии инициирования (в 5 раз). В PETN с наночастицами ванадия на второй гармонике рассеяние света назад не приводит к эффективной возможности повышения коэффициента освещенности. Однако капсюль оптического детонатора является плоскопараллельной пластиной [15], в которой появляется дополнительная возможность повышения коэффициента освещенности за счет рассеяния света в диапазоне углов полного внутреннего отражения (Kb). В пятом столбце таблицы представлены рассчитанные значения части энергии, которая после однократного рассеяния не может выйти из прессованной таблетки с коэффициентом преломления из первого столбца таблицы. Для матриц с малыми значениями ma Kb меньше 0.5, но для актуальных матриц этот параметр быстро растет, достигая 88%. Эта часть рассеянной энергии формирует компоненту света, которая может либо поглотиться на наночастице и нагреть ее, либо рассеяться на другой наночастице. Выйти из образца этот свет не в состоянии. Проведенное исследование позволяет сделать вывод о перспективности наночастиц ванадия в качестве добавки в таблетки PETN для создания капсюля оптического детонатора. Автор выражает благодарность научному руководителю аспиранту кафедры ХТТ и ХМ Галкиной Е.В.