Почему нельзя решить уравнение шредингера

Уравнение Шрёдингера

Дуальная корпускулярно-волновая природа квантовых частиц описывается дифференциальным уравнением.

Согласно фольклору, столь распространенному среди физиков, случилось это так: в 1926 году физик-теоретик по имени Эрвин Шрёдингер выступал на научном семинаре в Цюрихском университете. Он рассказывал о странных новых идеях, витающих в воздухе, о том, что объекты микромира часто ведут себя скорее как волны, нежели как частицы. Тут слова попросил пожилой преподаватель и сказал: «Шрёдингер, вы что, не видите, что всё это чушь? Или мы тут все не знаем, что волны — они на то и волны, чтобы описываться волновыми уравнениями?» Шрёдингер воспринял это как личную обиду и задался целью разработать волновое уравнение для описания частиц в рамках квантовой механики — и с блеском справился с этой задачей.

Тут необходимо сделать пояснение. В нашем обыденном мире энергия переносится двумя способами: материей при движении с места на место (например, едущим локомотивом или ветром) — в такой передаче энергии участвуют частицы — или волнами (например, радиоволнами, которые передаются мощными передатчиками и ловятся антеннами наших телевизоров). То есть в макромире, где живём мы с вами, все носители энергии строго подразделяются на два типа — корпускулярные (состоящие из материальных частиц) или волновые. При этом любая волна описывается особым типом уравнений — волновыми уравнениями. Все без исключения волны — волны океана, сейсмические волны горных пород, радиоволны из далеких галактик — описываются однотипными волновыми уравнениями. Это пояснение нужно для того, чтобы было понятно, что если мы хотим представить явления субатомного мира в терминах волн распределения вероятности (см. Квантовая механика), эти волны также должны описываться соответствующим волновым уравнением.

Шрёдингер применил к понятию волн вероятности классическое дифференциальное уравнение волновой функции и получил знаменитое уравнение, носящее его имя. Подобно тому как обычное уравнение волновой функции описывает распространение, например, ряби по поверхности воды, уравнение Шрёдингера описывает распространение волны вероятности нахождения частицы в заданной точке пространства. Пики этой волны (точки максимальной вероятности) показывают, в каком месте пространства скорее всего окажется частица. Хотя уравнение Шрёдингера относится к области высшей математики, оно настолько важно для понимания современной физики, что я его все-таки здесь приведу — в самой простой форме (так называемое «одномерное стационарное уравнение Шрёдингера»). Вышеупомянутая волновая функция распределения вероятности, обозначаемая греческой буквой ψ («пси»), является решением следующего дифференциального уравнения (ничего страшного, если оно вам не понятно; главное — примите на веру, что это уравнение свидетельствует о том, что вероятность ведёт себя как волна):

где x — расстояние, h — постоянная Планка, а m, E и U — соответственно масса, полная энергия и потенциальная энергия частицы.

Картина квантовых событий, которую дает нам уравнение Шрёдингера, заключается в том, что электроны и другие элементарные частицы ведут себя подобно волнам на поверхности океана. С течением времени пик волны (соответствующий месту, в котором скорее всего будет находиться электрон) смещается в пространстве в соответствии с описывающим эту волну уравнением. То есть то, что мы традиционно считали частицей, в квантовом мире ведёт себя во многом подобно волне.

Когда Шрёдингер впервые опубликовал свои результаты, в мире теоретической физики разразилась буря в стакане воды. Дело в том, что практически в то же время появилась работа современника Шрёдингера — Вернера Гейзенберга (см. Принцип неопределенности Гейзенберга), в которой автор выдвинул концепцию «матричной механики», где те же задачи квантовой механики решались в другой, более сложной с математической точки зрения матричной форме. Переполох был вызван тем, что ученые попросту испугались, не противоречат ли друг другу два в равной мере убедительных подхода к описанию микромира. Волнения были напрасны. Сам Шрёдингер в том же году доказал полную эквивалентность двух теорий — то есть из волнового уравнения следует матричное, и наоборот; результаты же получаются идентичными. Сегодня используется в основном версия Шрёдингера (иногда его теорию называют «волновой механикой»), так как его уравнение менее громоздкое и его легче преподавать.

