НЕЛИНЕЙНАЯ ОПТИКА
  
Содержание: 
Введение 
1.                 Поляризация диэлектрика в постоянном электрическом поле 
2.                 Поляризация диэлектрика в световом поле 
3.                 Нелинейное взаимодействие электромагнитных волн 
4.                 Генерация второй гармоники (ГВГ) 
5.                 Фазовый синхронизм при генерации второй гармоники Вынужденное комбинационное рассеяние (ВКР) Макроскопическая теория ВКР. Стоксово рассеяние. 
8.                 Антистоксово рассеяние Самофокусировка света 
Введение 
Появление мощных источников когерентного светового излучения (лазеров) привело к возникновению нового направления исследований, связанного с изучением взаимодействия мощных когерентных потоков электромагнитного излучения с веществом, получившего название "Нелинейная оптика" 
Оптические эффекты, характер которых зависит от интенсивности излучения, называют нелинейными, а область оптики, изучающая нелинейные оптические эффекты (оптика мощных световых потоков) - нелинейной оптикой 
До появления лазеров число нелинейных оптических эффектов можно было перечислить по пальцам. Существовавшие до лазеров источники давали световые волны слишком малой интенсивности и, как следствие, этого большинство наблюдаемых оптических эффектов не зависело от интенсивностей волн. Только после появления лазеров - источников света, позволяющих получить световые волны с напряженностями полей 10 7 ... 10 9 В/см и выше, т.е. поля, сравнимые с внутриатомными - нелинейные явления в оптике стали предметом пристального изучения 
Поляризация диэлектрика в постоянном электрическом поле 
Любой электромагнитный процесс в среде описывается уравнениями Максвелла: 
	(1) 	
где E ( r ,t ), H ( r ,t ) - вектора напряженностей электрического и магнитного полей в точке r в момент t ; r - плотность зарядов; j - плотность тока; D и B - вектора электрической и магнитной индукции 
Однако этих уравнений недостаточно для решения электромагнитной задачи, необходимы материальные уравнения, устанавливающие дополнительные связи между указанными векторами 
D = e 0 E + P' ,  B = m 0 H + M' ,  j = s E 	(2) 	
Уравнения (2) устанавливают связь между вектором макроскопической поляризации среды P' , вектором макроскопической намагниченности среды M' и векторами D и B , а также между плотностью тока j и напряженностью электрического поля E . Далее не будем учитывать магнитные свойства среды 
В изотропном случае макроскопическая поляризация среды зависит от напряженности электрического поля E . Коэффициентом пропорциональности в такой зависимости является диэлектрическая восприимчивость среды c ( E ), которая в общем случае тоже зависит от E . Если учесть эту зависимость, то для уравнения (2) получим: 
P' = e 0 c ( E ) E    (2) :   D =(1+ c ( E )) e 0 E = e ( E ) E 	(3) 	
Величина e ( E )=1+ c ( E ) называется диэлектрической проницаемостью 
В слабых полях восприимчивость среды (и диэл . проницаемость) - константа, не зависящая от напряженности электрического поля. Следовательно, реакция среды на внешнее поле - линейная: 
P' = e 0 c 0 E ,  D =(1+ c 0 ) e 0 E = e 0 e E ,  e =1+ c 0 	(3a) 	
Нелинейные эффекты проявляются лишь тогда, когда поля достаточно сильны и величины c и e уже нельзя считать не зависимыми от напряженности поля. Чтобы проиллюстрировать появление нелинейной зависимости величин c и e , вычислим их в рамках простой классической задачи. Рассмотрим газ, состоящий из атомов (два точечных заряда: ядро и электрон) без постоянного электрического дипольного момента. В отсутствие внешнего поля положение точечных зарядов совпадает. Поместим его в постоянное электрическое поле. Заряды в каждом атоме сместятся на некоторое расстояние. Для простоты будем считать, что смещение электрона совпадает с направлением внешнего электрического поля. Тогда можно не учитывать векторного характера величин, входящих в задачу, и оперировать скалярами. Таким образом, атомы приобретут дипольный момент d = e r 
Если было N атомов, то макроскопическая поляризация 
P' = N d = N e r 	(4) 	
На электрон действуют две силы: одна - действие электрического поля - F E = e E , а вторая - упругая - возвращает электрон в прежнее положение F У = - k r - q r 3 (эта сила в общем случае нелинейно зависит от смещения электрона). Приравняем их и получим уравнение для определения смещения электрона во внешнем поле 
e E = k r + q r 3 ; 	(5) 	
Из (4) выражаем r и подставляем в (5) и получаем нелинейное уравнение для поляризации: 
	(6) 	
