LASER POWER STABILIZATION METHOD

Abstract

Laser power stabilization method, based on thermal lens effect, is proposed. Output power vs. input power and stabilization ratio was computed.

Full Text

Введение. Постановка задачи Для задач контроля качества оптических волокон и метрологии требуются лазеры с высокой долговременной стабильностью мощности излучения. Существующие методы затрагивают только стабилизацию электрической мощности блока питания. Между тем изменение параметров окружающей среды (прежде всего температуры) приводит к изменению характеристик излучателя (например, за счет частичной разъюстировки), что ведет к уменьшению выходной мощности. В данной статье предлагается способ стабилизации мощности излучения лазера, основанный на явлении самодефокусировки лазерного пучка, проходящего через среду с отрицательной величиной dn/dT. Рис. 1. Принципиальная схема установки «Инфокоммуникационные технологии» Том 11, № 2, 2013 66 Петропавловский В.М., Топоркова Л.В. Принципиальная схема установки приведена на рис. 1. Лазерный пучок фокусируется линзой 2 в кювету 3, содержащую жидкость (например, ацетон). Частичное поглощение излучения приводит к выделению тепла и росту температуры жидкости. Происходит изменение показателя преломления жидкости Δη = дп/дТ-ΔΤ. Так как распределение интенсивности по сечению лазерного пучка неоднородно, то и изменение показателя преломления в разных точках будет различным. В случае одномодового лазера вблизи оси пучка изменение показателя преломления происходит по квадратичному закону. Такая структура аналогична рассеивающей линзе (так как дп/дТ < 0). Расходимость пучка, проходящего через эту линзу, увеличивается, поэтому увеличивается и ширина распределения интенсивности (диаметр) пучка в дальней зоне. Подобрав диаметр пучка в перетяжке после линзы 2 - w() , можно добиться постоянства интенсивности на оси пучка. Линза 4 коллимирует пучок, диафрагма 5 вырезает центральную часть пучка с постоянной мощностью. Фокусное расстояние тепловой линзы Разложим показатель преломления, распределение которого цилиндрически-симметрично, в ряд Макларена по радиальному смещению от оси: ґ*пЛ + ... (1) п(г) = и(0) + ri— \àr)r=о 1 + — 2 дг Л=о Для цилиндрически-симметричного распределения (дп дг)г=0 = 0. Луч, входящий в среду параллельно оси симметрии, будет идти по кривой, радиус кривизны которой рассчитан в [1] как - V · grad [ln η ], (2) где V - единичный вектор главной нормали к кривой. Предполагая, что градиент направлен по радиусу и что среда тонкая, (2) можно преобразовать к виду (3) R J _nL_r F ~ rR~ ҐЯ2 \ on dr (4) sr=0 Радиальное распределение показателя преломления можно выразить через радиальное распределение температуры, получаем 1 _ dn J~~ ~dT d2T dr2 (5) Λ=о то есть оптическая сила тепловой линзы зависит от распределения температуры в среде. В подавляющем большинстве сред (dn/dT) и (д2п/дг2) r = 0 отрицательны, так что тепловая линза рассеивает свет (F < 0). Распределение температуры в среде Функция, описывающая температуру в бесконечной цилиндрически-симметричной среде, в которой в момент времени t = 0 в бесконечномалом элементе цилиндра радиуса r мгновенно выделилась тепловая энергия [2]: G(r, r',t) = Q ΑπΌί exp 2 tl\ г + г ADt ґ , λ г г 2Dt (6) D = k/(pcp) - коэффициент температуропроводности; к - коэффициент теплопроводности; cp - удельная теплоемкость; p - плотность; /0(·) - модифицированная функция Бесселя; Q - интенсивность теплового импульса на единицу длины (скачок температуры, возникающий за счет выделения тепла). В случае, если источник тепла действует непрерывно (непрерывный лазер), температура на расстоянии г от оси в момент времени t вычисляется как результат интегрирования (6) по временам t', лежащим в интервале 0 < t' < t. В этом случае вместо количества мгновенно выделившейся тепловой энергии вводится скорость выделения тепла на единицу длины Q = od. Поэтому (6) можно переписать в виде: . , , 2nr'Qlрсп G(r,r',t') =- р x ΑπΌί' n dr xexp Ґ (г2 + /2^ h f V ADt' ) (7) Эффективное фокусное расстояние среды длиной L связано с радиусом кривизны луча R, вошедшего в среду на расстоянии r от оси симметрии В случае одномодового лазера с гауссовым профилем распределения интенсивности по сечению пучка рост температуры может быть найден как «Инфокоммуникационные технологии» Том 11, № 2, 2013 Петропавловский В.М., Топоркова Л.В. 67 T(r,t)= \df \dr' exp ' J AnDt (r2 + r ADt' '2 Λ (8) В этой формуле предполагается постоянство радиуса лазерного пучка w вдоль длины кюветы, то есть длина кюветы L должна быть меньше конфокального параметра пучка Ь0 = , λ - длина волны излучения. ^ Найдем теперь вторую производную температуры по радиусу при r = 0 и вычислим, используя (5) фокусное расстояние тепловой линзы. Интегрирование (8) по r дает: с» fd , 2ж г' г-ехр АжйГ (г2 + / ADt' ґ , λ rr 2D? X аР X—г exp / π w аР j2 \ V D w ! pcv = ( . (9) як w + 8Dt‘ 7 exp w2 +8Dt' Вычислить вторую производную по r проще до интегрирования по времени. Получим !L( дг2 аР D --Ï-Г ехр лк w +8Dt AaPD W w2 + 8 Dt' JJ Kk(w2+8Dt')2 ' Вычисление интеграла по t' дает d2T f AaPD (10) dr2 = f ι\2 лк{м> +8 Dt') аР dt'= 2якм> (l + tc/2t) (И) Параметр tc = w2/4D характеризует временной отклик среды на выделение тепла (время установления теплового равновесия). Подстановка (11) в (5) дает 2лк^{\.Л-іс!2і) F if) = при условии t >> t : aPL dn dT F = 2ккм? aPL dn dT. —Ç.C- P ’ 2 nkw aPL dn dT. .(12) Изменение диаметра пучка в дальней зоне. Мощность излучения на выходе стабилизатора Для гауссового пучка с радиусом перетяжки w0 радиус пучка W на экране в дальней зоне 0 2 7Щ λ (расстояние до перетяжки z » —= Ь0 ) может быть найден как λζ W = — (13) Если такой пучок с исходной мощностью Р пропустить через диафрагму радиусом а, то мощность пучка после диафрагмы будет РХ=Р ” ( 2 Λ ~ 1-ехр а 2 1 w ) (14) Прохождение пучка через линзу приводит к изменению радиуса кривизны волнового фронта R. Размер и положение перетяжки нового пучка меняются, следовательно, изменятся W и Р Для расчета параметров нового пучка воспользуемся формализмом комплексного параметра пуч- 1 1 а λ В - . В перетяжке q - мнимый ка q q R π\ν q = mw, . При распространении в среде длиной X q меняется по закону q2 = q1 + x. При прохождении через линзу радиус кривизны волнового фронта меняется по закону 1 1 1 То yR r Fj есть если отрицательная линза помещена в перетяжку пучка, то после линзы будет расходящийся пучок с R = F. В результате получим уравнение 1 іЛ гЛ πνη - +JC, (15) где X - расстояние до новой перетяжки; v0 - ее размер. Решая (15) получим Vo2=wn 1 + ґи Λ К 2 > Λ f F 1 + U'J \ b0 J 2 - (16) Размер пучка на экране станет W1: «Инфокоммуникационные технологии» Том 11, № 2, 2013 68 Петропавловский В.М., Топоркова Л.В. щ = λ(ζ + х) πνη (17) Подстановка (12), (15) и (16) в (14) дает ;Ç=f ; (18) Рг(Р) = Р-ü - exp w\P) = 1 + {P/Cj Построенный по формуле (18) график зависимости мощности излучения на выходе из стабилизатора от мощности излучения лазера приведен на рис. 2. Рис. 2. Зависимость мощности излучения на выходе из стабилизатора от мощности излучения лазера Как видно из рис. 2, функция РДР) имеет максимум вблизи Р = 14 мВт. При изменении Р на величину ±1,5 мВт выходная мощность Р1 меняется не более чем на 0,17 мВт, что означает повыше ние стабильности в 17 раз. Увеличение стабильности мощности входного излучения приведет к еще большим значениям коэффициента стабилизации. Таким образом, мы показали, что на основе эффекта тепловой линзы может быть создано устройство, стабилизирующее мощность излучения непрерывного лазера. К недостаткам метода можно отнести уменьшение мощности (в рассмотренном примере почти в 1,5 раза) и невысокое быстродействие (~1 мс). Увеличение размера диафрагмы а приведет к увеличению выходной мощности, но и к снижению коэффициента стабилизации установки ΔΡ/ΔΡ1. При необходимости повышения коэффициента стабилизации нужно уменьшить диаметр диафрагмы. К преимуществам метода относится простота конструкции стабилизатора и отсутствие внешних управляющих устройств. Подбор рабочей точки для требуемой величины выходной мощности и (или) коэффициента стабилизации осуществляется путем изменения диаметра диафрагмы, поглощения жидкости, длины кюветы или размера перетяжки.