Однако представить себе и принять, что нечто вроде электрона ведёт себя как волна, не так-то просто. В повседневной жизни мы сталкиваемся либо с частицей, либо с волной. Мяч — это частица, звук — это волна, и всё тут. В мире квантовой механики всё не так однозначно. На самом деле — и эксперименты это вскоре показали — в квантовом мире сущности отличаются от привычных нам объектов и обладают другими свойствами. Свет, который мы привыкли считать волной, иногда ведёт себя как частица (которая называется фотон), а частицы вроде электрона и протона могут вести себя как волны (см. Принцип дополнительности).

Эту проблему обычно называют двойственной или дуальной корпускулярно-волновой природой квантовых частиц, причем свойственна она, судя по всему, всем объектам субатомного мира (см. Теорема Белла). Мы должны понять, что в микромире наши обыденные интуитивные представления о том, какие формы может принимать материя и как она себя может вести, просто неприменимы. Сам факт, что мы используем волновое уравнение для описания движения того, что привыкли считать частицами, — яркое тому доказательство. Как уже отмечалось во Введении, в этом нет особого противоречия. Ведь у нас нет никаких веских оснований полагать, будто то, что мы наблюдаем в макромире, должно с точностью воспроизводиться на уровне микромира. И тем не менее дуальная природа элементарных частиц остается одним из самых непонятных и тревожащих аспектов квантовой механики для многих людей, и не будет преувеличением сказать, что все беды начались с Эрвина Шрёдингера.

Гидродинамика Шрёдингера на пальцах

В этой статье в качестве эксперимента я постараюсь максимально доступно рассказать, как работает новый метод расчёта гидродинамики, основанный на решении уравнения Шрёдингера.

Всем привет. В этой статье я хотел бы рассказать о новом методе расчёта гидродинамики, основанном на решении уравнения Шрёдингера вместо уравнений, типично используемых для гидродинамики вроде Навье-Стокса. Сам метод очень подробно и полно раскрыт в диссертации Albert Chern’а, названной «Fluid Dynamics with Incompressible Schrödinger Flow». Однако, статья Chern’а кому-то может показаться написанной на не самом доступном языке, поэтому своей статьёй я бы хотел в первую очередь если не объяснить в деталях, как работает этот метод, то хотя бы объяснить, какими интересными свойствами он обладает, и что же именно скрывается за его математикой. Попутно я кратко расскажу о том, как устроены классические методы расчёта гидродинаимики и как новый подход от них отличается. В качестве эксперимента я бы хотел попробовать написать статью так, чтобы каждый, кто отдалённо интересуется программированием физики, нашёл в ней что-то интересное, понятное, и новое для себя — от начинающего программиста до бывалых расчётчиков.

Вступление

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

Следующее очень важное следствие эквивалентности уравнения Шрёдингера и Навье-Стокса — это что решение одного из них эквивалентно решению другого. Так вот уравнение Навье-Стокса — нелиненое, его очень неудобно и неэффективно в общем случае решать, в то время как уравнение Шрёдингера — линеное и его решать гораздо проще. Чтобы составить представление, насколько же неудобным по сей день считается уравнение Навье-Стокса, могу сообщить, что существует целый международный фонд грантов для исследователей, которым хоть какую-то базу под них подстроит, так как(цитата):

Even basic properties of the solutions to Navier–Stokes have never been proven.

Уравнение Шрёдингера же, хоть и описывает мутную квантовую физику, поддаётся решению гораздо легче и эффективнее. Короче, я могу очень долго гудеть про то, как это невероятно и офигенно, но давайте уже перейдём к чему-то более конкретному.