Решим его относительно P', считая член с P' 3 малым. Пусть P' = P' 0 + P' 1 (два порядка малости), тогда, подставив их в (6), получим два уравнения (одно для членов нулевого порядка малости, другое - для членов первого порядка малости) и решим их 
	(7) 	
Сравнив полученное решение с (3), получаем 
Т.е. восприимчивость является нелинейной функцией напряженности поля. Если же поле достаточно слабое (значительно меньше внутриатомного), то вторым членом можно пренебречь (это означает, что смещение r мало и в выражении для F У мы пренебрегаем членом qr 3 ) и восприимчивость становится постоянной величиной 
До этого мы рассматривали случай изотропной среды. Когда среда анизотропна , восприимчивость и проницаемость вместо скаляров становятся тензорами второго ранга, а связь между векторами P' , D , E имеет вид 
	(8) 	
d ij - единичный тензор 
Для декартовой системы координат: 
	(8a) 	
Зная параметры внешнего электрического поля и тензор восприимчивости для данного кристалла, обычно определяемый экспериментальными методами, можно всегда рассчитать его поляризацию 
Поляризация диэлектрика в световом поле 
Рассмотрим поляризацию диэлектрика в высокочастотном поле на той же простейшей модели газа. Поскольку напряженность электрического поля теперь зависит от времени, необходимо решать динамическую, а не статическую задачу для движения электрона. Уравнение движения электрона запишется в виде 
	(9) 	
где F T - сила трения пропорциональна скорости движения электрона (так мы учитываем возможные потери энергии электроном); F E - сила, действующая со стороны внешнего электрического поля; F У - упругая сила. Упругую силу возьмем в линейном приближении (для случая слабого поля): F У = - k r . Подставив в (9), получим: 
	(10) 	
Последнее уравнение получено заменой r на выражение через поляризацию из (4), за w 0 2 принято k / m 
Пусть поле меняется по гармоническому закону E( t ) = E 0 cos w t , тогда решение для поляризации будем искать в виде P' = P' 0 cos ( w t+ j ). Дифференцируя это выражение нужное число раз, подставим его в (10): 
( w 0 2 - w 2 ) P 0 ( cos w t cos j - sin w t sin j ) - 
- g з w P 0 ( cos w t sin j + sin w t cos j ) = e²N / m E 0 cos w t 	(10a) 	
Приравняем по отдельности члены при cos w t и sin w t нулю: 
- ( w 0 2 - w 2 ) sin j - g з w cos j = 0 
( w 0 2 - w 2 ) P 0 cos j - g з w P 0 sin j = e 2 N E 0 / m 	(11) 	
Из первого равенства определяем фазу поляризации 
	(12) 	
и подставив во второе, получим 
	(11a) 	
Очевидно, решение для поляризации имеет вид 
	(13) 	
Выводы: 
Поляризация меняется с той же частотой w , что и внешнее поле 
Амплитуда поляризации существенно зависит от соотношения частот w и w 0 
Если w = w 0 (резонанс), амплитуда максимальна; 
Вдали от резонанса | w - w 0 | >> g з 
В этом случае фаза поляризации близка к нулю (см. (12)). Тогда поляризация 
	(13a) 	
т.е. восприимчивость зависит от частоты 
В предельном случае постоянного поля для восприимчивости получаем вновь формулу аналогичную (7 а ) ( w =0 (13a)): 
До сих пор предполагалось, что на электрон действует поле малой напряженности. Мы брали F У = - k r (линейное приближение, пригодное для случая малого смещения электрона). Теперь будем считать, что напряженность светового поля и смещение электрона могут быть достаточно большими, и для упругой силы возьмем F У = - k r - q r 3 : 
	(14) 	
	(15) 	
