Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Таммовские состояния

.pdf
Скачиваний:
59
Добавлен:
27.03.2016
Размер:
328.59 Кб
Скачать

Таммовские состояния.

В 1932 году была опубликована статья И. Е. Тамма «О возможной связи электронов на поверхности кристалла». В ней было показано, что граница кристалла при наличии на ней возмущений может служить источником особых состояний электронов. Они в основном локализованы вблизи этой границы, а уровни энергий попадают в запрещенную зону, т. е. выше потолка зоны проводимости либо ниже ее дна. В настоящее время такие состояния называются «таммовскими», иногда говорят об «уровнях Тамма». Впоследствии эти состояния были обнаружены экспериментально на многих поверхностях металлов и полупроводников. За это достижение шведская академия наук присудила И. Е. Тамму Нобелевскую премию. Тамм в своей работе рассматривал электроны, сильно связанные с ионным остовом атомов, т.е. достаточно локализованные. В рамках такой модели возможны перескоки электронов лишь между соседними атомами, которые происходят из-за перекрытия волновых функций соседних атомов. Чуть позже Шокли решил эту же задачу для почти свободных коллективизированных электронов и также продемонстрировал возможность появления уровней, лежащих вне зоны объемных уровней металла. По этой причине уровни, отщепляющиеся от более узких d-зон, называют уровнями Тамма, а уровни, отщепляющиеся от широких s,p-зон, – уровнями Шокли. Заметим, впрочем, что в экспериментах часто бывает трудно классифицировать наблюдаемые поверхностные уровни в соответствии с этой простой схемой.

Задача настоящей лекции состоит в том, чтобы продемонстрировать возникновение таммовских состояний на простой модели металла, в качестве которой мы рассмотрим одномерную цепочку из N атомов, и дать качественную интерпретацию полученного результата. В соответствии со сказанным выше, будем рассматривать такой металл в модели сильной связи, который в данном случае аналогичен методу, широко применяемому в квантовой химии, – методу линейной комбинации атомных орбиталей (ЛКАО). Далее, здесь мы не будем учитывать спин электронов, ограничиваясь, таким образом, рассмотрением немагнитных металлов. В ферромагнитных металлах, таких, как Fe, Co, Ni, таммовские уровни тоже наблюдаются, но там они расщепляются по спину, что усложняет картину энергетической структуры этих металлов. Нам же надо лишь продемонстрировать механизм появления таммовских состояний.

Для рассматриваемого одномерного твердого тела уравнение Шредингера в дираковских обозначениях имеет стандартный вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

E

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть

 

n

 

 

волновые функции атомов, составляющих цепочку из N атомов.

Разложим

 

 

 

волновую функцию всей цепочки

 

 

по волновым функциям отдельных атомов

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cn

 

n

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ч

 

 

 

 

 

 

 

n

отвечает состоянию электрона в поле атома с номером n, cn – коэффициент ЛКАО.

 

Подставим

 

 

 

в виде (2) в уравнение Шредингера (1) и умножим его слева на

m

 

 

, тогда оно

 

 

 

принимает следующий вид

cn (Hmn ESmn ) 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

где

Hmn m

 

H

 

n ,

Smn m

 

n

. Далее воспользуемся упрощающим предположением:

 

 

 

будем считать, что волновые функции, отвечающие различным атомам, ортогональны

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

n Smn mn

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1

загромождаются, но ничего нового, в свете поставленной нами задачи, не обнаруживается. По этой причине ниже мы будем пользоваться соотношением (4).

Предположим также, что матричные элементы гамильтониана отличны от нуля только для m n, n 1:

 

 

 

 

 

 

,

m n

 

 

m

 

n

Hmn

 

,

m n 1

 

 

 

 

 

H

 

 

(5)

 

 

 

 

 

 

m n, n 1

 

 

 

 

 

 

0,

 

 

Здесь – это энергия одного электрона в изолированном атоме,

 

– амплитуда перескока

электрона с одного атома на соседний атом. Поскольку мы рассматриваем сильно связанные электроны, то перескоки электрона с одного атома на атомы, следующие за ближайшим атомом, запрещены, – их волновые функции не перекрываются. Как следствие, матричные элементы гамильтониана Hmn в (5) с m n, n 1 положены равными нулю.