Решение классической гидродинамики на пальцах

Что вообще такое — уравнение гидродинамики? Что такое уравнение Навье-Стокса и как его понять? С ответом на этот вопрос гораздо лучше меня справились миллионы авторов статей по этому делу, например, классическая статья от нвидии, по которой многие начинали: https://developer.download.nvidia.com/books/HTML/gpugems/gpugems_ch38.html Однако, я попробую написать очень сжато и на пальцах, что это всё значит и что с этим обычно делают.
Уравнение Навье-Стокса описывает закон, которому обязана подчиняться скорость каждой точки пространства, заполненного равномерной несжимаемой жидкостью. Представьте себе, например, бассейн с водой, в котором выделили некоторый куб, достаточно далеко от стенок, поверхности и дна, в котором нет ничего кроме воды. Вода в нём может как угодно течь, но не может ни образовывать пузырей, ни с чем-то сталкиваться (мы для простоты опустим эти эффекты). Тогда само уравнение Навье-Стокса описывает закон, которому будет подчиняться скорость каждой точки воды в этом кубе:

\(\frac<\partial \vec u><\partial t>=-(\vec u \cdot \vec \nabla)\vec u-\frac<1><\rho>\vec \nabla p + \nu <\vec \nabla>^2 \vec u + \vec F\)
\(\vec \nabla \vec u = 0\)

прежде чем вообще смотреть на это уравнение, предлагаю сразу из него выбросить ненужное — то, что нам всё равно не пригодится для понимания и только место занимает. Это член, отвечающий за диффузию \(\nu <\vec \nabla>^2 \vec u\) (у идеальной жидкости один фиг диффузии нет), и за внешнюю силу \(\vec F\) (так как мы обойдёмся без неё). Остаётся система:
\(\frac<\partial \vec u><\partial t>=-(\vec u \cdot \vec \nabla)\vec u-\frac<1><\rho>\vec \nabla p\)
\(\vec \nabla \vec u = 0\)
Здесь перевёрнутый треугольник называется оператором Набла, который обозначает дифференцирование. Причём смысл этого оператора меняется в зависимости от того, где именно он стоит (например, перед вектором или скаляром). Я постараюсь объяснить смысл каждого его вхождения по порядку. На пальцах смысл всей формулы в следующем. \(\vec u(\vec x)\) — это значение скорости жидкости, которое определяется в каждой точке пространства \(\vec x\) . Уравнение описывает закономерности, которым обязана подчиняться эта величина, если она описывает поведение несжимаемой жидкости. Работает хоть для двумерного, хоть для трёхмерного случая. В левой части первого уравнения стоит \(\frac<\partial \vec u(\vec x)><\partial t>\) — это величина называется производной по времени и показывает, как быстро и куда(это вектор) изменится скорость в точке \(\vec x\) в момент времени \(t\) .

Нулевой вектор производной по времени обозначает, что скорость в этой точке сейчас не меняется, а, например, вектор (10, 0)[м/c 2 ] обозначает, что за следующую секунду скорость вырастет на 10[м/с] по оси x(если сама производная не поменяется).

Слагаемое вида \(-(\vec v \cdot \vec \nabla)\vec u\) называется адвекцией и говорит, что поле скоростей \(\vec u\) в этой точке утекает в направлении \(\vec v\) . В нашем же случае \(\vec u = \vec v\) , то есть поле скоростей сносит само себя. Это, кстати, и называется нелинейностью и из-за этого возникает миллион проблем при решении этого уравнения.

В принципе, смысл этого члена достаточно интуитивно можно представить именно как утекание каждой точки воды по вектору её скорости. Однако, в общем случае производная векторного поля \(\vec u\) по направлению \(\vec v\) обозначается как \((\vec v \cdot \vec \nabla)\vec u\) и обозначает, как меняется функция \(\vec u\) в направлении \(\vec v\) для этой точки.

