Косинский Ю.И., «Микроструктурная модель взаимодействия

advertisement
Косинский Ю.И.
Микроструктурная модель взаимодействия
электрического поля с оптически прозрачными средами
или принцип переизлучения на диполях в оптических
явлениях
Рассмотрим взаимодействие плоской электрической волны со средой,
представленной в виде произвольного числа p переизлучающих плоскостей –
бесконечно плоских решеток, в узлах которых расположены атомы. В атомах
такой решетки (переизлучающей плоскости) под действием внешнего
электрического поля с плоским фронтом волны наводятся диполи, которые
являются источником переизлучения плоской электрической волны по обе
стороны от плоскости. Для амплитуд электрических волн переизлучения
плоскостями в работе [1] была выведена система уравнений
zp
i k z z
i k zl
E l (z l )  r E 0 
r
' E j (z j )  l j
,
(1)
z j  z1
z j  zl

где
r
– эффективный коэффициент
между плоскостями, равный

r
r

r
1 r

r
r 
 2
1 r


,
r  i 2 N x y  ,




1
,   arc tg     ,


 r 


r   Cos( ) ,
Nx y
i
диполь–дипольного взаимодействия
Sin( )  1  r .
2
– поверхностная плотность расположения диполей на плоскости,
– поляризуемость атомов,
E 0 – амплитуда внешней электрической плоской волны,
(2)
(3)
(4)
E l ( zl )
–
амплитуда
электрической
волны,
источником
которой
l
является переизлучающая плоскость , расположенная в точке zl .
Решение системы уравнений (1) дает значение амплитуды электрической волны
для любого плоского излучателя, расположенного в пределах рассматриваемой
системы z1  zl  z p , где p – количество условных переизлучающих
плоскостей. Парциальные волны от плоского переизлучающего источника
распространяются в пространстве со скоростью света в вакууме в двух
противоположных направлениях от источника. При этом в работах [1], [2] было
показано выполнение закона сохранения потоков энергии для одной
переизлучающей плоскости и для резонатора, состоящего из двух
переизлучающих плоскостей. В связи с этим, оптическую среду можно
представить в виде многорезонаторной системы или полирезонатора, состоящего
из p плоских зеркал (переизлучающих плоскостей).
Модель взаимодействия электрического поля с веществом названа в данной
работе микроструктурной в связи с тем, что основная система уравнений (1)
записана для амплитуд волн переизлучающих плоскостей (атомов или диполей,
объединенных в структуру плоскости), а не для амплитуд волн отдельных
переизлучающих диполей, как это сделано в работе (3), название которой –
микроскопическая модель взаимодействия.
Для наглядности выпишем систему уравнений, связывающую амплитуды
парциальных волн отдельных источников переизлучения, для четырех
плоскостей, считая, что расстояние между любыми соседними плоскостями равно
z j  z j 1  z  mk . Из соотношения (1) имеем:

 r E
 r E
 r E

E1  r E o  E 2  i  E 3  2i  E 4  3i ,
E2
E3
E4

0
 i  E1  i  E 3  i  E 4  2i ,
0
 2i  E1  2i  E 2  i  E 4  i ,
0
 3i  E1  3i  E 2  2i  E 3  i ,


(5)
где   k z .
Решить систему уравнений (5) можно прямой подстановкой. Систему уравнений
для произвольно большого числа плоскостей решить прямой подстановкой
невозможно.
Систему уравнений (5) в случае произвольного числа плоскостей можно
записать в следующей форме
E l  rE 0 
где
( l 1) i
p l
 r  El j 
ji
j 1
p
– количество плоскостей переизлучения,
j
– индекс суммирования,
l 1
 r  E l  j  ji ,
j 1
(6)
l
– номер плоскости переизлучения.
Формы записи (1) и (6) эквивалентны. В правой части (6) выделены две
суммы, ответственные за переизлучение в обратном направлении и прямом.
Решить систему уравнений (6) можно несколькими методами. При условии,
что p  1 целесообразно обратиться к дифференциальному и интегральному
исчислению. Кратко изложим последовательность действий при решении системы
уравнений этим методом. Из системы уравнений (6) путем замены индексов
l  l  1 и l  l  1 находим следующие приращения.
El 1 (zl 1 )  El (zl )   El (zl ),
El (zl )  El 1 (zl 1 )   El (zl ),
1
El (zl )   El (zl )   El (zl ) ,
2
2 El (zl )   El (zl )   El (zl ).


С использованием преобразований
p l 1

j 1
E l 1 j 
p l 1

j 1
E l 1 j 
l
 E l 1 j 
