Óptica numérica · Hoja de problemas

La hoja de problemas · Óptica numérica

Seis problemas: presupuestar una FDTD antes de lanzarla, multiplicar de verdad una matriz característica, anular una reflexión con dos capas y contar lo que gana un objetivo de nueve lentes, elegir entre DDA y BEM para un nanocubo de plata, extrapolar una resonancia a rejilla cero y fabricar el índice 1,23 que no existe en ningún catálogo.

Dificultad ●○○–●●●, pistas sin soluciones cerradas. La brújula de este curso es corta y cara. Primero, y es el error que más simulaciones ha echado a perder: la longitud de onda que fija la rejilla es la del medio, no la del vacío. A 1550 nm dentro de silicio (nn = 3,48) la onda mide 445 nm, así que la regla de «10 a 20 celdas por λ» pide un Δx\Delta x 3,48 veces menor de lo que parece — y en 3D eso son 3,48⁴ ≈ 150 veces más coste. Un cálculo con la λ del vacío no avisa: converge, sale bonito y está mal. Segundo: la condición CFL la fija la velocidad máxima del dominio, que es la del vacío, no la del material más lento. Dividir cc por nn da un Δt\Delta t demasiado grande y la simulación explota en el aire, no en el silicio. Tercero, y es el más sutil: convergencia no es exactitud. Refinar la malla te lleva a la respuesta de tu modelo, con su geometría escalonada, su material tabulado y su dominio truncado; que dos rejillas coincidan no dice nada sobre si el modelo es el problema. Por eso todo lo de este curso se valida contra Mie o contra la Transfer Matrix antes de creérselo.

Necesitas: el curso 02 completo. Los problemas 2 y 3 usan la matriz característica del artículo 02 —y hay que multiplicarla, no mirarla—; el 1 y el 5 la condición CFL y el coste de FDTD del artículo 04; el 4 combina el criterio de discretización de DDA (artículo 05) con el recuento de incógnitas de BEM (artículo 07); el 5 invoca además la convergencia de DGTD del artículo 08; el 6 encadena el medio efectivo del artículo 03 con el antirreflejante del 02 y con la ecuación de la rejilla del curso 01, artículo 02.
Problema 1 ●○○ · Presupuestar una simulación antes de lanzarla

Quieres simular con FDTD un tramo de guía de silicio (nn = 3,48) sobre óxido, a 1550 nm, en un dominio 2D de 6 µm × 4 µm, dejando correr el campo 2 ps. (a) ¿Qué Δx\Delta x te pide la regla de las 20 celdas por longitud de onda, y cuánto te habría salido si tomas la λ del vacío? (b) Con Δx\Delta x = 20 nm, calcula el número de celdas, el paso temporal por CFL, los pasos necesarios y las actualizaciones de celda totales. ¿Cuánta memoria come? (c) Repite el presupuesto en 3D con un dominio de 6 µm × 4 µm × 4 µm. ¿Por cuánto se multiplica el coste, y por qué el factor no es simplemente el número de celdas nuevas?

Pista

Primero la λ\lambda dentro del material, λ0/n\lambda_0/n, y de ahí Δx\Delta x; después celdas por eje, Δt\Delta t por CFL y pasos.

Ojo: la CFL lleva la velocidad máxima del dominio, la del vacío, no c/nc/n; con c/nc/n el paso sale nn veces mayor y revienta en el aire.

Comprueba el paso: la CFL es más restrictiva en 3D, así que si Δt\Delta t no baja al añadir el tercer eje, has puesto 2\sqrt2 donde iba 3\sqrt3.

Problema 2 ●○○ · Multiplicar la matriz que todo el mundo mira y nadie usa

El artículo 02 dedujo la matriz característica de una capa y la fórmula del coeficiente de reflexión del apilamiento, y luego se contentó con dejar que el simulador las evaluara. Hazlo a mano. (a) Escribe MM para una capa de espesor óptico λ0/4\lambda_0/4 —es decir δ=π/2\delta = \pi/2— y sustitúyela en la fórmula de rr para demostrar que R=[(nf2n0ns)/(nf2+n0ns)]2R = \left[(n_f^2 - n_0 n_s)/(n_f^2 + n_0 n_s)\right]^2. (b) Evalúala para MgF₂ (nfn_f = 1,38) sobre vidrio (nsn_s = 1,52) y compárala con el 1,3 % que da el ejercicio 1 de aquel artículo. (c) ¿Qué dice la fórmula sobre el índice ideal, y cuánto vale RR sin recubrir?