Будем, как и раньше считать, что поле E( t ) меняется по гармоническому закону, рассматривая нерезонансный случай (| w - w 0 | >> g з ). Членами при g з и P' 3 пренебрегаем. Решение опять ищем в виде P'=P' 0 +P' 1 (два порядка малости), подставляем его в (15) и собирая отдельно члены нулевого и первого порядков малости, получаем: 
	(16) 	
	(17) 	
Первое уравнение мы уже решали, это решение вдали от резонанса - (13a). Подставляем его в (17): 
	(18) 	
Т.к. напряженность поля меняется по гармоническому закону, то 
E 3 (t) = 1 / 4 E 0 3 (3 cos w t + cos 3 w t) 	(19) 	
Уравнение (18) - это уравнение гармонического осциллятора, на который действует внешняя сила (правая часть уравнения), состоящая из двух компонент, одна из которых меняется с частотой w , а другая - с частотой 3 w . Поэтому решение будем искать в виде P' 1 =P' 1, w cos w t + P' 1,3 w cos 3 w t. Подставляя его в (18) и получаем: 
	(20) 	
	(21) 	
Объединяем (20-21) и получаем общее решение: 
P'= P' 0 + P' 1 = c(w ,E 0 ) E 0 cos w t + c(3w ,E 0 ) E 0 cos 3 w t 	(22) 	
где 
	(23) 	
Выводы: 
Поляризация в сильном световом поле является функцией не только частоты падающего излучения, но и его третьей гармоники. Известно, что заряд, совершающий гармоническое колебание с некоторой частотой, излучает монохроматическую электромагнитную волну той же частоты. Поэтому в рассмотренной задаче появляются две волны: одна с частотой w , другая - с частотой 3 w 
Таким образом, в рамках простейшей модели мы показали, каким образом из-за нелинейных свойств среды в сильном световом поле возникают высшие гармоники 
Нелинейное взаимодействие электромагнитных волн 
Тензор нелинейной восприимчивости 
Рассмотрим нелинейное взаимодействие двух электромагнитных полей. Одно из них, поляризованное вдоль j , описывается выражением: 
E j w 1 ( t ) = Re (E j w 1 exp i w 1 t) = 1 / 2 (E j w 1 exp i w 1 t + к.с. ), 	(1) 	
а второе, поляризованное в направлении k , - выражением 
E k w 2 (t) = Re(E k w 2 exp i w 2 t) 	
Если среда нелинейная, наличие этих двух полей может привести к появлению поляризации на частотах n w 1 +m w 2 , где n и m - целые числа. Записав i-компоненту поляризации на частоте w 3 = w 1 + w 2 в виде 
P i w 3 = w 1 + w 2 (t) = Re(P i w 3 exp i w 3 t), 	
определим тензор нелинейной восприимчивости (раньше мы использовали c ijk - тензор линейной восприимчивости) d ijk w 3 = w 1 + w 2 с помощью следующего соотношения для комплексных амплитуд 
	(2) 	
Подобным же образом вводим тензор восприимчивости на разностной частоте d ijk w 3 = w 1 - w 2 
	(3) 	
где согласно (1) E k - w 2 =(E k w 2 )* 
Рассмотрение взаимодействия электромагнитных полей начнем с записи уравнения Максвелла, выделив в явном виде поляризацию P : 

	(4) 	
Примечание: 
rot rot E = grad div E - С 2 E 	
Представив поляризацию в виде суммы линейного и нелинейного членов, перепишем первое уравнение 
	(5) 	
Возьмем ротор от обеих частей второго уравнения (4) и подставим rot H из (5) (см. тж. примечание), учитывая, что div E =0: 
	(6) 	