Слагаемое же \(-\frac<1><\rho>\vec \nabla p\) является ускорением, которое получает жидкость в точке из-за градиента давления.

Оператор \(\vec \nabla\) , действующий на скалярное поле(например, давление), называется градиентом. Если слева от некоторой точки давление больше, чем справа, то градиент в ней будет направлен вправо и будет увлекать за собой жидкость в этом направлении. Например, ветер всегда дует в направлении, обратном градиенту давления воздуха (отсюда и минус). Электрический ток течёт в направлении градиента электрического потенциала:
\(E=\vec \nabla \phi\)

Второе уравнение \(\vec \nabla \vec u = 0\) называется уравнение непрерывности, а оператор \(\vec \nabla\) здесь действует на вектор и называется дивергенцией.

Оператор дифференцирования, действующий на вектор, называется дивергенцией. Дивергенция, равная нулю, говорит, что для каждого маленького кубика сколько в него жидкости втекает, столько и вытекает. А так как любой объём можно разбить на маленькие кубики, то свойство будет справедливо и для объёма любой формы. Это свойство называют также условием несжимаемости, так как если бы в какой-то объём втекало больше жидкости, чем вытекало, это бы означало, что жидкость в объёме накапливается, сжимаясь. Другой случай применения дивергенции, который может помочь её представить — это теорема Гаусса:
\(\vec \nabla E=\rho\)
Эта теорема говорит, что напряжённость электрического поля, которая «вытекает» из некоторого объёма, всегда вызвана электрическим зарядом плотности \(\rho\) внутри этого объёма. Если в объёме заряда нет, то и дивергенция нулевая.

То есть, одним предложением уравнение Навье-Стокса можно описать так: темп изменения скорости определяется течением и градиентом давления, но жидкость при этом не может сжиматься.

Классическое решения уравнения Навье-Стокса

Посмотрим теперь, как это уравнение можно программно решить. Для этого можно использовать подход, который называется расщеплением — разбить сложный физический процесс, состоящий из нескольких элементарных, на отдельные чередующиеся стадии и считать, что на каждой стадии работает только один элементарный процесс, а остальные выключены. Как ни странно, можно доказать (см. статью выше), что это — на самом деле математически обоснованная стратегия. Поэтому будем считать, что состояние скоростей для каждой точки в текущий момент времени \(\vec u(\vec x, t)\) нам известно. А для расчёта состояния в следующий момент времени \(t+dt\) , разобьём сложный процесс гидродинамической эволюции на простые стадии:
1) снесём поле скоростей по течению. это может немного «сжать» жидкость.
2) найдём такое давление, чтобы жидкость «расжалась».
Первый шаг называется адвекцией, второй — проекцией.

Адвекция

Адвекция, или течение, можно приближённо посчитать достаточно легко — если известно, что в точке \(\vec x\) , в момент времени \(t\) скорость равна \(\vec u(\vec x, t)\) , то в момент времени \(t+dt\) скорость в неё притечёт жидкость из точки \(\vec x — \vec u(x, t)\cdot dt\) .
\(\vec u^*(\vec x, t+dt)=\vec u(\vec x — \vec u(\vec x, t), t)\)
То есть мы получили промежуточное значение скорости, котороже уже утекло по течению, но теперь в нём нарушено условие непрерывности.

Это особенно удобно программируется на GPU, так как это можно посчитать, если хранить скорость в текстуре и её обновлять, просто читая тексели со смещением \(- \vec u(x, t)\cdot dt\) и используя стандартную аппаратную линейную интерполяцию.

Проекция