В рамках сделанных предположений (4) и (5) уравнение (3) записывается в виде

 

N

 

 

 

 

 

cn 0

 

 

 

m,n 1 m,n m,n 1

E m,n

 

 

(6.1)

n 1

 

 

 

 

 

 

 

 

 

Это уравнение можно также записать в виде

 

 

 

 

 

 

 

 

 

 

 

Ac 0

 

 

 

 

 

(6.2)

Здесь c (c1, c2 ,...., cN ) – вектор с

компонентами,

являющимися коэффициентами ЛКАО, а

матрица А определяется следующим выражением

 

 

 

 

 

 

 

E

 

0

...

0

0

0

 

 

 

 

 

E

 

0

...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

E

 

...

...

0

 

 

A

 

0

0

 

...

 

...

...

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

...

...

 

E

 

0

 

 

 

 

0

0

...

0

 

E

 

 

 

 

 

 

 

 

 

0

0

0

...

0

 

 

 

 

 

 

E

 

Это так называемая симметричная трехдиагональная матрица Якоби размерности (N x N). Еще ее называют осцилляционной матрицей.

Уравнения (6.1), (6.2) – это уравнение на коэффициента ЛКАО cn . Понятно, что оно всегда имеет тривиальное решение c (0, 0, 0,...., 0) . Нетривиальное решение c 0 , как это видно из (6.2), возможно лишь при условии det A 0 , которое есть уравнение на собственные значения энергии En . Для того, чтобы не загромождать текст формулами, вывод определителя

трехдиагональной матрицы Якоби А приведен в приложении. Здесь же мы воспользуемся полученными там результатами. Для этого разделим в матрице все элементы на 2 и введем

обозначение

E

. Последнее – это обезразмеренная энергия цепочки атомов металла,

2

 

 

отсчитываемая от одноэлектронного уровня отдельного атома . Оказывается, что выражение для определителя det A существенно зависит от того, в каком интервале находятся значения параметра :

1) 1 1. Это условие означает, что мы ищем решение уравнения det A 0 в объемной зоне металла, т.е. 2 E 2 . В этом случае уравнение det A 0 можно записать как

det A

sin (N 1)

0,

arccos ,

(8)

2N sin

 

 

 

 

Решение уравнения (8) относительно имеет вид

2

 

n

 

,

n 1, 2,...., N.

(9)

N

 

 

1

 

 

 

 

Тогда для уровней энергии En имеем окончательно следующее выражение

 

E 2 cos

n

,

n 1, 2,..., N

(10)

 

n

 

 

N 1

 

 

 

 

 

 

 

 

При указанных значениях индекса n , нумерующего энергетические уровни, удовлетворяется сделанное выше предположение 1 1. Это означает, что все N уровней цепочки из N атомов находятся в зоне, то есть 2 En 2 . А поскольку всего уровней должно быть

N, то заниматься поиском уровней вне зоны, т.е. выше ее потолка ( E 2 ,

1) или

ниже дна ( E 2 ,

1) бессмысленно – и так понятно, что их там нет. Таким образом,

сам факт наличия поверхности, то есть обрыв цепочки при n 1,

n N еще не приводит к

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