Дальнейший анализ проведем для одномерного случая ( / x= / y=0). За направление распространения берем ось Z. Ограничимся рассмотрением взаимодействия колебаний трех частот и соответствующие поля возьмем в виде бегущих плоских волн: 
E i w 1 ( z,t ) = 1 / 2 [E 1i ( z ) exp i ( w 1 t-k 1 z) + к.с. ], 
E k w 2 ( z,t ) = 1 / 2 [E 2k ( z ) exp i ( w 2 t-k 2 z) + к.с. ], 
E j w 3 ( z,t ) = 1 / 2 [E 3j ( z ) exp i ( w 3 t-k 3 z) + к.с. ], 	(7) 	
где ijk - декартовы координаты. Заметим, что при P нел =0 решение уравнения (6) дается выражениями (7) с амплитудами, не зависящими от z . В качестве примера запишем i-компоненту нелинейной поляризации на частоте w 1 = w 3 - w 2 . Согласно (3) и (7) она имеет вид 
	(7a) 	
Вернемся к уравнению (6). В одномерном случае 
	(8) 	
Дифференцируем и полагаем, что изменение комплексных амплитуд полей достаточно медленное, т.е 
	(9) 	
Аналогичные выражения можно вывести для С 2 E j w 3 ( z,t ) и С 2 E k w 2 ( z,t ). Подставляя (9) в (6) и используя соотношение / t=i w 1 получим волновое уравнение для E i w 1 ( z,t ): 
	(10) 	
Предполагаем, что при взаимодействии конечного числа полей уравнение (6) должно удовлетворяться по отдельности для компонент с различными частотами. Поставив (7а) и заметив, что w 1 2 m 0 e =k 1 2 , получим 
	(11) 	
или (считая s функцией частоты) 
	(11a) 	
и аналогично 
	(11b) 	
  
	(11c) 	
Эти уравнения мы применим в дальнейшем при рассмотрении ряда конкретных случаев 
Генерация второй гармоники (ГВГ) 
Первый эксперимент по генерации второй гармоники света был выполнен Франкеном в 1961 году. Луч рубинового лазера с l = 694,3 нм фокусировался на поверхность пластины из кристаллического кварца. Выходящее излучение анализировалось спектрометром. Было найдено, что в нем содержится компонента с удвоенной частотой (т.е. с l = 347,15 нм). Эффективность преобразования в первых экспериментах была порядка 10 -8 . Использование более эффективных материалов, увеличение мощности лазера, обеспечение условий фазового синхронизма позволили в последние годы довести коэффициент преобразования почти до единицы 

Рис.1. Схема первых экспериментов по ГВГ. 
1 - рубиновый лазер, 2 - фокусирующая линза, 3 - кварцевая пластинка, 
4 - коллиматорные линзы, 5 - призма, 6 - фотопластинка (экран). 
Цвета показаны условно. 	
Применим уравнения (11a-11c) для рассмотрения ГВГ. Это частный случай взаимодействия полей трех частот, когда две частоты w 1 и w 2 одинаковы, а w 3 = 2 w 1 . Следовательно, необходимо анализировать только два уравнения: первое (или второе) и последнее. В целях упрощения будем считать, что потери мощности входного луча ( w 1 ) за счет преобразования во вторую гармонику малы, т.е. dE 1i / dz » 0. Следовательно, можно рассматривать только последнее уравнение (11c). Если среда прозрачна на частоте w 3 , то s 3 =0 и 
	(12) 	
где w = w 1 = 1 / 2 w 3 ,  D k = k 3 ( j ) - k 1 ( i ) - k 1 ( k ) ,  а k 1 ( i ) - волновое число волны с частотой w 1 , поляризованной по оси i . Если E 3j (0) = 0, т.е вторая гармоника на входе отсутствует, и кристалл имеет длину l , решением (12) будет 
	(13) 	
или 
	(14) 	
где e¦e 3 . Чтобы получить выражение для мощности второй гармоники P 2 w на выходе, воспользуемся соотношением 
	(15) 	
где S - площадь поперечного сечения пучка. Приняв e 1 »e 3 »e 0 n 2 приходим к коэффициенту преобразования 
	(16) 	
