Расчет радиационных параметров электронного перехода А–Х-димера кальция

Язык труда и переводы:
УДК:
539.194
Дата публикации:
06 декабря 2022, 00:19
Категория:
Перспективные направления исследования необратимых физических процессов
Авторы
Смирнов Александр Давыдович
МГТУ им. Н.Э. Баумана
Федоров Владимир Витальевич
МГТУ им. Н.Э. Баумана
Аннотация:
Для синтеза холодных молекул из холодных атомов необходимы данные о спектроскопических постоянных и радиационных параметрах. Проведен расчет радиационных параметров (коэффициентов Эйнштейна, сил осциллятора, факторов Франка — Кондона, волновых чисел электронно-колебательных переходов) электронного перехода димера кальция. Расчеты выполнены на основе потенциальных кривых возмущенного осциллятора Морзе, построенных в данной работе.
Ключевые слова:
потенциальная кривая, волновое уравнение, коэффициент Эйнштейна, сила осциллятора, фактор Франка — Кондона
Основной текст труда

Введение

За последние годы выросло число экспериментальных и теоретических публикаций, связанных с получением с помощью методики фотоассоциации при низких температурах холодных двухатомных молекул. Важность этих исследований  вызвана изучением динамики столкновения холодных молекул, проведением химических реакций при низких температурах. Для синтеза холодных молекул из холодных атомов необходимы данные о спектроскопических постоянных и радиационных  параметрах.

В работах [1–5] были рассчитаны радиационные параметры для электронных переходов A^{1}\Sigma _{u\bullet }^{+}-\bullet X^{1}\Sigma _{g\bullet }^{+}  и B^{1}\Sigma ^{+}-\bullet X^{1}\Sigma ^{+}  димеров щелочных металлов (Li2, Na2, K2, Rb2, Cs2,). Расчеты проводились квантово-химическим методом на основе полуэмпирических потенциальных кривых. Сравнение результатов расчета с экспериментальными данными показало эффективность метода расчета. В настоящей работе проведены расчетырадиационных параметров A^{1}\Sigma _{u\bullet }^{+}-\bullet X^{1}\Sigma _{g\bullet }^{+}  электронного перехода димера кальция. Спектральный анализ электронного перехода A^{1}\Sigma _{u\bullet }^{+}-\bullet X^{1}\Sigma _{g\bullet }^{+}   молекулы Ca2 был проведен с использованием лазерных методик и техники высокого разрешения в работе  [6]. Были получены молекулярные постоянные для основного и возбужденного электронных состояний. В отличие от димеров щелочных металлов основные электронные состояния димеров щелочноземельных металлов являются слабосвязанными. В работе [6] отмечается, что димеры атомов щелочноземельных металлов являются перспективными кандидатами для создания газовых эксимерных лазеров. 

Построение потенциальных кривых

Для построения потенциальных кривых использовалась модель потенциальной кривой возмущенного осциллятора Морзе (ВМ) [7]:

U\,(r)\;=V_{e}\,(\,y^{2}\;+\;\sum \limits _{n=4}^{\infty }{b_{n}\;y^{n}\;)},                                                              (1)

где y=1-exp[-\alpha (r-r_{e})] , r_{e}  —равновесное межъядерное расстояние; V_{e},\alpha ,b_{n} —параметры потенциальной функции ВМ.

Для расчета параметров ВМ димера кальция использовались  экспери-ментальные  спектроскопические постоянные, полученные в работе [6].

На рисунке приведены построенные потенциальные кривые для основного и возбужденного электронных состояний димера кальция.

Потенциальные кривые основного и возбужденного состояний электронного перехода А–Х-димера кальция

Расчет молекулярных постоянных

На основе построенных потенци- альных кривых были рассчитаны колебательные энергии, а также вращательные и центробежные постоянные. Колебательные энергии рассчитаны в  результате численного решения  волнового уравнения Шредингера методом Нумерова – Кулея [8]. Уравнение Шредингера для колебательного движения имеет вид

d^{2}\psi _{v}(R)/dR^{2}+8\pi ^{2\;}\mu /h^{2}[E_{v}-U(R)]\;\psi _{v}(R)\;=0.

где h — постоянная Планка; \mu — приведенная  масса  молекулы. В табл. 1 приведены  результаты расчета.

 