Теперь учтем наличие возмущения на поверхности. Для этого положим, что уровень энергии электрона в атоме на краю цепочки n 1 отличается от уровня энергии внутренних атомов в цепочке. Основанием для введения такого возмущения является следующие физические факторы. В приповерхностном атомном слое происходит перераспределение электронной плотности – имеют место так называемые фриделевские осцилляции. Это приводит, в частности, к тому, что в приповерхностной области атомы оказываются заряженными, причем их заряд меняется от атома к атому. В силу этого и взаимодействие этих атомов отличается от того, как они взаимодействуют друг с другом в объеме металла. Как следствие, имеет место явление межплоскостной релаксации, когда параметр решетки, то есть расстояние между соседними атомными плоскостями в приповерхностной области отличается от своего объемного значения, что в свою очередь также приводит к дополнительному перераспределению электронной плотности между атомами вблизи поверхности. В конечном счете, атомы поверхности оказываются в эффективном поле, которое, в частности, формирует так называемый дипольный барьер, наличие которого приводит к повышению работы выхода металла. Здесь мы учтем влияние этих физических факторов в рамках простейшего приближения, полагая, что электронный уровень только лишь поверхностного атома отличается от уровней энергии электронов всех остальных атомов в цепочке. Как будет показано ниже, если это поверхностное возмущение превышает некоторое критическое значение, то появляется таммовский уровень. В этом случае говорят, что имеется собственное таммовское состояние. Если же это поверхностное возмущение по тем или иным причинам недостаточно велико для появления таммовского состояния, то его появления можно добиться за счет адсорбции инородных атомов на поверхности, то есть за счет дополнительного взаимодействия атомов поверхности с адсорбатом. В этом случае говорят о несобственном таммовском состоянии. Итак, полагая, что у первого атома в цепочке энергия электронного уровня , перепишем матрицу А (7) в виде.

 

E

 

0

...

0

0

0

 

 

 

 

 

E

 

0 ...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

E

...

...

0

 

 

A

 

0

0

 

...

 

...

...

 

(11)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

...

...

 

E

 

0

 

 

 

 

0

0

...

0

 

E

 

 

 

 

 

 

 

 

 

0

0

0

...

0

 

 

 

 

 

 

E

 

3

В принципе, аналогичное возмущение должно присутствовать и на другом конце цепочки, то есть в нижнем правом углу матрицы А (11) мы также должны поставить E вместо E . Однако в дальнейшем мы перейдем к пределу N , и влиянием второй поверхности на электронную структуру вблизи первой поверхности можно пренебречь.

Теперь проанализируем уравнение на уровни энергии det A 0 при наличии поверхностного возмущения. Для этого в уравнении det A 0 , где матрица А определяется

формулой (11) разделим все матричные элементы на 2

и обозначим через поверхностное

возмущение как

 

. Обезразмеренная энергия как и раньше обозначена как

E

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Тогда уравнение на собственные значения энергии принимает следующий вид

 

 

 

 

 

 

 

 

1

0

...

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

... ...

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

0

(12)

 

det

0

 

 

0

 

...

 

...

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

0

 

 

... ...

 

1

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

0

 

 

0

...

0

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0

 

 

0

0

...

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

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

 

 

 

 

1

 

0 или, что то же

dN

 

 

(13)

 

dN 1

 

dN 2

 

dN 1

4

2

 

 

2

 

 

 

 

 

 

Теперь нам предстоит проанализировать это уравнение для трех различных случаев: (А)

1 1, (Б) 1 , (В) 1.

А) 1 1. Это условие означает, что мы ищем уровни энергии в зоне 2 En 2 , которую уже рассматривали раньше при отсутствии поверхностного возмущения. В этом случае выражение для определителя dN имеет вид (см. Приложение)

d

 

 

sin (N 1)

(14)

N

2N sin

 

 

 

 

и уравнение (13) принимает вид

 

 

 

 

 

 

 

sin (N 1)

 

(15)

 

 

 

 

 

sin N

 

График правой части этого уравнения представлен на Рис. 1. Простейший анализ показывает,

что если поверхностное возмущение мало, то есть N 1 N 1 , то уравнение (15) имеет

N N

N корней, то есть имеется N пересечений графика правой части уравнения с горизонтальной прямой. Это означает, что все N энергетических уровней находятся в зоне, и таммовских состояний нет. Ситуация меняется при увеличении поверхностного возмущения: как видно из

рисунка, при

 

N 1

1

исчезает самый левый корень, и, как следствие, корней в зоне

N

 

 

 

 

 

остается уже N 1. Естественно, возникает вопрос, где он находится, поскольку в системе из

4

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

есть N 1 1.

N

Рис.1

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

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

 

Б) 1 . Это условие означает, что мы ищем уровни энергии выше потолка зоны,

то есть

2 E . В этом случае определитель dN имеет вид (см. Приложение)

 

d

 

 

sh (N 1)

,

ch , 0

(16)

N

 

 

 

 