Фазовый синхронизм при генерации второй гармоники 
Из (16) следует, что предпосылкой для эффективной ГВГ является выполнение условия D k = 0, или, поскольку w 3 = 2 w , а w 1 = w 2 = w , 
D k = k 2 w - 2 k w = 0    k 2 w = 2 k w 	(17) 	
Если D k ¦ 0, то волна удвоенной частоты, генерируемая в некоторой плоскости (z 1 ), дойдя до другой плоскости (z 2 ), окажется не в фазе с волной удвоенной частоты, генерируемой в этой плоскости. Результат интерференции таких волн представлен в (16) множителем ( 1 / 2 D k l ) -2 sin 2 ( 1 / 2 D k l ). Два соседних максимума этой интерференции удалены на расстояние, называемое "когерентной длиной": 
	(18) 	
Она является в сущности максимальной длиной кристалла, которую можно использовать для ГВГ. Показатель преломления, как правило, растет с увеличением частоты, так что 
D k = k 2 w - 2 k w = (2 w /c)(n 2 w - n w ) 	(19) 	
Здесь использовано k= w n / c . Когерентная длина выражается формулой 
	(20) 	
в которой l - длина волны падающего света 
Пример 
Если l = 1 мкм и n 2 w - n w = 0,01 , то l c = 100 мкм 
Увеличение l c от 100мкм до 2см согласно (16) влечет за собой возрастание мощности второй гармоники в 4·10 4 раз 
Способ, который широко применяется для обеспечения условий фазового синхронизма, заключается в использовании анизотропных кристаллов, обладающих естественным двулучепреломлением . Используя связь k w = w Ц me 0 n w , вместо условия (17) получим условие n 2 w = n w , т.е. коэффициенты преломления на основной частоте и на удвоенной должны совпадать. В материалах с нормальной дисперсией показатель преломления обыкновенной и необыкновенной волн, распространяющихся в данном направлении, растет с частотой. Т.е. удовлетворить условию равенства коэффициентов преломления невозможно, если волны частот w и 2 w принадлежат одному типу (обыкновенные или необыкновенные). Однако фазовый синхронизм может осуществляться благодаря использованию волн разных типов 
В качестве примера рассмотрим зависимость показателя преломления необыкновенной волны в одноосном кристалле от угла q между направлением распространения и оптической осью (осью Z) кристалла. Эта зависимость имеет вид 
	(21) 	
Если n e 2 w < n o w , то существует угол q синх , при котором n e 2 w ( q синх ) = n o w . Таким образом, если волна частоты w распространяется под углом q синх к оси и имеет поляризацию, отвечающую обыкновенному лучу, то волна удвоенной частоты, возбуждаясь в том же направлении, будут обладать поляризацией необыкновенного луча. (См. рис.2) 

Рис.2. Поверхности показателей преломления для обыкновенного и необыкновенного лучей в отрицательном одноосном кристалле. 	
Угол q определяется пересечением сферы, представляющей собой поверхность показателей преломления для обыкновенного луча частоты w (желтая сфера) с эллипсоидом показателей преломления необыкновенного луча частоты 2 w (розовый эллипсоид). В случае отрицательного одноосного кристалла ( n e w < n o w ), угол, удовлетворяющий условию n e 2 w ( q синх ) = n o w , определяется так 
	(22) 	
откуда 
	(23) 	
Пример 
Генерация второй гармоники в кристалле KDP. Исходное излучение - рубиновый лазер ( l = 694,3 нм). Значения показателей преломления: n e w = 1,466, n e 2 w = 1,487, n o w = 1,506, n o 2 w = 1,534. Угол синхронизма, вычисленный по формуле (23), равен q синх = 50,4° 

Вынужденное комбинационное рассеяние (ВКР) 
Комбинационное или рамановское рассеяние света давно используется для изучения колебательных спектров молекул и оптической ветви колебаний кристаллических решеток. Ячейка, содержащая исследуемое вещество (жидкость, газ или кристалл), облучается светом с узкой спектральной линией. Спектральный анализ рассеянного излучения обнаруживает присутствие линий, смещенных вниз по частоте на величину, равную колебательным частотам облучаемого образца. Этот тип рассеяния называется стоксовым рассеянием 
В спектре рассеянного излучения присутствуют также частоты, равные сумме частоты падающего излучения и колебательных частот вещества. Это так называемое антистоксово рассеяние , интенсивность которого на несколько порядков меньше интенсивности стоксовой компоненты 
Указанные два типа рассеяния поясняются на рис.3 
( a ) 		Стоксово рассеяние, при котором поглощается лазерный фотон и вместе со стоксовым фотоном на частоте w c = w л - w u возникает квант колебаний молекулы ( u = 1) 	
( b ) 		Антистоксово рассеяние, при котором поглощаются лазерный фотон и колебательный квант, а испускается фотон на частоте w ac = w л + w u 	
( c ) 		Процесс поглощения фотонов частоты w c = w л - w u , стимулированный наличием лазерного излучения частоты w л 	
  	  	  	
