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

книги из ГПНТБ / Татевский В.М. Классическая теория строения молекул и квантовая механика

.pdf
Скачиваний:
101
Добавлен:
27.10.2023
Размер:
22.71 Mб
Скачать

ГЛАВА XXIX

МЕТОД ВАЛЕНТНЫХ СХЕМ

§1. Общая характеристика метода

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

выбирается

набор

заранее

заданных известных

функций

%Р(Х, у, Z),

зависящих каждая только от трех пространственных ко-

ординат: %\{х,у, г),

%N(x, у, z).

 

 

 

 

 

 

Каждая

из этих функций умножается на спиновую функцию г| (о),

т.е. а (а)

или

Р(а), и

составляются

все возможные

определители

вида:

 

 

 

 

 

 

 

 

 

 

 

ф А

=

-£=r

І Чі ( 0

г) Хі ( 2 ) ^2 (а2) Ъ

(2)

• • • Л * (<*2) %N

(2)

і

(XXIX, I)

 

 

Ут

 

 

 

 

 

 

 

 

 

где в качестве

функции г) (а)

в любом

столбце может

быть

постав­

лена либо а(о),

либо

Р(о).

 

 

 

 

 

 

 

При заданных функциях %, очевидно, можно составить 2N

таких

определителей *.

 

 

 

 

 

 

 

 

Поскольку функции Хр заранее заданы, то функции tyh, по­ строенные каждая совершенно определенным образом из функций Хр, также заранее заданы.

В принципе семейство функций Ф, из которого выбирается оптимальная приближенная функция основного состояния моле­ кулы, может быть определено как совокупность функций, представ­

ляющих каждая некоторую линейную

комбинацию функций %

 

 

2N

 

 

 

 

Ф = 2 С А

 

(XXIX, 2)

 

 

й=1

 

 

 

* Здесь мы изложим только такой вариант метода валентных схем, в котором

все

функции хь

• •. XN различны. Можно рассматривать

и такие определители

типа

(XXIX, 1),

в которых в некоторых парах

столбцов,

отличающихся спино­

выми функциями, стоят одинаковые координатные функции. Этого более общего варианта метода валентных схем мы рассматривать не будем.

Все возможные вариации коэффициентов С\ в этой линейной

комбинации определяют

семейство функций, из которого

выби­

рается приближенная волновая функция основного состояния.

Семейство

функций

Ф,

определенное на

функциях

tyk по

(XXIX, 2), по

существу,

есть

семейство функций,

рассматриваемое

в линейном варианте вариационного метода Ритца. Поэтому по­

строение приближенной функции

Ф

вида (XXIX, 2), в

принципе,

представляет собой непосредственное

применение линейного ва­

рианта прямого вариационного

метода

Ритца для' нахождения

функции Ф, приближенно

апроксимирующей истинную

функцию

ф 0 основного электронного

состояния

молекулы.

 

Из вариационного условия

(XXVII, 17)

 

 

б

j "

Ф ' Я Ф dv da =

О

(XXIX, 3)

и условия нормированности

функции Ф

 

 

 

J Ф*Ф dvda

= l

'

(XXIX, 4)

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

Именно, если ввести обозначения

Hkl

=

|

%H^t dv da

(XXIX, 5)

и

 

 

 

 

SH

=

j

dv da

(XXIX, 6)

и учесть, что функции % нормированы, т. е.

s f t f e = \ ^ k d v d a = \

(XXIX, 7)

то лучшие значения коэффициентов Ck в выражении (XXIX, 2) определяются из п = 2N уравнений

и - е) С, + и + eS1 2 ) С 2 + . . . +

(Я,„ -

е5,„) Сп = О

 

(XXIX, 8)

