Metrología Óptica de Frecuencias PDF
Metrología Óptica de Frecuencias PDF
Metrología Óptica de Frecuencias PDF
TESIS DOCTORAL
por
Juan Galindo Santos
Licenciado en Físicas
Directores:
Pedro Corredera Guillén
Aitor Villafranca Velasco
Madrid, 2017
Para Alicia, Claudia
Y el pequeño Rafael
Agradecimientos
En el año 2008 recalé en el antiguo Departamento de metrología del Instituto de
Física Aplicada, como técnico de laboratorio de calibración de Fibras ópticas, bajo la di-
rección del Dr. Pedro Corredera. Con él comencé a aprender los fundamentos y entresi-
jos de una tecnología llena de retos excitantes que abordar como las fibras ópticas. Con
el tiempo me propuso realizar la tesis doctoral compatibilizándolo con mi trabajo como
técnico de laboratorio, siendo mis directores de tesis, él mismo y el Dr. Aitor V. Velas-
co. Y después de este periplo, me encuentro escribiendo finalmente los agradecimientos.
A ellos les quiero mostrar mi primer agradecimiento.
Al Dr. Pedro Corredera por todas las lecciones aprendidas sobre las fibras ópti-
cas, la metrología y la vida. También le agradezco la confianza que ha depositado en mí
para llevar a cabo el trabajo de la presente tesis, así como la posibilidad que me ha brin-
dado de asistir a congresos tanto nacionales como internacionales.
Al Dr. Aitor. V. Velasco, por toda la ayuda que me ha prestado tanto en el labo-
ratorio como a la hora de la escritura, aplicando su buena pluma a la hora de la redac-
ción de la presente tesis. También quiero agradecerle todo el apoyo moral que me ha
prestado durante este tiempo.
Deseo agradecer a Joaquin Campos y Juan Diego Ania, que han sido directores
del Instituto de Óptica durante el tiempo de gestación de la tesis, por los momentos
compartidos departiendo sobre lo divino y lo humano en las horas de café.
Resumen .......................................................................................................................... 1
Introducción y objetivos................................................................................................. 5
1. El reto de medir el tiempo.................................................................................. 5
2. Objetivos ............................................................................................................ 6
3. Organización del presente documento ............................................................... 7
Capítulo1: Peines de frecuencias ópticas .................................................................... 11
1. Introducción histórica ...................................................................................... 11
2. Modelo teórico del peine de frecuencias ......................................................... 14
2.1. Fase portadora-envolvente........................................................................ 15
2.2. Espectro de frecuencias de un tren de pulsos ........................................... 18
2.3. Estabilización de las componentes frecuencias ........................................ 21
3. Aplicaciones de los peines de frecuencia en espectroscopía ........................... 24
4. Aplicaciones de los peines de frecuencia en el dominio del tiempo................ 26
5. Aplicaciones de los peines de frecuencia en metrología de frecuencias ......... 27
5.1. Medida de frecuencias ópticas absolutas .................................................. 27
5.2. Relojes atómicos ópticos .......................................................................... 29
Bloque I:
Capítulo 2: Desarrollo de relojes ópticos basados en celdas de absorción .............. 33
1. Introducción ..................................................................................................... 33
2. Desarrollo y caracterización de un reloj atómico de Rubidio .......................... 34
2.1. Estructura atómica del Rubidio ................................................................ 34
2.2. Diseño del reloj atómico de Rubidio ........................................................ 36
2.3. Caracterización espectroscópica de la celda ............................................. 37
2.4. Estabilidad del reloj atómico desarrollado ............................................... 41
3. Desarrollo y caracterización de un reloj óptico de cianuro de hidrógeno........ 42
3.1. Referencias frecuenciales en ventanas de comunicaciones ópticas.......... 42
3.2. Calibración experimental de la celda de absorción .................................. 44
3.3. Resultados experimentales y cálculo de incertidumbres .......................... 44
3.4. Corrección por presión y temperatura ...................................................... 50
3.5. Montaje para la caracterización del reloj óptico....................................... 54
3.6. Caracterización de la señal del reloj óptico .............................................. 55
4. Conclusiones .................................................................................................... 57
Capítulo 3: Síntesis de referencias ópticas monocromáticas mediante filtrado
Brillouin de peines de frecuencia ................................................................................ 61
1. Introducción ..................................................................................................... 61
2. Principio de operación ..................................................................................... 62
3. Montaje experimental ...................................................................................... 64
3.1. Generación ................................................................................................ 64
3.2. Filtrado ..................................................................................................... 65
3.3. Detección .................................................................................................. 66
4. Resultados ........................................................................................................ 67
4.1. Ancho de línea y estabilidad..................................................................... 67
4.2. Sintonización y filtrado ............................................................................ 69
5. Conclusiones .................................................................................................... 73
Bloque II:
Capítulo 4: Calibración de medidores de longitud de onda ..................................... 77
1. Introducción ..................................................................................................... 77
2. Funcionamiento del medidor de longitud de onda ........................................... 78
3. Técnica de calibración y fuentes de incertidumbre.......................................... 81
4. Montaje experimental ...................................................................................... 84
5. Resultados y análisis de incertidumbres .......................................................... 85
5.1. Banda C (1531 nm y 1541 nm) ................................................................ 85
5.1.1. Medidas preliminares ............................................................................ 85
5.1.2. Calibración de los láseres de referencia. ............................................... 86
5.1.3. Calibración de los medidores de longitud de onda ............................... 87
5.2. Calibración a 780 nm................................................................................ 89
5.3. Operación sostenida.................................................................................. 90
6. Conclusiones .................................................................................................... 92
Capítulo 5: Calibración de interrogadores de redes de Bragg ................................. 97
1. Introducción ..................................................................................................... 97
2. Principio de operación de los sensores basados en redes de Bragg en fibra.... 98
2.1. Introducción .............................................................................................. 98
2.2. Funcionamiento de las redes de Bragg como sensores............................. 99
2.3. Esquemas de detección de la longitud de onda de Bragg ....................... 102
3. Esquemas propuestos para la calibración de interrogadores.......................... 103
3.1. Método 1: Red de Bragg simulada ......................................................... 103
3.2. Método 2: Celdas de absorción .............................................................. 105
4. Demostración experimental ........................................................................... 106
4.1. Dispositivos bajo análisis ....................................................................... 106
4.2. Método 1: Red de Bragg simulada ......................................................... 106
4.3. Método 2: Celdas de absorción .............................................................. 107
5. Resultados y discusión ................................................................................... 108
5.1. Método 1: Red de Bragg simulada ......................................................... 108
5.2. Método 2: Celdas de absorción .............................................................. 110
6. Conclusiones .................................................................................................. 112
Conclusiones ................................................................................................................ 117
Anexo I: Caracterización de fuentes láser con un peine de frecuencias ................ 121
1. Introducción ................................................................................................... 121
2. Determinación de los signos de la frecuencia de offset y de batido .............. 122
3. Determinación del orden N del modo ............................................................ 125
Anexo II: El peine de frecuencias del Instituto de Óptica “Daza de Valdés” ....... 131
1. Introducción ................................................................................................... 131
2. El láser de femtosegundos ............................................................................. 132
3. Estabilización del peine de frecuencias ópticas ............................................. 136
3.1. Introducción ............................................................................................ 136
3.2. Estabilización de la frep del peine de frecuencias ópticas ....................... 136
3.2.1. El módulo de detección de fase (PHD110) ......................................... 137
3.2.2. El módulo de medida (MB110) .......................................................... 137
3.2.3. El módulo regulador (PID110) ........................................................... 138
3.2.4. El módulo de control del motor (MC110) .......................................... 139
3.2.5. El sintetizador de frecuencias (DDS120)............................................ 140
3.3. Estabilización de la f0 del peine de frecuencias ...................................... 141
3.3.1. Unidad electrónica de amplificación (OFD100) ................................. 143
3.3.2. El detector de fase digital (DXD200) ................................................. 144
3.3.3. La unidad de realimentación (PIC210) ............................................... 144
4. El contador de frecuencias (FXM50)............................................................. 145
5. Unidad de suministro de señales de radiofrecuencia (RFD10) ...................... 146
6. Las unidades de batido ................................................................................... 147
6.1. La unidad de detección de batido en aire BDU-FS ................................ 147
6.2. La unidad de detección de batido todo en fibra BDU-OF ...................... 149
6.3. La unidad electrónica de amplificación BFD100 ................................... 149
Referencias .................................................................................................................. 151
Acrónimos ................................................................................................................... 171
Publicaciones ............................................................................................................... 175
Resumen
La continua demanda de mejores prestaciones en las áreas de comunicaciones y
sensado ópticos ha llevado al desarrollo de equipos progresivamente más precisos y
complejos. Como consecuencia de estos avances, se hacen necesarias técnicas de cali-
bración de alta calidad, capaces de proporcionar precisión y trazabilidad a equipamiento
óptico cada vez más exigente. La presente tesis contribuye a cubrir estos requerimientos
crecientes proporcionando patrones de frecuencias ópticas de alta calidad, y demostran-
do su viabilidad al aplicar dichos patrones a la calibración de equipos comerciales para
comunicaciones y sensado.
1
Introducción y objetivos
3
Introducción y objetivos
Introducción y objetivos
1. El reto de medir el tiempo
5
Introducción y objetivos
go, el área científica más afectada por los OFCs ha sido sin duda la metrología de fre-
cuencias ópticas, convirtiéndose en una herramienta fundamental para la medida absolu-
ta de frecuencias [Margolis’03], la síntesis de frecuencias ópticas [Cundiff’01, Jost’02,
Schibli’05] o el desarrollo de relojes atómicos ópticos [Ye’01, Diddams’04, Stalna-
ker’07].
2. Objetivos
Dentro del marco descrito, podemos agrupar los objetivos de la tesis en dos blo-
ques principales:
6
Introducción y objetivos
El primer capítulo de esta tesis proporciona una visión general del desarrollo
histórico de los peines de frecuencias, así como sus aplicaciones más relevantes, con
particular énfasis en el área de la metrología de frecuencias.
7
Capítulo 1:
Peines de frecuencias ópticas
1. Introducción histórica
2. Modelo teórico del peine de frecuencias
2.1 Fase portadore-envolvente
2.2 Espectro de frecuencias de un tren de pulsos
2.3 Estabilización de las componentes frecuencias
3. Aplicaciones de los peines de frecuencia en espectroscopía
4. Aplicaciones de los peines de frecuencia en el dominio del tiempo
5. Aplicaciones de los peines de frecuencia en metrología de frecuencias
5.1 Medida de frecuencias ópticas absolutas
5.2 Relojes atómicos ópticos
9
Peines de frecuencias ópticas
11
Capítulo 1
Los OFCs solucionaron todas estas limitaciones, conectando el rango de las mi-
croondas y el rango de frecuencias ópticas en un sólo paso [Diddams’00]. Esto se con-
siguió mediante conjuntos equi-espaciados de líneas espectrales, generados por láseres
de modo anclado. Sin embargo, para lograrlo tuvieron que superarse numerosas barreras
tecnológicas. El desarrollo de los láseres pulsados empieza con la aparición del primer
láser de estado sólido de rubí realizado por Ted Mainman [Mainman’60]. A partir de es-
te punto surgieron dos corrientes de investigación con objetivos diferentes. Por una par-
te, diversos grupos de investigación centraron su interés en mejorar la estabilidad de la
frecuencia de emisión de láseres de baja potencia (por debajo del milivatio), desarro-
llando técnicas de estabilización que han permitido alcanzar estabilidades por debajo del
hertzio [Drever’83, Young’99, Träger’07]. Por otra parte, se empezaron a buscar méto-
dos para aumentar la potencia de emisión. Pero dado que la energía almacenada como
radiación en el medio láser está inevitablemente limitada, este aumento de potencia sólo
era posible reduciendo la duración de pulso.
El primer láser de modos anclados desarrollado fue un láser de He-Ne con an-
clado activo [Hargrove’64]. Este método de anclado se basa en la modulación de las
pérdidas dentro de la cavidad láser mediante un modulador (por ejemplo electro-óptico)
en sincronía con el tiempo de vuelo dentro de la cavidad [Paschotta’08]. La máxima po-
tencia del pulso atraviesa el modulador en el momento en el que las pérdidas son míni-
mas, aumentando las pérdidas en el momento que pasan las colas del pulso. La ventaja
de este método es que la tasa de repetición de pulsos puede ser elegida con cierta liber-
tad dentro de los límites impuestos por la geometría particular de diseño, pero presenta
un límite en la duración de los pulsos emitidos de decenas de picosegundos.
A principios de los años 70, Ippen desarrolló el primer láser de modos anclados
pasivo mediante la introducción en la cavidad láser de un absorbente saturable [Ip-
pen’72a], proporcionando pulsos de 1,5 ps de duración. En el método de anclado de
modos pasivo, la modulación de las pérdidas dentro de la cavidad es realizada por un
absorbente saturable (real o efectivo) [Sutherland’03]. Este mecanismo permite la gene-
ración de pulsos más cortos que en el caso del método de anclado activo, debido a que
cuanto menor es la duración de los pulsos, más rápida es la modulación de las pérdidas.
12
Peines de frecuencias ópticas
Sin embargo, esta drástica reducción en la duración de los pulsos emitidos por
los láseres pulsados seguía sin proporcionar la anchura espectral de una octava requeri-
da por la estabilización propuesta por Telle [Telle’99], fundamental para permitir el uso
metrológico de los peines (este requisito se explica en detalle en el epígrafe 2.3). El
fenómeno de autoensanchamiento espectral (también conocido como generación de su-
percontinuo) era un efecto bien conocido, producido fundamentalmente por un efecto no
lineal denominado automodulación de fase (SPM, del inglés Self-Phase Modulation)
[Boyd’08]. Sin embargo, no existía ninguna evidencia experimental de que durante la
generación del supercontinuo, todas las frecuencias que constituían el peine mantuvie-
ran la misma relación de frecuencias que en el pulso original. Esta evidencia experimen-
tal fue obtenida por Bellini y Hänsh mediante división de trenes de pulsos generados
por un láser Ti:zafiro en dos haces, posteriormente enfocados en un cristal no lineal de
CaF2 [Bellini’00].
13
Capítulo 1
Birks demostró que el supercontinuo podía ser obtenido utilizando fibras estranguladas
[Birks’00]. Partiendo de este espectro ensanchado, Jones implementó el método de au-
toreferencia discutido por Telle, permitiendo referenciar las componentes frecuenciales
del peine a un reloj atómico [Jones’00]. La estabilización de la fase portadora-
envolvente, permitió finalmente ser usados los OFCs como una “regla” para la medida
absoluta de frecuencias ópticas. Este desarrollo fue inmediatemente verificado con la
medida de la frecuencia de la emisión a 282 THz de un láser de Nd:YAG estabilizado
en Yodo [Diddams’00]. Desde entonces, las distintas variantes de los OFCs, incluyendo
aquellos basados en láseres de Ti:Zafiro [Diddams’01] y en resonadores de fibra dopada
con erbio (Er) e iterbio (Yb) [Rauschenberger’02, Adler’04] se han convertido en una
herramienta imprescindible en numerosas aplicaciones de metrología y espectroscopía.
Los OFCs consisten en trenes de pulsos ópticos ultracortos (en el rango de pico-
segundos a femtosegundos) emitidos por un láser de modos anclados, con una relación
de fase fija entre todos sus modos longitudinales [Ye’05]. Los parámetros que caracteri-
zan el tren de pulsos en el dominio temporal son la potencia de pico (Ppico); la duración
del pulso (P); la forma, determinada por la dispersión del medio láser o posibles efectos
no lineales; y el periodo de emisión (T) relacionado con la longitud de la cavidad del
láser. Si analizamos esta misma estructura en el dominio frecuencial, se obtiene un con-
junto discreto de líneas espaciadas regularmente a lo largo del espectro de emisión
(Figura 1), también denomidas “dientes”.
La frecuencia absoluta de cada uno de los dientes está definida por el valor de
dos frecuencias en el rango de las microondas: la frecuencia de repetición (frep) y la fre-
cuencia de offset (f0). La frecuencia de repetición define la separación entre dos modos
adyacentes del peine, y es inversamente proporcional al periodo de emisión. La frecuen-
cia de offset define la diferencia entre el primer modo del peine y el origen de frecuen-
cias, y está relacionada con los efectos dispersivos del material que produce la radiación
láser. Por lo tanto, la frecuencia de cada uno de los modos vendrá descrita por la ecua-
ción siguiente:
f N N f rep f 0 (1.1)
donde N es el índice de modo, que indica su posición dentro de la estructura del peine y
fN representa la frecuencia del modo N. La anchura del espectro en la que se encuentran
estos dientes se relaciona de forma inversa con la anchura de pulso P, de forma que a
14
Peines de frecuencias ópticas
Figura 1: (a) Sucesión de pulsos ópticos emitidos por un láser de modos an-
clados. (b) Representación en el espacio de la frecuencia del mismo tren de
pulsos.
15
Capítulo 1
Con esta descomposición del pulso, se puede definir para cada una de esas dos ondas
una velocidad de propagación en un medio: la velocidad de fase vp, asociada a la onda
portadora, y la velocidad de grupo vg, asociada a la envolvente. En medios no dispersi-
vos (medios en los cuales el índice de refracción es independiente de la longitud de on-
da) se cumple que ambas velocidades son iguales. Por contra, en los medios dispersivos,
como es el caso del material láser, estas dos velocidades son diferentes. Esto provoca
que la fase portadora-envolvente varíe de un pulso al siguiente.
CE
Las frecuencias permitidas en una cavidad láser vienen dados por la fórmula de
Fabry-Perot:
2 c
M M (1.3)
nL
k 2k
c 2 c
2
k kc (1.4)
c c
16
Peines de frecuencias ópticas
decir, que las derivadas de órdenes segundo y superiores en el desarrollo (1.4) son nu-
las, y que la derivada primera del vector de onda respecto a la frecuencia angular es
constante. La inversa de esta primera derivada define la velocidad de grupo del pulso:
k 1
cte
vg
(1.5)
c
k
i
i 0 i 1
c
1
k M kc M c k0 1 M (1.6)
vg vg
1
k0 kc c (1.7)
vg
ck 1 k
n c 0 (1.8)
M
vg M
2 vg
M M k0 vg (1.9)
L
2
M M k0vg Mrep k0vg (1.10)
T
O reescrito en frecuencia:
17
Capítulo 1
k0vg
f M M f rep (1.11)
2
k0 vg
f0 (1.12)
2
1 1
CE Lc (1.13)
v v p
g
Por otra parte, dado que la velocidad de fase de la onda portadora es igual a c/kc, susti-
tuyendo en (1.7), se puede escribir k0 de la siguiente forma:
1 1
k0 c CE (1.14)
v v L
p g
18
Peines de frecuencias ópticas
chos pulsos idénticos entre sí. Este tren de pulsos E(t) se puede representar matemáti-
camente como:
E t E t jT
j
p (1.16)
donde T es el tiempo entre pulsos y j es un número entero que índica el número del pul-
so. Dado la periodicidad de la ecuación (1.16) podemos representarla en forma de serie
de Fourier, utilizando un sumatorio de Poisson [Kammler’08].
1 ~
E t E t jT T E N T exp2iNt T
j
p
N
p
1 ~ (1.17)
E p N exp2iN t
N T
donde 0 es la fase portadora-envolvente inicial. Así pues, para el pulso j-ésimo, la ex-
presión del campo eléctrico será:
19
Capítulo 1
y por tanto, el campo eléctrico para el tren de pulsos se obtendrá sustituyendo (1.19) en
(1.16), resultando:
E t E t jT expi t j
j
e c CE cT 0 (1.20)
Una vez que tenemos la función que define nuestro tren de pulsos, podemos cal-
cular la transformada de Fourier para obtener el espectro correspondiente. De esta for-
ma, E() vendrá dado por la siguiente integral:
E Ee t jT expict j CE cT 0 exp i t dt
j
expi j CE cT 0 Ee t jT exp i c t dt
j
(1.21)
~
E t jT exp i t dt exp i jT E
e c c e c
(1.22)
Sustituyendo (1.22) en (1.21), la transformada de Fourier del campo eléctrico del tren de
pulsos es:
~
E expi j
j
CE cT 0 exp i j c T Ee c
(1.23)
~
expi0 Ee c expi j ce j T
j
~
E expi0 Ee c CE cT 2M (1.24)
M
20
Peines de frecuencias ópticas
2 CE
M M (1.25)
T T
CE f rep
f M M f rep M f rep f 0 (1.26)
2
ce f rep
f0 (1.27)
2
21
Capítulo 1
Consideremos dos modos del peine, N y M con M > N. Las frecuencias de estos
dos modos estarán definidas según las ecuaciones:
f N N f rep f 0
(1.28)
f M M f rep f 0
M M
f N f M 1 f 0 (1.29)
N N
22
Peines de frecuencias ópticas
23
Capítulo 1
Aplicando un enfoque alternativo, la salida del OFC puede utilizarse para ilumi-
nar directamente la muestra que se desea estudiar. Los trabajos de Baklanov [Bakla-
nov’77] y Eckstein [Eckstein’78] son los precursores de este tipo de espectrocopía. Una
ventaja interesante de la espectroscopía con OFC sobre la técnica de espectroscopía
láser es que las altas potencias de pico asociadas a los láseres de femtosegundos permi-
ten la generación de luz a longitudes de onda difícilmente obtenibles con láseres conti-
nuos. Este es el caso de las regiones ultravioleta en el rango de 40 nm a 200 nm [Jo-
24
Peines de frecuencias ópticas
25
Capítulo 1
26
Peines de frecuencias ópticas
tro métodos para llevar a cabo esta sincronicación: estabilización activa mediante im-
plementación electrónica [Ma’01, Shelton’2a]; modulación de fase cruzada para la sin-
cronización pasiva de dos láseres que comparten medio de ganancia [Leitenstorfer’95,
Wei’01]; sincronización de la repetición de ambos láseres al mismo oscilador óptico de
referencia [Bartels’03]; y correlación óptica cruzada entre los trenes de pulsos [Shi-
bli’03]. El nivel de sincronización obtenido permite aplicarse a la generación de fre-
cuencias suma o diferencia con alta potencia en cristales no lineales [Kaindl’00], la ge-
neración y modelado de pulsos [Shelton’01, Foreman’03] o la microscopía CARS
[Potma’02]. Asimismo, la capacidad de anclar en fase las ondas portadoras de dos láse-
res es aplicable a la generación de fuentes de femtosegundos sintonizables en regiones
espectrales que de otra forma no podrían ser alcanzadas [Rauschenberger’02, Joh-
son’10] o la distribucion de frecuencias ópticas de referencia a través de las redes de fi-
bras ópticas [Jones’03, Holman’03, Holman’04].
Finalmente, el desarrollo del láser alentó la idea de usar la luz para controlar la
evolución de las reacciones químicas mediante interferencia cuántica [Brumer’86, Ra-
bitz’00]. Los primeros estudios han demostrado el fenómeno de interferencia cuántica
usando dos pulsos controlados en su fase con los cuales se ionizada Rb [Yin’92] o se
controlaban las corrientes eléctricas en semiconductores GaAs en volumen [Haché’97],
observándose en ambos casos una conexión entre la dirección espacial y la fase relativa.
Recientemente, se ha observado que el GaAs es sensible a la fase portadora-envolvente
[Fortier’04], proporcionando un medio interesante para conseguir un detector de estado
sólido simple para la medida de la fase portadora-envolvente de forma directa sin la uti-
lización de técnicas interferométricas.
27
Capítulo 1
A partir del 2000 se produjo una explosión en el uso de OFCs basados en femto-
segundos para la medida de frecuencias absolutas. Las medidas más precisas que se rea-
lizaron tenían el objetivo de medir los patrones ópticos basados en transiciones con muy
altos factores de calidad. Una de las medidas más precisas fue la determinación de la
frecuencia de la transición 1S-2S del hidrógeno, obteniendo un nivel de precisión de 1,8
10-14 [Niering’00]. Investigadores de algunos laboratorios nacionales pusieron su obje-
tivo en medir la frecuencia de determinadas transiciones de trampas de átomos. Por
ejemplo, la transición cuadrupolar eléctrica del ión Hg+ fue determinada con un nivel de
precisión por debajo de 1 parte en 1014 [Udem’01]. Una precisión similar se obtuvo en
la medida de la transición cuadrupolar del ión Yb+ [Stenger’01]. Este incremento en la
precisión de la medida permitió corroborar postulados físicos fundamentales y determi-
nar constantes fundamentales hasta el siguiente decimal.
f X N f rep f 0 fbeat
(1.35)
Este cálculo requiere no obstante la determinación del número del modo más
próximo N con el que bate el láser desconocido, y de los signos de f0 y fbeat, cuya meto-
dología se presenta en detalle en el anexo I.
28
Peines de frecuencias ópticas
CE CE
TRANSFORMADA
DE FOURIER
f
f0
frep f
fbeat
29
Capítulo 1
I2, obteniendo una estabilidad mejor de 4,6·10–13 en un período de un año. Dado que la
estabilidad del sistema del OFC es mayor que la estabilidad de los patrones de frecuen-
cias ópticas actuales, es de esperar que la estabilidad de la señal de radiofrecuencia ob-
tenida sea básicamente la del patrón óptico.
30
Capítulo 2:
Desarrollo de relojes ópticos basados en cel-
das de absorción.
1. Introducción
2. Desarrollo y caracterización de un reloj atómico de rubidio
2.1 Estructura atómica del rubidio
2.2 Diseño del reloj atómico de rubidio
2.3 Caracterización espectroscópica de la celda
2.4 Estabilidad del reloj atómico desarrollado
3. Desarrollo y caracterización de un reloj óptico de cianuro de hidrógeno
3.1 Referencias frecuenciales en ventanas de comunicaciones ópticas
3.2 Calibración experimental de la celda de absorción
3.3 Resultados experimentales y cálculo de incertidumbres
3.4 Corrección por presión y temperatura
3.5 Montaje para la caracterización del reloj óptico
3.6 Caracterización de la señal del reloj óptico
4. Conclusiones
31
Desarrollo de relojes ópticos basados en celdas de absorción
33
Capítulo 2
cias de las transiciones en relojes ópticos, basta con estabilizar una fuente láser sintoni-
zable utilizando como referencia la absorción óptica producida a esa frecuencia [Bar-
ger’73, Cordiale’00, Madej’04, Träger’07]. Adicionalmente, es posible reducir el ancho
de línea de la referencia utilizando técnicas de saturación de la absorción [Dem-
trodër’15b].
34
Desarrollo de relojes ópticos basados en celdas de absorción
35
Capítulo 2
Las transiciones hiperfinas con origen en el mismo nivel fundamental dan como
resultado una única absorción debido al ensanchamiento Doppler, por lo que el espectro
del Rb estará compuesto por cuatro líneas de absorción [Steck’13, Steck’15, Dem-
trodër’10]. La diferencia de energía de los estados fundamentales desdoblados es muy
grande (alrededor de 3 GHz en el caso del Rb85 y de casi 7 GHz en el caso del Rb87),
pudiendo ser fácilmente resueltas mediante espectroscopía láser [Demtröder’15a]. Para
resolver las transiciones hiperfinas es necesario realizar la espectroscopía sub-Doppler o
espectroscopía de absorción saturada [Demtröder’15b].
36
Desarrollo de relojes ópticos basados en celdas de absorción
senta una mezcla de los isótopos Rb85 y Rb87 en sus abundancias naturales, utilizándose
como referencia para la estabilización una de sus líneas de absorción saturada o espec-
troscopía Sub-Doppler [Demtröder’15b].
Figura 10: Esquema experimental del reloj de Rb con bombeo colineal con-
trapropagante.
37
Capítulo 2
0,05
0,00
Transmitancia (u.a.)
-0,05
c
a
-0,10
-0,15
b
-0,20
-5 -4 -3 -2 -1 0 1 2 3 4 5
Frecuencia (GHz)
Línea Transiciones
a Rb87 (52S252P1,2,3)
b Rb85 (52S352P2,3,4)
c Rb85 (52S252P1,2,3)
d Rb87 (52S152P0,1,2)
0,1
0,0
-0,1
Potencia (u.a.)
-0,2
-0,3
-0,4
-0,5
-5 -4 -3 -2 -1 0 1 2 3 4 5
Desintonia en Frecuencia (GHz)
38
Desarrollo de relojes ópticos basados en celdas de absorción
transiciones, indicándose las frecuencias de cada uno de los picos en las tablas desde la
Tabla 4 a la Tabla 7.
1,2
4
1 3
Potencia (u.a.)
1,1 2
1,0
-3,2 -3,1 -3,0 -2,9 -2,8 -2,7 -2,6 -2,5 -2,4 -2,3 -2,2
Desintonizacion en Frecuencia (GHz)
39
Capítulo 2
1,8
1,7
3 4
Potencia (u.a.) 1,6
1,5
1,4 2
1 5
1,3
1,2
-2,2 -2,1 -2,0 -1,9 -1,8 -1,7 -1,6 -1,5 -1,4 -1,3
Desintonia en Frecuencia (GHz)
1,8
Potencia (u.a.)
1,7 2
1
1,6
1,5
-2,4 -2,3 -2,2 -2,1 -2,0 -1,9 -1,8 -1,7 -1,6 -1,5 -1,4 -1,3 -1,2 -1,1
Desintonía en frecuencia (GHz)
40
Desarrollo de relojes ópticos basados en celdas de absorción
1,76
1,74
Potencia (u.a.)
1,72 4
3
12
1,70 5
1,68
1,66
-2,4 -2,3 -2,2 -2,1 -2,0 -1,9 -1,8 -1,7 -1,6 -1,5 -1,4 -1,3 -1,2 -1,1
Desintonía en frecuencia (GHz)
En las líneas representadas en la Figura 13, Figura 14 y Figura 16, se pueden ob-
servar la existencia de cinco absorciones saturadas. De estas cinco absorciones satura-
das, únicamente tres corresponden con las transiciones hiperfinas del Rb. Las otras dos
corresponden con lo que se conoce como Crossover [Demtröder,15b], ubicadas exacta-
mente en la frecuencia promedio de las frecuencias de las transiciones hiperfinas adya-
centes. Para el diseño del reloj atómico óptico cualquiera de las cinco absorciones serían
adecuadas.
41
Capítulo 2
Ambas absorciones presentan una estabilidad similar, siendo válidas para ser
usadas como referencia frecuencial en la calibración de equipos ópticos. En particular,
presentan un valor inicial de 2,42·10–7 para t = 1 s, seguidos de una tendencia decrecien-
te con razón –1/2 para tiempos menores de 64 segundos. El valor mínimo de varianza,
3,15·10–8, se obtiene para 64 segundos. A partir de 64, la absorción 5 presenta una lige-
ra mejora en la estabilidad respecto a la absorción 4, siendo aquella por lo tanto preferi-
blemente seleccionada para labores de calibración de larga duración.
1E-6
()
1E-7
1E-8
1 10 100 1000
(s)
42
Desarrollo de relojes ópticos basados en celdas de absorción
1,1
Rama R Rama P
1,0
0,9
Transmitancia normalizada
0,8
0,7
0,6
0,5
0,4
1525 1530 1535 1540 1545 1550 1555 1560 1565
Longitud de onda (nm)
En este trabajo se ha adquirido una celda de absorción con una mezcla de tres
gases: cianuro de hidrógeno H12C14N, monóxido de carbono 12C16O y monóxido de car-
bono 13C16O a presiones de 5 Torr, 150 Torr y 150 Torr respectivamente (Wavelength
References, modelo Tri-H(80)-5/150/150-FCAPC). La celda tiene un tubo de vidrio de
longitud de 18,92 cm, protegido exteriormente por una carcasa de plástico con conecto-
res de fibra óptica SC/APC. Mediante cinco reflexiones internas, se consigue un camino
óptico total cercano a 80 cm. La combinación de gases cubre el rango espectral desde
los 1520 nm hasta los 1640 nm (ventanas C y L). El espectro de cada una de las espe-
cies químicas de la mezcla presenta dos ramas bien diferenciadas, denominadas rama R
y P. Este fenómeno se observa en la Figura 18 para el cianuro de hidrógeno (H12C14N),
43
Capítulo 2
44
Desarrollo de relojes ópticos basados en celdas de absorción
0,5
0,0
-0,5
Transmitancia (dB)
-1,0
-1,5
-2,0
-2,5
-3,0
1520 1530 1540 1550 1560 1570 1580 1590 1600 1610 1620 1630
Longitud de onda (nm)
0
Potencia (dB)
-1
-2
-3
1520 1530 1540 1550 1560
Longitud de onda (nm)
45
Capítulo 2
0,0
Potencia (dB)
-0,3
-0,6
0,0
Potencia (dB)
-0,3
-0,6
46
Desarrollo de relojes ópticos basados en celdas de absorción
H12C14N
Rama R Rama P
med. ref. Incer. med. ref. Incer.
/nm /nm /pm /pm /nm /nm /pm /pm
R23 1520,28112 1520,28055 0,57 0,59 P1 1534,52975 1534,52984 -0,09 0,59
R22 1520,72655 1520,72553 1,02 0,59 P2 1535,23775 1535,23629 1,46 0,59
R21 1521,1818 1521,18088 0,92 0,59 P3 1535,9536 1535,95333 0,27 0,59
R20 1521,6471 1521,64655 0,55 0,59 P4 1536,68225 1536,68095 1,30 0,59
R19 1522,12438 1522,12261 1,77 0,59 P5 1537,42025 1537,41922 1,03 0,59
R18 1522,61005 1522,60902 1,03 0,59 P6 1538,1688 1538,16813 0,67 0,59
R17 1523,1059 1523,10582 0,08 0,59 P7 1538,9279 1538,92767 0,23 0,59
R16 1523,61373 1523,61302 0,72 0,59 P8 1539,6979 1539,69784 0,06 0,59
R15 1524,13113 1524,13061 0,53 0,59 P9 1540,47845 1540,47871 -0,26 0,59
R14 1524,6591 1524,65859 0,51 0,59 P10 1541,2711 1541,27025 0,85 0,59
R13 1525,19785 1525,19702 0,83 0,59 P11 1542,0735 1542,07248 1,02 0,59
R12 1525,7467 1525,74586 0,84 0,59 P12 1542,8854 1542,88542 -0,02 0,59
R11 1526,30555 1526,30513 0,42 0,59 P13 1543,70997 1543,70908 0,89 0,59
R10 1526,8751 1526,87484 0,26 0,59 P14 1544,54325 1544,54348 -0,23 0,59
R9 1527,4553 1527,455 0,30 0,59 P15 1545,3892 1545,38859 0,61 0,59
R8 1528,0462 1528,04565 0,55 0,59 P16 1546,24565 1546,24448 1,17 0,59
R7 1528,6476 1528,64676 0,84 0,59 P17 1547,1127 1547,11115 1,55 0,59
R6 1529,2594 1529,25834 1,06 0,59 P18 1547,9879 1547,9886 -0,70 0,59
R5 1529,88073 1529,88044 0,29 0,59 P19 1548,8785 1548,87683 1,67 0,59
R4 1530,51343 1530,51304 0,39 0,59 P20 1549,77713 1549,7759 1,22 0,59
R3 1531,15707 1531,15613 0,93 0,59 P21 1550,68687 1550,68579 1,08 0,59
R2 1531,81027 1531,80978 0,48 0,59 P22 1551,6083 1551,60651 1,79 0,59
R1 1532,47377 1532,47395 -0,18 0,59 P23 1552,538 1552,53809 -0,09 0,59
R0 1533,14828 1533,14868 -0,40 0,59 P24 1553,4795 1553,48054 -1,04 0,59
P25 1554,43205 1554,43387 -1,82 0,59
P26 1555,3982 1555,3981 0,10 0,59
47
Capítulo 2
=2) de la calibración.
3 3
(a) (b)
2
2
1
Diferencia (pm)
Diferencia (pm)
1
0
-1
-1
-2
-2 -3
1520 1525 1530 1535 1530 1535 1540 1545 1550 1555 1560
Longitud de onda (nm) Longitud de onda (nm)
Figura 24: Diferencias entre las longitudes de onda medidas y las de refe-
rencia para las ramas (a) R y (b) P del H12C14N. En rojo se muestra el valor
promedio de las diferencias.
48
Desarrollo de relojes ópticos basados en celdas de absorción
0 3
-1 1
0
Diferencia (pm)
Diferencia (pm)
-2 -1
-2
-3 (a) -3
-4 (b)
-4 -5
1560 1565 1570 1575 1575 1580 1585 1590 1595
Longitud de onda (nm) Longitud de onda (nm)
Figura 25: Diferencias entre las longitudes de onda medidas y las de refe-
rencia para las ramas (a) R y (b) P del 12C16O. En rojo se muestra el valor
promedio de las diferencias.
Finalmente, la Tabla 11y Figura 23 muestran el caso del 13C16O, también usando
los valores de referencia del NIST [Gilbert’02]. Las diferencias encontradas son (-1.84
± 0.24) pm para la rama R y (-2.04 ± 0.27) pm para la rama P.
13
C16O
Rama R Rama P
med. ref. Incer. P med ref. Incer.
/nm /nm /pm /pm ./nm /nm /pm /pm
R19 1596,15653 1596,1595 -2,98 0,59 P1 1610,657 1610,6596 -2,60 0,59
R18 1596,5873 1596,5895 -2,20 0,59 P2 1611,63675 1611,6393 -2,55 0,59
R17 1597,0429 1597,0454 -2,50 0,59 P3 1612,64463 1612,6457 -1,07 0,59
R16 1597,5255 1597,5271 -1,60 0,59 P4 1613,6762 1613,6788 -2,60 0,59
R15 1598,03325 1598,0349 -1,65 0,59 P5 1614,737 1614,7388 -1,80 0,59
R14 1598,5671 1598,5686 -1,50 0,59 P6 1615,8236 1615,8255 -1,90 0,59
R13 1599,1269 1599,1284 -1,50 0,59 P7 1616,9378 1616,9392 -1,40 0,59
R12 1599,7123 1599,7141 -1,80 0,59 P8 1618,07745 1618,0797 -2,25 0,59
R11 1600,3241 1600,3258 -1,70 0,59 P9 1619,245 1619,2473 -2,30 0,59
R10 1600,9615 1600,9636 -2,10 0,59 P10 1620,4404 1620,442 -1,60 0,59
R9 1601,62585 1601,6274 -1,55 0,59 P11 1621,66165 1621,6641 -2,45 0,59
R8 1602,3161 1602,3174 -1,30 0,59 P12 1622,9105 1622,9132 -2,70 0,59
R7 1603,0323 1603,0334 -1,10 0,59 P13 1624,1886 1624,1898 -1,20 0,59
R6 1603,7738 1603,7756 -1,80 0,59 P14 1625,4916 1625,494 -2,40 0,59
R5 1604,5426 1604,5439 -1,30 0,59 P15 1626,8233 1626,8257 -2,40 0,59
R4 1605,3368 1605,3385 -1,70 0,59 P16 1628,1836 1628,1851 -1,50 0,59
R3 1606,15777 1606,1593 -1,53 0,59
R2 1607,0042 1607,0064 -2,20 0,59
R1 1607,877 1607,8799 -2,90 0,59
49
Capítulo 2
0 0
(a) (b)
-1 -1
Diferencia (pm)
Diferencia (pm)
-2 -2
-3 -3
-4 -4
1595 1600 1605 1610 1610 1615 1620 1625 1630
Longitud de onda (nm) Longitud de onda (nm)
Figura 26: Diferencias entre las longitudes de onda medidas y las de refe-
rencia para las ramas (a) R y (b) P del 13C16O. En rojo se muestra el valor
promedio de las diferencias.
H12C14N 12
C16O 13
C16O
Promedio Incert. Promedio Incert. Promedio Incert.
/pm /pm /pm /pm /pm /pm
Rama R 0,60 0,18 –1,60 0,25 –1,84 0,24
Rama P 0,49 0,35 –1,74 0,52 –2,04 0,27
50
Desarrollo de relojes ópticos basados en celdas de absorción
Tabla 13: Comparación entre las presiones parciales en las celdas de refe-
rencia usadas (incluida su referencia bibliográfica) y la célula bajo calibra-
ción.
12
C16O
Rama R Rama P
20 -0,67 1,38 1 -1,27 0
19 -3,15 -1,16 2 -1,5 -0,2
18 -1,27 0,63 3 -1,73 -0,3
17 -1,7 0,14 4 -2,4 -0,83
16 -1,6 0,17 5 -2,4 -0,75
15 -1,9 -0,14 6 -1,9 -0,14
14 -1,7 0,06 7 -1,6 0,2
13 -1,7 0,05 8 -1,7 0,13
12 -2,7 -0,96 9 -1,7 0,18
11 -1,85 -0,12 10 -1,85 0,06
10 -1,65 0,06 11 -1,3 0,65
9 -1,7 0,01 12 -2,35 -0,38
8 -1,4 0,27 13 -1,73 0,27
7 -1,7 -0,04 14 -2,77 -0,74
6 -1,9 -0,3 15 -3,25 -1,19
5 -1,3 0,24 16 1,6 3,68
4 -1,4 0,08
3 -1,07 0,33
2 -0,85 0,46
1 -1,4 -0,17
0 -0,95 0,15
51
Capítulo 2
13
C16O
Rama R Rama P
19 -2,98 -0,81 1 -2,6 -1,34
18 -2,2 -0,09 2 -2,55 -1,22
17 -2,5 -0,49 3 -1,07 0,37
16 -1,6 0,31 4 -2,6 -1,06
15 -1,65 0,16 5 -1,8 -0,15
14 -1,5 0,31 6 -1,9 -0,12
13 -1,5 0,31 7 -1,4 0,49
12 -1,8 0,01 8 -2,25 -0,34
11 -1,7 0,06 9 -2,3 -0,37
10 -2,1 -0,38 10 -1,6 0,41
9 -1,55 0,17 11 -2,45 -0,37
8 -1,3 0,42 12 -2,7 -0,66
7 -1,1 0,62 13 -1,2 0,78
6 -1,8 -0,15 14 -2,4 -0,36
5 -1,3 0,28 15 -2,4 -0,36
4 -1,7 -0,18 16 -1,5 0,66
3 -1,53 -0,12
2 -2,2 -0,9
1 -2,9 -1,7
2 5
(a) 4 (b)
1
3
Diferencia (pm)
Diferencia (pm)
2
0
1
0
-1
-1
-2 -2
1560 1565 1570 1575 1575 1580 1585 1590 1595
Longitud de onda (nm) Longitud de onda (nm)
52
Desarrollo de relojes ópticos basados en celdas de absorción
2 2
(b)
1 (a)
1
0
Diferencia (pm)
Diferencia (pm)
0
-1
-1
-2
-3 -2
1595 1600 1605 1610 1610 1615 1620 1625 1630
Longitud de onda (nm) Longitud de onda (nm)
corregida
Gas Rama. /pm /pm
Prom. Incer. Prom. Incer.
12 R –1,60 0,25 0,06 0,22
C16O
P –1,74 0,52 0,04 0,54
13 R –1,84 0,24 –0,11 0,26
C16O
P –2,04 0,25 –0,23 0,32
Por otra parte, las líneas de absorción pueden ser sensibles a cambios de tempe-
ratura y campos electromagnéticos. De forma general, esta sensibilidad es pequeña en
comparación con el cambio de presión. Debido al hecho de que el cambio de presión
depende de la tasa de colisión molecular, hay una pequeña dependencia del cambio de
la presión con la temperatura. Esta dependencia puede estudiarse con la fórmula [Gil-
bert’02]:
Δ Δ
53
Capítulo 2
54
Desarrollo de relojes ópticos basados en celdas de absorción
Lock 110 fabricado por Toptica, el cual genera una señal de control que es enviada a la
fuente de alimentación ILX a través de un BNC que permite la modulación de la co-
rriente de alimentación del diodo láser. El 5 % de la segunda rama en el segundo aco-
plador es introducido en la unidad de batido todo-en-fibra para poder medir la frecuen-
cia de emisión con el OFC (Figura 29).
0
Potencia (dBm)
-1
FWHM= 15,8 pm
-2
Figura 30: Línea P10 del H12C14N medida durante la calibración con el WM
(negro) y su ajuste gaussiano (rojo).
55
Capítulo 2
-30
-40
-50
Potencia (dBm)
-60
-70
-80
-90
-100
5 10 15 20 25 30 35 40 45 50 55
Frecuencia (MHz)
Figura 31: Batido filtrado del laser DFB con un diente del IO-FFC.
Tabla 17: Frecuencia del laser anclado en la línea P10 del H12C14N medida
con el OFC.
56
Desarrollo de relojes ópticos basados en celdas de absorción
Tabla 18: Frecuencia del laser anclado en la línea P10 del H12C14N medida
con el WM
4. Conclusiones
57
Capítulo 3:
Síntesis de referencias ópticas monocromáti-
cas mediante filtrado Brillouin de peines de
frecuencia.
1. Introducción
2. Principio de operación
3. Montaje experimental
3.1 Generación
3.2 Filtrado
3.3 Detección
4. Resultados
4.1 Ancho de línea y estabilidad
4.2 Sintonización y filtrado
5. Conclusiones
59
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
Sin embargo, cuando se desea extraer una componente frecuencial arbitraria del
61
Capítulo 3
OFC fuera de las referencias disponibles, se hace necesario el filtrado óptico del peine.
Para esta tarea se han propuesto AWGs (del inglés Array-Waveguide Gratings), redes
de Bragg en fibra o anclado mediante inyección óptica, [Cai'95, Fukushima'03]. No obs-
tante, todas las alternativas mencionadas presentan graves limitaciones en términos de
velocidad de respuesta, capacidad de sintonización y/o potencia de salida. En 2009,
Subías y su grupo demostraron la amplificación Brillouin selectiva de un OFC para
propósitos metrológicos [Subías'09]. Las líneas obtenidas en la banda C de comunica-
ciones ópticas mantenían una precisión en la frecuencia relativa en el orden del Hz, de-
mostrando que el filtrado no deterioraba la estabilidad intrínseca del modo filtrado. A
pesar de ello, el OFC de partida se basaba en modulación cruzada de fase (XPM)
[Agrawal’07], por lo que el sistema no podía proporcionar medidas frecuenciales abso-
lutas ni era utilizable para caracterización durante periodos largos de tiempo.
2. Principio de operación
f N N · f rep f 0 (3.1)
Cada uno de los dientes que componen el OFC puede considerarse un láser con-
tinuo con una frecuencia absoluta conocida con muy baja incertidumbre y por tanto
adecuado para ser utilizado como referencia óptica. Para extraer una única de estas refe-
rencias, proponemos un filtrado basado en el efecto de dispersión Brillouin estimulada
62
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
2nva
B (3.2)
p
1
B (3.3)
2TB
estando este valor en el rango de los MHz y siendo inversamente proporcional al cua-
drado de p [Nikles’97].
En nuestro caso, puesto que el ancho de banda del SBS está en el orden de la de-
cena de MHz, es posible amplificar un único diente del OFC, generando así referencias
monocromáticas. Esta amplificación da como resultado una señal óptica que hereda las
características de anchura de línea y estabilidad en frecuencia del peine original, pero
cuya longitud de onda puede seleccionarse de manera gruesa modificando las condicio-
63
Capítulo 3
3. Montaje experimental
Para demostrar la viabilidad del sistema propuesto, se ha filtrado una única línea
espectral del OFC auto-referenciado disponible en el laboratorio del IO-CSIC. Tal como
se muestra en la Figura 32, podemos dividir el montaje experimental en tres secciones:
la generación del OFC, el filtrado SBS y caracterización de la señal resultante.
WM
FILTRADO
Bombeo
Brillouin Aislador PC1
90/10 VOA
EDFA
1
Aislador
2
FOTF 3
50/50
Contador
FOP
90/10 PC3 MZM PC2
D
OSA o WM
ESA DETECCIÓN
3.1. Generación
La sección de generación comprende el OFC auto-referenciado y un filtro óptico
sintonizable que nos permite seleccionar el rango de frecuencias del que posteriormente
extraeremos el diente particular. El filtro óptico utilizado es un filtro en fibra modelo
FOTF-025121333 de Agiltron, con rango de sintonización entre 1510 nm y 1580 nm, y
un ancho de banda de 1 nm. El peine filtrado es dividido en un acoplador 50/50 cuyas
salidas se dirigen a las secciones de filtrado y detección.
64
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
ra 33, y está descrito en detalle en el Anexo II. Consiste en un láser de modos anclados
basado en un oscilador de anillo en fibra óptica dopada con Er [Adler’04]. El oscilador
genera pulsos de anchura de femtosegundos centrados en 1560 nm, con una tasa de re-
petición sintonizable entre 98 MHz y 102 MHz. La frecuencia de repetición se fija con
un sintetizador de frecuencias referenciado a un reloj de Rb, sincronizado a su vez a la
referencia internacional de tiempo mediante GPS (del inglés Global Positioning Sys-
tem). La estabilidad obtenida es del orden de 10–13 s a tiempos cortos. Para este experi-
mento, se utilizó la salida auxiliar en fibra, que extrae el 6% de la potencia del oscila-
dor.
Figura 33: Esquema del OFC. WDM (Multiplexor por división de onda), PC
(Control de polarización), M (Espejo), HWP (Lámina media onda), PBS
(Cubo divisor de haz por polarización), G (Red de difracción), D (Detector),
A (Apertura), L (Lente), PPLN (Cristal no lineal de Niobatio de Litio)
3.2. Filtrado
El filtrado del diente del peine se lleva a cabo en 1 km de fibra monomodo (Pire-
lli SMF28). Como láseres de bombeo Brillouin contrapropagante, se utilizaron secuen-
cialmente dos láseres DFB conectorizados con fibra óptica, comunes en el área de co-
municaciones ópticas, de longitudes de onda nominales 1531,90 nm (Láser 1) y 1542,14
nm (Láser 2). Los láseres se estabilizaron en las pendientes de las absorciones P11
(1531,5879 nm) y P25 (1540,82744 nm) de una celda de acetileno (12C2H2), y se ampli-
ficaron con un EDFA (EDFA-TV-24-FC/APC-B-11 de Accelink) con +23 dBm de po-
tencia óptica de salida. Las condiciones de bombeo se optimizaron mediante un atenua-
dor óptico variable en fibra (VOA) y un controlador de polarización (PC1).
65
Capítulo 3
3.3. Detección
El diente filtrado, extraído a través del tercer puerto del circulador, se divide
mediante un acoplador 50/50. Un 50% de la señal se monitoriza en un analizador de es-
pectros ópticos (OSA), con el fin de asegurar la superposición entre la banda de ganan-
cia Brillouin y el OFC. El 50% de la señal restante es utilizada para caracterizar la esta-
bilidad y la anchura de la señal obtenida del filtrado, mediante batido con el propio OFC
auto-referenciado en una unidad de batido todo en fibra. Puesto que el diente filtrado ha
sido obtenido del mismo OFC auto-referenciado con el que se pretende medir, es nece-
sario inducir un desplazamiento en frecuencia conocido para que el batido resultante no
sea nulo. Esto se consigue mediante un modulador óptico Mach-Zehnder de banda late-
ral única (MZM) operando en condiciones de polarización V para la supresión de la
portadora. Antes del MZM, se incluye un segundo control de polarización (PC2) para
optimizar la modulación. La frecuencia de modulación elegida fue de 430 MHz, de tal
forma que la frecuencia de batido se encontrara dentro de la banda de la electrónica de
detección. A la salida del MZM, el tercer control de polarización (PC3) permite optimi-
zar la potencia de la señal de batido. La frecuencia de batido es medida con un fotode-
tector PIN de alta sensibilidad modelo FPD510 (Menlo Systems). El valor instantáneo
de la frecuencia se midió con un contador de frecuencias de cuatro canales con 1 mHz
de resolución modelo FXM50 (Menlo Systems), mientras que el ancho de línea de la
señal filtrada se midió con un Analizador de espectros eléctricos (ESA) modelo
HM5014-2 (Hameg) con 9 kHz de resolución.
66
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
4. Resultados
0
Potencia normalizada (dB)
-10
-20
-30
-6 -4 -2 0 2 4 6
Frecuencia (MHz)
Para cada uno de los láseres seleccionados se midió la frecuencia absoluta del
bombeo Brillouin, el diente del OFC original y del diente una vez filtrado, durante un
período continuo de 8 horas. Existen diversos estudios estadísticos posibles para anali-
zar la estabilidad de las señales y sus fuentes de ruido, incluyendo la varianza de Allan
[Allan’66] y sus modificaciones posteriores [Riley’08]. En particular, nuestro estudio se
basa en la Desviación de Allan modificada (MVAR), regida por la siguiente ecuación:
2
yk m yk
M 3 m 2 j m 1 i m 1
1
MVAR ( ) 4
2 m M 3m 2 j 1 I j k i
(3.4)
67
Capítulo 3
f t f nom
y t (3.5)
f nom
1E-9
Desviación de Allan modificada
1E-10
1E-11
1E-12
1E-13
1 10 100 1000 10000
(s)
68
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
1E-8
1E-9
1E-11
1E-12
1E-13
1 10 100 1000 10000
(s)
69
Capítulo 3
-10
Potencia normalizada (dB)
-20
-30
-40
1470 1500 1530 1560 1590 1620 1650
Longitud de onda (nm)
11,05
11,00
Desplazamiento Brillouin (GHz)
10,95
10,90
10,85
10,80
10,75
10,70
1525 1530 1535 1540 1545 1550 1555 1560 1565 1570 1575
Longitud de onda (nm)
70
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
0
6 dB
-10
-20
46 dB
-30
Potencia (dBm)
-40
-50
-60
-70
-80
-90
-100
1540,0 1540,2 1540,4 1540,6 1540,8 1541,0 1541,2 1541,4 1541,6 1541,8
Longitud de onda (nm)
71
Capítulo 3
-5
Potencia normalizada (dB)
-10
-15
-20
-25
-30
-35
170 180 190 200 210 220 230
Frecuencia (MHz)
-1
Potencia normalizada (dB)
-2
-3
-4
-5
-6
-7
-20 -15 -10 -5 0 5 10 15 20
Frecuencia (MHz)
72
Síntesis de referencias ópticas monocromáticas mediante filtrado Brillouin de pei-
nes de frecuencias
5. Conclusiones
73
Capítulo 4:
Calibración de medidores de longitud de on-
da.
1. Introducción
2. Funcionamiento del medidor de longitud de onda
3. Técnica de calibración y fuentes de incertidumbre
4. Montaje experimental
5. Resultados y análisis de incertidumbres
5.1 Banda C (1531 nm y 1541 nm)
5.1.1 Medidas preliminares
5.1.2 Calibración de los láseres de referencia
5.1.3 Calibración de los medidores de longitud de onda
5.2 Calibración a 780 nm
5.3 Operación sostenida
6. Conclusiones
75
Calibración de medidores de longitud de onda
77
Capítulo 4
Los WM calibrados en este trabajo, al igual que la mayoría de los equipos co-
merciales, están basados en el interferómetro de Michelson [Derikson’98]. La Figura 42
muestra esquemáticamente el funcionamiento de este tipo de interferómetros. La radia-
ción coherente procedente de la fuente desconocida es colimada y dividida en dos haces.
Uno de los haces se refleja en un espejo fijo, mientras que el otro haz se refleja en un
espejo móvil. La interferencia resultante depende de la diferencia de camino óptico en-
tre los brazos, obteniéndose máximos de interferencia en:
m 2n d (4.1)
78
Calibración de medidores de longitud de onda
n L N (4.2)
n L N 0
0 0
(4.3)
N n
0
0 (4.4)
N n0
79
Capítulo 4
La determinación del número de ciclos que componen cada uno de los respecti-
vos interferogramas requiere de un dispositivo electrónico denominado contador de cru-
ces por cero (Zero-crossing counter), que produce un pulso eléctrico cada vez que la in-
tensidad del interferograma cambia de signo en el flanco de subida o de bajada. Deno-
minaremos C1 al contador asociado a la longitud de onda de referencia y C2 al contador
asociado a la longitud de onda desconocida será C2. Ambos contadores se inician si-
multáneamente, estando el instante de inicio de medida definido por el disparo del de-
tector asociado a la longitud de onda desconocida (D2). Cuando el contador C2 alcanza
un valor de cuentas predeterminado, N0, ambos contadores se detienen, obteniéndose en
el contador C1 un número entero de ciclos Nr más una fracción de ciclo, :
N N0
(4.5)
N Nr
0
N r n
0
0 (4.6)
N0 n
MNr n N M n
0
0 r 0
0
(4.7)
MN0 n N0 n
80
Calibración de medidores de longitud de onda
[Owens’67].
2406030 15997
n0 1 108 8342,13 (4.8)
130 1 38,9 1 2
2
donde el superíndice cero indica que el índice de refracción está calculado en las condi-
ciones estándar y está expresado en micras. Esta ecuación fue corregida posteriormen-
te por el propio Edlén [Edlén’66] para el caso más genérico de aire seco a temperatura y
presión no estándar, resultando la siguiente expresión:
nT , P 1
n
0
1 0,00138823 P
(4.9)
1 0,003671T
0,0457 8
nT , P, h nT , P h 5,722 10 (4.10)
2
81
Capítulo 4
donde fOFC es la frecuencia del láser de referencia, N es el índice del diente con el que
bate el láser, frep es la frecuencia de repetición del peine, f0 es la frecuencia de cero y fbeat
es la frecuencia de batido entre el láser y el diente N. Por lo tanto, la incertidumbre en la
medida de la frecuencia ( u f en k = 1) depende de la incertidumbre en la medida de las
OFC
frecuencias frep, f0 y fbeat. Las incertidumbres asociadas a la frep y f0 están dadas por la ca-
lidad del OFC (estabilidad del reloj atómico usado para sincronizar, los bucles de reali-
mentación, etc...). En cambio, la incertidumbre de la fbeat depende de la estabilidad en
frecuencia y de la anchura espectral del láser.
82
Calibración de medidores de longitud de onda
83
Capítulo 4
uK WM
u 2f u 2f
WM OFC
(4.13)
4. Montaje experimental
84
Calibración de medidores de longitud de onda
las absorciones P11 (1531,5879 nm) y P25 (1540,82744 nm) de una celda de acetileno
12
C2H2, trazable al NIST [Gilbert’01a].
1
Sbatido 2 (4.14)
batido
1
laser 2
donde Sbatido es la señal del espectro de batido adquirida por el ESA, y laser es la an-
chura de línea espectral del láser. En la Tabla 21 se presentan las anchuras de línea pro-
porcionadas por el fabricante y las anchuras de línea de estos mismos láseres medidas
en el laboratorio con el método descrito. En la Figura 46 se muestran los datos obteni-
dos de la señal del batido para ambos láseres, junto sus respectivos ajustes.
85
Capítulo 4
Longitud de onda
Anchura de
nominal Anchura de línea
Modelo línea (fabricante)
(nm)/Frecuencia ITU (medida) (MHz)
(MHz)
(THz)
EP1550-NLW-BBI-001 1542.14/194.40 0.1 0.66 ±0.01
EP1550-DM-VAD-001 1531.90/195.80 <1 2.96 ±0.01
Tabla 21: Anchura de líneas nominales y medidas de los láseres usados co-
mo referencias.
-35 -40
-40 a
-45
b
-45
Potencia (dBm)
Potencia (dBm)
-50
-50
-55
-55
-60
-65
-60
-70
-75 -65
5 10 15 20 25 30 35 40 45 50 55 5 10 15 20 25 30 35 40 45 50 55
Figura 46: Batidos de los láseres con el OFC obtenidas del ESA. (a) Láser
centrado en 1542 nm, (b) láser centrado en 1531 nm.
86
Calibración de medidores de longitud de onda
Tabla 23: Frecuencia medida con el OFC de los láseres anclados en las líne-
as P25 y P11 de la célula de 12C2H2.
Tabla 24: Cálculo de las incertidumbres de los láseres inducidos por la dis-
persión del índice de refracción del aire.
Cada uno de los láseres de referencia (Láser 1 y Láser 2) se han medido cinco
veces con cada uno de los WM. La calibración de ambos WM se realizó a temperatura
87
Capítulo 4
Tabla 25: Frecuencias medidas de los dos láseres de referencia con los dos
equipos a calibrar.
La corrección de los WM está dada por la ecuación (4.12), mientras que la incer-
88
Calibración de medidores de longitud de onda
Comparados estos valores con la Tabla 24, se observa que a 780 nm, las incerti-
dumbres debidas a la resolución en el sensor temperatura y presión se duplican. Esto es
debido a que la dispersión del aire a esta longitud de onda presenta una variación con la
longitud de onda mucho más rápida que para las longitudes en la 3ª y 4ª ventana de co-
municaciones.
89
Capítulo 4
dad relativa. Estos factores se evalúan en la Tabla 29, junto con la constante de correc-
ción correspondiente:
90
Calibración de medidores de longitud de onda
45
40
35
30
Correción (MHz)
25
20
15
10
5
0
-5
-10
-15
-20
5 6 6
0,0 5,0x10 1,0x10 2,5x10
Tiempo (s)
50
45
40
35
Corrección (MHz)
30
25
20
15
10
-5
33 34 35 36 37 38
Temperatura (oC)
91
Capítulo 4
WM (Figura 48), se ponen en evidencia las carencias de la corrección interna del equi-
po. Las aproximaciones introducidas en el algoritmo de compensación de la dispersión
del aire en los brazos del interferómetro contribuyen de manera sistemática a la correc-
ción de la medida.
45
40
35
30
25
Corrección (MHz)
20
15
10
5
0
-5
-10
-15
-20
5 6 6
0,0 5,0x10 1,0x10 2,5x10
Tiempo (s)
6. Conclusiones
92
Calibración de medidores de longitud de onda
93
Capítulo 5:
Calibración de interrogadores de redes de
Bragg
1. Introducción
2. Principio de operación de los sensores basados en redes de Bragg en fibra
2.1 Introducción
2.2 Funcionamiento de las redes de Bragg como sensores de tensión y tempera-
tura
2.3 Esquemas de detección de la longitud de onda de Bragg
3. Esquemas propuestos para la calibración de interrogadores
3.1 Método 1: Red de Bragg simulada
3.2 Método 2: Celdas de absorción
4. Demostración experimental
4.1 Dispositivos bajo análisis
4.2 Método 1: Red de Bragg simulada
4.3 Método 2: Celdas de absorción
5. Resultados y discusión
5.1 Método 1: Red de Bragg simulada
5.2 Método 2: Celdas de absorción
6. Conclusiones
95
Calibración de interrogadores de redes de Bragg
Las redes de Bragg en fibra óptica (Fiber Bragg Grating o FBG) se emplean ca-
da vez con mayor frecuencia en el sensado cuasi-distribuido de temperatura, deforma-
ciones, vibraciones, etc., ya que superan a otras tecnologías tradicionales en términos de
ventajas mecánicas, posibilidades de multiplexación y operación remota [Kashyap’09].
Las FBG se fabrican inscribiendo en el núcleo de la fibra óptica una modulación perió-
dica del índice de refracción, comportándose en reflexión como un filtro de banda estre-
cha [Krohn’14]. La longitud de onda central de la banda reflejada, conocida también
como longitud de onda de Bragg, depende del período de la red inscrita en el núcleo y
del índice de refracción efectivo. Cualquier modificación en dicho índice o en el período
causado por efectos externos, tales como deformación o temperatura, resulta en un cam-
bio de la longitud de onda de Bragg, permitiendo su uso para aplicaciones de sensado.
En particular, las aplicaciones de las FBG abarcan desde sensores in situ en la industria
médica, para monitorizar funciones biológicas, detectar y tratar tumores y para el cuida-
do post-traumatológico [Pinet’08, Poeggel’15]; hasta la supervisión de la salud estructu-
97
Capítulo 5
Para superar estas limitaciones, en este trabajo presentamos dos técnicas para la
calibración de interrogadores de FBG. Estas técnicas propuestas son la aplicación de los
resultados previamente obtenidos durante la tesis, los cuales permiten la calibración de
estos equipos de forma absoluta y con trazabilidad a la unidad de tiempo del SI. El pri-
mer método utiliza una red de Bragg simulada, donde la longitud de onda reflejada es
medida con el WM patrón calibrado directamente con el OFC en el capítulo 4 de esta
tesis; en el segundo método se realiza la medida de la longitud de onda de absorciones
de dos celdas de gases: una calibrada por el NIST y la segunda calibrada tal como se
presentó en el capítulo 2 de esta tesis. En este capítulo se demuestran experimentalmen-
te estos métodos sobre cuatro sistemas comerciales.
2.1. Introducción
Una FBG consiste en la modulación del índice de refracción de su núcleo en una
región que puede extenderse típicamente desde 1 mm hasta 20 mm [Krohn’14]. La am-
plitud de esta modulación del índice de refracción (n) en la región de la red de Bragg
se encuentra típicamente en el rango de 10–5 hasta 10–3. Esta estructura periódica, mos-
trada esquemáticamente en la Figura 50, produce la reflexión de toda radiación cuya
longitud de onda cumpla la condición de Bragg [Kersey’97]:
B 2neff (5.1)
donde B es la longitud de onda de Bragg reflejada por la red, neff es el índice de refrac-
ción efectivo del núcleo de la fibra, y es el periodo de la red.
98
Calibración de interrogadores de redes de Bragg
Figura 50: Red de Bragg inscrita en el núcleo de una fibra óptica con un es-
paciado .
Cuando una FBG es iluminada con una fuente de radiación de espectro ancho,
ésta refleja una componente espectral estrecha y centrada en B, mientras que el resto de
componentes espectrales atraviesan la red, como se observa en la Figura 51. Mediante la
selección de un periodo distinto para cada red, es posible multiplexar varias FBG cen-
tradas en distintas longitudes de onda en una única fibra óptica, lo que permite desarro-
llar esquemas de sensado quasi-distribuido.
I FBG
I
B
99
Capítulo 5
dneff
n 2 T
P12 P11 P12
dT
B 2neff 1 (5.2)
eff
2
neff
Para el caso de fibras cuyo núcleo es sílice el término que involucra a los coefi-
cientes piezoeléctricos de la ecuación (5.2) es aproximadamente igual a:
n2
1
eff
P12 P11 P12 0,78·10 6 ( 1 ) (5.3)
2
Así pues, para el caso de temperatura constante, la ecuación (5.2) se reduce a la siguien-
te expresión:
Esta relación implica que para una longitud de onda de 1550 nm, 1 nm de variación en
100
Calibración de interrogadores de redes de Bragg
1 dneff
B 2neff T (5.5)
n dT
eff
Dado que:
1 dneff
(5.6)
neff dT
1 dneff dneff
B 2neff T B 1 T (5.7)
n dT n dT
eff eff
1 dneff
6,67·10 6 (º C -1 ) (5.8)
n dT
eff
101
Capítulo 5
Otros posibles métodos para compensar la influencia cruzada utilizan dos FBG
colocadas en un mismo punto pero con diferentes respuestas a temperatura y tensión
[Kersey’97, Xu’94]. Si denotamos por K1 y K2 las respuestas a la tensión de cada una
de las redes, y por KT1 y KT2 sus respuestas a la temperatura, se puede escribir la si-
guiente ecuación matricial:
1 K 1 KT 1
(5.10)
2 K 2 KT 2 T
Este conjunto de dos ecuaciones admiten solución siempre que la matriz sea invertible,
lo que ocurre si la razón entre las respuesta de los dos sensores a la tensión es diferente
de la razón entre las respuestas a la temperatura. Esta diferencia en responsividad se
puede obtener reduciendo el diámetro del recubrimiento de una de las redes de forma
que la sensibilidad a la tensión aumenta. Una variante de este método utiliza dos FBG
con B muy separados, por ejemplo 850 nm y 1300 nm [Xu’94], permitiendo la deter-
minación de tensión y temperatura con precisión de ± 10 y ± 5 ºC respectivamente.
La combinación con redes de Bragg de largo período también se puede utilizar para la
discriminación de temperatura y tensión [Patrick’96, Bhatia’97].
Para la elección de uno u otro método, se debe tener en cuenta el mayor ancho
de banda de los filtros ópticos respecto a los láseres sintonizables, la dependencia de la
red con la polarización de la radiación óptica y la posibilidad de la aparición de efectos
de interferencia producidos cuando la fuente utilizada es altamente coherente. Si la apli-
cación particular no requiere de la medida absoluta de la longitud de onda B, sino de la
102
Calibración de interrogadores de redes de Bragg
103
Capítulo 5
Con respecto al uso de FBG calibradas como referencia, nuestro método presen-
ta varias ventajas: la forma del espectro se mantiene constante en todo el rango de me-
dida y no se ve afectada por la deformación y la temperatura, y la forma gaussiana re-
sultante permite una determinación unívoca del pico, ya que el pico en longitud de onda
y la longitud de onda central del espectro reflejado coinciden. Esta característica hace
que nuestro método de calibración sea adecuado para todos los tipos de interrogadores
104
Calibración de interrogadores de redes de Bragg
INTERROGADOR
FBG
ESPEJO
CÉLULA DE ABSORCIÓN
105
Capítulo 5
4. Demostración experimental
106
Calibración de interrogadores de redes de Bragg
Con este método se han calibrado únicamente los FBGI estáticos SM125 fabri-
cados por Micron Optics, ya que el software del equipo fabricado por HSM no permite
detección de valles.
107
Capítulo 5
5. Resultados y discusión
108
Calibración de interrogadores de redes de Bragg
K U (k = 2)
/pm /pm
SM125 (unidad 1) –0,9 1,1
FS42 –0,6 1,1
SM130 16,4 1,1
Tabla 32: Corrección media e incertidumbre para cada uno de los tres equi-
pos comerciales calibrados.
5 5
3 3
Corrección (pm)
2 2
Corrección (pm)
1 1
0 0
-1 -1
-2 -2
-3 -3
-4 -4
-5 -5
1510 1520 1530 1540 1550 1560 1570 1580 1590 1600 1500 1510 1520 1530 1540 1550 1560 1570 1580 1590 1600
Longitud de onda (nm) Longitud de onda (nm)
109
Capítulo 5
22
Corrección (pm)
18
16
14
12
1510 1520 1530 1540 1550 1560 1570 1580 1590 1600
110
Calibración de interrogadores de redes de Bragg
111
Capítulo 5
5 5
4 (a) SM125_unidad 1 4
(b) SM125_unidad 2
3 3
2 2
Corrección (pm)
Corrección (pm)
1 1
0 0
-1 -1
-2 -2
-3 -3
-4 -4
-5 -5
1510 1520 1530 1540 1550 1560 1570 1580 1590 1510 1520 1530 1540 1550 1560 1570 1580 1590
Longitud de onda (nm ) Longitud de onda (nm)
6. Conclusiones
112
Calibración de interrogadores de redes de Bragg
tiempo del SI: a través del WM calibrado en el primer método y a través de las celdas
calibradas con trazabilidad al patrón del segundo en el último método.
En el primer método se utiliza una FBG simulada compuesta por un filtro óptico
sintonizable y un espejo. La longitud de onda de esta FBG simulada se calibra frente a
un WM de referencia previamente calibrado con un OFC. Este primer método es aplica-
ble a todos los sistemas comerciales, independientemente de sus prestaciones software,
pero su incertidumbre (±1,1 pm) está limitada por el ancho espectral del filtro sintoniza-
ble y la estabilidad de la longitud de onda del filtro óptico. Los resultados del primer
método ponen de manifiesto la importancia de incorporar en los FBGI una celda de re-
ferencia para la autocalibración durante la operación del equipo, lo que permite obtener
correcciones por debajo de la incertidumbre de calibración.
113
Conclusiones
115
Conclusiones
Conclusiones
La presente tesis ha proporcionado nuevas referencias frecuenciales de alta cali-
dad para su uso en la calibración en tiempo y frecuencia de equipamiento óptico, de es-
pecial relevancia en las áreas de comunicaciones ópticas y sensado distribuido.
117
Conclusiones
ancho de banda de emisión del peine, utilizando como bombeo dos fuentes láser estabi-
lizadas en la pendiente de dos absorciones moleculares de la celda de acetileno de refe-
rencia (12C2H2). Las anchuras de línea de las referencias ópticas sintetizadas estaban por
debajo de los 75 kHz, mientras que la estabilidad de las señales sintetizadas resultaron
ser de 4·10–11 para tiempos de observación cortos.
118
Anexo I:
Procedimiento de caracterización de fuentes
láser con un peine de frecuencias
1. Introducción
2. Determinación de los signos de la frecuencia de offset y de batido
3. Deternimación del orden del modo N en la medida de una frecuencia láser
119
Procedimiento de caracterización de fuentes láser con un peine de frecuencias
Consideremos fX la frecuencia del láser bajo test, y N el orden del modo del OFC
más próximo a fX. Como puede observarse en la Figura 58, podemos calcular fX a partir
de las frecuencias principales del OFC (la frecuencia de repetición frep y la frecuencia de
offset f0) y la frecuencia de batido entre el láser y el modo N:
Las frecuencias de repetición y offset son conocidas por la propia configuración del
OFC, mientras que la frecuencia de batido se obtiene de manera directa en el fotodetec-
tor en el que se hacen interferir la emisión del OFC y la señal bajo test. No obstante, pa-
121
Anexo I
ra poder calcular fX sigue siendo necesario determinar el valor de N y los dos signos de
la ecuación I.1. El proceso para resolver estas incógnitas se detalla en los epígrafes 2 y 3
del presente anexo.
Los signos de f 0 y f beat en la ecuación (I.1) deben ser determinados para cada
medida, analizando para ello el comportamiento de la frecuencia f beat cuando se aplica
un cambio en la frecuencia de repetición o en la frecuencia de offset. Un cambio en la
frecuencia f rep induce un cambio en la separación de los modos del OFC, semejante al
fuelle de un acordeón, donde uno de los extremos está fijo en la frecuencia f0 (Figura
122
Procedimiento de caracterización de fuentes láser con un peine de frecuencias
59).
En caso de modificar f0, el espaciado entre los modos del OFC permanece inalte-
rado y es todo el OFC el que se mueve en conjunto hacia frecuencias mayores o meno-
res, dependiendo del signo de f0 y la variación introducida. Por ejemplo, si f 0 es positi-
va, un aumento f 0 también positivo inducirá un desplazamiento del OFC hacia fre-
cuencias mayores (Figura 60).
123
Anexo I
Estas dependencias pueden aprovecharse para determinar el signo de fbeat y f0,
observando el comportamiento de fbeat ante pequeñas variaciones de los parámetros de
configuración del OFC. En primer lugar, se parte de un valor de frep utilizado para la
medida y se aumenta ligeramente. Si la frecuencia láser se encuentra a la derecha del
modo con el cual se está batiendo (Figura 61a), el aumento de la frep producirá que el
modo en cuestión se acerque al láser, con lo que la frecuencia de batido disminuirá. En
este caso el signo de la frecuencia de batido será positivo. En el caso de que el láser esté
a la izquierda del diente más cercano (Figura 61b), el aumento de la frep producirá que el
modo en cuestión se aleje del láser, con lo que la frecuencia de batido aumentará.
Caso 1: Supongamos que fbeat es positiva (luego el láser está a la derecha del
modo) y aumentemos la f0 en una pequeña cantidad. Si la fbeat disminuye, implica que el
modo en cuestión está acercándose a la frecuencia láser desconocida, luego el OFC está
siendo desplazado hacia frecuencias mayores. En este caso la f0 tendrá el signo positivo.
Caso 2: Supongamos que fbeat es de nuevo positiva (luego el láser está a la dere-
cha del modo) y aumentemos f0 en una pequeña cantidad. Si fbeat aumenta, implica que
el modo en cuestión está alejándose de la frecuencia láser desconocida, luego el OFC
está siendo desplazado hacia frecuencias menores. En este caso f0 tendrá el signo nega-
tivo.
Caso 3: Supongamos que fbeat es negativa (luego el láser está a la izquierda del
modo) y aumentemos f0 en una pequeña cantidad. Si fbeat disminuye, implica que el mo-
do en cuestión está acercándose a la frecuencia láser desconocida, luego el OFC está
siendo desplazado hacia frecuencias menores. En este caso f0 tendrá el signo negativo.
Caso 4: Supongamos que fbeat es de nuevo negativa (luego el láser está a la iz-
quierda del modo) y aumentamos f0 en una pequeña cantidad. Si fbeat aumenta, implica
que el modo en cuestión está alejándose de la frecuencia láser desconocida, luego el
OFC está siendo desplazado hacia frecuencias mayores. En este caso f0 tendrá el signo
positivo.
Todos los casos posibles que se pueden dar en la determinación de los signos de
fbeat y f0, se encuentran resumidos en la Tabla 36.
124
Procedimiento de caracterización de fuentes láser con un peine de frecuencias
fbeat f0
+ Si f0 aumenta y fbeat disminuye +
Si frep aumenta y fbeat disminuye
Si f0 aumenta y fbeat aumenta –
– Si f0 aumenta y fbeat aumenta +
Si frep aumenta y fbeat aumenta
Si f0 aumenta y fbeat disminuye –
El valor del modo, N, es un número entero muy elevado, cuyo orden puede esti-
marse como:
fX
N (I.2)
f rep
Por ejemplo, para una longitud de onda en torno a 1550 nm, y una frecuencia de repeti-
ción de 100 MHz, el orden del modo será del orden de 2·106. Para determinar su valor
exacto, existen dos alternativas. Si tenemos un conocimiento previo aproximado de la
frecuencia del láser bajo test con una resolución menor que la mitad de la frecuencia de
repetición (por ejemplo mediante un WM de alta precisión, o si el láser está anclado en
alguna absorción conocida), es suficiente con redondear el valor obtenido con la si-
guiente ecuación:
fWM f 0 f beat
N (I.3)
f rep
donde fWM es la frecuencia aproximada del láser obtenida con un WM, por ejemplo.
Asimismo, para una segunda frecuencia de repetición, frep2, el orden del modo del OFC
será (N + m) si frep1 es mayor que frep2, o (N – m) si es menor, obteniéndose una nueva
frecuencia de batido fbeat2. Para el caso en que frep1 es mayor que frep2, obtenemos:
125
Anexo I
f X N m f rep2 f 0 f beat2 (I.5)
2
uN (I.7)
f rep1 f rep 2
Esta ecuación implica que cuanto mayor sea la incertidumbre en la medida de fbeat, ma-
yor debe ser la diferencia entre frep1 y frep2 para que el valor de la incertidumbre (I.7) se
mantenga por debajo de 0,5, asegurando así la determinación unívoca del modo:
2
f rep1 f rep 2 4 (I.8)
0,5
126
Procedimiento de caracterización de fuentes láser con un peine de frecuencias
127
Anexo II:
El peine de frecuencias del Instituto de ópti-
ca “Daza de Valdés”
1. Introducción
2. El láser de femtosegundos
3. Estabilización del peine de frecuencias ópticas
3.1 Introducción
3.2 Estabilización de la frep del peine de frecuencias ópticas
3.2.1 El módulo de detección de fase (PHD110)
3.2.2 El módulo de medida (MB110)
3.2.3 El módulo regulador (PID110)
3.2.4 El módulo de control del motor (MC110)
3.2.5 El sintetizador de frecuencias (DDS120)
3.3 Estabilización de la f0 del peine de frecuencias
3.3.1 Unidad electrónica de amplificación (OFD100)
3.3.2 El detector de fase digital (DXD200)
3.3.3 La unidad de realimentación (PIC210)
4. El contador de frecuencias (FXM50)
5. Unidad de suministro de señales de radiofrecuencia (RFD10)
6. Las unidades de batido
6.1 La unidad de detección de batido en aire BDU-FS
6.2. La unidad de detección de batido todo en fibra BDU-OF
6.3 La unidad electrónica de amplificación BFD100
129
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
131
Anexo II
EDFA Y INTERFEROMETRO
GENERADOR f-2f
SUPERCONTINUO
BOMBEO OSCILADOR
EDFA BDU-FS
EN FIBRA
LASER
DESCONOCIDO
Figura 62: Esquema global del OFC auto-referenciado del IO-CSIC. Las
líneas azules corresponden a conexiones internas de fibra óptica, las líneas
rojas a radiación óptica en aire y las líneas negras a conexiones eléctricas.
2. El láser de femtosegundos
Los láseres de fibra óptica se basan en osciladores con fibras ópticas dopadas,
principalmente mediante tierras raras como Er e Yb [Dong’16, Agrawal’08]. La genera-
ción de pulsos en láseres de fibra óptica se consigue mediante técnicas de anclado de
modos activos con moduladores de ganancia [Hanna’89, Takada’90]; o pasivos median-
te diversos absorbentes saturables [Loh’93, Abraham’93, Yamashita’05], interferóme-
132
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
Diodos de Diodos de
Generador de Control de bombeo bombeo Control de
supercontinuo dispersión dispersión
HNLF variable variable
WDM WDM
Amp.1 Er:Fiber Er:Fiber
Amp.2
Divisor 50/50
Diodo de
bombeo Acoplador de salida
WDM MRO Rejection
Monitor
Oscilador de Er Er:Fiber
Elemento
Mode-Locking
Port
El primer amplificador (Amp. 1), está diseñado para proporcionar máxima po-
tencia y mínima duración de pulso a fin de generar el espectro de una octava necesario
para la detección de la frecuencia f0. Está bombeado por dos diodos láser a 980 nm con
una potencia total máxima de 980 mW. A la salida del amplificador, la potencia óptica
media es de 250 mW, correspondiente a 2,5 nJ de energía por pulso. Presenta un ancho
de banda de 60 nm centrado en 1550 nm. El ensanchamiento de pulso introducido por el
amplificador se compensa mediante el sistema de prismas mostrado en la Figura 64
[Fork’84]. Esta alternativa presenta pérdidas más bajas que las soluciones basadas en
redes de difracción [Treacy’69], y permite obtener pulsos con una duración de alrededor
133
Anexo II
de 85 fs.
prisma 2
Espejo
prisma 1
1
Potencia normalizada (dB)
0,1
0,01
1000 1100 1200 1300 1400 1500 1600 1700
Longitud de onda (nm)
134
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
0,1
0,01
1480 1500 1520 1540 1560 1580 1600 1620 1640
Longitud de onda (nm)
El OFC presenta una tercera salida en fibra, denominada MRO rejection port.
Esta salida recibe alrededor del 6% de la potencia óptica total existente en el oscilador,
extraída mediante el correspondiente acoplador. Su espectro de salida se muestra en la
Figura 67.
-5
-10
Potencia normalizada (dB)
-15
-20
-25
-30
-35
-40
1480 1500 1520 1540 1560 1580 1600 1620 1640
Longitud de onda (nm)
135
Anexo II
3.1. Introducción
La estabilización del OFC requiere de la estabilización de las dos radiofrecuen-
cias que definen su estructura global de modos: la frecuencia de repetición y la frecuen-
cia de offset. Ambas frecuencias se estabilizan utilizando como referencia un reloj ató-
mico de Rb. La estabilización de la frep es la más sencilla de llevar a cabo, bastando con
comparar en un detector de fase la diferencia entre la frep del láser de femtosegundos con
la frecuencia objetivo establecida por un sintetizador de radiofrecuencias. La señal de
error servirá para variar la longitud total de la cavidad resonante de forma que la dife-
rencia entre ambas frecuencias sea nula. Sin embargo, la estabilización de f0 precisa de
su caracterización mediante un interferómetro f-2f [Jones’00].
Phase out
RF1RF2
Ordenador
Offset
136
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
PHD 110
Detector de fase
RF1
RF2
in
Busy out
137
Anexo II
provenientes del módulo PHD110, así como el voltaje suministrado al actuador piezo-
eléctrico de la cavidad láser. Todas estas señales son digitalizadas y enviadas al ordena-
dor que controla el sistema mediante la conexión USB.
MB 110
Módulo de medida
Output
USB
Input
En el caso del módulo regulador PID110 (Figura 71), la señal que se pretende
minimizar es la diferencia de fase relativa entre la frecuencia de repetición del OFC y la
frecuencia de referencia proporcionada por el sintetizador. De las cuatro conexiones
presentes en el frontal del módulo, únicamente se utilizan en este caso la entrada del fo-
todetector y la salida. La entrada del fotodetector se conecta con la salida de fase del
módulo PHD110 mientras que la salida del módulo se conecta con la entrada de alimen-
tación del actuador piezoeléctrico que controla el movimiento de uno de los colimadores
fibra-aire de la cavidad resonante.
138
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
PID 110
Regulador PID
ON min. int.
inten. set
OFF
HV/ inp.
P 100
Out err.
I mod. reg. reg. inp.
D monitor
ON ON
OFF OFF
reg. mod. in
3 4 5
outp. range 2 6
1 7
set point
0 8
pos
gain
neg
mon. out
output
139
Anexo II
MC 110
Módulo de control
de motor
USB
motor
140
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
ON
OFF
DDS SYNTHESIZER
FREQUENCY MAX 120 [MHz]
100.000.000,000,00 Hz
PHASE AMPLITUDE
OUT
141
Anexo II
en donde las partes fundamentales son el cristal (C) de Niobatio de Litio periódicamente
pulido (PPLN), una red de difracción (G) y un fotodetector PIN (D) modelo FPD150-F
fabricado por Menlo Systems. La señal es introducida a través de la apertura de entrada,
y guiada hacia el cristal no lineal mediante un primer espejo (M1). La lámina retardado-
ra de media onda (W) permite ajustar la distribución de potencia entre el detector y la
apertura de salida auxiliar. A ambos lados del PPLN se sitúan dos lentes asféricas de 25
mm de focal (L1, L2), de tal forma que la ubicación del cristal coincide con los puntos
focales de las lentes. La primera lente enfoca la radiación óptica sobre el cristal a fin de
conseguir la máxima eficiencia en la generación del espectro doblado. La segunda lente
tiene como función colimar la radiación de salida del cristal. Un cubo divisor por pola-
rización (PBS) divide el haz en dos: uno de los haces es guiado por un espejo M3 hacia
el detector, mientras que el segundo haz es dirigido por un segundo espejo (M2) hacia la
apertura de salida. La señal de batido se enfoca sobre el fotodetector mediante una ter-
cera lente (L3). La red de difracción G selecciona el rango espectral alrededor de 1050
nm donde la zona inferior del espectro original se superpone con la zona superior del
espectro doblado. Sobre el fotodetector se produce el batido heterodino, dando como re-
sultado una frecuencia de batido igual a f0.
OSCILOSCOPIO
ESA
Contador
FAST IN
MONITOR Ext. 10 MHz
Offset
ON
MONITOR Invert 20 MHz
Dither OFF MONITOR
PHASE DOWN
Offset out IN 20 MHz
SLOW In OUT
COUNTER UP
142
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
Atenuador APD
50 MHz variable
1:100 MONITOR
PHASE
1:3 COUNTER
20 MHz
DC Detector
RF
143
Anexo II
UP
Contador
DOWN
DAC
SLOW
20 kHz
Buffer FAST
MONITOR
50 kHz
Los tres puertos de salida de la unidad proporcionan señales con voltajes entre
–3 V y +3 V. La salida Fast suministra una señal sin filtrar que es introducida en el bu-
cle de realimentación PIC210 como señal de error. La salida Monitor está filtrada por
un filtro electrónico pasa-baja con frecuencia de corte en 50 MHz y es utilizada para
monitorizar f0 en el contador de frecuencias.
144
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
beo.
CONTROLES
CONECTORES
Offset On/off
Ext. Offset
Offset
Gain coarse In
1 Gain
Invert On/off Coarse Monitor
HF Gain
HF Gain Dither
Gain Fine
Gain Fine
Out
Invert On/off
Offset Out
145
Anexo II
señales introducidas simultáneamente por sus cuatro canales (Figura 79), utilizando los
pulsos de contaje generados en la unidad RFD10 e introducidos por el puerto MP. Las
lecturas de frecuencias de todos los canales son transmitidas al ordenador, junto con la
descripción de las características del pulso de disparo, por una conexión RS232. La re-
solución en la medida de la frecuencia en cada uno de los cuatro canales del contador es
de 1 mHz.
B
A 1- (II.1)
1000
146
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
Indicador
de estado
+ 20 dB
Ref. In
10 MHz 10 MHz
1:3 10 MHz
+ 20 dB
x2 20 MHz
20 MHz
+ 20 dB
x2 20 MHz
20 MHz
El OFC cuenta con dos unidades de batido para metrología de frecuencias: una
unidad de batido en aire (BDU-FS, del inglés Beat Detection Unit Free-Space) y una
unidad de batido todo en fibra (BDU-OF, del inglés Beat Detection Unit Optical-Fiber).
En ambos casos, la señal eléctrica resultante es introducida a un módulo electrónico de
amplificación modelo BFD100 fabricado por Menlo Systems.
147
Anexo II
HWP2
L
D
Entrada del
M2 láser
M1 desconocido
La señal obtenida en el detector contiene información del batido del láser bajo
test con todos los dientes del peine, dentro de la zona espectral seleccionada por la red
de difracción G. A fin de clarificar las diferentes frecuencias involucradas en la medida,
en la Figura 83 se presenta un ejemplo experimental de señal de batido observada en un
ESA. La señal comprende dos batidos muy intensos a 100 MHz y 200 MHz, correspon-
dientes a los batidos entre modos del peine. Entre estos dos picos, se observan dos bati-
dos menos intensos, a 30 MHz y 70 MHz respectivamente del primer batido máximo. El
primero de éstos es la medida deseada, es decir, el batido del láser bajo test con el modo
N del peine (fbeat), mientras que el segundo resulta del batido del láser con el modo de
índice (N + 1) del peine.
-10
-30
fbeat,2=((n +1)frep+ f0) -flaser
-40
-50
-60
-70
100 110 120 130 140 150 160 170 180 190 200
Frecuencia (MHz)
148
El peine de frecuencias del Instituto de Óptica “Daza de Valdés”
Aislador
Atenuador APD
50 MHz variable
1:100 MONITOR
PHASE
1:3 COUNTER
30 MHz
Detector
DC
RF
149
Referencias
Referencias
A
151
Referencias
[Birch’93] K. P. Birch, M. J. Downs, “An Updated Edlén Equation for the Refractive
Index of Air”, Metrologia 30, 155 (1993).
[Brabec’00] Th. Brabec, F. Krausz, “Intense few-cycle laser fields: Frontiers of nonli-
near optics”, Rev. Mod. Phys. 72, 545 (2000).
152
Referencias
[CMCs] http://kcdb.bipm.org/appendixC/default.asp.
153
Referencias
[Derikson’98] D. Derikson, “Fiber Optic Test and Measurement”, Prentice Hall PTR
(1998).
154
Referencias
[Di Sante’15] R. Di Sante, “Fibre Optic Sensors for Structural Health Monitoring of
Aircraft Composite Structures: Recent Advances and Applications”, Sensors (Basel) 15,
18666 (2015).
[Edlén’53] B. Edlen, “Dispersion of standard air”, J. Opt. Soc. Am. 43, 339 (1953).
155
Referencias
standard by use of a combined frequency chain and femtosecond comb”, Opt. Lett. 29,
566 (2004).
[Essen’55] L. Essen, J. V. L. Parry, “An atomic standard of frequency and time interval:
a cæsium resonator”, Nature 176, 280 (1955).
156
Referencias
[Gardner’05] F. M. Gardner, “Phaselock Techniques (3rd Edition)”, John Wiley & Sons
(2005).
[Hall’76] J. L. Hall, S. A. Lee, “Interferometric real time display of CW dye laser wave-
length with subDoppler accuracy”, Appl. Phys. Lett. 29, 367 (1976).
[Hall’06] J. L. Hall, “Nobel Lecture: Defining and measuring optical frequencies”, Rev.
Mod. Phys. 78, 1279 (2006).
157
Referencias
[Hänsch’06] T. W. Hänsch, “Nobel Lecture: Passion for precision”, Rev. Mod. Phys.
78, 1297 (2006).
158
Referencias
[Jones’97] J. D. Jones, “Review of fiber sensor techniques for temperature strain dis-
crimination,” 12th annual conference in Optical Fiber Sensors, OSA Technical Digest
Series 16, paper OTuC1 (1997).
159
Referencias
[Kersey’92a] A. D. Kersey, “Multiplexed fiber optic sensors”, Proc. SPIE 1797, 161
(1992).
160
Referencias
EFTF and the IEEE International Frequency Control Symposium, Besançon (France) 2,
651 (1999).
[Krohn’14] D. A. Krohn, “Fiber optic sensors: fundamentals and applications (4th ed)”,
Society of Photo-Optical Instrumentation Engineers (2014).
161
Referencias
88
Sr+ clock transition using a femtosecond optical frequency comb”, Phys. Rev. A 67,
032501 (2003).
162
Referencias
[Proctor’92] B. Proctor, F. Wise, “Quartz prism sequence for reduction of cubic phase
in a mode-locked Ti.Al2O3 laser”, Opt. Lett. 17, 1295 (1992).
163
Referencias
[Russell’03] P. St. J. Russell, “Photonic crystal fibers”, Science 299, 358 (2003).
164
Referencias
[Shiller’02] S. Schiller, “Spectrometry with frequency combs”, Opt. Lett. 27, 766
(2002).
165
Referencias
[Takada’90] A. Takada, H. Miyazawa, “30 GHz picosecond pulse generation from ac-
tively mode-locked erbium-doped fibre laser”, Electron. Lett. 26, 216 (1990).
166
Referencias
167
Referencias
[Ye’01] J. Ye, L. S. Ma, J. L. Hall, “Molecular Iodine clock”, Phys. Rev. Lett. 87,
270801 (2001).
168
Referencias
169
Acrónimos
Acrónimos
ADC: Analog-to-Digital Converter
171
Acrónimos
HARPS: High Accuracy Radial Velocity Planet Searcher.
NIR: Near-infrared.
PID: Proportional-Integral-Derivative.
PIN: P-type/Intrinsic/N-type.
172
Acrónimos
SMA: SubMiniature version A.
173
Publicaciones
Publicaciones
Presentaciones en congresos
Presentaciones orales:
(6) J. Galindo-Santos*, M. Alcon-Camas, S. Martin-Lopez, A. Carrasco-Sanz, P.
Corredera, “Application of Brillouin scattering to optical frequency combs”,
Photonics Europe. Bruselas (Bélgica), 2012.
PUBLICACIÓN: Proc. of SPIE Vol. 8434, 843410-1.
175
Publicaciones
*
(7) J. Galindo-Santos , M. Alcon-Camas, A. Carrasco-Sanz, S. Martin-Lopez, P.
Corredera, “Filtering of an Optical Frequency Comb with Stimulated Brillouin
Scattering”, EOS Annual Meeting. Aberdeen, Scotland, (Gran Bretaña), 2012.
PUBLICACIÓN: ISBN 978-3-9815022-4-4.
Posters:
(13) J. Galindo-Santos*, A. Carrasco-Sanz, S. Martin-López, P. Corredera, “Estabi-
lización de un láser de diodo en absorciones saturadas de Rb como patrón de
176
Publicaciones
frecuencia”, VII Reunión Española de Optoelectrónica OPTOEL 2011. San-
tander, (España).
PUBLICACIÓN: Libro de comunicaciones de VII Reunión Española de Op-
toelectrónica, OPTOEL’11 en formato electrónico (ISBN: 978-84-86116-31-
6), S1-7 (2011).
177
Publicaciones
PUBLICACIÓN: Libro de abstracts de Reunión de jóvenes investigadores
SINFOTON: 2ª Feria de Otoño en formato electrónico (https://www.sinfoton-
cm.es/difusión/eventos/), 3 (2016).
178