Рис.3. Переходы при вынужденном комбинационном рассеянии. 	
Т.к. антистоксово излучение определяется молекулами, находящимися в возбужденном состоянии, то его интенсивность ниже интенсивности стоксова излучения на величину множителя exp (-- w u / kT ). На рис.3 ( c ) представлен также обратный процесс, при котором фотон стоксовой частоты поглощается 
До недавнего времени в спектроскопии комбинационного рассеяния применялись интенсивные источники некогерентного излучения (например, ртутные лампы). В последнее время когерентные лазерные источники вытеснили ртутную лампу 
Характерные частоты колебаний атомных групп в молекулах ([2] с.370) 
Частота, см -1 	Колеблющаяся 
атомная группа 	Тип соединения 	
445-550 	S-S 	Алифатические дисульфиды 	
510-594 	C-Br 	Алифатические соединения 	
750-850 		Парапроизводные бензола 	
884-899 		Циклопентан и монопроизводные 	
939-1005 		Циклобутан и производные 	
990-1050 		Бензол и одно- и трехзамещенные бензолы 	
» 1340 		Ароматические соединения 	
» 1380 		Нафталин и производные 	
» 1630 	C=N 	Ароматические соединения 	
1654-1670 	C=N 	Алифатические соединения 	
1974-2260 	C ¦ C 	Алифатические соединения 	
2150-2245 	C ¦ N 	Нитрилы 	
4160 	H-H 	H 2 	
Характерно, что частоты мало меняются от соединения к соединению 
  
Если среду, способную к комбинационному рассеянию, поместить в оптический резонатор, то при наличии поля лазерной накачки усиление стоксовой компоненты способно скомпенсировать потери, и на частоте w c возникает генерация. Генерация при ВКР представляет собой практический способ преобразования излучения импульсных лазеров (например, лазера на неодимовом стекле) в когерентное излучение, сдвинутое по частоте на колебательную частоту вещества 
Эксперименты по исследованию влияния интенсивности лазерной накачки на интенсивность стоксовой компоненты показали, что по достижении некоторой критической интенсивности накачки интенсивность стоксовой компоненты резко возрастает, а затем идет насыщение (см. рис.4) 

Рис.4. Зависимость интенсивности генерации стоксовой компоненты от интенсивности накачки лазера 	

Макроскопическая теория ВКР. Стоксово рассеяние. 
В экспериментах по ВКР было обнаружено, что выходное излучение содержит несколько стоксовых ( w л - w u ), ( w л - 2 w u ), ... и антистоксовых ( w л + w u ), ( w л + 2 w u ), ... компонент. Из рис. 3 видно, что процесс излучения стоксовой компоненты приводит к увеличению населенности колебательного уровня (u=1) , поэтому становится возможным излучение на антистоксовой частоте. Стоксова ( w c ) и антистоксова ( w ас ) компоненты могут, в свою очередь, служить исходным излучением, генерирующим частоты w с - w u = w л - 2 w u и w ас + w u = w л + 2 w u . Аналогично можно объяснить появление комбинационных частот более высоких порядков 
Чтобы пояснить основные особенности возникновения ВКР, получим условие усиления или генерации на первой стоксовой частоте w с = w л - w u , т.к. первоначально может усиливаться только эта компонента. Для возникновения других спектральных компонент требуется либо наличие молекул в возбужденном состоянии, либо присутствие стоксовой компоненты первого порядка 
Для анализа используется такая модель: рассеивающая среда состоит из N независимых осцилляторов (т.е. ансамбль осцилляторов не поддерживает волновое движение с отличной от нуля групповой скоростью), каждый характеризуется своим положением z (одномерный случай / x= / y=0) и нормальной колебательной координатой X( z,t ). Уравнение движения для осциллятора имеет вид 
	(1) 	