Таблица 1. Результаты расчета колебательных энергий E_{v} , вращательных B_{v} и центробежных D_{v} молекулярных постоянных для основного и возбужденного состояний  электронного  перехода    A^{1}\Sigma _{u\bullet }^{+}-\bullet X^{1}\Sigma _{g\bullet }^{+} молекулы  Ca2

v

Ev

Bv

Dv

Ev

Bv

Dv

 

Х

A

0

32,23

0,0456

9,38(–8)

68,08

0,0593

4,2608(–8)

1

95,11

0,0449

9,71(–8)

203,23

0,0585

4,2601(–8)

2

155,88

0,0441

10,08(–8)

336,90

0,0580

4,2612(–8)

3

214,55

0,0433

10,47(–8)

469,07

0,0576

4,2621(–8)

4

271,13

0,0426

10,89(–8)

599,71

0,0572

4,2631(–8)

5

325,60

0,0418

11,35(–8)

728,79

0,0568

4,2641(–8)

6

377,87

0,0409

11,85(–8)

856,29

0,0565

4,2651(–8)

7

426,24

0,0401

12,40(–8)

982,18

0,0561

4,2661(–8)

8

476,40

0,0393

13,00(–8)

1106,42

0,0557

4,2670(–8)

9

522,47

0,0384

13,66(–8)

1228,98

0,0552

4,2680(–8)

10

566,43

0,0375

14,34(–8)

1349,84

0,0548

4,2689(–8)

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

B_{v}=\beta \langle v\left|R^{-2}\right|v\rangle ;                                                               (2)

D_{v}=\beta ^{2}\sum _{u\neq v}{\frac {\langle v\left|R^{-2}\right|v\rangle ^{2}}{\left(E_{u}-E_{v}\right)}}.                                                          (3)

Рассчитанные молекулярные постоянные приведены в табл.2 и сравниваются с экспериментальными  значениями.

 

Таблица 2. Сравнение результатов расчета молекулярных постоянных (см–1)  с экспериментальными данными

Молекулярные

постоянные

Х

А

\omega _{e}

\omega _{e}

 

64,98

65,07

136,61

136,65

\omega _{e}x_{e}

\omega _{e}x_{e}

 

–1,061

–1,087

–0,723

–0,721

\omega _{e}y_{e}

\omega _{e}y_{e}

 

3,26(–3)

3,46(–3)

–2,62(–3)

–2,85(–3)

B_{e}

B_{e}

 

4,60(–2)

4,61(–2)

5,81(–2)

5,82(–2)

a_{e}

a_{e}

 

–7,37(–4)

–7,41(–4)

–2,53(–4)

–2,97(–4)

 

D_{e}

D_{e}

9,3(–8)

9,1(–8)

4,4(–8)

4,3(–8)

Примечание. 1-я строка — результаты расчета по уравнению Шредингера; 2-я строка — экспериментальные данные [6].

Как видно из хорошо согласующихся результатов сравнения молекулярных постоянных, рассчитанных по уравнению Шредингера на основе построенных потенциальных кривых и экспериментальных данных, что характеризует адекватность построенных потенциальных кривых.   

Расчет радиационных параметров

К радиационным параметрам относятся коэффициенты Эйнштейна для спонтанного излучения A_{\nu '\nu ''} , cилы осциллятора f_{\nu '\nu ''} , факторы Франка — Кондона (ФФК)  q_{\nu '\nu ''} , волновые числа  электронно-колебательных переходов \nu _{\nu '\nu ''} .

Коэффициент Эйнштейна (с–1) и сила осциллятора для поглощения (безразмерная величина) связаны с функцией дипольного момента  электронного перехода от межъядерного расстояния R_{e}(r)  (в ат. ед.) соотношениями [10]:

A_{\nu '\nu ''}=2,062\times 10^{-6}{\frac {\nu _{\nu '\nu ''}^{3}(2-\delta _{0,\Lambda '+\Lambda ''})}{(2-\delta _{0,\Lambda '})}}\times \left[\langle \psi _{\nu '}(r)|R_{e}(r)|\psi _{\nu ''}(r)\rangle \right]^{2};                     (4)

f_{\nu '\nu ''}^{\text{погл}}=3,0376\times 10^{-6}{\frac {\nu _{\nu '\nu ''}(2-\delta _{0,\Lambda '+\Lambda ''})}{(2-\delta _{0,\Lambda ''})}}\times \left[\langle \psi _{\nu '}(r)|R_{e}(r)|\psi _{\nu ''}(r)\rangle \right]^{2},                  (5)

где \Lambda — проекция орбитального момента количества движения электронов на межъядерную ось ( \Lambda =0,1,2,\dots для состояний \Sigma ,\Pi ,\Delta ,\dots ); \delta _{0,\Lambda } — символ Кронекера, равный 1, если \Lambda =0 и нулю   для других   значений  \Lambda ; \nu _{\nu '\nu ''} — волновое число  электронно-колебательного перехода; \psi _{\nu '}(r),\psi _{\nu ''}(r)   — колебательные  волновые функции возбужденного и основного электронных состояний.

ФФК характеризуют относительное распределение интенсивностей электронно-колебательных полос и представляют собой квадраты интеграла перекрывания колебательных волновых функций комбинирующих электронных состояний

q_{\nu '\nu ''}=\langle \nu '|\nu ''\rangle =\left|\int \limits _{0}^{\infty }\Psi _{\nu '}(r)\Psi _{\nu ''}(r)dr\right|^{2}.                               (6)

Для расчета волновых чисел  колебательных полос электронного  перехода A^{1}\Sigma _{u\bullet }^{+}-\bullet \bullet X^{1}\Sigma _{g}^{+}  использовалось выражение

    \nu _{\nu '\nu ''}=T_{e}'+E'(\nu )-E''(\nu ),                                              (7)

где  T_{e}' — электронная энергия возбужденного состояния; E'(\nu ),E''(\nu ) — колебательные энергии для возбужденного и основного электронных состояний.

Были рассчитаны радиационные параметры для электронного перехода A^{1}\Sigma _{u}^{+}-\bullet X^{1}\Sigma _{g}^{+}  молекулы Ca2 (0\leq v'\leq 8,0\leq v''\leq 5) .

В табл. 3 приведены результаты расчета.

 

Таблица 3. Результаты расчета радиационных параметров A^{1}\Sigma _{u}^{+}-\bullet X^{1}\Sigma _{g}^{+}     электронного перехода молекулы Сa2

v''

v’ = 0

v’ = 1

v' = 2

v' = 3

v' = 4

v' = 5

v' = 6

v' = 7

v' = 8

 

0

0,00151

1,13(+5)

8,31(–4)

14286,7

0,00852

6,52(+6)

4,70(–3)

14422,6

0,02515

1,96(+6)

1,39(–2)

14558,0

0,05152

4,10(+6)

2,85(–2)

14693,0

0,08241

6,69(+6)

4,56(–2)

14827,6

0,10970

9,08(+6)

6,08(–2)

14961,7

0,12650

1,07(+7)

7,03(–2)

15095,4

0,12984

1,12(+7)

7,24(–2)

15228,6

0,12096

1,06(+7)

6,76(–2)

15361,5

 

1

0,00754

5,61(+5)

4,61(+5)

14233,9

0,03272

2,48(+6)

1,80(–2)

14359,8

0,02515

5,47(+6)

3,90(–2)

14495,3

0,09978

7,86(+6)

5,50(–2)

14693,0

0,10010

8,04(+6)

5,53(–2)

14764,7

0,07181

5,88(+6)

3,97(–2)

14898,8

0,03323

2,78(+6)

1,84(–2)

15032,5

0,00585

4,99(+5)

3,25(–3)

15165,8

0,00102

8,91(+4)

5,71(–4)

15298,6

 

2

0,02020

1,49(+6)

1,11(–2)

14163,2

0,06524

4,90(+6)

3,59(–2)

142991,1

0,09676

7,41(+6)

5,33(–2)

14434,5

0,08075

6,30(+6)

4,45(–2)

14569,5

0,03413

2,71(+6)

1,88(–2)

14704,0

0,00217

1,75(+5)

1,19(–3)

14838,1

0,00808

6,71(+5)

4,49(–3)

14971,8

0,03680

3,12(+6)

2,05(–2)

15105,1

0,05829

5,03(+6)

3,25(–2)

15237,9

 

3

0,03853

2,82(+6)

2,12(–2)

14104,7

0,08906

6,64(+6)

4,91(–2)

14240,6

0,08105

6,15(+6)

4,46(–2)

14375,9

0,02685

2,07(+6)

1,48(–2)

14510,9

0,00002

1,32(+3)

1,34(–5)

14645,5

0,02284

1,84(+6)

1,26(–2)

14779,6

0,05132

4,27(+6)

2,84(–2)

14913,3

0,04710

3,95(+6)

2,61(–2)

15046,6

0,01928

1,65(+6)

1,07(–2)

15179,4

 

 

4

0,05882

4,27(+6)

3,24(–2)

14048,3

0,09240

6,93(+6)

5,09(–2)

14184,2

0,04137

3,11(+6)

2,27(–2)

14319,6

0,00018

1,37(+4)

9,83(–5)

14454,6

0,02406

1,88(+6)

1,33(–2)

14589,2

0,05068

4,05(+6)

2,80(–2)

14723,3

0,03305

2,69(+6)

1,83(–2)

14856,9

0,00394

3,26(+5)

2,18(–3)

14990,2

0,00523

4,45(+5)

2,92(–3)

15123,0

 

 

5

0,07668

5,52(+6)

4,23(–2)

13994,2

0,07604

5,58(+6)

4,19(–2)

14129,6

0,00902

6,71(+5)

4,94(–3)

14265,5

0,01323

1,01(+6)

7,31(–3)

14400,4

0,04768

3,70(+6)

2,63(–2)

14534,9

0,03073

2,43(+6)

1,69(–2)

14669,1

0,00136

1,09(+6)

7,45(–4)

14802,8

0,01180

9,75(+5)

6,55(–3)

14936,0

0,03601

3,03(+6)

2,00(–2)

15068,9

Примечание.  Первая строка —  q_{\nu '\nu ''} ;  вторая  строка — A_{\nu '\nu ''} , с–1; третья строка — f_{\nu '\nu ''}^{\text{погл}} ;  четвертая строка – \nu _{\nu '\nu ''} , см–1.

При расчетах A_{\nu '\nu ''}  и f_{\nu '\nu ''}^{\text{погл}}  использована зависимость дипольного момента электронного перехода от межъядерного расстояния R_{e}(r) для электронного перехода A^{1}\Sigma _{u\bullet }^{+}-\bullet \bullet X^{1}\Sigma _{g}^{+}  молекулы Ca2 [11]. Радиационные параметры рассчитаны впервые.

Литература
  1. Смирнов А.Д. Расчет спектроскопических постоянных и радиационных параметров для электронных переходов А–Х и В–Х-димера лития. Оптика и спектроскопия, 2012, т. 113, № 4, с. 387–394.
  2. Смирнов А.Д. Расчет спектроскопических постоянных и радиационных параметров для электронных переходов А–Х и В–Х-димера натрия. Оптика и спектроскопия, 2010, т. 109, № 5, с. 739–745.
  3. Смирнов А.Д. Расчет спектроскопических постоянных и радиационных параметров для электронных переходов А–Х и В–Х-димера калия. Вестник МГТУ им. Н.Э. Баумана. Естественные науки, 2013, № 2, с. 67–85.
  4. Смирнов А.Д. Расчет спектроскопических постоянных и радиационных параметров для электронного перехода А–Х-димера цезия. Журнал прикладной спектроскопии, 2010, т. 77, № 5, с. 661–667.
  5. Смирнов А.Д. Энергетические и радиационные свойства электронного перехода В–Х-димеров цезия и рубидия. Инженерный журнал: наука и инновации. 2013, вып. 6, 53 с.
  6. Vidal C.R. The molecular constants and potential energy curves of the Ca2 system from laser induced fluorescence. J Chem Phys, 1980, vol. 72, no. 3, pp. 1864–1874.
  7. Huffaker J.N. Diatomic molecules as perturbed Morse oscillators. J Chem Phys, 1976, vol. 64, pp. 3175–3181.
  8. Cooley J.W. An improved eigenvalue corrector formula for solving the schrödinger equation for central fields. Math Comput, 1961, vol.15, pp. 363–374.
  9. Цюлике Л. Квантовая химия. Т. 1. Основы и общие методы. Москва, Мир, 1976, 512 с.
  10. Laher R.R., Khakoo M.A., Antic–Jovanovic A. Radiative transition parameters for the А–Х band system of Ag2. J Mol Spect, 2008, vol. 248, pp. 111–121. DOI: https://doi.org/10.1016/j.jms.2007.12.003
  11. Bussary B.Theoretical spectre of Ca2 molecule. J Chem Phys, 2006, vol. 125, pp. 114315–114330.
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.