2N sh

 

 

 

и уравнение (14) принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

sh (N 1)

 

(17)

 

 

 

 

 

 

 

 

 

sh N

 

5

Рис. 2

Графическое решение этого уравнения представлено на Рис.2. Как следует из этого рисунка, уравнение (17) имеет решение, только если поверхностное возмущение превышает некоторое

критическое значение

 

N 1

. Это как раз то условие, при котором исчезает самый левый

N

 

 

 

 

корень уравнения (15), то есть это то условие, при котором из зоны исчезает один из уровней. Как мы видим, он появляется выше потолка зоны. Это и есть таммовский уровень.

В) 1. Это условие означает,

что мы ищем уровни энергии

ниже дна зоны, то есть

E 2 . В этом случае определитель dN имеет вид

 

 

d

 

(1)N

sh (N 1)

, ch ,

0

(18)

N

 

 

 

 

 

2N sh

 

 

и уравнение (13) принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

sh (N 1)

 

 

(19)

 

 

 

 

 

 

 

 

 

 

 

 

sh N

 

 

6

Рис. 3

Графическое решение этого уравнения представлено на Рис.3. Как следует из этого рисунка, уравнение (19) имеет решение, только если абсолютное значение поверхностного возмущения

превышает некоторое критическое значение, точнее

 

N 1

 

. Это как раз то условие, при

N

 

 

 

котором исчезает самый правый корень уравнения (15), то есть это то условие, при котором из зоны исчезает один из уровней в зоне. Как мы видим, он появляется ниже дна зоны. Это и есть таммовский уровень

Важно отметить, что при достаточно большой величине числа атомов в цепочке (N 1)

уровни в зоне расположены квазинепрерывно – расстояние между ними порядка 4 . В то же

N

время, как это следует, например, из рассмотрения случая (Б), расстояние от уровня, отщепившегося от потолка зоны до ближайшего к нему уровня на потолке зоны вообще не зависит от N . Это же следует из рассмотрения случая (В), где расстояние от уровня, отщепившегося от дна зоны до ближайшего к нему уровня на дне зоны также не зависит от N . Следовательно, таммовский уровень не есть один из зонных уровней, он действительно отщепляется от потолка либо от дна зоны и находится от нее на достаточно большом расстоянии.

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

7

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

Продемонстрировав, таким образом, механизм возникновения таммовского уровня, вернемся к задаче нахождения волновой функции электронов, отвечающей этому уровню. Эта задача, как следует из (2) и (6.1), сводится к определению коэффициентов ЛКАО из системы рекуррентных уравнений (6.1), которую в более развернутом виде можно записать следующим образом

 

 

 

 

1

 

0

 

c1

 

c2

2

 

 

2

 

 

 

 

 

 

1

c

c

 

1

c

 

0

(20)

 

 

 

 

 

 

 

2

1

2

 

 

2

 

3

 

 

 

 

 

 

………..

 

 

 

 

 

 

 

1

c

 

c

1

c

 

0

 

 

 

 

 

 

 

2 n 1

 

n

2

 

n 1

 

 

 

 

Здесь мы рассмотрим частный предельный случай бесконечного числа атомов в цепочке N . Тогда уравнение (17) принимает вид e , из чего следует, что условие появления таммовского уровня выше потолка зоны металла записывается как 1.

Как видно из (20), рекуррентные уравнения на коэффициенты cn , начиная со второго, совпадают друг с другом по форме. Будем искать решение в виде cn cn . Подстановка этого выражения в третье уравнение в (20) приводит к квадратному уравнению на параметр c

1

c2

c

1

0 ,

(21)

2

2

 

 

 

 

Решение уравнения (21) с учетом выражения для ch имеет следующий вид c e . Тогда в силу линейности рекуррентного уравнения ответ для cn можно записать как

c Aen Be n .

(22)

n

 

Из условия нормировки волновой функции (2) на единицу следует следующее требование на коэффициенты ЛКАО

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

cn

 

1 ,

 

 

(23)

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

которое может быть выполнено только при выполнении условия

A 0 . Тогда, подставляя

выражение для c

n

в виде