ji
ji

 E l 1 j 
j 1
p l
 E l  j  ji
j 1
ji

i
p l
 E l  j  ji
j 1

j 1
l2
 i
i
l 1
 E l  j  ji
j 1
ji

 i
l 1
j 1
(8)
 E l 1 ,
 E l  i ,
 E l  i ,
 E l  j  ji
(7)
 E l 1 .
(9)
средние значения (8) приращений амплитуд парциальных волн, вычисленные для
произвольной плоскости zl , имеют следующие функциональные зависимости.
r
r
E l (z l )  iSin( )
E 0  (l 1)i  iSin( )
1 r
1 r
r
 iSin( )
1 r
l 1
 El j 
ji
 E l  j  ji

j 1
(10)
,
j 1
2 El (zl )  2(1  Cos( ))El (zl )  2( i  1)
где zl  (l  1)z 
p l
l 1
L,
p 1
r
El (zl ),
1 r
L –толщина слоя cреды,
  k z .
(11)
(12)
kz  1, или   z . Поделим соотношение (11)
i
и, ограничиваясь первыми членами разложения Sin( ), Cos( ) и  ,
В оптическом диапазоне
на (z )2
когда
1  Cos( )
2
2

k
,
2
(z )
Sin( )
i  1
 k,
 ik ,
z
z
(13)
получим
d2
E (z )  (k 2  2kk ' )El (zl )  0 ,
2 l l
(dz )
(14)

r 1
i r
1
2k

 2kN xy

 2kN .
где k '  i
1  r z
z
z x y z
(15)
N – объемная плотность атомов оптической среды.
Уравнения (6), (10) после ряда преобразований можно привести к
интегральному виду.
zp
E l (zl )  iAE 0  ikzl  ik '   ikzl  E j (z j ) 
ikz j
dz j 
zl
 ik ' 
ikzl
zl
 E j (z j ) 
 ikz j
(16)
dz j ,
z1
z
p
d
ikz
ikzl
 ikzl
E l ( zl )   AkE 0   kk ' 
E j ( z j )  j dz j 

dzl
z
l
 kk ' 
где
ikzl
zl
 E j (z j ) 
 ikz j
(17)
dz j ,
z1
A  k ' z , z  zl 1  zl .
(18)
Переход от сумм к интегралам и дифференциалам является всего лишь
математическим приемом, облегчающим решение поставленной задачи.
Сущность рассматриваемых явлений остается неизменной.
Интегральные уравнения (14), (16), (17) также можно получить
непосредственно из уравнения (1) следующим образом. Уравнение (1) можно
записать и преобразовать так
(1  r ) E l ( zl )  rE 0 
ikzl
 r
zl 1
ikzl
 E j (z j ) 
 ikz j
z j  z1
r
 ikzl
zp
 E j (z j ) 
ikz j
 rE l ( zl ) 
,
z j  z1
r
r ikzl zl
 ikz
ikzl
E l ( zl ) 
E0  

E j (z j )  j 

1 r
1 r
z j  z1
z
r ikzl p
ikz


E j (z j )  j .

1 r
z j  zl 1
(19)
Из уравнения (19) путем непосредственного перехода от сумм к интегралам и
используя соотношение (15) легко получается уравнение (16).
Уравнение (17) легко получить прямым дифференцированием уравнения
(16). Уравнение (14) получается прямым дифференцированием уравнения (17).
Результат решения дифференциального уравнения (14) следующий
E l ( zl ) 
k'
 ik 1 2 z l
k
C1
 C2 
ik 1
k'
zl
k
,
(20)
Для сокращения последующих записей введем обозначение
k'
u  1  2  1  4N .
k
(21)
Константы C1 и C2 непосредственно следуют из решения уравнения (14) и
обычно определяются из дополнительных начальных условий. В качестве этих
условий можно использовать уравнения (16) и (17). Последнее уравнение
выводится из первого путем его дифференцирования, поэтому достаточно
применить одно из них. Подставив решение (20) в уравнение (16), после
вычисления интегралов получим четыре группы слагаемых при различных
экспонентах
интегралы
L

zl

C1 
 ikuz j
ikz l ,   ikzl , ikuzl ,   ikuzl . Вычисляем
 C2 
ikuz j


ikz j
dz j 
следующие




C1
 ik (1 u ) L   ik (1 u )zl 
ik (1  u)

C2
 ik (1 u ) L   ik (1 u )zl ,
ik (1  u)
(22)
zl

0

C1 
 ikuz j
 C2 
ikuz j


 ikz j
dz j 





C1
  ik (1 u )zl  1 
ik (1  u)
C2
  ik (1 u )zl  1 .
ik (1  u)
Подставив соотношения (20), (22) в уравнение (16), получим
C1 
 iukz
 C2 
iukz
 iAE 0 
  ik (1 u )z )  ik '   ikz