( # „ , - eS n l ) С, + . . . + Пп - е) Сп

= О

 

Эта система п однородных уравнений имеет решения, отличные от нуля только если определитель системы равен нулю, т. е. если

Нц — г

Hl2 — e S 1 2 . . . H m — eSln

(XXIX, 9)

 

 

= 0

fni — £Sr a i

• • • Нпп

є

 

Последнее уравнение определяет п значений е т , которые ему удов­ летворяют

Є ї , в 2

е т > . . . . е „

(XXIX, 10)

Каждому такому 'значению ет, если его подставить в уравнения (XXIX, 8), соответствует свой набор коэффициентов Си, именно

С\т), С{2т), . . . . C( f e m ) С{] (XXIX, I I )

определяющий одну из функций Ф

(XXIX, 2),

удовлетворяющую

поставленным условиям (XXIX, 3)

и

(XXIX, 4).

 

 

 

 

Таким

образом, в результате

получим п функций

 

 

 

 

 

 

 

 

Ф Т

= 2 Л

 

 

 

 

(XXIX,

12)

 

 

 

 

 

т =

ft

 

п

 

 

 

 

 

 

 

 

 

 

 

1,

2

 

 

 

 

 

 

 

удовлетворяющих условиям

(XXIX, 3)

и

(XXIX, 4).

 

 

 

Из этих функций надо выбрать ту, для которой значение энер­

гий, вычисленное с этой

функцией:

 

 

 

 

 

 

 

 

 

 

£<"» =

j "

ФтРФт

dv da =

2

С^'С^Ны

 

 

(XXIX, 13)

будет наименьшим.

 

 

 

 

 

 

 

ft.

і.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Легко показать, что. энергия Ёт

(XXIX, 13), вычисленная с функ­

цией Ф т

(XXIX, 12),

 

совпадает

 

с

соответствующим

корнем

е т

уравнения

 

(XXIX, 9).

Действительно,

подставив

в

уравнения

(XXIX, 8)

Вт вместо е, получим

уравнения для

коэффициентов

СІт).

Умножим к-ое уравнение на С{кт)" и просуммируем

по к. Тогда

по­

лучим'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

С« т ) * 2

С(Ґ

(НЫ

-

-

С

 

 

(XXIX, 14)

Отсюда

 

ft

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«« 2

СҐ'СҐЗкІ

 

-

2

Ф'Ф^Ы

 

 

 

(XXIX, 15)

 

 

ft,

/

 

 

ft,

 

 

 

г

 

 

 

 

 

 

Из условия

(XXIX, 4)

и

выражения

(XXIX, 15)

получим

 

 

 

 

е т

=

2

С(кЩт)Ны

= £ < т )

 

 

 

(XXIX, 16)

 

 

 

 

 

ft,

г

 

 

 

 

 

 

 

 

 

 

Поэтому,

выбрав

Ф т ,

соответствующую

наименьшему корню

уравнения

(XXIX, 9), получим приближенную

функцию

наинизшего

по энергии электронного состояния системы (при выбранной кон­ фигурации ядер).

Обозначив наименьший из корней уравнения (XXIX, 9) через

е0 , а соответствующую ему функцию (XXIX, 2) через Ф0 ,

установим

таким образом приближенную функцию и приближенное

значение

энергии наинизшего по энергии электронного состояния

(при вы­

бранной конфигурации ядер):

 

§ 2. Упрощение метода при использовании функций г|^» собственных для операторов Sz и S2

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

шать уравнение

(XXIX, 9)

степени п =

2N, т. е. очень высокой

степени. Во-вторых, при вычислении матричных элементов

Ны и

Ski оказалось бы, что многие из них равны

нулю и вследствие

этого

матрица \\Ны

BSMW распадалась бы на диагонально расположен­

ные блоки, а определитель

(XXIX, 9) — н а

произведение определи­

телей более низкого порядка. Иными словами, задача свелась бы все равно к нескольким более простым. Такое сведение описанной задачи к более простым можно предусмотреть прежде, чем пере­ ходить к ее решению, если вместо линейной комбинации всех 2N функций tpfc рассматривать отдельно линейные комбинации, состав­

ленные каждый раз из выбранных определенным образом

функций

фи или их определенных комбинаций.

 

 

 

Под методом валентных схем обычно и понимается такой ва­

риант построения приближенной волновой функции Ф

основного

электронного состояния

системы, при котором эта функция

строится

не непосредственно

как линейная комбинация 2N

определителей

вида

(XXIX, 1),

а

как линейная

комбинация некоторых

функций,

каждая из которых

является определенной линейной комбинацией

некоторых

из

2N

определителей

вида (XXIX, 1),

построенной по

соответствующему

рецепту.

 

 

 

 

Сущность этого метода состоит в следующем. Если в какой-

либо

функции

Uph содержится

iVa* столбцов, в которых

стоят спи­

новые функции

a (а), и N$k)

