Аннотация: Приведены методы генерации структурных формул и изображений для описания формы ядер в невозбуждённом состоянии во всём диапазоне массовых чисел, включая изотопы и изомеры.
Ключевые слова: форма ядра, структурная формула, 3D-графика.
Фiзико-математичнi науки
УДК 539.143.5
Жижко Володимир Абрамович
незалежний експерт,
Жижко Владимир Абрамович
независимый эксперт
Zhyzhko V.A.
independent expert
ПОБУДОВА АТЛАСУ ФОРМ АТОМНИХ ЯДЕР
ПОСТРОЕНИЕ АТЛАСА ФОРМ АТОМНЫХ ЯДЕР
CONSTRUCTION ATLAS FORMS OF ATOMIC NUCLEUSES
Анотацiя: Наведено методи генерацii структурних формул та зображень для опису форми ядер у незбудженому станi для всього дiапазону масових чисел, включаючи iзотопи та iзомери.
Ключовi слова: форма ядра, структурна формула, 3D-графiка.
Аннотация: Приведены методы генерации структурных формул и изображений для описания формы ядер в невозбуждённом состоянии во всём диапазоне массовых чисел, включая изотопы и изомеры.
Ключевые слова: форма ядра, структурная формула, 3D-графика.
Summary: The methods of generating structural formulas and images offered to describe the shape of the nuclei in the unexcited state over the entire range of mass numbers, including isotopes and isomers.
Keywords: shape of nucleus, structural formula, 3D-graphics.
В [1] предложен подход к описанию форм ядер, основанный на тезисе о минимальном количестве сильных связей между нуклонами. Элементарной ячейкой ядерного пространства оказался тетраэдр, а формы всех ядер можно рассматривать как конфигурации из тетраэдров. Это накладывает ограничения на форму конфигураций, что позволяет выявить некоторые свойства ядер, обусловленные их формой. К примеру, удалось объяснить химический состав продуктов деления U235 исходя только из формы этого ядра. Этот подход по сути есть исследование свойств симметрии структур из тетраэдров. Подобный подход оказался продуктивным в физике твёрдого тела, структурной химии. Данная статья – это продолжение и детализация работы [1].
Для названий основных понятий примем следующие сокращения:
ПИФ – половина интерфейса (участок поверхности блока, пригодный для соединения с другим блоком),
КФ – конфигурация (взаимное расположение тетраэдров в блоке),
МУГ – матрица углов между гранями поверхности КФ.
Будут рассмотрены КФ для всей периодической таблицы элементов, т.е. для ядер со значениями массовых чисел А от 1 до 240. Следует также учесть наличие изотопов (около 2000) и изомеров (несколько сотен) в указанном диапазоне А. В среднем каждому изобару должно соответствовать около 10 разных КФ, для лёгких ядер меньше, для тяжёлых – больше. В таблице 1 приведены основные блоки для формирования всех КФ. Основным блокам присвоены имена из одной буквы. Формулы некоторых лёгких и средних ядер даны в табл.2 для сплошного заполнения всего диапазона А, что станет понятно ниже.
Таблица 1. Базовые блоки для построения всех ядер Таблица 2. Формулы некоторых ядер
Рис.1 Блоки из 5 тетраэдров: а – блок S (ПИФы 5+), б – ПИФы – 2-
Рис.2 КФ из 6 тетраэдров Рис.3 КФ из 7 тетраэдров
Рис.4 Блок K ( ПИФ 6+) Рис.12 Блок Y = O3G (ПИФы 6=)
Рис.5 Блок G (ПИФ 6=) Рис.6 Блок W (ПИФы 5-)
Рис. 7 Блок М (3 ПИФа 6+) Рис.8 Блок Н (ПИФ 6-)
Рис. 9 Блок О (4 ПИФа 6=) Рис.10 “Сферическая” КФ,модель ядра Ne(20,10)
а б
Рис.11 а – блок J (ПИФы 5- и 6-), б - изображение того же блока после
операции Smooth (см. п.4)
2. Генерация формул полного набора КФ
Как показано в [1], допустимыми сочетаниями ПИФов являются (2+,2-), (5+,5), (6+,6-), (6=,6=). Между ПИФами есть зависимости. Например, если на часть ПИФ (6=) некоторой КФ поместить блок D, то получится ПИФ (6-). Это можно записать в виде
6=(2-,2+)D => 6- или проще (6=) + D => (6-).
Подобные трансформационные соотношения можно написать для всех ПИФов
(6=) + D => (6-)
(6-) + T => (5-)
(5-) + S => (5+)
(5+) + T => (6+)
(6+) + D => (6=)
Видно, что набор ПИФов является замкнутым, с помощью нескольких небольших добавок можно перейти от исходного ПИФа к желаемому. Это значит, что если при построении множества КФ из простых блоков взять любую разрешённую пару ПИФов, то это позволит построить большую часть КФ. Иначе говоря, множества КФ, построенные на разных парах ПИФ, будут в значительной мере пересекаться. Это позволяет ограничиться для начала комбинациями из блоков с любой парой ПИФов, проще всего (6=,6=). Это блоки O, F, G, X,Y. Массовое число очередной КФ есть сумма массовых чисел блоков, образующих КФ, A(G)=10, A(O) =16, A(F)=17, A(X)=46, A(Y)=47. Будем учитывать только те варианты КФ для одного значения А, которые имеют разные значения к-ва граней поверхности КФ. Количество граней несложно подсчитать – оно состоит из суммы граней составляющих КФ блоков за вычетом граней в местах соединения блоков. При соединении двух блоков через ПИФ 6= общая поверхность уменьшается на 12 граней. Общая поверхность имеет
Sum( Qi *Si ) - 12 * ( Sum (Qi) - 1) граней,
де Qi – к-во блоков i-того вида, входящих в формулу КФ, Si - к-во граней блоков: S(X) = S(Y) =36, S(O)=S(F)= 24, S(G)= 16. Поверхность блока O (и F) состоит из 4 ПИФов 6=, у блоков X и Y по 2 ПИФа 6=, у блока G – один.
Добавление блока O или F в КФ увеличивает к-во вариантов продолжения КФ на 3 ПИФа, добавление блока X или Y не меняет к-во вариантов, добавление блока G уменьшает это к-во (т.е. блок G - это “заглушка” КФ).
Ниже приведена простая программа вычисления всевозможных комбинаций блоков O, F, G, X, Y . Программа – на языке Clipper 5 для работы с таблицами локальных баз данных. Таблица formula в программе состоит из 7 целочисленных полей:
qA - массовое число полученной КФ, qS – к-во граней поверхности КФ,
qO, qF, qG, qX, qY - к-во соответствующих блоков в КФ,
а также вспомогательного символьного поля для упорядочения результата.
USE formula ; ZAP ; nnpp := 1
FOR o := 0 TO 15 ; FOR f := 0 TO 15
FOR x := 1 TO 5 ; FOR y := 1 TO 5
maxg := 2*(o+f+1) // макс.к-во незанятых ПИФов
FOR g := 0 TO maxg
ac := 16*o + 17*f + 47*y + 46*x +10*g // массовое число
s := 24*(o+f) +16*g +36*(x+y) - 12*(o+f+g+x+y - 1) // граней поверхности
acs := LTRIM(STR(ac)); IF ac < 100 ; acs := "0"+acs ; ENDIF
ss := LTRIM(STR(s)); IF s < 100 ; ss := "0"+ ss ; ENDIF
IF ac < 241 ; APPEND BLANK; REPLACE ;
formula->qo WITH o, formula->qf WITH f,;
formula->qg WITH g, formula->qa WITH ac,;
formula->qx WITH x, formula->qy WITH y,;
formula->qs WITH s, formula->ikey WITH acs+ss
END IF
END FOR //g
END FOR ; END FOR
END FOR ; END FOR
INDEX ON ikey TO fa UNIQUE// упорядочение таблицы по значениям массового числа с // с выбором неповторяющихся к-в граней поверхности КФ
GO TOP
// печать результата
? "Формулы конфигураций из блоков с ПИФ 6= “
DO WHILE ! eof()
IF qo > 0 ; ? qo,”O” ; ELSE; ? " " ; END IF
IF qf > 0; ?? qf, “F” ; ELSE; ?? " " ; END IF
IF qx > 0; ?? qx, “X” ; ELSE; ?? " " ; END IF
IF qy > 0; ?? qy, “Y” ; ELSE; ?? " " ; END IF
IF qg > 0 ; ?? qG,“G” ; ELSE; ?? " " ; END IF
?? qa, qs, nnpp
nnpp++ ; SKIP
END DO
Результат данной программы – почти 1000 формул КФ, ниже приводится фрагменты файла печати
Ф о р м у л а Масс Граней NN
………………..
1O 2G 36 32 8
1F 2G 37 32 9
2O 1G 42 40 10
1O 1F 1G 43 40 11
…………………………..
1F 1Y 2G 84 56 60
4O 2G 84 68 61
1O 4F 84 72 62
3O 1F 2G 85 68 63
5F 85 72 64
…………………………………….
3F 2X 1G 153 100 281
5O 1F 1X 1G 153 112 282
3O 5F 2G 153 116 283
9F 153 120 284
1O 3X 154 96 285
3F 1X 1Y 1G 154 100 286
4O 2F 1X 1G 154 112 287
2O 6F 2G 154 116 288
9O 1G 154 124 289
……………………………………….
5Y 235 132 896
3F 4X 235 144 897
5F 1X 2Y1G 235 148 898
5O 1F 3X 235 156 899
3O 5F 2X 1G 235 160 900
1O 9F 1X 2G 235 164 901
8O 3F 1X 1G 235 172 902
6O 7F 2G 235 176 903
3O11F 235 180 904
13O 1F 1G 235 184 905
……………………………………….
Чтение результата простое, например, ядро с массовым числом А=153 может быть образовано в виде КФ, которые можно описать одной из следующих формул: 3F2XG, 5OFXG, 3O5F2G, 9F. Для набора блоков (O, F, X, Y, G) количество K вариантов КФ для одного значения массового числа A возрастает с ростом массового числа почти линейно
K = int ((A- 60)/15).
Естественно предположить, что наиболее стабильному ядру из изобаров с массовым числом А соответствует наиболее компактная КФ, т.е. КФ с наименьшим к-вом S граней поверхности. Такие КФ для более чем 200 значений А, полученные из блоков с ПИФ 6=, приведены в табл. 3.
Таблица 3. Компактные КФ из блоков с ПИФ (6=).
Содержание табл.3 понятно. Например, ядру с А=237 соответствует формула 4XFO2G. Форма ядра с А=46 и формулой O3G показана на рис.12. Из табл.1 видно, что O3G = X, у блока Х два ПИФа (6=). Поэтому прогнозируемыми являются формулы КФ со значениями А, кратными 46, а именно:
А=92 (Zr) ф-ла 2Х; А=138 (Ba) – 3Х: A=184(W) – 4X; A=230(Th) – 5X.
Таким образом, таблицы 1, 2 и 3 содержат описание форм ядер всей периодической системы элементов. Некоторые замечания по таблице 1.
В отличие от табл.3, таблица 1 – результат ручного подбора на макете вариантов КФ, исходя из соображений компактности и симметрии в той мере, в какой это позволяет набор тетраэдров. С другой стороны, не все КФ в табл.1 можно отождествлять с формами ядер, некоторые КФ – это кластеры, конструктивные блоки. Примеры: блок S (А=5) один из важных для “сборки” многих КФ, но стабильных ядер с А=5 не существует. Для А=10 приведено 3 блока, два из которых W и G важны для “сборки” многих КФ.
Получить компактную и симметричную КФ для изобара с А=16, к которому принадлежит и дважды магический изотоп О(16,8), оказалось возможным только предположив отсутствие одного тетраэдра в центре КФ (чего можно достичь, убрав один тетраэдр из центра КФ с А=17).
Детальный учёт различных блоков со всеми комбинациями ПИФов позволяет построить к-во КФ, необходимое для отображения форм всех изотопов и изомеров (около 2.5 тысяч). Алгоритм подобен тому, что для (6=).
3. Операции на уровне поверхности конфигураций
Описание с помощью формул не является полным в том смысле, что формула описывает перечень блоков, образующих КФ, но не указывает порядок соединения блоков (в основном это относится к блокам с ПИФ 6=). Это означает, в частности, что одна формула может описывать несколько КФ. Пример – формула 3F2XG (A=153, S=100) может быть реализована в следующих вариантах
F F F G
FXXG XFXG GFXX FFXF и т.д.
F F F X
(Кстати, это обстоятельство решает вопрос о тысячах вариантов КФ).
В [1] намечен способ более точного описания КФ, чем с помощью формул.
Это описание взаимного примыкания граней тетраэдров с помощью матрицы смежности. Возможен и другой точный способ описания КФ, основанный на описании взаимного примыкания рёбрами граней поверхности КФ. Такой способ описания несколько проще. Поверхность КФ состоит из граней тетраэдров – правильных треугольников, расположенных под некоторыми углами друг к другу. С учётом оговорок [1], можно считать, что таких углов на поверхности КФ всего может быть 4. Приняв за меру угла двугранный угол между гранями тетраэдра, получим допустимые значения углов на поверхности КФ: 1,2,3,4. Под мерой угла между гранями на поверхности КФ будем понимать – сколько тетраэдров можно уложить между этими гранями. Значение 0 и 5 означают прилегание, но такие грани уже не принадлежат поверхности (а находятся в объёме КФ). В таком представлении КФ можно описать как матрицу углов между гранями (МУГ) поверхности КФ. МУГ – квадратная симметричная матрица размером A*A с нулевой диагональю.
В каждой строке и каждом столбце МУГ есть 3 ненулевых члена (каждая грань имеет трёх соседей). Номера строк и столбцов соответствуют номерам граней поверхности. Из-за низкой плотности заполнения МУГ, удобной программной моделью МУГ может служить таблица side с полями:
ncf - номер КФ (<массовое число> . <номер варианта изобара> )
nss - номер грани поверхности
nsn1 - номер грани-соседа 1, ang1 - угол 1 (между nss и nsn1)
nsn2 - номер грани-соседа 2, ang2 - угол 2
nsn3 - номер грани-соседа 3, ang3 - угол 3
Оперирование с описанием поверхности предполагает ручное кодирование данных (которое заключается в нумерации граней поверхности и заполнении таблицы side) для младших блоков (табл.1) и дальнейшую программную обработку, которая включает:
(далёкий дискретный аналог теоремы Остроградского-Гаусса).
Операция “сборки” формирует МУГ блока большего размера из двух МУГ блоков меньшего размера и дописывает новый МУГ в таблицу side.
Из-за большого (для статьи) размера программы ограничимся короткими выдержками и комментариями. Описание последовательности углов в ПИФах имеет вид
pify := { { 3,2,3,2,3,2,"6=",6}, { 3,3,2,3,3,2,"6+",6}, { 3,2,2,3,2,2,"6-",6},;
{ 3,3,3,3,3,0,"5+",5}, { 2,2,2,2,2,0,"5-",5} }
Это описание есть точное определение ПИФов. Исполнительная часть программы начинается с проверки формальной правильности описаний блоков. Процедура распознавание ПИФов блока выполняет поиск последовательностей соседних граней, чередование углов между которыми соответствует одному из описанных в массиве pify.
Процедура соединения двух блоков состоит из ряда рутинных операций, наиболее содержательными из которых являются:
От МУГ сравнительно несложно перейти к развёртке поверхности КФ. Развёртка – это расположение всех граней поверхности в одной плоскости путём изменения углов между гранями и, как следствие, разведения некоторых граней (разрыв по рёбрам). На рис.13 показан вариант развёртки блока W, цифры на рёбрах взяты из МУГ и обозначают углы для складывания развёртки в объёмную фигуру.
Рис.13 Вариант развёртки блока W Рис.14 Потенциал ядерных сил
Размер текста описанной программы составляет около 500 строк.
Как результат построения, МУГ не обладает достаточной наглядностью, поэтому следующим шагом является визуализация формы КФ.
4. Получение изображений трёхмерных объектов
“Объёмное” изображение на основе таблиц МУГ можно получить в рамках подхода, принятого в пакете OpenSCAD. Этот подход основан на том, что пользователь составляет задание не на уровне команд меню, а пишет программный код на объектно-ориентированном языке, а основные объекты уже определены. Подобный подход возможен и в пакете SketchUp, где есть окно для кода на языке Ruby.
Далее излагается более простой путь - работа с командами меню. Рисунки в п.1 получены с использованием пакетов Wings3D и SketchUp (оба пакета freeware). SketchUp имеет обширную библиотеку плагинов, в частности, объект ‘тетраэдр” создаёт плагин Geodesic (с параметром subdivision=0). Для пакета Wings3D этот объект – встроенный. Работа с этими пакетами - это ручное составление сложных КФ из простых блоков, в рассматриваемой ситуации - по формулам из таблиц 1 и 2.
Сборка КФ из тетраэдров в пакетах 3D-графики – это сложный процесс из нескольких операций. Объявленные в пакетах свойства “прилипания” успешно проявляются в основном на ортогональных телах. Добавление очередного тетраэдра в SketchUp можно выполнять с помощью достаточно трудоёмкой операции “копирование тела поворотом вокруг ребра” (folding on selected axis) с установкой величины угла поворота в поле числовых значений. Существенно проще добавлять тетраэдры в Wings3D. Последовательность операций такова:
С помощью этого средства можно создать набор “лёгких” блоки с А < 20, все остальные КФ можно получить соединением меньших блоков с помощью команд начального сближения Move и Rotate и команды “прищепки” Tweak.
Оба пакета содержат средства всестороннего обзора 3D-объектов на экране (рендеринга), а также средства портирования файлов для обмена с другими пакетами. Wings3D v 1.5.4 не имеет средств печати, поэтому можно воспользоваться печатью из SketchUp или Windows. Путь к печати:
Wings3D – экспорт .wings => .3ds
SketchUp – импорт .3ds => .skp, экспорт .skp => 2D-графика .bmp
Затем можно открыть файл .bmp программой Windows Paint и Сохранить как в формате:
Коротко о сферах применения каждого из рассмотренных инструментов.
1. С точки зрения практичности наиболее удобной является работа с макетом, например, с намагниченными тетраэдрами. Скорость работы с макетом примерно на два порядка выше скорости работы с программой 3D графики из-за того, что в программах операции перемещения, поворота, совмещения, кругового обзора (как и все остальные) можно выполнять только последовательно. Макет – удобное средство, чтобы быстро получить предметное осязаемое представление о форме КФ, заданной некоторой формулой.
2. Представление КФ в виде МУГ или её объёмного аналога нужно для численных расчётов свойств моделируемых ядер (спектры, электрические и магнитные моменты, характеристики рассеяния).
3. Графика удобна для наращивания коллекции изображений КФ, в том числе с элементами анимации (напр., формат gif). Можно указать и преимущества графики как информационного объекта над “материальным” макетом.
5. Физические аспекты модели формы ядер
Всё вышеизложенное – это описание техник оперирования с моделью, в которой элементарной ячейкой пространства ядра является тетраэдр.
Остаётся выяснить, как соотносятся геометрические и топологические характеристики модели со свойствами других моделей и реальных ядер.
1. Предложенная модель формы ядер позволяет описать несколько тысяч различных КФ, что даёт возможность сопоставить индивидуальную форму каждому изотопу, которых имеется до десятка и более у каждого из ядер (в зависимости от массового числа). Основной результат данной работы содержится в таблице 3 и алгоритме получения этой таблицы.
2. Принято считать, что магические ядра имеют форму, близкую к сфере. Более точно было бы утверждать, что магическим ядрам соответствуют локальные минимумы внешнего квадрупольного момента (ВКМ). Модель характерна наличием многих вариантов КФ для одного значения А. Сделав некоторое предположение о распределении протонов и нейтронов в КФ, можно вычислить его ВКМ. Сравнив с известными значениями ВКМ для разных КФ с одинаковым значением А (изобаров), можно сопоставить форму КФ каждому изобару, в том числе магическим.
3. В [1] сделана попытка объяснить величину энергии фотоделения U235 6.2 МэВ особенностями строения этого ядра, а именно, как резонанс продольных колебаний линейной цепочки кластеров YYYYY (5Y). При этом увеличение частоты облучения в n= 2,3,4… раз должно привести к изменению состава продуктов деления (увеличению доли Ar и K), а также увеличению выхода энергии деления и к-ва свободных нейтронов. Энергия в 12.5 МэВ малоинтересна из-за наложения эффектов фотоделения и гигантского дипольного резонанса (ГДР). Данные по выходу нейтронов из мишени U235 при обстреле электронами, протонами и дейтонами с энергиями 25, 50, 75 и 100 МэВ приведены в [2] . Рост выхода нейтронов действительно наблюдается, по составу продуктов деления и выходу энергии деления данные не приводятся, механизмы реакции также не обсуждаются. Можно считать, что фотоделение вызывают гамма-кванты тормозного излучения электронов или дебройлевские волны протонов и дейтонов. В [1] показано, что ядра с массовыми числами от 230 до 235 имеют практически одинаковую структуру и частоты фотоделения. В диапазон этих значений А попадают 26 изотопов от актиния до плутония.. Одно из следствий – для целей генерации нейтронов можно эффективно использовать смеси указанных изотопов.
4. Представление о форме ядер в виде вытянутых или сплюснутых эллипсоидов [3] является упрощённым, на что указывает и различие вычисленных и измеренных значений магнитных моментов ядер. Модель формы ядер из тетраэдров позволяет говорить о более разнообразных формах, чем даже трёхосные эллипсоиды.
5. Какие из КФ с фиксированным значением А (изобаров) отнести к изотопам (N+Z=A, Z – различны), какие – к изомерам, - этот вопрос может быть решён на основании расчётов электрических и магнитных свойств КФ или характеристик рассеяния. Поскольку КФ - синоним формы ядра, это даёт основание для численного расчёта спектра энергий. КФ – это конкретизация распределённого объёмного потенциала взаимодействия нуклонов. В данной модели потенциал состоит из (2A - S/2) однотипных слагаемых ядерных сил (рис.14) и Z(Z-1)/2 кулоновских слагаемых.
6. Недавние работы [4] по структуре ядерной материи километровых масштабов (“ядерные макароны” коры нейтронных звёзд) дают основание считать тему элементарной ячейки ядерного вещества и формы ядер такой, которая сохраняет актуальность.
7. Рассмотренные в п.2 ПИФы составлены из 5 или 6 граней и соответственно имеют такое же к-во внешних вершин. Такое же к-во вершин характерно для полиэдрических кластеров – фуллеренов, полых молекул сферической или эллипсоидной формы [5], первых наноструктур. Такое совпадение неслучайно, поскольку Р.Б.Фуллер проектировал свои прочные “геодезические” поверхности из треугольников и объёмных конструкций с треугольными гранями 0 – тетраэдров.
Выводы
Для предложенной модели ядерного пространства с элементарной ячейкой в виде тетраэдра приведен алгоритм построения набора формул для описания формы ядер химических элементов во всём диапазоне массовых чисел, что позволяет создать полный атлас индивидуальных форм для всех изотопов и изомеров.
Предложены разные подходы к формированию изображения форм ядер в привязке к пакетам 3D графики. Затронуты вопросы использования модели как источника информации о пространственном распределении потенциалов взаимодействия нуклонов для расчёта различных ядерно-физических величин.
Литература
Примiтка. Усi таблицi та малюнки розробленi автором.