Pista

Con δ=π/2\delta = \pi/2 la matriz se queda antidiagonal; mete sus cuatro elementos en la fórmula de rr y multiplica arriba y abajo por nfn_f.

Ojo: el espesor óptico es nfd=λ0/4n_f d = \lambda_0/4, no d=λ0/4d = \lambda_0/4; fuera de la longitud de diseño la matriz recupera su diagonal.

Contrasta (b) con el ejercicio 1 del artículo 02: fórmula cerrada y simulador tienen que dar el mismo número.

Problema 3 ●●○ · Dos capas, cero reflexión, y lo que se juega un objetivo

Este problema retoma el antirreflejante del problema anterior y el objetivo fotográfico que el artículo 02 menciona de pasada. (a) Multiplica las matrices de dos capas de cuarto de onda, de índices n1n_1 (la exterior) y n2n_2, y demuestra que la reflexión se anula exactamente cuando n2/n1=ns/n0n_2/n_1 = \sqrt{n_s/n_0}. (b) Con n1n_1 = 1,38 (MgF₂), ¿qué material haría falta debajo? ¿Cuánto vale RR con el par 1,38/1,70? (c) Un objetivo de nueve lentes tiene 18 superficies vidrio-aire. Calcula la transmisión total sin recubrir, con monocapa de MgF₂ y con una multicapa moderna de 0,2 % por cara. (d) ¿Adónde va la luz que no se transmite, y por qué eso importa más que la pérdida?

Pista

El producto de dos matrices antidiagonales es diagonal, así que MM sale en una línea; para (c), eleva la transmisión de una cara a la 18.

Ojo: la transmisión es (1R)18(1-R)^{18}, no 118R1-18R; con seis superficies más, la segunda cuenta llega a dar negativa.

Comprueba (a) en el límite: con n1=n0n_1 = n_0 la condición tiene que degenerar en la de una sola capa del problema 2.

Problema 4 ●●○ · Un nanocubo de plata: dónde falla el criterio que te han dado

Hay que calcular el espectro de extinción y el campo cercano de un nanocubo de plata de 60 nm de arista a 400 nm, donde su índice complejo es mm = 0,17 + 1,95ii. (a) Aplica el criterio mkd<1|m|kd < 1 del artículo 05: ¿qué espaciado máximo entre dipolos permite, y cuántos dipolos por arista son? (b) Mira el número que te ha salido y decide si el criterio basta. ¿Qué otra cosa fija en realidad la discretización de este objeto? (c) Con dd = 2 nm, cuenta las incógnitas y la memoria de la matriz densa de DDA. (d) Haz el mismo recuento con BEM, mallando la superficie con triángulos de 3 nm de lado. ¿Cuál eliges y por qué?

Pista

Aplica el criterio con m|m| y k=2π/λk = 2\pi/\lambda, mira cuántos dipolos por arista deja y decide si un cubo cabe ahí; luego cuenta incógnitas y memoria.

Ojo: m|m| es el módulo del índice complejo, no su parte real; con 0,17 el espaciado sale once veces mayor.

Comprueba BEM contra el ejercicio 1 del artículo 07: si no cuadra, revisa el área del triángulo equilátero, 3l2/4\sqrt3\,l^2/4 y no l2/2l^2/2.

Problema 5 ●●● · Tres rejillas y la respuesta que ninguna da

Se simula con FDTD la resonancia de un microdisco de silicio y, refinando la rejilla, se obtiene λres\lambda_{\text{res}} = 1548,2 nm con Δx\Delta x = 20 nm, 1551,0 nm con 10 nm y 1551,7 nm con 5 nm. (a) Deduce del propio conjunto de datos el orden de convergencia observado. ¿Coincide con el que promete FDTD? (b) Usa la extrapolación de Richardson para estimar el resultado a rejilla cero, y di cuánto error arrastra la corrida más fina. (c) ¿Qué Δx\Delta x haría falta para bajar el error a 0,05 nm por fuerza bruta, y cuánto costaría en 3D? ¿Y con DGTD de orden pp = 3? (d) ¿Qué habría que sospechar si el cociente de diferencias sucesivas no hubiera salido tan limpio?