где Г - постоянная затухания, выбранная так, что наблюдаемая ширина линии спонтанного комбинационного рассеяния равна Dn= G / 2p ; w u - резонансная частота колебаний молекулы в отсутствие затухания; m - масса; F( z,t ) - возбуждающая сила 
Возбуждающую силу можно получить, рассматривая электромагнитную энергию в молекулярной среде. Плотность энергии, запасенной в электрическом поле E = 1 / 2 e E 2 при использовании равенства 
e = e 0 (1 + N a) = e 0 {1 + N [ a 0 + ( a/ X) 0 X]} 	(2) 	
может быть записана в виде 
E = 1 / 2 e 0 {1 + N [ a 0 + ( a/ X) 0 X]} E 2 	(3) 	
Сила, действующая на единицу объема поляризуемой среды, равна e / X, откуда делением на N получаем силу, действующую на один осциллятор 
F( z,t )= 1 / 2 e 0 ( a/ X) 0 <E( z,t )> 2 	(4) 	
<E( z,t )> означает усреднение за несколько колебаний, предпринимаемое потому, что молекула неспособна реагировать на эти колебания. Из (4) видно, что при отличной от нуля дифференциальной поляризуемости ( a / X) 0 колебания молекул могут возбуждаться электрическим полем 
Дальнейшая задача - показать, как колебания молекул воздействуют на электромагнитное поле. В соответствии с (2) колебания молекул с частотой w u вызывают модуляцию диэлектрической проницаемости с той же частотой. Это приводит к фазовой модуляции поля излучения (появляются боковые составляющие, смещенные на w u друг от друга). Т.е. происходит обмен энергией между электромагнитными полями различных частот., разделенных интервалами, кратными w u 
Полное поле является суммой лазерного ( w 2 ) и стоксова ( w 1 ) полей: 
E( z,t ) = 1 / 2 E 1 ( z ) exp i w 1 t + 1 / 2 E 2 ( z ) exp i w 2 t + к.с. 	(5) 	
<E( z,t )> 2 = 1 / 4 E 2 (z) E 1 * (z) exp i ( w 2 -w 1 )t + к . с . 	(6) 	
(6) (4) (1) 	  	
	(7) 	
Здесь использованы соотношения / t=i w и 
X( z,t ) = 1 / 2 X( z ) exp i w t + к.с. 	(8) 	
Из (7) следует, что на частоте w=w 2 -w 1 молекулярные колебания имеют комплексную амплитуду 
	(9) 	
Поляризация, наведенная полем частоты w 1 , имеет вид 
P = e 0 N a ( z,t ) E( z,t ) = e 0 N [ a 0 + ( a/ X) 0 X( z,t )] E( z,t ). 	(10) 	
Используя (5), (9) и (10) для нелинейная поляризации (второй член поляризации, пропорциональный X E) получаем: 
	(11) 	
Примечание: 	
Осуществив умножение в формуле (11) (см. тж. примечание), получим составляющие поляризации, осциллирующие с частотами w 1 , w 2 , 2w 1 -w 2 и 2w 2 -w 1 . Рассмотрим сначала составляющую нелинейной поляризации, имеющую частоту w 1 : 

P нел w 1 ( z,t ) = 1 / 2 P нел w 1 ( z ) exp i w 1 t + к.с. , 	(12) 	
где 
	(13) 	
Коэффициент пропорциональности между полем и поляризацией представляет собой восприимчивость. Нелинейная комбинационная восприимчивость подобно линейной восприимчивости имеет лоренцеву форму линии 
Форма линии стоксова рассеяния имеет вид 
	(13a) 	
