Página tipo blog en el que voy a publicar mis notas de aprendizaje, en especial de temas como matemáticas, física y quizá algo de programación
Comunicaciones ópticas: la fibra y el LED
Todo lo que hace falta para el primer examen de un curso de comunicaciones ópticas, de la ley de Snell al LED, con un simulador debajo de cada fórmula. La idea no es memorizar: es mover un valor y ver qué le pasa al resultado. Cada simulador muestra además el paso a paso con tus números reemplazados, tal como se escribe en un examen.
Los parámetros marcados con el ícono de enlace (n1, n2, D, λ, L) son una sola fibra compartida: si cambias n1 en el primer simulador, cambian también los modos, la dispersión y el acoplamiento de más abajo. Los botones «Cargar» ponen los datos de los problemas resueltos en clase.
1. La luz dentro del vidrio
El vidrio es casi un dieléctrico perfecto: no conduce y absorbe muy poco, por eso se usa para guiar luz. Dentro de un medio la luz va más lenta, y cuando cruza de un medio a otro se dobla según la ley de Snell (los ángulos se miden desde la normal):
n1sinθ1=n2sinθ2v=nc
Qué es cada cosa:n1 es el índice del medio de donde viene la luz y n2 el del medio al que pasa. θ1 (incidencia) y θ2 (refracción) se miden desde la normal: la recta perpendicular a la frontera en el punto donde llega el rayo, nunca desde la frontera. Si la luz pasa a un medio más denso (n2>n1) se acerca a la normal; si pasa a uno menos denso, se aleja.
Ojo con los nombres: en Snell «1» y «2» significan antes y después de la frontera. En la fibra, en cambio, n1 es siempre el núcleo y n2 el revestimiento. Coinciden solo cuando la luz va del núcleo al revestimiento.
Simulador 0 · La ley de Snell en una sola frontera
Cargar:
Pasa a un medio MÁS denso (n₂ > n₁): se ACERCA a la normal (θ₂ < θ₁)
θ1(incidencia, desde la normal)°
n1(medio de donde viene)
n2(medio al que pasa)
θ225.74°
θcno existe (n₁ ≤ n₂)
potencia reflejada / refractada4.3 % / 95.7 %
¿Qué es cada variable y dónde está en el dibujo? (clic para abrir la tabla)
Variable
Qué es
Se mide desde
Dónde está en el dibujo
Valor
n1
Índice del medio de donde viene la luz
—
Zona de arriba (azul)
1.000
n2
Índice del medio al que pasa la luz
—
Zona de abajo
1.480
normal
Recta perpendicular a la frontera en el punto donde llega el rayo
—
Línea vertical punteada
90° con la frontera
θ1
Ángulo de incidencia
la normal
Arco ámbar, arriba a la izquierda
40.00°
θ1′
Ángulo de reflexión (siempre = θ₁)
la normal
Arco magenta, arriba a la derecha
40.00°
θ2
Ángulo de refracción (el que despejas con Snell)
la normal
Arco verde, abajo a la derecha
25.74°
θc
Ángulo crítico: el θ₁ con el que θ₂ = 90°. Solo existe si n₁ > n₂
la normal
Línea roja punteada arriba a la izquierda
no existe
90∘−θ1
Ángulo con la FRONTERA. Es el error típico: NO se usa en Snell
la frontera
Arco gris pequeño
50.00°
Paso a paso (con tus valores)
n1sinθ1=n2sinθ2⇒sinθ2=n2n1sinθ1=1.4801.000sin(40.0∘)=0.4343θ2=sin−1(0.4343)=25.74∘n1≤n2⇒no hay aˊngulo crıˊtico: la luz siempre puede pasarv1=n1c=3.000×108smv2=n2c=2.027×108sm
La fibra tiene un núcleo de índice n1 rodeado de un revestimiento de índice n2<n1. Si el rayo llega a la frontera con un ángulo mayor que el ángulo crítico, no sale: se refleja entero (reflexión total interna, RTI) y queda atrapado.
n1sinθc=n2sin90∘⇒sinθc=n1n2(solo si n1>n2)
¿Con qué ángulo puede entrar la luz desde el aire para que luego haya RTI? Eso define el cono de aceptación (semiángulo θa) y la apertura numérica, que mide cuánta luz es capaz de captar la fibra. La deducción es pregunta clásica de examen:
Cuidado: θc se mide desde la normal y φ=90∘−θc desde el eje. Y θa es el ángulo máximo de entrada, no el mínimo.
Simulador 1 · ¿Entra o no entra la luz?
Cargar:
Mueve θ y cruza el borde del cono verde lima.
θ(entrada, desde el eje)°
n1(nuˊcleo)
n2(revestimiento)
n0(medio exterior)
θ1=φ(dentro, desde el eje)6.74°
θi=90∘−θ1(en la interfaz)83.26°
θc80.57°
φmax=90∘−θc9.43°
NA0.2425
θa(cono=2θa)14.03° (28.07°)
¿Qué es cada variable y dónde está en el dibujo? (clic para abrir la tabla)
Variable
Qué es
Se mide desde
Dónde está en el dibujo
Valor
n0
Índice del medio exterior (aire = 1)
—
Zona izquierda, antes de la cara
1.00
n1
Índice del núcleo (por donde viaja la luz)
—
Franja azul central
1.480
n2
Índice del revestimiento (n₂ < n₁)
—
Franjas grises arriba y abajo
1.460
normal ①
Perpendicular a la cara de entrada: es paralela al eje
—
Línea punteada horizontal en la entrada
—
normal ②
Perpendicular a la frontera núcleo–revestimiento
—
Línea punteada vertical en el primer rebote
—
θ
Ángulo con que llega la luz desde el aire (Snell ①: n₀ sen θ = n₁ sen θ₁)
normal ① (el eje)
Arco ámbar a la izquierda de la entrada
10.00°
θ1=φ
Ángulo del rayo ya dentro del núcleo
normal ① (el eje)
Arco celeste a la derecha de la entrada
6.74°
θi=90∘−θ1
Ángulo con que el rayo golpea el revestimiento (el que se compara con θc)
normal ②
Arco magenta en el primer rebote
83.26°
θc
Ángulo crítico: el mínimo θᵢ para que haya RTI (Snell ②: n₁ sen θc = n₂ sen 90°)
normal ②
Línea roja punteada en el primer rebote
80.57°
θ2
Ángulo de lo que escapa al revestimiento (solo si θᵢ < θc)
normal ②
Rayos verdes que salen al revestimiento
no hay (RTI)
θa
Máximo θ que todavía se guía: semiángulo del cono de aceptación. NA = sen θa
normal ① (el eje)
Cono verde lima a la izquierda
14.03°
φmax=90∘−θc
Máximo ángulo con el eje que puede tener un rayo guiado
el eje
Es el θ₁ que sale cuando θ = θa
9.43°
Paso a paso (con tus valores)
n0sinθ=n1sinθ1⇒θ1=sin−1(1.4801.00⋅sin10.0∘)=6.74∘en la interfaz: θi=90∘−θ1=83.26∘vsθc=sin−1(n1n2)=sin−1(1.4801.460)=80.57∘83.26∘≥80.57∘⇒reflexioˊn total interna: se guıˊaNA=n12−n22=2.1904−2.1316=0.2425θa=sin−1(n0NA)=14.03∘
2. Tipos de fibra y número de modos
Por el modo de propagación la fibra es monomodo (núcleo de unos 8–10 µm, un solo camino posible; la de las redes troncales y largas distancias) o multimodo (núcleo de 50 o 62.5 µm, muchos caminos). En casi todas, el revestimiento mide 125 µm. Por el perfil de índice, la multimodo puede ser de índice escalonado (SI, step index) o de índice gradual (GI, graded index):
Índice escalonado (SI): n₁ constante, el rayo rebota en zig-zag y recorre más camino a la misma velocidad → llega tarde.
Índice gradual (GI): n baja hacia el borde; el rayo que se aleja del eje va por vidrio más «rápido» (v = c/n) y compensa su camino largo → todos llegan casi juntos. Más caro.
Cuántos modos caben lo dice la frecuencia normalizadaV: crece con el diámetro y la NA, y baja con la longitud de onda. Por debajo de 2.405 solo cabe uno.
V=λπDNAN≈2V2monomodo⟺V<2.405
Simulador 2 · ¿Cuántos modos caben?
Cargar:
D(diaˊmetro del nuˊcleo)µm
λµm
n1
n2
NA0.2425
V63.48
N≈V2/22015
Dmaxmonomodo1.89 µm
A escala. Cada punto ámbar ≈ un modo (se dibujan 1500 de 2015).
Paso a paso (con tus valores)
V=λπDNA=0.6000μmπ(50.00μm)(0.2425)=63.483N≈2V2=2(63.483)2≈2015modosmonomodo si V<2.405⇒D<πNA2.405λ=π(0.2425)2.405(0.6000)=1.894μm
3. Dispersión modal y ancho de banda
En una fibra de índice escalonado todos los modos viajan a la misma velocidad c/n1, pero no recorren lo mismo: el modo axial va recto y llega primero; el que rebota en zig-zag con el ángulo crítico recorre más y llega después. El pulso que sale del extremo lejano llega ensanchado.
En clase se calcula contando triángulos: cada cruce del núcleo es un triángulo con un cateto de avance, el diámetro D y una hipotenusa. Como hay tantos catetos como hipotenusas, Lzz=(L/cateto)⋅h. El atajo Ln1/n2 da lo mismo y muestra algo importante: el camino extra no depende del diámetro.
BW=Δt1(algunos libros: 2Δt1)
Simulador 3 · Dispersión modal: el zig-zag llega tarde
Cargar:
Lkm
n1
n2
Dµm
B(tasa de bits)Mb/s
taxial=Ln1/c4.9333 µs
tzig-zag5.0009 µs
Lzz−L(camino extra)13.699 m
Δt67.58 ns
Δtpor km67.58 ns/km
BW14.80 MHz
BvsBWB ≤ BW ✓
Sube la tasa de bits o la longitud: cuando Δt se acerca a la duración del bit, los pulsos se montan unos sobre otros y los ceros dejan de verse como ceros. Eso es lo que limita el ancho de banda. Cambia a índice gradual y mira cuánto más aguanta.
4. Pérdidas de acoplamiento: meter la luz en la fibra
Antes de recorrer un solo metro ya se pierde luz en la unión fuente–fibra. Son tres pérdidas, y su suma es la pérdida de acoplamiento de entrada (ICL, input coupling loss):
No intercepción: si el área que emite la fuente As es mayor que el núcleo Ac, lo que cae fuera se pierde. Si la fuente es más pequeña, esta pérdida no existe.
Apertura numérica: la fibra solo acepta los rayos dentro de su cono. Para un LED de superficie (lambertiano) la fracción aceptada es exactamente NA2: va al cuadrado.
Reflexión (Fresnel): una parte rebota en la cara del vidrio.
Todas son negativas en dB, así que se suman: −12.04 dB y −0.2 dB dan −12.24 dB, no −12.02.
Simulador 4 · Pérdidas al meter la luz en la fibra (ICL)
Cargar:
Lambertiana (cos θ) · ancho a media potencia 120°
Dentro del cono verde cae el 5.88 % de la luz = NA² = 0.0588.
Ø fuente (aˊrea de emisioˊn)µm
Ø nuˊcleo Dµm
NA:
NA = √(n₁² − n₂²) = 0.2425 con n₁ = 1.480, n₂ = 1.460 (se cambian en los simuladores de arriba).
Reflexión:
PfuentemW
UIlossno existe
NAloss-12.306 dB
Rloss-0.166 dB
ICL-12.472 dB
Pacoplada0.1132 mW
Paso a paso (con tus valores)
As=1256.6μm2≤Ac=1963.5μm2⇒UIlossno existe (la fuente cabe entera)NAloss=10log(NA2)=20log(0.2425)=−12.306dBρ=(n1+1n1−1)2=(2.4800.480)2=0.03746⇒Rloss=10log(1−ρ)=−0.166dBICL=UIloss+NAloss+Rloss=(0.00)+(−12.31)+(−0.17)=−12.47dBPacoplada=Pfuente⋅10ICL/10=(2.000mW)(0.05660)=0.1132mW
5. Presupuesto de enlace
Del transmisor al detector la potencia va bajando: la ICL a la entrada, cada conector, cada empalme (mecánico o por fusión), los acopladores, la atenuación de la propia fibra α (en dB/km) y la pérdida de salida (OCL). Además se reserva un margen de seguridad. Si lo que llega es muy poco, la SNR cae y la tasa de errores (BER) sube. En dB todo se suma y se resta:
Absorción UV:αUV=AeB/λ; cae rápido con λ y desde ~1 µm se desprecia.
Impurezas: iones metálicos (Co, V, Ni, Cr, Cu, Fe) y grupos OH⁻ (agua). Hoy existen fibras «secas» sin el pico de agua.
Absorción IR: vibración de los enlaces Si–O del SiO₂, con picos en 9.1, 12.5 y 21 µm; su cola sube después de 1.6 µm.
Curvaturas: macro y microcurvaturas (y cuántas hay).
Dispersión de Rayleigh: fluctuaciones microscópicas del índice congeladas al solidificar el vidrio; αR∝1/λ4.
Defectos geométricos: núcleo no concéntrico o no circular.
La suma deja tres «ventanas» de baja pérdida: 0.85 µm (la primera, barata), 1.31 µm y 1.55 µm (la mínima, ~0.2 dB/km). Los picos de agua están en ≈ 0.95, 1.24 y 1.38 µm.
Simulador 6 · Atenuación de la fibra según λ
Cargar:
Total Rayleigh ∝ 1/λ⁴ Absorción UV Absorción IR Iones OH⁻ (agua)· Toca o arrastra sobre la gráfica para mover λ.
λµm
contaminacioˊn OH (0 = fibra seca)
Ade RayleighdB·µm⁴/km
αRayleigh0.1386 dB/km
αUV0.0292 dB/km
αIR0.0204 dB/km
αOH≈ 0 dB/km
αtotal0.188 dB/km
Paso a paso (con tus valores)
αR=λ4A=(1.550)40.80=0.1386kmdBαUV=AeB/λ=0.0292αIR≈0.0204αOH≈0αtotal=αR+αUV+αIR+αOH=0.188dB/kmen L=1.000km:αL=0.19dB⇒queda el 95.8%de la potencia
Curvas UV e IR con el modelo de Keiser para sílice con 6 % de GeO₂; los picos de OH son ilustrativos (su altura real depende de cuánta agua quedó en el vidrio).
7. Semiconductores y el LED
En un semiconductor los electrones viven en dos bandas separadas por una brecha de bandaWg: la banda de valencia (abajo) y la banda de conducción (arriba). Un conductor no tiene brecha; un semiconductor tiene ≈ 1 eV; un aislante, varios eV. El silicio puro comparte sus 4 electrones de valencia con sus vecinos; al doparlo:
Tipo N: el fósforo (5 e⁻) es donador: cuatro enlazan, el quinto queda libre.
Tipo P: el boro (3 e⁻) es aceptor: falta un enlace y queda un hueco.
Al unir P y N y polarizar en directa (+ al lado P), los electrones de la banda de conducción caen sobre los huecos de la banda de valencia (recombinación) y la energía Wg sale como un fotón. En inversa no circula corriente y no hay luz. La longitud de onda la fija la brecha:
Potencia vs corriente: la potencia óptica crece linealmente con la corriente; la pendiente es la eficienciaη=ΔIΔP[AW].
Patrón de radiación: potencia vs ángulo; su ancho se mide entre los puntos de media potencia.
Tipos: el de superficie emite en un cono ancho (~120°); el de borde concentra el haz (~30°) y acopla mucho mejor a la fibra.
Corriente de modulación: en digital, pulsos de corriente dan pulsos de luz. En analógica hace falta una corriente DC de polarización para que la señal no baje de cero y se recorte:
i(t)=IDC+I^sinωt⇒P(t)=PDC+P^sinωt
Simulador 8 · LED: curva P–I y corriente de modulación
Cargar:
η(pendiente P–I)W/A
IDCmA
I^mA
fMHz
Pmax(1000 = sin lıˊmite)mW
IDC500.0 mA
I^100.0 mA
PDC10.00 mW
P^2.000 mW
Pentre8.000 y 12.00 mW
T=1/f1.000 µs
Paso a paso (con tus valores)
η=ΔIΔP=0.0200AW=20.0AmWPDC=ηIDC=(0.0200)(0.5000)=10.00mWP^=ηI^=2.000mWi(t)=IDC+I^sinωt=500.0+100.0sin(2π⋅1.00×106t)mAP(t)=PDC+P^sinωt=10.00+2.000sin(2π⋅1.00×106t)mWI^≤IDC⇒la corriente nunca baja de cero: sin recorte
La animación va lenta a propósito: el tiempo está en «períodos», no en segundos reales.
Simulador 9 · Patrón de radiación del LED
Superficie · ancho a media potencia 120°Borde · ancho a media potencia 30°
ancho del LED de borde (a media potencia)°
n1
n2
θade la fibra14.03°
acopla el LED de superficie5.88 %
acopla el LED de borde47.07 %
Paso a paso (con tus valores)
superficie (lambertiana): P(θ)=P0cosθ⇒η=1−cos2θa=sin2θa=NA2=(0.2425)2=0.0588borde: P(θ)=P0cos20.0θ⇒η=1−cos21.0θa=0.4707NAloss:10log(0.0588)=−12.31dB (superficie)vs−3.27dB (borde)
Formulario
Velocidad en el medio
v=c/n
Ángulo crítico
sinθc=n2/n1
Ángulo axial máximo
φ=90∘−θc
Apertura numérica
NA=sinθa=n12−n22
Frecuencia normalizada
V=πDNA/λ
Número de modos
N≈V2/2
Monomodo
V<2.405
Retardo modal (SI)
Δt=cLn1(n2n1−1)
Ancho de banda
BW=1/Δt
No intercepción
10log(Ac/As)
Pérdida por NA
10log(NA2)=20logNA
Reflexión
10log(1−ρ),ρ=(n1+1n1−1)2
Potencia en dBm
10log(P/1mW)
Presupuesto
PTx−PRx=∑+M+αL
LED
λ(μm)=1.24/Wg(eV)
Eficiencia del LED
η=ΔP/ΔI
Trampas frecuentes
La pérdida por NA va al cuadrado: 10log(NA2), no 10log(NA). Con NA = 0.2425 son −12.31 dB, no −6.15 dB.
Monomodo es V<2.405, no V=1.
Las pérdidas en dB se suman con su signo: una pérdida nunca «devuelve» potencia.
Si el problema da un ICL, entra en el presupuesto del enlace.
Redondear θc a 80.56° en vez de 80.57° cambia Δt de 67.58 a 67.72 ns: arrastra decimales hasta el final.
El germanio es un semiconductor (≈ 0.66 eV), no un aislante.