Математические
структуры и моделирование 2014. №4(32). С. 6-12
УДК 621.315
МОДИФИЦИРОВАННЫЕ ФОРМУЛЫ ОПТИЧЕСКОГО ПОКАЗАТЕЛЯ ПРЕЛОМЛЕНИЯ
С.Ю. Ланина
к.ф.-м.н., сотрудник лаборатории моделирования, обработки и управления информацией, e-mail: [email protected] В.В. Еремина
к.ф.-м.н., доцент факультета математики и компьютерных наук, e-mail: [email protected] О.С. Косолапова аспирант, e-mail: [email protected]
Федеральное государственное бюджетное образовательное учреждение «Амурский государственный университет», г. Благовещенск
Аннотация. Рассматривается вывод формул для вычисления оптического показателя преломления и оптического показателя поглощения, проверяется адекватность полученных зависимостей реальным спектрам рассматриваемых характеристик.
Ключевые слова: оптический показатель преломления, оптический показатель поглощения.
Взаимосвязь химического состава диэлектрика, его структуры с его поляризационными свойствами обычно выявляются за счёт анализа частотных спектров, отображающих преломление или поглощение им светового излучения. Кроме того, оптические постоянные в области собственного поглощения вещества являются фундаментальными константами, характеризующими целый ряд физико-химических параметров вещества. Измерение оптических характеристик материалов в высокочастотной области электромагнитного спектра позволяет рассчитать его характерные диэлектрические параметры, экспериментальное определение которых другими способами в настоящий момент технически не реализуемо. В свою очередь, актуальна и обратная задача, связанная с прогнозированием оптических характеристик перспективных материалов, проводимым на базе имитационного моделирования их поляризационных спектров.
Как известно [1], при рассмотрении взаимодействия вещества с переменным электромагнитным полем выражения комплексной диэлектрической проницаемости образца и его комплексного показателя преломления n(jw) принимают вид подобающих комплексных функций:
n{ju) = п(ш)+ jxM; £{ju) = £'{и)+ j£"{u),
где н(ш) и х(^) — соответственно, частотные зависимости показателя преломления и показателя поглощения.
В свою очередь, принимая во внимание соотношение п2 = е, частотные зависимости разбираемых функциональных зависимостей общепринято связывать между собой на базе непосредственно вытекающего выражения:
е"М + зе"(ш) = п2М - Х2(^) + 32п(ш)х(ш). (2)
Или
е'М = п2М - х2М; (3)
е"(ш) = 2п(ш)х(ш).
Рассмотрим математическое описание (3) в качестве системы линейных алгебраических уравнений, отвечающих частному виду функций п(ш) и х(^) при скалярном значении частоты, характерному для действия статического внешнего поля. В этом случае, исключение переменной х, даёт уравнение
п4 - е'п2 - 0, 25е"2 = 0, (4)
решение которого даёт два вещественных и два мнимых корня:
^ - Те'2 + е"2 , е' + Те'2 + е''2 п1,2 = ±\1 -2-; п3,4 = -2-. (5)
Поскольку исследуемая характеристика, по своему определению, описывает численную функцию вещественной части комплексной характеристики п(]ш), то из полученных решений следует отбросить мнимые корни п1,2. Рассматривая оставшуюся пару вещественных корней и принимая во внимание физическую сущность рассматриваемой функции, получим, что практической ценностью обладает единственный (положительный) корень.
Следовательно, частотная зависимость оптического показателя преломления может быть представлена в виде:
( ) Уе'2Н+ е''2Н+ е'Н
пМ = \1 --2-. (6)
Действуя аналогичным образом, из исходной системы уравнений исключим переменную п:
X4 + е'х2 - 0, 25е''2 = 0, ^
/-т-vёfТP'f 1-е' + те^Ге2 (7)
х1,2 = ±у-2-' хз,4 = ±у-2-.
Рассуждая аналогично, получим, что частотная зависимость оптического показателя поглощения имеет вид:
х Уе'2М + е'2М- е'М. (8)
Учитывая, что уравнения (6) и (8), описывающие частотные характеристики комплексного показателя преломления, включают в себя выражение амплитудной частотной характеристики £дт(^) комплексной диэлектрической проницаемости, то они могут быть эквивалентно представлены в виде:
= у^Щ+Щ (9)
хМ = у . (10)
Таким образом, принимая во внимание структуры частотных зависимостей оптических показателей преломления и поглощения представленного типа, можно отметить, что они характеризуют взаимосвязь частотных спектров реальных составляющих комплексной диэлектрической проницаемости материала.
С целью проведения качественной оценки практической эффективности рассматриваемых уравнений был выполнен ряд вычислительных экспериментов, общим объектом которых послужил гипотетический диэлектрик, состоящий из частиц одной и той же разновидности. Разберём выявленные с их помощью особенности внешнего вида имитационных графиков и(ш) и х(ш), достаточно объективно характеризующих общие особенности изучаемых спектральных зависимостей применительно к раздельному рассмотрению упругих и релаксационных колебаний частиц.
Первая часть проведённого исследования заключалась в расчётах спектра и(ш) для условия существования упругих колебаний частиц, коэффициент затухания которых полагался в десять раз меньшим частоты их собственных колебаний ш0. Вторая часть эксперимента представляла моделирование спектра и(ш) для релаксационных колебаний тех же самых частиц, коэффициент затухания которых задавался в десять раз больше частот их собственных колебаний. Полученные имитационные спектры изображены на рис. 1.
Анализ графиков показал, что графический образ частотного спектра и(ш), имитируемого для упругих колебаний частиц, существенно отличается от его реальных аналогов, которые чаще всего представляют собой резонансный всплеск с острым углом в нижней части. В свою очередь, спектр, рассчитанный на базе использования формулы (6), обладает явно выраженным горизонтальным участком, находящимся вблизи частоты собственных колебаний (рис. 1, а). При этом для моделирования спектра и(ш) в случае релаксационных колебаний тех же самых частиц подобная вычислительная аномалия отсутствует, а полученные графики хорошо соответствуют данным физического эксперимента в рамках их оценки с качественной точки зрения (рис. 1, б).
Причина полученного несовпадения может заключаться в следующем: с одной стороны, вещественная частотная характеристика комплексной диэлектрической проницаемости разбираемого образца, детализируемая в области электромагнитного резонанса его частиц для случая их упругих колебаний, обладает как положительными, так и отрицательными значениями; с другой -
Рис. 1. Расчетные спектры оптического показателя преломления: а) — упругие колебания
частиц; б) — релаксационные колебания
его мнимая частотная характеристика имеет практически нулевые значения во всем диапазоне частот за исключением узкого энергетического выброса с максимумом на резонансной частоте. Следовательно, присутствие горизонтального участка, отмеченного выше на спектре и(ш), оказывается вызванным тем, что при совместном выполнении условий е (ш) ~ 0 и е (ш) < 0 зависимость вида (6) даёт приблизительно нулевые значения, обусловленные взаимным подавлением текущих величин е (ш):
л/е'2(ш) + 0 - е'(ш) = 0. (11)
Сложившаяся ситуация может быть исправлена за счёт использования только абсолютных величин вещественной частотной характеристики е(]ш), т.е. видоизменения структуры традиционного выражения и(ш) к его альтернативному представлению в виде:
и(ш) = М Уе,2(ш)+ ^ш' + КМ. (12)
Результаты имитационного моделирования частотной зависимости оптического показателя преломления, проведённого на базе использования модифицированной формулы (12), изображены на рис. 2. Объективная оценка вытекающих изменений внешнего вида исследуемого спектрального графика показывает, что принятие предлагаемых авторами мер, направленных на его качественное улучшение, оказывается весьма эффективным. Следовательно, применение уравнения (12) для проведения практических расчётов частотных зависимостей и(ш) позволяет имитировать оптические спектры материалов, наиболее адекватные их физически измеряемым свойствам.
Аналогичным образом были построены имитационные спектры, полученные с помощью уравнения (8) и изображённые на рис. 3. Необходимо отметить, что физические измерения разбираемых зависимостей показывают, что внешний
Рис. 2. Модифицированный расчётный спектр п(ш)
Рис. 3. Расчётные спектры оптического показателя поглощения: а) упругие колебания частиц;
б) релаксационные колебания
вид их графиков представляет совокупность энергетических выбросов, симметричных относительно своих вертикальных осей.
Анализ представленных графиков показывает, что расчёт оптических характеристик релаксационных колебаний даёт изображение полосы поглощения, качественно соответствующее их реальной картине. В свою очередь, при рассмотрении упругих колебаний частиц имитационный спектр опять принимает неестественный вид.
Причина очередной вычислительной аномалии, также как и в предыдущем случае, может быть связана с фактом существования отрицательных значений е (ш), характерных для упругих поляризационных процессов, и устранена за счёт использования только её абсолютных величин. Для этого уравнение частотной зависимости х(ш), обладающее традиционной структурой (8), должно быть модифицировано к виду:
Рис. 4. Модифицированный расчётный спектр х(ш)
хм = 1/ ^<*" + . (13)
Результаты имитационного моделирования частотной зависимости оптического показателя поглощения, проведённого на базе предлагаемой формулы (13), изображены на рис. 4.
Оценка полученных изменений исследуемого графика показывает качественное улучшение его вида. Следовательно, модифицированное уравнение (13) позволяет имитировать оптические спектры показателя поглощения материалов, наиболее адекватные их физически измеряемым свойствам.
В заключении отметим, что полученные результаты позволяют обоснованно сформулировать следующие выводы: частотные зависимости оптических показателей преломления и поглощения светового излучения диэлектрическим материалом целиком определяются взаимосвязью величин вещественной и амплитудной частотных характеристик его комплексной диэлектрической проницаемости; максимумы полос поглощения соответствуют резонансным частотам, величины которых меньше значений частот собственных колебаний частиц. При этом их физически наблюдаемая разница особенно ярко выражается в релаксационных процессах.
Литература
1. Юхневич Г.В. Инфракрасная спектроскопия воды. М.: Наука, 1973.
2. Золотарев В.М., Морозов В.Н., Смирнова Е.В. Оптические постоянные природных и технических сред. Справочник. Л.: Химия. 1984.
3. Еремина В.В. Особенности поляризационных спектров конденсированных диэлектрических сред. I // Информатика и системы управления. 2009. № 2(20). С. 60-65.
12 C.fö. flaHHHa u gp. MogH$H^poBaHHbie fyopMyxu onTunecKoro.
MODIFIED FORMULA OPTICAL REFRACTIVE INDEX
S.Y. Lanina
Researcher of the Laboratory of modeling, information processing and management, candidate of physical and mathematical sciences, e-mail: [email protected]
V.V. Eremina
Associate Professor of Information and Control Systems of the Faculty of Mathematics and Computer Science, PhD in Physics and Mathematics, e-mail: [email protected]
O.S. Kosolapova
post-graduate student, Department of Information and councils-Barking Systems Faculty of Mathematics and Informatics, e-mail: [email protected]
Amur State University, Blagoveshchensk
Abstract. We consider the derivation of the formulas for the calculation of the optical refractive index and optical absorption coefficient, and verify their adequacy to the real spectra of the characteristics.
Keywords: optical refractive index, optical absorption coefficient.