Проекция берёт скорость, для которой нарушено условие непрерывности \(\vec u^*\) и ищет такое давление, которое её «выправит» до нормальной скорости \(\vec u\) . Умные мужики доказали, что такое поле можно найти единственным образом и оно всегда будет градиентом некоторого скалярного поля (давления, в нашем случае):
\(\vec u(\vec x, t+dt)=\vec u^*(\vec x, t+dt) + \vec \nabla p\)
Помножим обе стороны этого равенства на оператор дифференцирования:
\(\vec \nabla \vec u(\vec x, t+dt)=\vec \nabla \vec u^*(\vec x, t+dt) + \vec \nabla^2 p\)
«ПОГОДИ-КА СУСЕЛ, ЭТО ЕЩЁ ЧТО» — можете меня спросить вы. Всё по порядку, но на самом деле отсюда для общего понимания достаточно знать, что если \(\vec u^*(\vec x)\) известно(а оно известно), то отсюда можно найти давление \(p(\vec x)\) . Если вспомнить, что в нашем случае дивергенция скорости равна нулю, то остаётся вот такое выражение.
\(\vec \nabla^2 p=-\vec \nabla u^*\)

В правой части этого равенства стоит дивергенция скорости, которую можно легко приблизительно посчитать, если известна скорость \(\vec u^*\) (а она известна). В левой части стоит штука, которая называется лапласианом давления.

Лапласиан — это оператор дифференцирования (ещё называется оператор набла) в квадрате, то есть применённый дважды к скалярному полю. Первый раз применяем оператор дифференцирования — получаем градиент. Второй раз — получаем дивергенцию. Таким образом оператор лапласа — это дивергенция градиента скалярного поля. Его можно представить как изменение потока скорости через маленький кубик, которое будет вызвано давлением в точке. Ещё одна аналогия — как поменяется дивергенция электрического поля в объёмчике, если в него положить заряд плотностью \(\rho\) (опять же, теорема Гаусса):
\(\vec nabla \vec E = \rho\) , \(\vec nabla \phi=\vec E\) => \(\vec nabla^2 \vec phi = \rho\)

Уравнение вида «лапласиан чего-то неизвестного равен чему-то известному» называется уравнением Пуассона. Что бы это ни значило, существует стандартный итеративный алгоритм, который позволяет его решить, то есть найти такое давление, чтобы его лапласиан был равен чему угодно. «Что угодно» мы знаем — это дивергенция промежуточной скорости, поэтому считаем по ней давление. Далее для давления считаем градиент и вычитаем результат из промежуточной скорости, чтобы получить окончательную скорость для следующего шага по времени:
\(\vec u=\vec u^* + \vec \nabla p\)

Шаги адвекции и проекции повторяем до посинения, рассчитывая всё дальше и дальше эволюцию поля течений по времени. Для визуализации можно, например, напускать частиц, которые могу сноситься этим полем скоростей. Результат выглядит так:

Важно понять, что в этом видосе, равно как и во всех остальных гифках этой статьи, жидкость на самом деле находится в большом кубе (границы которого не показаны), а не только там, где видны частицы. Частицы только уносятся полем скоростей, как, например, частицы дыма уносятся полем скоростей воздухе. Сами частицы никакой роли в физике процесса не играют и только позволяют относительно наглядно его продемонстрировать. Частицы обычно добавляются заранее туда, где ожидаются какие-то интересные турбулентности.

Важные особенности классического подхода

«Всё здорово, сусел, но в названии статьи ты написал что-то там про Шрёдингера! Он вообще где? Зачем нам это всё?» — спросите вы. Вопрос резонный. Но всю крутость подхода со Шрёдингером можно осознать, только если иметь представление о слабых сторонах классического солвера, который мы рассмотрели в предыдущей главе. В чём же они заключаются? Давайте об этом поговорим.