c Be n в условие (23) и суммируя получающуюся геометрическую

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

прогрессию, получаем выражение для коэффициента B

e2 1 .

Окончательное выражение

для коэффициентов ЛКАО имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

e2 1e n

 

 

(24)

 

 

 

n

 

 

 

 

 

 

 

 

 

8

Подстановка этого выражения в первое уравнение (20) дает e , то есть уравнение (17) в

рассматриваемом здесь частном предельном случае N . Поставляя это

выражение в

формулу для волновой функции (2), получим для нее окончательное выражение

 

 

 

 

 

 

 

 

 

e2 1e n

 

n

(25)

 

 

 

n 1

 

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

таммовского состояния пропорциональна l

 

1

, то есть чем ближе таммовский уровень к зоне

 

 

 

 

1,

0 , тем больше глубина проникновения таммовского состояния вглубь кристалла.

Приложение

Поскольку нас интересует решение уравнения det A 0 , то удобно перейти к новым переменным, разделив все элементы в матрице А на 2 и обозначив диагональные элементы

 

E

. Тогда

 

– это просто обезразмеренная энергия,

которая теперь отсчитывается от

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

уровня . При таком обезразмеривании матрица А приобретает вид

 

 

 

 

 

 

 

 

 

1/ 2

0

...

0

0

0

 

 

 

 

 

 

 

1/ 2

 

1/ 2

0

...

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1/ 2

 

1/ 2

...

...

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

0

0

1/ 2

...

1/ 2

...

...

 

(А1)

 

 

 

 

 

0

...

...

1/ 2

 

1/ 2

0

 

 

 

 

 

 

 

0

0

...

0

1/ 2

 

1/ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

...

0

1/ 2

 

 

 

 

 

 

 

 

 

 

Раскладывая определитель этой матрицы размерности (N x N) по первой строке и далее по первому столбцу получим рекуррентное уравнение

d

 

d

 

 

1

d

 

(А2)

N

N 1

 

N 2

 

 

4

 

 

 

 

 

 

 

 

 

Будем искать решение этого уравнения в виде dN uN . Подставляя такое выражение для dN в уравнение (А2), получим квадратное уравнение на u ,

u2 u

1

0 ,

(А3)

4

 

 

 

вид решений которого существенно зависит от величины . Действительно,

9

 

 

 

1

 

 

 

 

u1,2

 

2 1

 

 

 

 

2

 

 

 

Теперь рассмотрим различные значения параметра :

1) 1 1. В этом случае удобно

воспользоваться

cos,

0 . Тогда (А4.1) можно записать в виде

u1,2 12 cos i sin

и решение рекуррентного уравнения (А2) записывается в виде

 

 

 

 

(А4.1)

 

следующей

параметризацией:

1

e

i

,

(А4.2)

2

 

 

 

 

 

d

N

Aei N Be i N

(А5)

 

 

 

Константы А и В определяются из граничных условий, а именно из очевидных значений определителя при частных значениях индекса N = 1 и 2.

d

Aei Be i cos

 

 

 

 

1

 

 

 

 

 

(А6)

 

 

Ae2i Be 2i 2

1

 

cos2

1

d2

 

 

 

 

4

 

4

 

 

 

 

 

 

 

Решая систему этих двух уравнений относительно констант А и В и подставляя полученный выражения в (А5) получим после некоторых преобразований следующее выражение для определителя dN

d

 

 

sin (N

1)

 

(А7)

N

2N sin

 

 

 

 

 

 

 

2) 1 . В этом случае удобно воспользоваться следующей параметризацией: ch,

0 .

Тогда (А4) можно записать в виде

 

 

 

 

 

 

u1,2 12 c h s h 12 e

и решение рекуррентного уравнения (А2) записывается в виде

d

N

Ae N Be N

(А5)

 

 

 

Константы А и В опять определяются из граничных условий, а именно из очевидных значений определителя при частных значениях индекса N = 1 и 2.

d Ae Be c h

 

 

 

 

1

 

 

 

 

(А6)

 

 

Ae2 Be 2 2

1

c h2

1

d2

 

 

 

4

4

 

 

 

 

 

 

10