столбцов, содержащих

спиновые

функ­

ции

Р(а),

то

каждая

из

таких

функций яр* соответствует

такому

спиновому состоянию всех электронов системы, для которого сум­

марная проекция спина всех электронов на какое-либо

выбранное

направление (например, ось OZ) равна

 

 

 

 

 

 

 

 

 

 

 

 

(XXIX, 17)

Очевидно, что максимальное значение Sz

=

N/2-h/2n, минимальное

значение Sz

=

—N/2-h/2n.

могут быть разбиты

на

группы

Таким образом, все функции

так, что в

каждую группу войдут

функции

с одним

и тем же

значением Sz,

г. е. с одним и тем же значением разности

М?'

NfK

Обозначая значение Sz

в единицах

h/2n

верхним индексом

а, полу­

чим

такие

возможные

группы функций

грь(а = —N/2,

 

1,0,

.+ 1

+N/2 ) (для четного числа N электронов)

 

 

(XXIX, 18)

Можно показать, что матричные элементы, входящие в уравне­ ние (XXIX, 9), определенные на функциях разных групп (XXIX, 18), равны нулю, т.е.

 

Н0'

°' =

J* ^ а ) * Щ ( а

, )

dvdo

=

0

( X X I X , 19)

 

 

 

а ф &

 

 

 

 

 

 

S 0 ,

°' =

J" ^ a ) ' ^ a , )

d v d a

= 0

 

( X X I X , 20)

Поэтому матрица коэффициентов НЫ

Е8Ы

уравнений

(XXIX, 9)

распадается

на диагональные блоки,

 

а уравнения (XXIX, 8 ) н а

отдельные системы уравнений, каждая

из которых относится только

к функциям

г|з(0' одной

группы. Следовательно, рассматриваемая

задача распадается на отдельные задачи, каждая из которых бу­ дет включать только функции i|/a> одной группы.

Поэтому имеет смысл с самого начала рассматривать не мно­ жество функций Ф, определенных на всех грь уравнением (XXIX, 2), а совокупность множеств Ф<а), каждое из которых определено на функциях ф только одной группы 1\)<-а\

Перенумеруем функции г|/ст>

в каждой

группе (XXIX, 18) ин­

дексом / ( / = 1 , 2 ,

fa). Тогда

множества

функции

ф№, каждое

из которых определено

на одной

группе функций i|)( f a ) ,

будут

 

ф № = 2

CJaWfa)

 

( X X I X , 21)

Если потребовать, чтобы функция ф(ст) удовлетворяла

усло­

виям (XXIX, 3) и (XXIX, 4),

то для определения

соответствующих

значений коэффициентов

в ее представлении

(ХХ1Х,21)

полу­

чим систему уравнений, аналогичную (XXIX, 8),

но

содержащую

меньше уравнений, именно

столько, сколько функций

i j 3 ( f a )

в дан­

ной группе. Соответствующим образом понизится и порядок урав­

нения,

анолигичного

(XXIX, 9).

 

 

 

 

 

 

 

 

В

методе валентных

