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

книги из ГПНТБ / Кизель В.А. Отражение света

.pdf
Скачиваний:
9
Добавлен:
24.10.2023
Размер:
11.45 Mб
Скачать

1 0 8 МОЛЕКУЛЯРНАЯ ТЕОРИЯ ОТРАЖЕНИЯ СВЕТА [ГЛ. 3

Р = Л^рЕэфф =

(Е + Евтор)

(10.5)

(где г), Ф, £ — углы Эйлера осей молекулы;

— число

молекул в единице объема), а затем интегрированием по всему объему среды 2, точнее, по всему пространству,

 

Cppüa'l

Z I

(Вакуум)

Рис. 44. Расположение взаимодействующих диполей.

0 —начало координат;

т , 0,

і ,

А — номера

диполей, влияющих на избранный

(находящийся в сфере

а);

X ,

Y,

Z — точка

вне среды 2; в среде I показан

 

фронт

первичной

волны.

в котором поле молекул отлично от нуля, за исключени­ ем «сферы молекулярного действия» (поверхность о с радиусом а), окружающей рассматриваемый диполь:

г

р(гѴ_4

)

dv' =

Еэфф (г, t) = Е (г, t) + I

rot rot —-----

£------

-

V '

 

 

 

 

 

= E (г, t) +

EBTOp (r, t). (10.6)

Позволительность усреднения (10.4) мотивируется предположением «идеальным беспорядком», т. е. равно­ вероятностью всех ориентаций и большой величиной N\. Замена суммирования интегрированием Р(г) означает как бы «размазывание» дипольного момента «точечного» диполя по объему пространства, приходящегося в сред­ нем на один такой диполь.

§ 10] СРЕДЫ МОЛЕКУЛЯРНОГО СТРОЕНИЯ 109

Замена дискретной среды континуальной обычно мо­ тивируется, кроме сказанного выше, линейностью урав­ нений колебаний поля; однако некоторое затруднение представляет учет особенностей, которыми обладает кон­ тинуальное поле в местах самих осцилляторов (точки R = 0). С этим связана необходимость выделения сферы о из объема интегрирования. Излагаемые ниже неболь­ шие расхождения в данных различных авторов вызваны неодинаковым выбором и учетом «сферы действия» мо­ лекулы.

При расчете межмолекулярные силы и эффекты «свя­ зи осцилляторов» не учитываются. Не учитывается также возможность изменения частоты в излучающих центрах.

Решение интегродифференциального уравнения (10.6)

ищется ') в виде

 

 

 

 

Еэфф (г, t) = Аэфф (г) еіаі_— Аэфф е'(“' - к*г),

(10.7)

где k3= n |k |s a;

п — некоторая

величина, подлежащая

определению; |к|

= Л . Иначе

говоря,

ищется

«волна

поляризации»

С

 

 

 

 

 

 

 

 

Р = іУхрАэфф е‘1“/_кэГ).

 

(10.8)

Возникновение волны со скоростью

<СС можно

объяснить следующим схематическим рассуждением. Пусть на поверхность раздела падает нормально к ней плоская волна Е (см. рис. 44). Эта волна возбуждает диполи, лежащие на поверхности раздела, с одинаковой амплитудой и фазой.

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

см. § 16).

') Очевидно, по аналогии с макроскопическим рассмотрением гл. 1.