Основа любого расчётного метода — это то, как в нём представлены моделируемые данные. В рассмотренном нами подходе мы храним значение скорости для каждой точки. Например, в текселях двумерной или трёхмерной текстуры. Этот способ здорово работает, если требуется описать ровное поле течений, в котором нет особенностей (так называются завихрения и разные другие неоднородности). Неоднородностей обычно нет в вязких жидкостях вроде мёда или майонеза, поэтому метод очень здорово подходит, чтобы моделировать майонез. Но более текучие среды (например, вода, воздух и дым) отличаются тем, что в них существенную роль играют злополучные турбулентные течения — мелкие завихрения, имеющие очень сложную и нерегулярную структуру, даже образующие фракталы, которые очень неудобно описывать просто их значениями в каждой точке текстуры/массива. Если попытаться их моделировать, то все мелкие особенности быстро смазываются и расплываются, что соответствует поведению вязкой жидкости. Такое поведение называется численной вязкостью — это вязкость жидкости, которая появляется не потому что она является частью уравнения, которое мы решаем, а это паразитная вязкость, всплывающая как паразитное следствие нашего метода решения. Более того, напомню, что первое, что мы сделали, не успев взглянуть на уравнение Навье-Стокса — выкинули из него вязкость, так в ней недостатка точно не будет.

А вот избавиться от вязкости гораздо труднее, чем случайно её посчитать. Один из способов — это измельчать расчётную сетку. Чтобы таким методом получить что-то хоть как-то похожее на дым, понадобится сетка минимум 1024x1024x1024, то есть как минимум гигабайт памяти, если хранить по 1 байту на узел. А хранить захочется как минимум трёхкомпонентную скорость, то есть, скорее всего, 32 гигабайта в сумме. Это не только не разумно с точки зрения затрат памяти, это ещё и очень медленно. Другой способ — это представлять скорость не её направлением в каждой точке, а как сумму маленьких элементарных вихрей. Этот метод называется также методом дискретных вихрей. В нём вообще всё не так просто с процессами порождения новых вихрей и удаления старых, с поддержанием нужной плотности (так как вихри друг друга уносят, как частицы) и ещё миллион проблем, можете сами почитать, если интересно. Другой подход основан на том, что в реальных течениях вихри имеют свойство образовывать вращающиеся нити. Представьте медленно движущийся жгут, вокруг которого быстро вращается жидкость. Если такой жгут замыкается в кольцо, получается тороидальный вихрь, образующий знакомое кольцо дыма:

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

Однако, у классического метода есть одно очень важное положительное свойство — в нём вообще нет параметров. Обратите внимание, что для расчёта используется только скорость и больше вообще ничего — ни вязкости, ни даже плотности. В уравнении Навье-Стокса без вязкости есть плотность, но её можно «спрятать» в нормировку давления, поэтому можно сказать, что в исходном уравнении параметров также нет. Забегая вперёд, замему, что в солвере на уравнении Шрёдинге будет параметр. Загадочный.

На следующей странице мы рассмотрим, как же применить уравнения Шрёдингера, чтобы смоделировать тот же самый процесс, и какой в этом профит. Будет много картинок.

Уравнение Шредингера

Вы будете перенаправлены на Автор24

Предпосылки вывода уравнения Шредингера

Основная идея волновой механики заключается в том, что для таких малых тел, как электрон, нельзя с определенностью сказать, где оно находится в данное время и куда направляется. Можно установить только относительную вероятность его нахождения в том или ином месте и наличие определенного количества движения в определенный момент времени.

В соответствии с волновой механикой какая-либо система – атом, молекула, электрон и т.д. – описывается функцией состояния или волновой функцией, обозначаемой $\psi$ («пси»), которая является функцией координат всех частиц, образующих эту систему. Следовательно, величина $\psi$ зависит только от положения всех частиц в пространстве.

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

где $\hbar$ – постоянная Планка ($6,626\cdot 1034 Дж\cdot с$ или $6,626\cdot 10-27 эрг\cdot с$), а $p$ – количество движения (импульс) электрона в пучке, т.е. его масса, умноженная на его скорость.

Физическое подтверждение волновой природы электрона было продемонстрировано в 1927 – 1928 гг. Дейвиссоном, Джермером и Томсоном, которые показали, что пучок электронов может испытывать дифракцию на подходящей решетке (атомы в кристалле золота), аналогичную дифракции пучка света.

