книги из ГПНТБ / Кизель В.А. Отражение света
.pdf1 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 = |
4я |
|
|
|
|
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, введе нием комплексного угла % при сохранении неизменного вида формул Френеля или же записью этих уравнений в явно комплексном виде. Поскольку в данном случае мнимая часть весьма мала, второй способ удобнее для