схем

при четном

числе электронов

часто

заранее

предполагают

(хотя

это и не может быть обосновано и

во

многих

случаях

противоречит

экспериментальным

данным),

что для основного электронного состояния

Sz

= 0. Тогда

рассмат­

ривают только функции i])( f a )

той группы,

для которой а =

0, т. е..

Na

= N$ и Sz

= 0. Мы рассмотрим

дальнейшие

преобразования

в

методе

валентных

схем

только

для

группы

функций aj)(f0) с

Sz

= 0,

хотя,

строго

говоря, всю описанную

ниже

процедуру

следовало

бы

производить

для функций каждой группы ip(fCT>

(XXIX, 18)

и затем

из

результатов

расчета

приближенной

функ­

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

Д ля каждой группы например группы ip( f 0 ) , можно по­ лучить дальнейшее упрощение задачи, если использовать для со­ ставления функции Ф не сами функции ip( f a ) , а некоторые их ли­ нейные комбинации ф(с т .s ). Мы рассмотрим такое преобразование только для группы функций ipj.°>. Для остальных групп они могут быть выполнены аналогично. Смысл этого преобразования со­

стоит

в том, что можно построить такие линейные комбинации

ф(°. 8)

функций группы гр<0), что каждая такая линейная комбина­

ция

будет соответствовать не только определенному значению

суммарной проекции спинов всех электронов, но и определенному значению квадрата вектора общего спина системы. Если строить затем искомую функцию Ф из таких линейных комбинаций ф(°-8 ), то задача еще упрощается, так как система уравнений, опреде­ ляющих коэффициенты разложения искомой функции Ф по этим линейным комбинациям функций гр<0), снова распадается на от­ дельные системы с меньшим числом неизвестных, а соответствую­ щее характеристическое уравнение распадается на несколько урав­ нений более низкой степени. Введение функций ф(а 8 > как линейных комбинаций гр^, соответствующих определенному значению Sz и определенному значению квадрата вектора общего спина си­ стемы, позволяет заранее учесть эти возможности. Не доказывая здесь сказанного, чтобы не загромождать изложения, ниже мы приведем только метод построения функций ф(°-8>, соответствую­ щих определенной проекции общего спина системы и определен­ ному значению квадрата общего вектора спина системы.

Функции ф(°'я ) строятся из функций яр<0> следующим образом*. Первой функцией ф(°>8), которую можно построить из i|)<f0', будет

(о.

Л

 

ф\

2 ) В 2 ^

(XXIX, 22)

 

f

квадрату век­

Можно показать, что эта функция соответствует

тора спина системы, равному

 

 

5 2 ( 4 + ! )

( Х Х 1 Х -2 3 )

т.е. вектору спина системы с квантовым числом S = JV/2. Другие функции Ф<0 '8 ) могут быть построены следующим путем.

Возьмем функции гр<°> и выделим во всех соответствующих определителях какие-либо два столбца, например столбцы с но­ мерами і и /.

*

Функции

Ф 1 0 ' *\ встречающиеся в этом параграфе, все

не нормированы.

Чтобы

получить

нормированные функции, нужно умножить

каждую функцию

Ф**7, ^ на соответствующий нормирующий множитель. Так как метод валентных

схем в дальнейшем не используется, на этом вопросе мы

останавливаться не-

будем. Условие нормированности приближенной функции Ф

(XXIX, 30), встре­

чающейся

ниже, учитывается при решении задачи вариационным методом.

14 Зак.

454

417

Построим линейную комбинацию функций ip(f0) вида

 

 

 

 

o ( 0 ' - T - 1

) =

2 Y

l . / ( f ) ^

 

 

 

(XXIX, 24)

 

 

 

 

 

 

 

f

 

 

 

 

 

 

Здесь Yij(/)

имеет следующие

значения:

 

 

 

 

 

 

Y(/(/) =

+ 1,

если

в /-ом

столбце

i j ^ 0 )

стоит а (о), в /-ом В (а)

 

 

 

— 1,

»

 

 

 

»

В (а),

— а

(а)

(XXIX, 25)

 

 

О,

»

 

 

 

»

а (а),

— а

(а)

_ —

 

 

 

 

 

 

 

0,

»

 

 

 

»

В (а),

- В (а)

 

 

 

 

 

 

 

 

 

 

 

/

л?

л

 

 

 

Можно

показать,

что

построенная функция

ф \

2

/

соответ­

ствует значению квадрата общего вектора спина

системы,

рав­

ному

 

 

^iMi-1 )^-

 

 

 

 

 

 

 

 

 

 

 

 

( x x i x , 2 6 )

т. е. для этой функции квантовое число спина S =

JV/2

1. Выде­

ляя

разные

пары столбцов в определителях ip( f 0 ) ,

можно построить

набор таких функций. Далее

можно

во всех определителях группы

i j 3 ( f 0 )

выделить

две пары

столбцов, например пару

столбцов

с но­

мерами і и / и пару столбцов с номерами т и п, и построить функции

 

 

 

 

 

 

S Y / / t f ) Y M n ( / > * } 0 )

 