Рисунок 1. Дифракция пучка электронов

На преграду с двумя узкими щелями направлен параллельный пучок моноэнергетических (т.е. обладающих одинаковой кинетической энергией) электронов (рис. 1. а). За преградой находится фотопластина $Фn$. При закрытии щели номер $2$ и экспонировании в течение времени $t$ почернение на проявленной фотопластине будет характеризоваться кривой $1$ (рис. 1. б). При закрытии щели номер $1$, соответственно, почернение на фотопластине будет соответствовать кривой $2$. Однако в случае, когда открыты обе щели картина почернения фотопластины (рис. 1. в) отнюдь не эквивалентна наложению двух первых картин. Зато она аналогична картине, получающейся при интерференции двух когерентных световых волн.

Готовые работы на аналогичную тему

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

Общий вид уравнения Шредингера

Рисунок 2. Эрвин Шрёдингер (1887 — 1961)

Волновое уравнение, применяемое для расчета стационарных состояний системы, можно записать в символическом виде:

где $H$ представляет собой определенный способ выражения общей энергии системы, а $E$ – числовое значение этой энергии. Для всех систем, которые обычно интересуют химиков, общая энергия представляет собой сумму кинетической энергии $Т$ и потенциальной энергии $V$:

Это соотношение было широко использовано физиком-теоретиком Гамильтоном, поэтому $H$ часто называют функцией Гамильтона, а $\mathcal H$ гамильтонианом системы.

Уравнение Шредингера на примере атома водорода

Рассмотрим модель атома водорода, предложенную Бором. Для простоты предположим, что тяжелое ядро закреплено (оно почти, но не совершенно неподвижно, когда электрон движется вокруг него). Тогда полная кинетическая энергия $Т$ системы представляет собой просто кинетическую энергию электрона

где $m$ – масса электрона и $\nu$ – его скорость. Потенциальная энергия системы есть просто энергия, возникающая вследствие электростатического взаимодействия (гравитационные силы приблизительно в $10^<18>$ раз меньше), и ее можно выразить как

где $e$ — заряд электрона, $r$ — радиус орбиты, знак минус появляется вследствие того, что заряд одной из частиц положителен $(+)$, а другой отрицателен $(-)$. Поэтому для атома водорода функция Гамильтона в классической (т.е. доквантовомеханической) физике равна:

Если использовать понятие количества движения электрона $p=m\nu$, данное уравнение запишется в следующем виде:

Теперь для перехода от классического описания этой или какой-либо другой системы к описанию при помощи волновой механики, необходимо взять функцию Гамильтона (уравнение 6) и произвести в ней определенные замены: в функции Гамильтона количество движения следует заменить выражением

Таким образом, гамильтониан для атома водорода в его квантовомеханической форме $<\mathcal H>$ следует записать в виде

Если теперь это выражение гамильтониана подставить в общее волновое уравнение (уравнение 1), то получим:

Это и есть волновое уравнение для атома водорода. Из уравнения 9 следует, что нужно вторые производные функции $\psi $ сложить и умножить на $-<<\hbar >^2>/<8<\pi >^2m>$, затем к этому добавить $\left(-/\right)\psi $, тогда получим величину, тождественную Е$\psi $. Если найдена функция $\psi $, то говорят, что она является решением волнового уравнения, и ее называют волновой функцией. Вообще, может быть несколько различных функций $\psi_1$, $\psi_2$, . , $\psi_n$, которые являются решениями уравнения 9, причем каждой соответствует свое значение энергии $Е_1$, $Е_2$, . , $Е_n$.


источники:

http://gamedev.ru/code/articles/shrodinger_hydrodynamics

http://spravochnick.ru/himiya/atomnye_i_molekulyarnye_orbitali/uravnenie_shredingera/