Pista

Al dividir hh por dos, las diferencias sucesivas van en razón 2p2^p: de ahí sale pp medido, no supuesto, y con él Richardson da el límite a rejilla cero.

Ojo: DGTD de orden pp converge como hp+1h^{p+1}, no como hph^p; confundirlo dobla el refinamiento que le pides.

Comprueba el signo: el valor extrapolado tiene que caer fuera de las tres corridas; si cae entre dos, has sumado la corrección al revés.

Problema 6 ●●● · Fabricar el índice que no existe: la córnea de la polilla

El problema 2 dejó pendiente algo incómodo: el antirreflejante perfecto aire-vidrio pide nfn_f = 1,233 y no hay ningún material sólido con ese índice. La naturaleza lo resolvió hace cien millones de años: el ojo de las polillas nocturnas está cubierto de conos submicrométricos. (a) Para que el conjunto se comporte como un medio homogéneo y no como una rejilla de difracción, ¿qué periodo máximo Λ\Lambda admite si el sustrato tiene nsn_s = 1,52 y hay que cubrir hasta 400 nm? (b) Modela un cono como cinco láminas de fracción de volumen de vidrio decreciente y calcula el índice efectivo de cada una con Bruggeman. ¿Qué índice sale a media altura, y con qué se parece? (c) ¿Qué altura mínima necesita el cono, y por qué? (d) Un fabricante escala el molde a Λ\Lambda = 280 nm para abaratar la litografía. ¿Qué le pasa a su antirreflejante?

Pista

Lee la ecuación de la rejilla al revés, como condición de no difracción: ningún orden salvo el cero puede propagarse. En (b), Bruggeman lámina a lámina.

Ojo: el medio que fija Λ\Lambda es el sustrato, no el aire: la condición lleva λ/ns\lambda/n_s.

Comprueba Bruggeman en los extremos: con fracción de vidrio 1 tiene que devolver nsn_s exacto, y con 0, exactamente 1.

Fin del curso 02. Con esto tienes lo que hace falta para no empezar una simulación a ciegas: sabes qué método discretiza qué, qué criterio fija la resolución de cada uno —y, más importante, cuándo ese criterio no es el que manda—, cuánta memoria y cuánto tiempo va a costar antes de lanzarlo, y cómo saber si el número que sale ha convergido a algo o sólo ha dejado de moverse. Sabes multiplicar una matriz característica, presupuestar una FDTD en 2D y en 3D, elegir entre volumen y superficie con un recuento de incógnitas, extrapolar tres rejillas a rejilla cero y reconocer el momento en que un medio efectivo se convierte en una rejilla de difracción.

Lo que queda fuera conviene decirlo: aquí no hay optimización inversa —el diseño fotónico moderno se hace con gradientes adjuntos, no probando geometrías—, ni métodos espectrales para estructuras periódicas (RCWA, el caballo de batalla de la litografía), ni nada sobre la incertidumbre de los datos de material, que en plasmónica suele dominar el error total por encima de cualquier refinamiento de malla.

Y hay una deuda concreta que el curso siguiente cobra. Aquí han aparecido tres veces la polarizabilidad de una esfera, el factor de Clausius-Mossotti y un denominador que puede anularse, y las tres veces se ha pasado de largo diciendo «eso es una resonancia plasmónica». El curso 03, Plasmónica, es lo que hay dentro de ese denominador: de dónde sale la permitividad negativa de un metal, por qué la resonancia concentra el campo en volúmenes cien veces menores que la longitud de onda, cuánto vale esa concentración y qué se puede hacer con ella —detectar una sola molécula, calentar un tumor, medir un nanómetro con color—. Todo lo que aprendiste aquí sobre BEM y campos cercanos era la herramienta; allí está el problema para el que se inventó.