Кандидат физико-математических наук, доцент
Нижний Новгород; Университет Людвига-Максимилиана, Мюнхен
Независимый эксперт; приглашенный ученый
УДК 539.19+544.18
Введение
Светочувствительные инициирующие взрывчатые вещества (ВВ) к настоящему времени довольно широко изучены и нашли применение в ряде конкретных технических устройств гражданского и военного назначения. Разработаны различные типы оптических детонаторов и воспламенителей, предложены методы тестирования прочности материалов и стойкости конструкций при одновременном нагружении их поверхности большой площади продуктами взрыва тонких слоев светочувствительного взрывчатого вещества при лазерном инициировании [1-5].
Наиболее перспективные светочувствительные взрывчатые вещества создаются на основе энергонасыщенных комплексов перхлоратов металлов с тетразолами [6, 7]. Исследование синтеза и изучение свойств новых перспективных светочувствительных взрывчатых веществ проводится постоянно в лабораториях ряда стран. Однако физико-химические свойства многих светочувствительных металлокомплексов еще не в полной мере изучены. Использование методов квантово-химических расчетов изолированных молекул подобных комплексов в газовой фазе может быть полезным для изучения свойств таких веществ и получения новой количественной и качественной информации по этому вопросу [8, 9].
Автором проводилось изучение свойств ряда подобных веществ с использованием методов квантово-химических расчетов изолированных молекул [10-17]. Изучались такие комплексные перхлораты, как перхлорат (5-гидразинотетразол)ртути(II), перхлорат бис(3-гидразино-4-амино-1,2,3-триазол)меди(II), перхлорат (5-цианотетразолато-N2)пентаамминокобальта(III), перхлорат (5-нитротетразолато-N2)пентаамминокобальта(III), перхлорат цис-бис(5-нитротетразолато-N2)тетраамминoкобальта(III) и перхлорат цис-бис(1-метил-5-аминотетразолато-N3,N4)тетраамминoкобальта(III). Расчеты выполнялись с использованием специализированной квантово-химической программы Gaussian версий 03 и 09. Первоначальное построение и оптимизация структур молекул проводились с использованием полуэмпирического метода РМЗ. Дальнейшее получение результатов проводилось в рамках метода теории функционала плотности (DFT) с использованием комбинированного функционала B3LYP и ряда базисных наборов электронных функций, таких как 3-21G, 3-21G**, LANL2DZ, 6-31G(d), 6-31+G(d) и 6-31+G(d,p).
Для указанных светочувствительных ВВ определялись оптимальные с энергетической точки зрения структуры, определялись геометрические характеристики молекул, такие как взаимное расположение различных групп, длины связей, величины углов между связями. Для некоторых ВВ рассматривались несколько различных молекул-изомеров, различающихся пространственным расположением входящих в их структуру перхлорат-анионов. Для дальнейшего более детального изучения использовалась молекула-изомер, имеющая минимальное значение рассчитанной общей энергии. Общая энергия в конкретном случае представляла собой сумму электронной энергии и энергии нулевых колебаний молекулы. Наряду с общей энергией молекул рассчитывались и другие энергетические характеристики, характеризующие различные возможные аспекты их разложения, в частности, энергетические зазоры HOMO-LUMO, энергии атомизации, энергии диссоциации некоторых групп и отдельных атомов. Определялись также инфракрасные колебательные спектры молекул, которые могут в какой-то степени характеризовать их чувствительность к внешнему, особенно импульсному облучательному, воздействию. Основное внимание здесь уделялось динамике связей, ответственных за возможные механизмы разложения молекул. Проводился совместный анализ результатов, полученных для всех исследуемых таким образом светочувствительных комплексных перхлоратов. Было проведено сопоставление полученных результатов с отдельными известными экспериментальными и расчетными данными, полученными другими исследователями для ВВ этого класса. Отмечена определенная взаимосвязь полученных в работах расчетных свойств молекул с известными экспериментальными результатами по их инфракрасным спектрам и порогам инициирования лазерным излучением.
Согласно существующим представлениям [7, 18] на первой стадии термического распада у изученных к настоящему времени амминатов кобальта(III) происходит диссоциация молекул аммиака и их выход из внутренней сферы комплекса. Это приводит к запуску последующих процессов перестройки комплекса вплоть до его окончательного разложения. В то же время на основании результатов выполненных расчетов в работе [15] было сделано предположение, что первой стадией термического распада этих молекул может также быть открытие тетразольного кольца с последующим отрывом молекулы азота.
Полученные в работах [13, 15] ограниченные и выполненные с использованием малого базисного набора расчетные результаты не позволяют делать окончательные выводы о предпочтительности того или иного первичного механизма распада рассматриваемых комплексных перхлоратов в каждом конкретном случае. Чтобы более детально разобраться с этим вопросом требуется проведение ряда дополнительных расчетов и использование в расчетах более высокого уровня теории. Для этого, в качестве предварительной меры, необходимо более тщательное рассмотрение полученных ранее результатов с учетом подробного рассмотрения процесса этих начальных стадий распада рассматриваемых перхлоратов амминотетразолатов кобальта (III). Это перхлораты (5-цианотетразолато-N2)пентаамминокобальта(III), (5-нитротетразолато-N2)пентаамминокобальта(III), цис-бис-(5-нитротетразолато-N2)тетраамминoкобальта(III) и цис-бис(1-метил-5-аминотетразолато-N3,N4)тетраамминoкобальта(III). В работе [17] краткие результаты подобного расчетного рассмотрения приведены для условий отрыва лигандов аммиака от молекул этих перхлоратов. В данной работе такие же краткие результаты представляются для условий разрушения тетразольных лигандов.
Результаты расчетов
Квантово-химические расчеты процесса растяжения связей между атомами азота в тетразольных кольцах выполняли с использованием программы Gaussian 09 [19] в приближении теории функционала плотности с использованием гибридного функционала B3LYP. Для проведения сравнительных расчетов разрушения тетразольных колец применяли малый базисный набор 3-21G. Для рассмотрения исходных состояний молекул и дополнительных уточняющих расчетов процесса растяжения связей в тетразольных кольцах использовали более серьезный базисный набор 6-31+G(d). Диаграммы растяжения связей между атомами азота в тетразольных кольцах получали путем последовательного увеличения расстояния между этими атомами и последующего расчета электронной энергии молекулы в каждой точке диаграммы с учетом релаксации ее структуры. В последующем будем иметь в виду, что результаты, полученные с использованием базисного набора 6-31+G(d), приводятся с указанием этого базиса, а результаты, полученные с использованием базисного набора 3-21G, приводятся без указания базиса.
Оптимизированная с использованием базисного набора 6-31+G(d) структура молекулы CP показана на рис. 1. Длины связей между атомами в тетразольном кольце молекулы приведены в табл. 1. На рис. 2 приведена диаграмма растяжения наиболее длинной связи между атомами азота N8-N13 в тетразольном кольце молекулы CP, а на рис. 3 показана структура молекулы в процессе разрыва тетразольного кольца. Видно, что достигнутая степень растяжения нагружаемой связи приводит к существенному повреждению тетразольного кольца, однако существенного необратимого разрушения не происходит, так как отсутствует переход к какой-то иной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 1. Структура молекулы CP в исходном состоянии (6-31+G(d)).
Табл. 1. Длины связей между атомами в тетразольном кольце молекулы CP (6-31+G(d))
Рис. 2. Диаграмма растяжения связи N8-N13 в тетразольном кольце молекулы CP.
Рис. 3. Структура молекулы CP в процессе разрыва тетразольного кольца по связи N8-N13 при x = 4.98 Å.
Оптимизированная с использованием базисного набора 6-31+G(d) структура молекулы NCP показана на рис. 4. Длины связей между атомами в тетразольном кольце молекулы приведены в табл. 2. На рис. 5 приведена диаграмма растяжения наиболее длинной связи между атомами азота N15-N16 в тетразольном кольце молекулы NCP, а на рис. 6 показана структура молекулы в результате разрыва тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца, при котором происходит отрыв молекулы азота. В результате разрушения тетразольного кольца система переходит к иной, практически стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 4. Структура молекулы NCP в исходном состоянии (6-31+G(d)).
Табл. 2. Длины связей между атомами в тетразольном кольце молекулы NCP (6-31+G(d))
Рис. 5. Диаграмма растяжения связи N15-N16 в тетразольном кольце молекулы NCP.
Рис. 6. Структура молекулы NCP после разрыва тетразольного кольца по связи N15-N16 при x = 5.29 Å.
Оптимизированная с использованием базисного набора 6-31+G(d) структура молекулы BNCP показана на рис. 7. Длины связей между атомами в обоих тетразольных кольцах молекулы приведены в табл. 3. На рис. 8 приведена диаграмма растяжения наиболее длинной связи между атомами азота N12-N19 в одном из тетразольных колец молекулы BNCP, а на рис. 9 показана структура молекулы в результате разрыва этого тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца, при котором происходит отрыв молекулы азота. В результате разрушения тетразольного кольца система переходит к иной, практически стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 7. Структура молекулы BNCP в исходном состоянии (6-31+G(d)).
Табл. 3. Длины связей между атомами в тетразольных кольцах молекулы BNCP (6-31+G(d))
Рис. 8. Диаграмма растяжения связи N12-N19 в тетразольном кольце молекулы BNCP.
Рис. 9. Структура молекулы BNCP после разрыва тетразольного кольца по связи N12-N19 при x = 5.59 Å.
Оптимизированная с использованием базисного набора 6-31+G(d) структура молекулы TCP показана на рис. 10. Длины связей между атомами в обоих тетразольных кольцах молекулы приведены в табл. 4. Для молекулы TCP рассмотрен разрыв четырех связей между молекулами азота, по две в каждом тетразольном кольце. В одном из этих случаев наряду с расчетом в базисном наборе 3-21 G выполнен также расчет в базисном наборе 6-31+G(d).
Рис. 10. Структура молекулы TCP в исходном состоянии (6-31+G(d)).
Табл. 4. Длины связей между атомами в тетразольных кольцах молекулы TCP (6-31+G(d))
На рис. 11 приведена диаграмма растяжения наиболее длинной связи между атомами азота N38-N50 в расположенном на рис. 10 справа тетразольном колеце молекулы TCP, а на рис. 12 показана структура молекулы в результате разрыва этого тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца, при котором происходит отрыв молекулы метилазида. В результате разрушения тетразольного кольца система переходит к иной стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 11. Диаграмма растяжения связи N38-N50 в тетразольном кольце молекулы TCP. Кружок – стабильное состояние после разрыва тетразольного кольца.
Рис. 12. Структура молекулы TCP после разрыва тетразольного кольца по связи N38-N50 при x = 3.22 Å.
На рис. 13 приведена диаграмма растяжения наиболее длинной связи между атомами азота N38-N50 в расположенном на рис. 10 справа тетразольном колеце молекулы TCP, а на рис. 14 показана структура молекулы в результате разрыва этого тетразольного кольца. Расчет в данном случае выполнен с использованием базисного набора 6-31+G(d). Видно, что подобная степень растяжения нагружаемой связи, как и в случае с базисным набором 3-21G, приводит к существенному разрушению тетразольного кольца, при котором происходит отрыв молекулы метилазида. В результате разрушения тетразольного кольца система переходит к иной стабильной структуре с соответствующей релаксацией растягивающего усилия. Следует отметить, что энергетический барьер реакции разрыва тетразольного кольца в этом случае составляет 241 кДж/моль, тогда как при использовании базисного набора 3-21G он был равен 196 кДж/моль.
Рис. 13. Диаграмма растяжения связи N38-N50 в тетразольном кольце молекулы TCP (6-31+G(d)). Кружок – стабильное состояние после разрыва тетразольного кольца.
Рис. 14. Структура молекулы TCP послее разрыва тетразольного кольца по связи N38-N50 при x = 3.53 Å (6-31+G(d)).
На рис. 15 приведена диаграмма растяжения связи между атомами азота N51-N52 в расположенном на рис. 10 справа тетразольном кольце молекулы TCP, а на рис. 16 показана структура молекулы в результате разрыва этого тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца, при котором происходит образование в структуре молекулы азидной группы. В результате разрушения тетразольного кольца молекула переходит к иной стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 15. Диаграмма растяжения связи N51-N52 в тетразольном кольце молекулы TCP. Кружок – стабильное состояние после разрыва тетразольного кольца.
Рис. 16. Структура молекулы TCP после разрыва тетразольного кольца по связи N51-N52 при x = 3.71 Å
На рис. 17 приведена диаграмма растяжения связи между атомами азота N33-N41 в расположенном на рис. 10 сверху тетразольном колеце молекулы TCP, а на рис. 18 показана структура молекулы в результате разрыва этого тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца, при котором происходит отрыв молекулы цианамида. В результате разрушения тетразольного кольца система переходит к иной стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 17. Диаграмма растяжения связи N33-N41 в тетразольном кольце молекулы TCP. Кружок – стабильное состояние после разрыва тетразольного кольца.
Рис. 18. Структура молекулы TCP после разрыва тетразольного кольца по связи N33-N41 при x = 3.69 Å
На рис. 19 приведена диаграмма растяжения наиболее длинной связи между атомами атомами азота N42-N53 в расположенном на рис. 10 сверху тетразольном колеце молекулы TCP, а на рис. 20 показана структура молекулы в результате разрыва этого тетразольного кольца. Видно, что подобная степень растяжения нагружаемой связи приводит к существенному разрушению тетразольного кольца с перестройкой его атомов. В результате разрушения тетразольного кольца система переходит к иной стабильной структуре с соответствующей релаксацией растягивающего усилия.
Рис. 19. Диаграмма растяжения связи N42-N53 в тетразольном кольце молекулы TCP. Кружок – стабильное состояние после разрыва тетразольного кольца.
Рис. 20. Структура молекулы TCP после разрыва тетразольного кольца по связи N42-N53 при x = 3.19 Å
В табл. 5 указаны значения энергетических барьеров реакций разрушения тетразольных лигандов и энергии стационарных состояний образующихся после разрушения систем для молекул NCP, BNCP и TCP. Сразу же отметим, что молекулярные расчеты, выполненные с использованием малого базисного набора 3-21G, могут давать только качественные, но не количественные результаты для энергетических характеристик молекул. Их результаты могут быть использованы только в сопоставительных целях для первоначальных оценок.
Табл. 5. Энергетические барьеры реакции разрушения тетразольных лигандов и энергии стационарных состояний после разрушения
* Basis set 6-31+G(d)
Так вот в этом плане, основываясь на полученных расчетных значениях энергетических барьеров реакции разрушения тетразольных лигандов, можно очень условно наметить ряд стойкости рассмотренных комплексных перхлоратов к внешнему энергетическому воздействию. Он будет иметь следующий вид, начиная от более стойкого CP, и далее BNCP, NCP и TCP. В случае же рассмотрения расчетных значений энергетических барьеров реакции отрыва молекул аммиака [17], условный ряд стойкости рассмотренных комплексных перхлоратов к внешнему энергетическому воздействию будет иметь совершенно иной вид, а именно: TCP, CP, NCP и BNCP. Разброс результатов в этом случае значительно меньше, чем в случае разрыва тетразольных колец. Таковы результаты предварительных расчетных оценок, выполненных с использованием малого базисного набора 3-21G.
Однако, как известно, определенное доверие к результатам квантово-химических расчетов может быть проявлено только при использовании базисных наборов, начиная с 6-31G(d) и выше. А вот тут для имеющихся еще более ограниченных результатов проявляется уже совсем иная картина. Так энергетические барьеры для реакции отрыва молекул аммиака [17] в молекулах NCP и BNCP при использовании базисного набора 6-31+G(d) существенно ниже таковых, полученных при использовании базисного набора 3-21G. В то же время в данной работе получено, что энергетический барьер для реакции разрыва тетразольного кольца в молекуле TCP при использовании базисного набора 6-31+G(d) ощутимо выше барьера, полученного при использовании базисного набора 3-21G. Таким образом, окончательные выводы о реальных первичных механизмах разложения и об относительной стойкости рассмотренных перхлоратов амминотетразолатов кобальта (III) можно будет делать только после выполнения более полных расчетов этих молекул с использованием базисного набора не ниже 6-31+G(d). Подобная задача и будет рассмотрена в следующей работе этого направления.
Заключение
Выполнено расчетное определение условий разрушения тетразольных лигандов в молекулах перхлоратов (5-цианотетразолато-N2)пентаамминокобальта(III) (вещество CP), (5-нитротетразолато-N2)пентаамминокобальта(III) (вещество NCP), цис-бис-(5-нитротетразолато-N2)тетраамминoкобальта(III) (вещество BNCP) и цис-бис-(1-метил-5-аминотетразолато-N3,N4)тетраамминoкобальта(III) (вещество TCP). Разрушение осуществлялось путем растяжения и разрыва одной из N-N связей тетразольного кольца. Расчеты проводились с использованием двух базисных наборов электронных функций, 3-21G и 6-31+G(d). Для первого из базисных наборов получены расчетные результаты для всех рассматриваемых молекул комплексных перхлоратов. Эти результаты позволяют выявить тенденцию снижения стойкости рассматриваемых молекул к разрыву тетразольного кольца, которая уменьшается от вещества CP к веществу TCP. Результат, полученный для одного из перхлоратов с использованием более полного и более точного базисного набора 6-31+G(d), довольно ощутимо отличается от результата, полученного с использованием малого базисного набора 3-21G. Это указывает на необходимость проведения всей серии выполненных и ряда дополнительных расчетов с использованием более точного базисного набора.
Рецензии:
20.11.2021, 1:49 Ашрапов Улугбек Товфикович
Рецензия: Голубев Владимир Константинович в статье выполнил расчетное определение условий разрушения тетразольных лигандов в молекулах перхлоратов. Насколько эти результаты расчетов соответствуют экспериментальным результатам в реальных условиях?
21.11.2021, 7:37 Ашрапов Улугбек Товфикович Отзыв: Аммиачные комплексы перхлората кобальта с производными тетразола (BNCP) являющийся средством инициирования повышает безопасность проведения взрывных работ и может использоваться в горнодобывающей отрасли. Поэтому, статья автора Голубева Владимира Константинович имеет актуальность. Статью рекомендую к публикации в журнале SCI-ARTICLE. |