Төмен өлшемді кванттық жүйелердегі шашырау параметрлері мен спектрлерін сандық зерттеу әдісі
Алдыңғы бөлімдерде (1.1-1.2) көрсетілген есептеу сызбасы гармоникалық тұзақтың шектеулі геометриясындағы атомдардың соқтығысуы кезіндегі шашырау параметрлерін есептеуге бағытталған. Ол сондай-ақ ангармоникалық тұзақтан туннелдеу эффекттерін зерттеуге мүмкіндік береді.
Жұмыстың басқа маңызды мәселесі – тұзақ ангармонизмінің тұзақтағы молекулалық (екіатомды молекулалар) және атомдық спектрлерге әсерін зерттеу. Қарастырылып отырған есептеу сызбасының модификациясы бұл мәселені зерттеуге мүмкіндік береді. Меншікті мәндеріне есепті шешу үшін (1.1-1.4) гамильтонианы үшін [15] жұмысында берілген спектрлік әдіс қолданылады. Бұл әдіс корреляциялық функциясын есептеуге негізделген, мұндағы бастапқы күйді спектрдің туынды тесттік функциясы ретінде таңдауға болады [14]. Ал стационар емес Шредингер теңдеуінің шешімі (1.1-1.4) гамильтонианмен жоғарыда қарастырылған сандық әдіс көмегімен анықталады.
(1.1-1.4) түрдегі гамильтонианмен ангармоникалық тұзақтағы кванттық жүйелердің төменөлшемді спектрін есептеу үшін төмендегі Шредингер түріндегі қарапайым дифференциалдық теңдеулер жүйесінің шешімін сандық зерттеуді жоспарлаймыз
(1.3.1)
ДАТЕТ көрінісінде жеке туындылардағы бастапқы төменөлшемді Шредингер теңдеуін жуықтайды. Мұнда
(1.3.2)
мұндағы
(1.3.2) теңдеуін сандық интегралдау үшін есептеу математикасының стандартты әдістерін қолдану жоспарлануда: радиал туындыларды жуықтау үшін жоғарғы ретті ақырғы айырым әдісі және ішкі кеңістіктегі кері итерация алгоритмі.
1.4. Осцилляторда өрнектелу (ОӨ) әдісімен азөлшемді кванттық жүйелер спектрлерін есептеу
(1.1) Шредингер теңдеуін ОӨ әдісімен шешуге кірісейік [16]. Бұл әдіс өрістің скалярлық кванттық теориясының идеяларына және әдістеріне негізделген. Әдістің толық сипаттамасы [12] жұмыста баяндалған. Әдістің артықшылығы, оның аясында алынған потенциалдардың кең класы үшін ауытқу теориялар қатары өте тез үйлеседі. Нөлдік жуықтаудағы байланысқан жүйе гамильтонианының негізгі күйі кейбір Ω эффективті жиілігі бар осцилляторлық толқындық функциямен бейнеленеді деп болжанады. Ω жиілігі байланысқан күй беретін потенциалдан негізгі кванттық үлесін ескеретіндей анықталады. Ол үшін бастапқы гамильтониандағы әзірге белгісіз Ω жиілігі бар осцилляторлық бөлігі бөлініп алынады, ал каноникалық айнымалыларды, координата және импульсті туу және жою операторлары арқылы өрнектеледі, содан кейін оларды қалыпты түрде реттейді. Ары қарай осылай түрленген гамильтонианда каноникалық айнымалылар бойынша сызықты және квадраттық қосылғыштар болмауы қажет. Бұл шарт Ω - осциллятор жиілігін анықтайды және осцилляторда өрнектелу шарты деп аталады. Осылайша, алынған ауытқу теорияларының қатары эффективті қосылады.
R1 бірөлшемді кеңістіктегі жүйе гамильтонианы:
(1.4.1)
мұндағы:
(1.4.2)
потенциалдық энергия. Таза осцилляторлық бөлігін жазамыз:
(1.4.3)
мұндағы Ω әзірге ерікті оң параметр болып табылады. каноникалық айнымалыларын туу және жою операторлары арқылы көрсетейік:
(1.4.4)
Ары қарай каноникалық айнымалыларға арналған өрнектерді (1.4.3)-ге қоямыз және туу және жою операторлары бойынша реттейміз. Сонда қалыпты түрде көрсетілген потенциал мына түрге ие болады ([12]-ден толығырақ қараңыз):
(1.4.5)
мұндағы ~ - нормалдану символы және
(1.4.6)
Осы қатынастарды қолдана отырып, (3.3) гамильтонианды қайта жазайық ([12]-ден толығырақ қараңыз):
(1.4.7)
мұндағы
ал - негізгі күй немесе (1.4.3) гамильтониан вакуумының энергиясы. Осцилляторда өрнектелу әдісінің шартынан алынған үшін теңдеу бойынша минимум энергияны анықтайтын теңдеумен сәйкес келеді:
(1.4.8)
осциллятор жиілігі потенциалының бастапқы параметрлерінің функциясы болып табылады.
Осылайша, негізгі күй энергиясын есептеу үшін, жиілігі үшін (1.4.8) теңдігі орындалу қажет. потенциалдың нақты түрін қолдана отырып, келесі оңтайландырылған есепке келеміз:
(1.4.9)
(1.4.9)-дегі үшін теңдеу жалғыз бір оң айғақты түбірге ие болады.
Гамильтонианның бастапқы есебіндегі әсерлесуінің белгіленуі ауытқу теориясының (АТ) аясында деңгей энергиясына жоғарғы түзетулерді есептеуге мүмкіндік береді. ОК-тегі АТ бірінші реті нөлге тең. Екінші реті мына формуламен анықталады ([12]-ден толығырақ қараңыз)
(1.4.10)
мұндағы
(1.4.11)
ОӨ әдісіндегі негізгі күйдің толқындық функциясы, ал - сәйкес гамма және бета функциялар,
(1.4.12)
Осылайша, (1.4.1) бастапқы гамильтонианның негізгі күйінің энергиясы
(1.4.13)
Әртүрлі әсерлесу потенциалының параметрлері кезіндегі ОӨ әдісімен есептеу нәтижелері 4-7 суреттерде көрсетілген. Графиктерде сондай-ақ ОӨ әдісінің ауытқу теориясының екінші ретіне сәйкес келетін қисықтар көрсетілген.
Тұтас сызықтар – энергияның сандық есептелуі. Штрих үзік сызықтар – ОӨ әдісінің нөлдік жуықтауындағы энергия. Үздік сызықтар – ОӨ әдісіндегі екінші жуықтаудағы энергия.
Сурет 4 – (1.1) теңдеуі үшін ОӨ әдісінің қолданылу аймағы
- [25]-да есептелген дельта-функция түріндегі қосылған әсерлесу жағдайындағы энергия, - ОӨ нөлдік жуықтаудағы энергия, - ОӨ екінші жуықтаудағы энергия.
Сурет 5 – Әртүрлі әдістердің есептерін салыстыру; - энергияның сандық есептелуі
4-суретте әрқайсысында үш қисығы бар графиктер көрсетілген. Екі қисық ОӨ әдісімен есептеуге сәйкес келеді, ал үшіншісі сандық әдіспен есептелген. Бірінші графикте үш қисық қосылады да екі әдістің жақсы сәйкес келуі бақыланады. Графиктерден көрініп тұрғандай, теріс мәндер кезінде ОӨ әдісі кезіндегі сандық мәндерден ауытқи бастайды. ОӨ әдісіндегі екінші ретті АТ дәл осы потенциал параметрлері кезінде аз ауытқиды және теріс аймақтарында өте жақсы дәлдік көрсетеді. Алайда, оң үшін нөлдік сияқты, ОӨ әдісіндегі екінші ретті АТ да кезінде сандық есептеуден айтарлықтай ауытқиды.
5-суретте ОӨ әдісімен есептеулер сандық есептеулермен және кері шашырау ұзындығына тәуелді дельта-функция түріндегі псевдопотенциалды қолдану арқылы [8]-де жүргізілген есептеулермен салыстырылады.
Әртүрлі потенциал параметрлері кезінде ОӨ әдісінің ауытқу себебі ОӨ-те толқындық функция, нөлдік жуықтауда (1.4.11) гауссианмен бейнеленеді. Осылайша ОӨ әдісінің ауытқуы осцилляторлық толқындық функцияның гаусстық түрінің шынайы ТФ түрінен ауытқуымен түсіндіріледі, оны 6- және 7-суреттерден көруге болады.
Тұтас сызықтар – сандық есептеудегі ТФ. Штрихпунктирлі сызықтар – ОӨ әдісіндегі нөлдік жуықтаудағы ТФ. - энергияның сандық есептелуі, - ОӨ әдісіндегі нөлдік жуықтаудағы энергия, - ОӨ әдісіндегі екінші жуықтаудағы энергия.
Сурет 6 – үшін әртүрлі кезінде сандық есептеудегі ТФ мен ОӨ әдісіндегі ТФ салыстыру.
Тұтас сызықтар – сандық есептеудегі ТФ. Штрихпунктирлі сызықтар – ОӨ нөлдік жуықтаудағы ТФ. - энергияның сандық есептелуі, - ОӨ нөлдік жуықтаудағы энергия, - ОӨ екінші жуықтаудағы энергия.
Сурет 7 – үшін әртүрлі кезінде сандық есептеудегі ТФ мен ОӨ әдісіндегі ТФ салыстыру
6,7-суреттерде көрініп тұрғандай, кезінде оң аймағында пайда болатын ОӨ әдісінің сандық есептеуден ауытқуы осы екі әдістің толқындық функцияларының ауырмашылығына сәйкес болады. Осылайша, (2.1) берілген гамильтониан үшін ОӨ әдісінің қолданылуы әсерлесу потенциалы параметрлерінің функциясы болып табылады.
Бұл жұмыста бірөлшемді гармоникалық тұзақтағы екібозонды байланысқан жүйені сипаттайтын (2.1) ШТ негізгі күйінің энергиясы кері шашырау ұзындығының функциясы ретінде есептелді. Алынған тәуелділік [21] жұмысындағы псевдопотенциалдық әдісте орындалған есептеуден бірқатар ерекшеленеді. Біздің жағдай Гаусстың реалистік потенциалын қолданумен ерекшеленеді, оны таңдау орынды болып табылады, өйткені ол атомдарда қандай да бір құрылымның бар екенін ескереді. Сондай-ақ соңғы радиус потенциалынан дельта-функция потенциалына ауысу айтарлықтай жуықтау болып табылады. Осыған ұқсас үшөлшемді реалистік потенциалмен және нөлдік радиустың псевдопотенциалымен орындалған есептеулердің салыстырылуы екіөлшемді гармоникалық тұзақ үшін [30] жұмысында келтірілген.
Осцилляторда өрнектелу әдісін әсерлесу потенциалы параметрлерінің функциясы ретінде қолданылу аймағы анықталды. ОӨ әдісінің сандық есептеуден ауытқуы осцилляторлық ТФ гаусстық түрінің шынайы ТФ түрінен ауытқуымен түсіндіріледі. ОӨ әдісінің қолданылу аймағын анықтау келешекте ангармоникалық тұзақтағы екіатомды жүйе спектрін есептеу кезінде қажет болады.
Ангармоникалық тұзақтағы екіатомдық жүйенің спектрін есептеу кезіндегі негізгі қиындық, бұл жағдайда салыстырмалы координата мен массалар орталығының координаталарының айнымалыларын бөлу мүмкін болмайды. Есептеуді альтернативті әдістермен жүргізу мүмкіндігі нәтижелердің сенімділігіне себепкер болады.