Антистоксово рассеяние 
Антистоксово излучение на частоте w 3 = w 2 + w u является результатом комбинационного рассеяния света молекулой, находящейся в возбужденном колебательном состоянии (u=1) . При классическом подходе к задаче мы должны найти поляризацию на w 3 , наведенную электрическим полем: 
E( z,t ) = 1 / 2 [E 1 ( z ) exp i w 1 t + E 2 ( z ) exp i w 2 t + E 3 ( z ) exp i w 3 t + к.с. ], 	(14) 	
где w 3 - w 2 = w 2 - w 1 
В выражении для поляризации по аналогии с (11) найдем член, соответствующий возбуждению молекулярных колебаний силой, пропорциональной E 3 E 2 * . Из (13) заменой частот и индексов у E получим 
	(15) 	
Важно, что мнимые части (13) и (15) имеют разные знаки, поэтому 
антистоксова волна, распространяясь в среде, активной в комбинационном отношении, в присутствии излучения лазера ( w 2 ), но в отсутствии стоксова излучения (w 1 = w 2 - w u ) будет затухать 
Существует, однако, еще одна компонента поляризации на частоте w 2 : 
P нел w 3 (z) ~ E 2 E 2 E* 1 exp [ i (2w 2 -w 1 ) t] 	(16) 	
Она не содержит E 3 и может рассматриваться как верхняя боковая частота [w 2 + (w 2 - w 1 )] в спектре модулированных колебаний диэлектрической проницаемости с несущей w 2 и модулирующей w 2 - w 1 частотами. Эта компонента является источником излучения с частотой w 3 
Если дополнить (16) пространственной зависимостью поляризации, то 
P нел w 3 (z) ~ E 2 E 2 E* 1 exp [- i (2 k 2 - k 1 ) r ] 	(17) 	
Этому члену соответствует поле E 3 exp(-i k 3 r ), причем 
k 3 = 2 k 2 - k 1 	(18) 	
Следовательно, антистоксова волна может излучаться только в направлениях, удовлетворяющих условию (18). См. рис.5. А так как | k i | = w i n i / c , то антистоксова компонента распространяется в направлениях, составляющих коническую поверхность с половинным углом b при вершине и осью лазерного луча 
Реальная ситуация сложнее. Помимо наличия стоксовых и антистоксовых компонент высоких порядков, имеет место отклонение от направлений, рассчитанных по формуле (18) из-за эффекта самофокусировки 

Рис.5. Диаграмма для определения направления распространения антистоксова излучения. 	

Рис.6. Схема эксперимента по изучению комбинационного рассеяния. 
1 - рубиновый лазер; 2 - линза; 3 - ячейка с бензолом; 4 - экран 
Цвета показаны условно. 	
Самофокусировка света 
Выше уже упоминалось, что ВКР в среде наступает только при превышении некоторого порога интенсивности электрического поля. Однако измеренная пороговая интенсивность часто оказывается ниже ожидаемой. Расхождения между теорией и экспериментом могут быть весьма значительными: в некоторых жидкостях соответствующие пороги отличаются в сотни и более раз, что обусловлено явлением самофокусировки. В таком случае диаметр пучка по мере распространения в среде уменьшается и на некотором расстоянии пучок собирается в "фокусе". В фокальной области плотности мощности лазерного излучения очень велики и могут привести к разрушению материала. Это явление имеет непосредственное отношение к импульсным лазерам с очень высокой мощностью излучения, поскольку разрушению может подвергаться и активный элемент лазера 
В первой лекции были выведены зависимости c (w ,E 0 ) и c (3w ,E 0 ) (формулы (23)), на основе которых можно записать: 
e общ = 1 + c + b E 2 , 	(19) 	
тогда 
n общ = Ц e общ ¦ n + n 2 E 2 , где n 2 = b / 2n 	(20) 	
Если n 2 >0, то в местах большой напряженности поля - показатель преломления больше. Т.е. в нелинейном материале сам пучок формирует положительную линзу. Это так называемая крупномасштабная самофокусировка . Существует также мелкомасштабная самофокусировка , обусловленная нарастанием возмущений в пучке в поле мощной световой волны 
  
На рисунках показано применение ВКР 

Рис.7a. КАРС спектроскопия. 	

Рис.7b. Многопроходные кюветы.