ik '  ikz
ik '  ikz
ikz
 ik ' 
 ikz
C1
( ik (1 u ) L 
ik (1  u)
C2
( ik (1 u ) L   ik (1 u )z ) 
ik (1  u)
C1
(  ik (1 u )z  1) 
ik (1  u)
C2
(  ik (1 u )z  1).
ik (1  u)
(23)
Из уравнения (23) можно выделить группу слагаемых при экспонентах. Каждая
группа независимо от других групп слагаемых должна равняться нулю, т.к.
уравнение (23) должно выполняться при любых значениях координаты плоскости
zl
в пределах
0  zl  L . В результате будет получено четыре уравнения.
k'
k'

)  0,
k (1  u) k (1  u)
k'
k'
ikuz C2 (1 

)  0,
k (1  u) k (1  u)
C1k '
C2 k '
ikz (iAE0 

)0
k (1  u) k (1  u)
C k ' ik (1 u) L
C2 k ' ik (1 u) L
 ikz ( 1



)0
k (1  u)
k (1  u)
  ikuz C1 (1 
(24)
(25)
(26)
Из первых двух равенств (24) находится соотношение, тождественное равенству
(21)
k ' u2  1

k
2
Из равенства (26) находится соотношение между константами
C2   C1
1  u 2ikuL

1 u
(27)
Подставив функциональную зависимость (27) в соотношение (25), получим
равенство для константы C1
(1  u) iAE 0 ikuL  C1
k '  1  u ikuL 1  u ikuL 


 0



k 1 u
1 u
(28)
Из соотношений (27), (28) находим значения констант
k
1  u ikuL
C1  iA E0

,
k'
D
k
1  u ikuL
C2  iA E0

,
k'
D
(29)
(30)
здесь введено обозначение для знаменателя
D
1  u ikuL 1  u ikuL


 .
1 u
1 u
(31)
Подставив значения констант (29), (30) в функциональную зависимость (20) и
воспользовавшись соотношениями (15), (18), (21) получим функциональную
зависимость для амплитуд парциальных волн. Источниками этих электрических
волн являются совокупность диполей, расположенных в одной плоскости с
координатой zl . Значения амплитуд являются функциями параметров, которые
характеризуют всю систему переизлучения в целом :
z  zl 1  zl ,
L  z p  z1 , u , k , E0 .
ikz
El (zl )  E0
(1  u ) iku ( L  zl )  (1  u ) iku ( L  zl )
D


. (32)
Волны от источника переизлучения распространяются в двух противоположных
направлениях со скоростью света в вакууме. Электрическое поле, создаваемое
парциальной волной от отдельного источника с координатой
во времени имеет следующую зависимость.
zl
в пространстве и
El  El (zl )   ikz l   i ( t  kz ) , z  zl ,
El  El ( zl ) ikzl  i ( t  kz ) , z  zl ,
(33)
(34)
где значения E l ( zl ) следует подставить из соотношения (32).
При соблюдении условий установившегося процесса, в котором учтены все
обратные связи между условными плоскими источниками переизлучения, вклад в
суммарное (измеряемое) электрическое поле волны прямого направления вносят
источники, расположенные слева от плоскости наблюдения, а вклад в суммарное
электрическое поле волны обратного направления вносят источники,
расположенные справа от плоскости наблюдения. Это утверждение является
аксиомой, не требующей доказательства.
zl
i k ( zl  z j )

i k zl
E l  E 0

E j (z j ) 
,
(35)
z j  z1

E l1 
zp
 E j (z j ) 
z j  zl 1
i k ( z j  zl 1 )
,
где значения амплитуд волн источников переизлучения
(36)
E j ( z j ) необходимо
брать из формулы (32).
С вычислением суммарных измеряемых полей (35), (36) будет решена
задача микроструктурной модели взаимодействия волн электрического поля со
средой, представленной в виде пластины толщиной L для случая нормального
падения электрической волны на пластину. Это будет сделано в следующей
статье.
Литература
1. Косинский Ю.И., Принцип переизлучения электрического поля на
диполях в стационарных электрических явлениях, а также магнитных
и
оптических явлениях, 11, (2002).
2. Косинский Ю.И., Резонаторная система из переизлучающих плоских
решеток, в узлах которых расположены диполи, 1–8, (2002).
3. М.Борн, Э. Вольф, Основы оптики, 125–126, (1970).
Download