(XXIX,

27)

 

 

 

 

 

 

 

f

 

 

 

 

 

ЗДЄСЬ Утп(ї) ПОДЧИНЯЮТСЯ ТЄМ ЖЄ СООТНОШЄНИЯМ,

ЧТО И Yij(f), т.

е.

Y m n ( f ) = ± i

или

0 в

зависимости от того, стоят ли

в столбцах

с номерами

тип

определителя ар(?0)

а и р или р

и а, или а и а,

или р и р соответственно. Это будут функции со

значением кван­

тового числа

S полного спина системы, равным

 

 

 

 

 

 

 

 

5

=

- j -

-

2

 

(XXIX,

28)

Выделяя

разные

пары

столбцов в

определителях

t|)( f 0 ) ,

можно по-

строить

набор таких функций

ф \

' 2

/.

 

 

N),

Действуя подобным образом, можно, наконец

(при

четном

разбить все N столбцов каждого определителя на пары, например

на пары

 

 

 

/),

(т,

n),

(s,

t), ...

 

 

 

 

 

 

 

 

 

 

и построить

функцию

 

 

 

 

 

 

 

 

 

 

 

Ф ( 0 ' 0 ) -

2

Ъ,

(/)

Утп (f) Y „ if> • • • +Г

 

(XXIX, 29)

со значением квантового числа S полного спина системы, рав­ ным 0.

Производя различные разбиения столбцов определителей на пары, можно построить набор функций ф(°>°) и перенумеровать их каким-либо индексом, например индексом k. Тогда набор функций ф(о, о) будет

 

 

 

Ф'»-«1,Ф(».»)

ф М ) . „

В

методе

валентных

схем обычно

принимают (хотя это в об­

щем

случае

не

всегда

верно), что для основного состояния си­

стемы с четным

числом

электронов

S =

0.

Поэтому приближенную функцию основного электронного со­

стояния

часто отыскивают как

линейную комбинацию только функ-

ции Фц'

':

 

 

 

Ф =

2 С Й Ф 6 0 , 0 )

(XXIX, 30)

 

 

к

 

Именно из этого множества функций (при вариациях Ch) и выби­ рают обычно приближенную функцию основного электронного со­ стояния при четном числе электронов.

Следует отметить,

что функции

Ф( £0 , 0 )

построены

по

опреде­

ленному

правилу

из

заранее заданных

функций

(XXIX, 1) и

поэтому

также являются вполне определенными для каждой кон­

кретной

задачи и

не

содержат никаких

варьируемых

параметров

или варьируемых

функций. Таким

образом,

функции

ФІ°'0 >

пред­