п о МОЛЕКУЛЯРНАЯ ТЕОРИЯ ОТРАЖЕНИЯ СВЕТА [ГЛ. 3

Выбрав решение

(10.7)

и (10.8),

можно

записать:

Евтор =

(‘

р (г' - ' - 4 )

 

 

 

 

 

\ rot rot —-

----------- du' =

 

 

 

 

V'

 

 

 

 

 

 

 

 

= N $

rot rot j-i- Аэфф (r') еш е

' е

dv',

Евтор =

Af1ße-tW 1J rot rot {Аэфф (r') В (R)} dv' = A BTOp eiat,

 

 

 

 

- -1R

 

(Ю.9)

 

 

 

 

 

( 10. 10)

 

 

 

 

 

 

 

Очевидно,

 

 

 

 

 

 

 

 

Ѵ*Аэфф +

п9к2А,фф =

0,

 

(10.11)

 

 

 

V2ß+k2ß =

0.

 

(10.12)

Это — волновые

уравнения

для

волн

со скоростью

с/іг

и с.

 

rot

берется по координате г точки на­

Операция rot

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

При некоторых добавочных предположениях о выборе «сферы действия» и величины а можно показать, что

при а-*- 0

 

 

 

 

 

 

 

 

J

rot rot (Аэфф (r') В (R)} du' =

 

 

 

 

 

 

 

= rot rot j Аэфф (r') В (R) du'

Аэфф (г). (10.13)

 

Пользуясь

теоремой

Грина

и

учитывая (10.11)

(10.13), можно получить

 

 

 

 

 

 

 

j

Аэфф (г' ) B ( R ) d v ' =

 

 

 

 

 

 

 

 

 

A

j ,

d B

__

R дАЭфф

ds'

 

= K - _

l ) k 2 { f Аэффалг

 

В —

 

 

 

 

д

дв_

р ДАзфф

ds', (10.14)

 

 

J ЛэФФ gN'

п

 

 

§ 10]

СРЕДЫ МОЛЕКУЛЯРНОГО СТРОЕНИЯ

111

где интеграл берется по поверхности 2, ограничивающей объем V' снаружи и по поверхности о, а щр— произ­

водные по нормалям к этим поверхностям. Далее,

А ”т° р = (»1= - f) k5 rot rot2I{ А »

Ф

Ф

*

1

 

 

 

 

(«2

 

 

 

 

1) к 2

Х

X rot rot J {Дэфф (r')

-

В (R)

ЗА.эфф ( г') ds'

а

 

 

dN'

 

 

 

 

 

 

 

---- Т ” АэФФ

 

гвтор +

а:втор I

(10.15)

 

 

 

 

 

и отсюда

 

 

 

 

 

Еэфф = Е ~Ь Евтор -{- Евтор.

 

(10.16)

Очевидно, что в среде 2 прошедшее поле Ed определится из соотношений

и

Dd =

rt2Ed =

Ed + 4яР

(10.16а)

Ed =

 

 

 

n2 —1Р,

 

а также из выражений (10.7) — (10.9)

 

Erf =

Л^Аэфф (г)

^

—tnlkls .г

Л^А°фф (г) в

d еш ,

где

 

 

 

(10.17)

 

 

 

 

sd~— s3-

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

Е^хор = Аэфф (г) - f - TViß

еш =

^ ФФ. ( г ) ^ - ^ ^ е ^ Г ‘кэГ. (10.18)

112 МОЛЕКУЛЯРНАЯ ТЕОРИЯ ОТРАЖЕНИЯ СВЕТА [ГЛ. 3

Если переписать (10.16) без временных множителей, используя (10.18), получим

4 ф ф ( г ) е- ВД5эГ = Е 0 ( г ) е- ^ І з г +

+ А^фф(г)Ч- Л”фф (г) ~fzt] N 1ße~‘”|k|S3r. (10.18а)

Поскольку экспоненты с разными показателями дол­

жны быть линейно независимы, отсюда следует

 

4л п* + 2дгй

1

(10.19)

3 «2— 1

іѵіР —

 

Ео (г) e—,'li{/sr= А^тор (г),

(10.20)

или иначе

 

 

(10.21)

Еэфф

Евтор»

 

Е -|- ЕВТОр =

0.

 

(10.22)

Соотношение (10.22) носит название теоремы ногашения Озеена. Физический смысл изложенного заключа­ ется в следующем. Эффективное поле в среде (10.6) состоит из внешнего поля Е и поля вторичного излучения частиц среды. Последнее поле состоит из излучения то­ чечных дипольных осцилляторов в вакууме. Оно может

быть разделено на две

составляющие — одна из них,

ЕІтор , имеет скорость

с и (внутри среды) гасит волну

Е; отсюда название теоремы (10.22). Вторая состав­ ляющая Еа, имеющая скорость с/п, дает преломленную волну; неизвестная константа п, введенная в (10.7), оп­ ределяется из (10.19).

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

в Автор (г)

имеет вид

 

е —[|k[R g —M |k |s ^ r',

Автор — это

сумма излучений точечных осцилляторов

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

§ 10] СРЕДЫ МОЛЕКУЛЯРНОГО СТРОЕНИЯ 113

дение о виде закона модуляции, необходимом для со­ блюдения (10.21) и (10.22).

Проведенное рассуждение показывает, что первичная волна гасится в основном осцилляторами, лежащими у поверхности (ибо вклад осциллятора в Еатор быстро спадает по мере удаления его от поверхности). Очевид­

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

Именно поэтому в среде 2 волна Е отсутствует на заметных расстояниях от поверхности раздела.

Эффективное поле в месте нахождения какого-либо осциллятора [определяемое (10.21)] создается компо­

нентой Евтор — в основном полем ближних осциллято­ ров. Эти качественные соображения в настоящее время подтверждены точным расчетом (см. стр. 117).

Из формулы погашения можно получить условия для амплитуд (формулы Френеля) и определить направление к0. Прямым вычислением Е2 для верхнего полупростран­ ства молено получить выражение для отраженной волны.

Путем некоторых несложных, но длительных преоб­ разований, ограничиваясь приблилеением кг^>1 и исполь­ зуя известную формулу

rot rot Ae£|k,rs= — к2 [s [sA]] e['lkfrs,

выралеения (10.18a) — (10.22) можно привести к виду

E0e -£lftlsr

2л дГ о sin (ф + ф)

Ы

м ”, ф ф ( г ) ] ] г ' |к |‘ 2 г ,

Л2 — 1

^ sin ф COS ф

 

 

(10.23)

где

 

 

 

 

 

 

ч

 

Ф =

sN', ф — ф =

 

SSd ,

отсюда

S s= s;

молено получить

Пsin lj)= sin ф.

Сопоставляя (10.17) и (10.23), можно получить формулы Френеля для амплитуд Ed.

Если точку наблюдений выбрать вне среды (в полу­ пространстве 1), то в этой точке D = E BHDm, Рвпеш=0, Е„пеш=Е+Ег. В этом случае выделять сферу о не нужно,

8 В. А. Кизель

114

МОЛЕКУЛЯРНАЯ ТЕОРИЯ ОТРАЖЕНИЯ СВЕТА

[ГЛ. 3

и операцию rot rot можно менять местами с интегриро ванием без поправок. Тогда получим

г

 

с

Ег = Аг (г) еш — rot rot

R

dv , Аг — Авторѣ

V'

 

 

 

аналогичными несложными, но длительными вычисле-

ниями получим фг = Я — ф, фг = s rs и формулы

Френеля для амплитуд отраженного света.

Подробности всех промежуточных вычислений см., на­ пример, в работах [02, 03].

Все сказанное можно повторить в несколько иной форме. После составления основного уравнения (10.1) ищутся условия, при которых внутри среды, в результа­ те интерференции первичной и всех вторичных волн, опи­ сываемой (10.1), формируется плоская световая волна (преломленная), вообще говоря, с иными значениями к и п, чем у падающей, а последняя в результате интер­ ференции в среде гасится (отсюда название — «теорема погашения»). Таких условий должно быть, очевидно, два — одно представляет собой закон преломления (ус­ ловие для ф, к), другое есть условие для п. Только при этих значениях п, ф, к устанавливается искомое волно­ вое поле1). Расчет поля вне среды при соблюдении ука­ занных условий дает отраженную волну. Иначе говоря, интерференционное поле содержит два типа волн— распространяющихся со скоростью о/л и со скоростью ц. Первые образуют отраженную волну и гасят первичную в среде, вторые образуют преломленную волну в среде. Эти утверждения справедливы лишь для установивше­ гося процесса.

§ 11. Дальнейшие уточнений

\

Формула Лоренц — Лорентца

(10.19)

 

п* + 2 _

3

_

 

п2 — 1

4яЛ/]ф

°

 

получается в качестве условия для л лишь в упрощен-

*) Таким образом, теорема, по существу, есть граничное условие.

§ п]

ДАЛЬНЕЙШИЕ УТОЧНЕНИЯ

1 1 5

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

В этих расчетах выбирается такой размер сферы а, чтобы в нем находилось какое-то число молекул, сум­ марное усредненное влияние которых на выделенную молекулу принимается равным нулю (ср. также [1]).

При более совершенном выводе соотношения Ло­ ренц— Лорентца следует вместо простого усреднения по распределению в пространстве и ориентациям, которые считаются случайными и некоррелированными [формула (10.4)], учесть функцию распределения, точнее, унарную и бинарную радиальные функции распределения и кор­ реляцию ориентаций (уже не считая осцилляторы изо­ тропными).

Попытка подобного учета проведена, например, в работе [2], где учитывалась и тернарная функция рас­ пределения. Было получено выражение

я * -1

<ß> <#i>

1+ ff

m

n

7i3 + 2 ~

3

1 +<ЛЛ>£(1 + RY

{

;

где D и R — некоторые функции, при термодинамическом равновесии зависящие лишь от (А^), Т, ю и учитывающие указанные распределения. Считается, что Объем а выделяется так же, как и выше, и последующая про­ цедура расчета интеграла остается без изменений; ис­ пользуется лишь значение и из (11.1) [3, 4].

В теории [5] выбиралось такое значение а, чтобы в нем помещалась только одна рассматриваемая молекула,

и я-с А, а решение искалось только

для точек, удален­

ных от поверхности на г^>Я; тогда

при вычислениях

возможно разложение (10.18) в ряды по степеням а/К. Соответственно при расчетах для п получено более слож­ ное условие

я2+ 2

п 2 + 1 0 / 2 л а \ 2 , 2 , [ 2 п а \ 3 _ п

( 11.2 )

п 2 — 1

10 ( X ) г 3 1 [ X )

U '

 

8*

1 1 6

МОЛЕКУЛЯРНАЯ ТЕОРИЯ ОТРАЖЕНИЯ CDETA

[ГЛ. 3

В этих расчетах формулы Френеля получались как

первые

члены некоторых разложений по степеням

а/К,

т. е. не были вполне точными и, строго говоря, были при­ менимы лишь на расстояниях от поверхности, много боль­ ших а. Последний вывод, проведенный по той же схеме,

но существенно

более точный1), дан в работе [6]; для

п получено условие

 

 

 

 

 

 

 

о

(7

sin ап\ /

1

1

~

i

\

I

 

3*—

[[cos а л -----— )

- д

а

д а )

+

 

 

 

 

,

1 -j-

іа sin art)

__Q

(11.3)

 

 

 

 

n2 — 1

an

J

 

 

 

 

 

 

 

Здесь формулы Френеля получены как строгий резуль­

тат2)* (см. также

[7, 8]). Ограничиваясь

членами до

третьего порядка по а, из (11.3)

получим

 

п2 + 2

, п2 + 10 а2 +

і -у- а3 = С.

(11.4)

п 2 — 1

 

 

 

Оценки [05] говорят, что формулы Френеля практически становятся применимыми уже на расстояниях порядка а от поверхности, т. е. во всех случаях, интересных для оптики.

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

кулы уравнения (10.6) — (Ю.8) логически неправомерно. По-видимому, изложенный вывод скорее дает моле­ кулярную интерпретацию явления и физического смысла

')

Уточнения касались опять-таки вычисления интегралов.

что

2)

Появление мнимых членов в (11.2) и (11.3)

связано с тем,

при вычислении

производных под интегралами (10.14)

и (10.15),

на­

пример,

 

 

 

 

 

 

 

 

3В [R) _

д

е

 

 

 

 

 

 

дN'

w

-

I

 

 

 

 

учитывается также и дробь

отбрасываемая

при

упрощенном

kR

расчете для kR^> 1.

§ И] д а л ь н е й ш и е у т о ч н е н и я 1 1 7

констант, нежели доказывает единственность решения (10.6). Этот вопрос до сего времени дискутируется (см., например, работу [9], где автор метко замечает, что «молекула в глубине среды не видит поля Е», а так­ же работу [10]). Здесь, по-видимому, следует искать и определять некоторое «самосогласованное поле».

Предположение о беспорядочном расположении ос­ цилляторов принципиально не существенно. Были сде­ ланы расчеты для простой кубической решетки, где по­ верхность среды образована одной из граней решетки и, таким образом, представляет собой квадратную сетку молекул; вместо интегрирования производилось сумми­ рование полей излучения указанных сеток. Этот метод предложен Эвальдом [11] и приводит к тем же резуль­ татам. Впоследствии он был усовершенствован Сивухиным [12] путем явного учета неоднородных компонент волнового поля внутри решетки и конечности расстояния, на котором нижележащие слои перестают давать свой вклад в поле у поверхности.

Можно показать, что влияние одних слоев на другие весьма быстро спадает по мере их взаимного удаления. В работе [10] получены следующие относительные зна­ чения энергии диполь-дипольной связи между избранным ^ диполем и смежным слоем — сеткой:

Э (0) — 0,359

(в скобках указан порядковый номер

Э (і) — 0,013

слоя, в котором находится избранный

Э (2) — 0,00002

диполь).

Э (3) — 0,00000004 и т. д.

 

Была показана возможность динамического равнове­ сия поля Е с полем вторичного излучения осциллято­ ров — диполей и их колебаниями. Отметим, что подобная же возможность равновесия при наличии неоднородных волн (для частного случая полного внутреннего отраже­ ния) была показана также в работе Сотского и Федорова [6]. Вытекающая из изложенных теорий и расчетов комплексность п может быть учтена, как в гл. 1, введе­ нием комплексного угла % при сохранении неизменного вида формул Френеля или же записью этих уравнений в явно комплексном виде. Поскольку в данном случае мнимая часть весьма мала, второй способ удобнее для

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