|
|
|
Основные понятия квантовой механикиКлассические механика и электродинамика при попытке применить их к объяснению атомных явлений приводят к результатам, находящимся в резком противоречии с опытом. Глубокое противоречие теории с экспериментом свидетельствует о том, что построение теории, применимой к атомным явлениям, требует фундаментального изменения в основных классических представлениях и законах. В данной главе мы рассмотрим несколько примеров использования новых базовых принципов квантовой механики, следуя первой главе книги Ланау и Лифшица. Напомним, что при переходе от квантовой к классической механике движение, описываемое волновой функцией, в общем случае отнюдь не переходит в движение по определенной траектории. Для того чтобы получить движение по определенной траектории, надо исходить из волновой функции особого вида, заметно отличной от нуля лишь на очень малом участке пространства (так называемый волновой пакет). В квазиклассическом пределе размеры этого участка надо уменьшать вместе с константой Планка. Тогда можно утверждать, что в квазиклассическом пределе волновой пакет будет перемещаться в пространстве по классической траектории частицы. Следуя задаче N17 из задачника З.Флюгге , изучим волновой пакет в случае свободного движения частицы. Рассмотрим одномерный волновой пакет наиболее общего вида
Для того, чтобы интеграл сходился, амплитуда C(k) должна стремиться к нулю при k-> Построим волновой пакет таким образом, чтобы в начальный момент Разложение этой функции по плоским волнам имеет вид
Вычислив коэффициент
Вычислим этот интеграл и посмотрим на его поведение во времени. Попытаемся проделать все дальнейшие вычисления из задачника вместе с Maple. Зададим аргумент экспоненты под интегралом > restart: Попробуем вычислить этот интеграл > otv:=int(exp(arg),k=-infinity..infinity); В данном случае ответ означает следующее утверждение: если > whattype(otv); Для данного типа выражения можно подсчитать число операндов в выражении и их последовательность > nops(otv); Таким образом, необходимый нам ответ - это второй операнд в выражении типа функция > op(2,otv); Естественно, что тот же самый ответ можно получить, следуя классическому подходу З. Флюгге. Далее мы повторим все вычисления из задачника для того, чтобы приобрести практический опыт работы с системами символьных вычислений. Начнем с упрощения аргумента экспоненты стоящей под интегралом > factor(arg); Прямой подход не приводит к необходимому ответу и, поэтому, используем ряд встроенных процедур из библиотеки (пакета) student > arg:=student[completesquare](arg,k); Аргумент разлагается на два слагаемых. Первое, квадратичное по > In:=int(exp(-b^2*z^2),z=-infinity..infinity); Мы опять получили условие на знак квадрата параметра > interface(showassumed=0); этот интеграл легко вычислить > In(c); Заметим, что в данном случае мы задали табличный интеграл в виде отображения, которое будет использоваться нами в дальнейшем. Таким образом, для вычисления необходимого нам интеграла остается научиться заменять переменную интегрирования. Итак, выделим квадратичное слагаемое в аргументе > arg1:=op(1,arg); и определим новую переменную > k0:=solve(arg1,k); После этого произвести замену переменных можно следующей встроенной командой из пакета student > In1:=student[changevar](k-k0[1]=y,Int(exp(arg1),k=-infinity..infinity),y); Как видим, изменилось не только подинтегральное выражение, но и пределы интегрирования. Тем не менее, теперь можно вычислить этот интеграл, накладывая дополнительные ограничения на значения параметров. Можно использовать и другой путь - просто выделить коэффициент > assume(a>0,h>0,t>0,m>0); и подставить его в табличный интеграл > In(b0); Теперь рассмотрим вторую часть аргумента, которая не зависит от > arg2:=op(2,arg); Выделим явно зависимость от переменной > arg2:=student[completesquare](arg2,x); Итак, после всех подготовительных работ, окончательный результат имеет вид > psi:=A*a/sqrt(2*Pi)*In(b0)*exp(arg2): После вычисления волновой функции волнового пакета, изучим поведение плотности волнового пакета > rho:=evalc(simplify(abs(psi)^2)); Заметим, что в данном случае мы вынуждены "заставить" систему Maple произвести комплексные вычисления, используя встроенную команду > rho:=simplify(rho); Детальное обсуждение полученного результата может быть найдено в задачнике Флюгге. Теперь мы рассмотрим плотность потока вероятности для волнового пакета, которое имеет вид
Перепишем это определение в системе Maple в абстрактном виде > s:=phi->(h/(2*I*m))*(conjugate(phi)*diff(phi,x)-phi*diff(conjugate(phi),x)); Подставляя в данное определение волновую функцию волнового пакета, мы можем вычислить плотность Однако, если рассматриваемого нами волнового пакета мы вычислим производные волновых функций > assume(q,real); assume(x,real); > simplify(diff(psi,x)/psi); > simplify(diff(conjugate(psi),x)/conjugate(psi)); то можно видеть, что производные пропорциональны самим функциям. Это означает, что в данном случае поток пропорционален плотности > v:=simplify(s(psi)/rho); Заметим так же, что без дополнительных предположений о переменных > v:=simplify(evalc(v)); Для того, чтобы посмотреть на дисперсию скорости на графике, подставим частные значения параметров задачи в данное выражение для величины > v1:=simplify(subs(h=1,m=1,a=1,q=2,v)); Далее нарисуем зависимость отношения плотности потока вероятности к вероятности, т.е. зависимость "скорости" от координаты в различные моменты времени > p0:=plot(subs(t=0,v1),x=-5..10,color=black, Данные прямые линии характеризуют линейную дисперсию скорости. Рассмотрим, как изменяется распределение скорости во времени (анимация) > plots[animate](v1,x=-5..10,t=0..5,frames=30,color=blue,title='Изменение скоросто во времени'); Достаточно сложное изменение скорости связано с тем, что различные части волнового пакета двигаются с разной скоростью, при этом часть волнового пакета даже движется в обратном направлении, т.е. с отрицательной скоростью. Рассмотрим на связанное с этим стандартное "расплывание" волнового пакета. Подставим частные значения параметров в выражение для плотности > rhoS:=simplify(subs(A=(Pi)^(-1/4),h=1,m=1,a=1,q=2,rho)); Посмотрим, как изменяется со временем распределение плотности по координате (анимация) > plots[animate](rhoS,x=-5..15,t=0..5,frames=30,color=red,numpoints=200,title='Изменение распределения плотности во времени'); Полная вероятность при этом не зависит от времени и равна единице > int(rhoS,x=-infinity..infinity); На наш взгляд более наглядно процесс "размазывания" волнового пакета виден на графике логарифма плотности > LogRho:=simplify(log(rhoS)): Достаточно наглядно видно, что различные части волнового пакета двигаются с различными скоростями (даже по направлению), что и приводит к возрастанию неопределенности координаты частицы после процесса измерения в начальный момент времени Проследим, каким образом происходит в уравнении Шредингера переход к классической механике, рассматривая для простоты всего одну частицу во внешнем поле. Подставим в уравнение Шредингера > restart: следующее выражение для волновой функции > phi:=A(x,t)*exp(I/h*S(x,t)); получим > eq1:=expand(subs({psi(x,t)=phi}, eq/phi)); В этом уравнении имеются чисто действительные и чисто мнимые члены. Мы подразумеваем, что S и A - действительны, надо не забыть только сказать об этом и системе > assume(A(x,t),real); assume(S(x,t),real); assume(m>0,h>0); Приравнивая обе части в отдельности к нулю, получим два уравнения > eq_r:=evalc(Re(rhs(eq1)))-evalc(Re(lhs(eq1)))=0; > eq_i:=evalc(Im(rhs(eq1)))-evalc(Im(lhs(eq1)))=0: Пренебрегая в первом уравнении членом, содержащим h , получим > eq_HJ:=coeff(lhs(eq_r),h,0)=0; известное классическое уравнение Гамильтона-Якоби для действия S . Мы видим, что при h->0 классическая механика справедлива с точностью до величин первого, а не нулевого порядка по h включительно. Введем плотность вероятности > eq_rho:=Diff(rho(x,t),t)+(1/m)*Diff(rho(x,t)*Diff(S(x,t),x),x)=0; Проверим это > expand(value(subs({rho(x,t)=A(x,t)^2},eq_rho) + eq_i*2*A(x,t)^2 )); Это уравнение имеет наглядный физический смысл, так как В 20-тые годы, т.е. в самом начале развития квантовой механики наряду с де-бройлевской волновой интерпретацией квантовой механики предпренимались так же попытки поместить квантовую механику в классическое прокрустово ложе. При этом первое полученное нами уравнение предлагалось рассматривать как классическое уравнение Гамильтона-Якоби с добавочным квантовым потенциалом > Uq:=psi->-diff(psi,x$2)/2; отвечающем некоторой квантовой силе > Fq:=psi->-diff(Uq(psi),x); Далее мы посмотрим как выглядит квантовый потенциал и квантовая сила в случае гауссовского волнового пакета, рассмотренного нами ранее > f_g:=(x,t,alpha)->(1/sqrt(alpha*(1+I*t))*sqrt(Pi))*exp(-x^2/alpha^2/(1+I*t)); В начальный момент времени данный пакет расположен в окрестности x=0 со средним моментом p=0 . Посмотрим на изменение данного пакета со временем > with(plots): Warning, the name changecoords has been redefined Здесь синим и красным цветом изображены вешественная и мнимая части данной функции. Посмотрим теперь как изменяется со временем плотность вероятности > animate(Re(f_g(x,t,1))^2+Im(f_g(x,t,1))^2,x=-5..5,t=1..10,frames=10,numpoints=50,color=magenta); Из предыдущих параграфов данной главы нам известно, что гауссовский пакет "расплывается" в координатном представлении, что мы и видели на графиках. При этом в импульсном представлении распределение не изменяется со временем. Напомним, что квантовый потенциал является мерой кинетической энергии, которая изменяется. Это изменение при постоянстве распределения импульсов должно быть связано с изменением фазы волновой функции S(x,t) . Определим модуль гауссовского волнового пакета > A_g:=(x,t,alpha)->simplify(sqrt(evalc(conjugate(f_g(x,t,alpha))*f_g(x,t,alpha)))); Посмотрим что такое квантовая сила > simplify(Fq(A_g(x,t,alpha))); и как она меняется со временем > animate(Fq( (A_g(x,t,1)) ),x=-5..5,t=0..5,frames=5,numpoints=100); Видно, что убывание со временем очень быстрое. Тем не менее, распределение импульсов не меняется и это можно понять изучив изменение фазы волнового пакета > S_g:=(x,t,alpha)->simplify(evalc(arctan(Im(f_g(x,t,alpha))/Re(f_g(x,t,alpha))))); > animate(S_g(x,t,1) ,x=-10..10,t=0..40,frames=40,numpoints=100); Нулевая в начальный момент времени фаза S(x,t) затем расплывается, быстро осциллируя при этом. Для того, чтобы понять изменение кинетической энергии, которая включает в себя и изменение модуля и изменение фазы мы введем плотность кинетической энергии > tau:=(x,t,alpha)->simplify(-evalc(conjugate(f_g(x,t,alpha))*diff(f_g(x,t,alpha),x$2)/2)); Посмотрим как изменяются со временем вещественная и мнимая части кинетической энергии > p1:=animate(Re(tau(x,t,1)),x=-5..5,t=1..5,frames=10,numpoints=50,color=blue): Докажем теперь, что распределение импульсов при этом все-таки не меняется. Для того, чтобы вычислить матричный элемент <p^2> , который пропорционален интегралу от плотности кинетической энергии > evalf(int(tau(x,1,1),x=-infinity..infinity)); > evalf(int(tau(x,2,1),x=-infinity..infinity)); > evalf(int(tau(x,4,1),x=-infinity..infinity)); Аналогично проверим сохранение следующего матричного элемента <p^4> > p4:=(t,alpha)->int(evalc(conjugate(f_g(x,t,alpha))*diff(f_g(x,t,alpha),x$4)),x=-infinity..infinity); > evalf(p4(0,1)); > evalf(p4(2,1)); > evalf(p4(4,1)); Таким образом, мы видим, что локальные изменения рассмотренных нами величин в координатном пространстве несут в себе слишком мало информации. Например, быстрые изменения фазы и амплитуды в координатном пространстве ничего не могут сказать о распределении моментов частиц в импульсном пространстве. Таким образом рассмотрение квантовой механики с помощью добавления квантового потенциала в классическую механику практически невозможно. Именно поэтому данная интерпретация квантовой механики была отвергнута в свое время. В сдедующем параграфе мы рассмотрим более адекватную интерпретацию квантовой механики и перехода ее к квазиклассическому пределу. В то время как в классической механике в каждый данный момент частица обладает определенными координатами и скоростью, в квантовой механике деле обстоит совершенно иным образом. Действительно, получая в любой момент времени информацию о значении координаты и скорости, например, электрона, означало бы наличие определенной траектории, каковой электрон не обладает. Таким образом, в квантовой механике координаты и скорость электрона являются величинами, которые не могут быть одновременно точно измерены, т.е. не могут одновременно иметь определенных значений. Если интервалы значений, в которых могут оказаться координата и импульс частицы обозначить
который был получен Гейзенбергом. Мы видим, что чем с большей точностью известна координата частицы (т.е. чем меньше Далее мы воспользуемся представлением Де Бройля, в котором частица с импульсом Здесь Рассмотрим свойства данного волнового пакета при следующей средней скорости электрона > restart; with(plots): Warning, the name changecoords has been redefined Итак, определим значения константы Планка и массы электрона > hbar:=1.055*10^(-34); Зададим наиболее вероятное значение волнового числа > v:=5*10^5; и интервал неопределенности > delta_k := k[0]/50; Зададим число волн в пакете > N:=40; Предположим, что соответствующие волновые числа > kRange:=4*delta_k: Построим теперь волновой пакет > c1:=hbar/(2*m): c2:=1/(2*delta_k)^2: Посмотрим, как выглядит данный волновой пакет в нулевой момент времени > y:=subs(t=0,psi): Используя анимацию, посмотрим как изменяется данный волновой пакет со временем, т.е. изучим движение волнового пакета > animate(psi,x=0..3*10^(-8),t=0..6*10^(-14),frames=20, Построим теперь другой волновой пакет, в котором соответствующие волновые числа > kRange2:=2*delta_k: В начальный момент времени > y:=subs(t=0,psi2): данный пакет занимает "больше" места в пространстве, т.е. затухает медленнее с ростом значения координаты Используя анимацию, посмотрим на движение данного волнового пакета > animate(psi2,x=0..3*10^(-8),t=0..6*10^(-14),frames=20, На этом мы закончим рассмотрение волновых пакетов в квантовой механике.
Базь А.И., Зельдович Я.Б. Рассеяние, реакции и распады в нерелятивистской квантовой механике, 2-изд. перераб. 1971 год. 544 стр. 4.56 Мб. Книга посвящена вопросам квантовой механики, связанным с ее приложениями к атомным и ядерным процессам. По характеру изложения книга заполняет разрыв между учебниками и оригинальной литературой. Наряду с конкретными задачами рассматриваются современные общие методы, на примере нерелятивистской теории разъясняется понятие перенормировки, важное для теории элементарных частиц. Подробно рассмотрены следующие вопросы. Свойства систем с малой энергией связи. Системы со случайным вырождением — атом водорода, трехмерный гармонический осциллятор. Аналитические свойства волновой функции и матрицы рассеяния. Функция Грина уравнения Шредингера. Точное решение задачи об осцилляторе с переменной частотой под действием внешней силы. Квазиклассические свойства вырожденного ферми-газа. Многомерная квазиклассика, квазиклассическое приближение в нестационарном случае. Свойства нестабильных систем. Свойства многоканальных систем. Пороговые явления. Описание системы из трех тел с помощью уравнений Фаддеева. Изучение теории перенормировок на примере нерелятивистской модели Ли. Рассчитана на научных работников, аспирантов и студентов старших курсов соответствующих специальностей. В первом томе (гл. 1-6) дано описание аппарата углового момента (орбитального н спинового), сложения угловых моментов, их геометрического, алгебраического и теоретико-группового аспектов. Второй том содержит широкий круг применений аппарата углового момента к таким задачам квантовой физики, как эффект Зеемана. спектроскопии многоэлектронных атомов и молекул, мультипольное электромагнитное излучение, угловые распределении в реакциях, структура ядра и т.д. В приложении даны таблицы коэффициентов и операторов Вигнера и Рака, сферических и тензорных функций и редуцированных матричных элементов Собственно теория групп изложена с учетом использования её в физических приложениях, причём наибольшее внимание уделено симметрической группе, группе вращении и важнейшему для приложений разделу - теории представлений. Перед тем как перейти к приложениям, автор кратко излагает основные положения и аппарат квантовой механики, и теорию атомных спектров. Развитая в книге общая теория применяется к атомным спектрам в форме, позволяющей использовать её для более широкого круга проблем - ядерных спектров, теории поля и элементарных частиц и т.п. В связи с этим изложены такие вопросы, как свойства коэффициентов векторной связи и коэффициентов Ракa, а также обращение времени. Книга рассчитана на научных работников и аспирантов физиков, особенно физиков-теоретиков, работающих в области атомной и ядерной спектроскопии, изучения структуры молекул, физики твёрдого тела, а также математиков, интересующихся физическими приложениями теории групп. Для студентов и аспирантов физико-технического и инженерно-физического профиля, а также научно-технических работников, занятых в различных областях современной физики. Цель автора состояла в том, чтобы изложить квантовую механику, с самого начала исходя из физических фактов и экспериментов, связанных с микромиром, а не путем постепенного перехода от классических понятий к квантовым (как это обычно делается). В книге последовательно изложен широкий круг проблем теории квантованных полей и физики элементарных частиц, теории многих тел и квантовой статистики, обсуждаются многие принципиальные вопросы и понятия современней теоретической физики (понятие состояния, понятие частицы, законы сохранения, операции симметрии и т. д.). Книга представляет интерес для широкого круга физиков — как специалистов теоретиков, так и экспериментаторов. Много полезного найдут здесь лекторы вузов, читающие курс квантовой механики. Кинга может являться пособием для аспирантов и студентов старших курсов физических факультетов, специализирующихся по теоретической физике. Для научных работников, занимающихся теорией атома, теоретической спектроскопией, теорией ядра, теорией излучения, а также для аспирантов и студентов старших курсов соответствующих специальностей. Содержание учебного пособия соответствует курсу лекций, который авторы читают в МГТУ имени Баумана Н.Э. Для студентов технических университетов и вузов. Кроме того, мы решили изложить основы вторичного квантования, без знания которого невозможно понять современную теорию излучения. Нам кажется, что это особенно важно для студентов физиков не теоретиков, которые вряд ли будут слушать специальные курсы по квантовой теории поля. Учитывая, что настоящий курс рассчитан, главным образом, на студентов физиков широкого профиля, мы решили остановиться преимущественно на основных вопросах квантовой механики, опуская различные детали, носящие узкоспециальный характер. Книга может рассматриваться как учебное пособие для студентов физических специальностей пединститутов, университетов, а также вузов, где читаются основы квантовой механики. Книга предназначена студентам и аспирантам, углубленно изучающим курс квантовой механики. Кроме того, авторами развит новый подход к исследованию свойств ВКБ-приближения, полезный при дальнейших приложениях метода (в частности, в случае комплексных коэффициентов дифференциального уравнения). Для научных работников в области математики и физики, а также аспирантов и студентов старших курсов соответствующих специальностей. |