ставляют собой для каждой задачи набор заранее заданных функ­ ций, на базе которых фактически строится линейный вариант пря­ мого вариационного метода Ритца при отыскании Ф в виде (XXIX, 30).

Дальнейшее решение задачи проводится точно так же, как это было описано в § 1 настоящей главы. Путем, описанным выше, получаются системы уравнений, аналогичные (XXIX, 8), но с зна­ чительно меньшим числом неизвестных коэффициентов С?[, харак­

теристическое

уравнение,

аналогичное

(XXIX, 9),

и определяются

его корни em 0 )

и функции

Ф т ' 0 ) ,

в том числе его наименьший ко­

рень е( 0 °-0 )

и

соответствующая

ему

функция

из множества

(XXIX, 30)

 

f ^ S c W 0 , 0 '

 

 

 

(XXIX, 31)

 

 

 

к

 

 

 

Если

заранее не делать предположения о значении квантового

числа

S

общего вектора спина системы

для основного состояния,

то все

описанные расчеты должны быть

проведены для каждого

набора функций ф(0 'в ) и таким путем получены на базе каждого набора этих функций функции Фo0 ''s , , отвечающие наинизшим

корням &o°'S) для

каждого

набора ф^8К Выбирая из этих корней

ej,°'s> наинизший,

получим

приближенное значение энергии для

электронного состояния системы, наинизшего по энергии при вы­

бранной конфигурации ядер, а соответствующая

ему функция

14*

419

Фо

будет найденной приближенной функцией

этого

состояния.

Строго говоря, подобные выкладки должны

быть

повторены

для

каждой группы функций

oj)(fa) (XXIX, 18). Из

результатов

вы­

брана функция Ф, отвечающая самому низкому

значению

є из

всех

полученных описанным

путем для каждой

группы, однако

часто всего этого не делают, а исходят из какого-либо предполо­ жения о значении Sz и S для основного состояния и работают только с одной группой функции ilpW, как было описано выше для

группы l|)( f 0 ) .

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

ветствующую ему приближенную

функцию

Ф и энергию Е~(°> как

функцию

параметров RU

. . . , R 3 K

- 6 ,

нужно

исследовать

величины

е (о,s) к а к

ф у Н К ц И И

RU

- T I >

R3K_6

на

минимум. Основному

элек­

тронному

состоянию

будет

соответствовать

величина е'

s>

имею-

 

 

 

 

 

 

 

ТП

9

 

щая наиболее

глубокий минимум, лежащий

ниже диссоциацион­

ных пределов,

и отвечающая этой величине

е'т 5 ) функция Фет s > .

ГЛАВА XXX

ПРИБЛИЖЕННОЕ КВАНТОВОМЕХАНИЧЕСКОЕ ОПИСАНИЕ РЯДОВ МОЛЕКУЛ. ПОСТАНОВКА ЗАДАЧИ. ИСХОДНЫЕ ПРИБЛИЖЕНИЯ

§1. Постановка задачи

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

Второй

вопрос, решение

которого непосредственно вытекает

как следствие рассмотрения

первого вопроса, — это

вопрос о по­

лучении и

пути использования квантовомеханических

выражений,

описывающих свойства не одной молекулы, а больших рядов мо­ лекул.

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

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

(или

измерив

экспериментально),

например,

энергии

основных

*

Эти закономерности изложены во I I части

книги в

гл. X V I и

X V I I .

** Именно, в

приближении

Фока — Рузана

достаточно предположить при­

ближенное равенство функций ре

в пределах

объемов эквивалентных

фрагментов

любых

молекул. Как следствие,

получается приближенное

равенство

функций ри

в пределах объемов эквивалентных фрагментов. Так что вместо двух предполо­ жений, которые необходимо было ввести в гл. XXV при использовании точного решения электронного уравнения, при использовании волновой функции в при­ ближении Фока Рузана требуется только одно предположение-

Соседние файлы в папке книги из ГПНТБ