Al trabajar con sensores de metasuperficies en COMSOL, surge un interés particular en la aplicación de la Transparencia Inducida Electromagnéticamente (EIT) y los Estados Ligados en el Continuo (BIC) para la detección de índice de refracción. Estos fenómenos ópticos permiten resonancias de alta precisión en metasuperficies, logrando sensibilidades a nanoescala.

La simulación de metasuperficies a menudo implica la definición de la geometría, las propiedades de los materiales y el análisis de resonancias. Por ejemplo, en COMSOL, el proceso podría comenzar definiendo un modelo de dominio de frecuencia:

% Inicializar el modelo
model = createpde("electromagnetic","frequency-domain");

% Definir la geometría (ej. una matriz de nanoestructuras)
% model.geometry.create
% model.selection.all;
% model.geometry.createRectangle('rect1', [-w/2 -h/2], [w, h]);
% ... definir la estructura periódica ...

% Especificar coeficientes del material
epsilon_r = 10; % Permitividad relativa
mu_r = 1;       % Permeabilidad relativa
specifyCoefficients(model,"RelativePermittivity",epsilon_r,"RelativePermurbility",mu_r);

% Aplicar condiciones de contorno (ej. campos eléctricos/magnéticos)
% applyBoundaryCondition(model,"edge", ...);

% Definir el rango de frecuencias a simular
freq_start = 100e12; % 100 THz
freq_end = 200e12;   % 200 THz
freq_vector = linspace(freq_start, freq_end, 100);

% Resolver la PDE
results = solvepde(model, freq_vector);

La clave reside en ajustar los parámetros estructurales y las propiedades del material para excitar resonancias EIT o BIC. Por ejemplo, al diseñar una matriz de nanobarras con un período específico en relación con la longitud de onda de operación, se puede observar una ventana de transmisión, un fenómeno análogo a EIT.

La EIT es ventajosa para la detección de índice de refracción debido a su curva de dispersión empinada. Pequeños cambios en el índice de refracción local (Δn) provocan desplazamientos significativos en el pico de resonancia. A diferencia de las cavidades F-P tradicionales que requieren Δn ≈ 0.1 para la detección, las metasuperficies basadas en EIT pueden detectar cambios de Δn = 0.005, manifestándose como un desplazamiento del pico de transmisión mayor que el ancho de línea.

Los BIC, por otro lado, ofrecen un enfoque distinto. En aplicaciones como la detección de glucosa, cuando la estructura satisface ciertas condiciones geométricas y de materiales, los canales de disipación se cierran debido a interferencia cuántica. Esto resulta en factores de calidad (Q) extremadamente altos, del orden de 105 o más. Durante las simulaciones de barrido paramétrico en COMSOL, los BIC se manifiestan como "acantilados de datos" donde la intensidad del campo aumenta drásticamente.

La optimización para la detección de índice de refracción implica un equilibrio entre alta sensibilidad y límite de detección. Un enfoque práctico es monitorizar la demora del grupo y la transmitancia durante los barridos paramétricos en COMSOL:

% Simulación paramétrica para detectar cambios en el índice de refracción del analito
material_index_values = 1.33:0.001:1.7; % Rango de índices de refracción
delay_times = zeros(size(material_index_values));
transmittances = zeros(size(material_index_values));

for i = 1:length(material_index_values)
    % Actualizar la permitividad relativa del material del analito
    model.param.set('n_analyte', num2str(material_index_values(i)));
    
    % Ejecutar el estudio
    results = model.study('std1').run;
    
    % Extraer la demora del grupo y la transmitancia
    % ... (código para extraer resultados, ej. de datos de campo o de transmisión) ...
    % delay_times(i) = extract_group_delay(results);
    % transmittances(i) = extract_transmittance(results);
end

% Analizar las curvas de demora del grupo y transmitancia para características EIT

La aparición de un patrón de "doble pico y un valle" en la curva de demora del grupo indica la apertura de la ventana EIT. En esta condición, el efecto de acumulación de fase inducido por los cambios en el índice de refracción se amplifica, lo que permite sensibilidades de fase muy altas, del orden de 2π × 104 rad/RIU, superando significativamente la detección basada en intensidad.

La fabricación de estructuras para BIC presenta desafíos debido a su baja tolerancia a errores. Desviaciones mínimas en las dimensiones de la estructura pueden degradar significativamente el factor Q. El análisis de perturbaciones en COMSOL puede revelar que asimetrías estructurales superiores a λ/20 pueden transformar un BIC en un cuasi-BIC. Por lo tanto, es aconsejable incorporar parámetros de rugosidad aleatoria en las simulaciones para prever el impacto de las imperfecciones de fabricación.

En resumen, EIT y BIC son mecanismos complementarios en la detección con metasuperficies. EIT crea resonancias agudas, mientras que BIC ofrece un confinamiento de energía extremo. La implementación de estas técnicas en sensores de índice de refracción puede mejorar drásticamente la sensibilidad y el rendimiento.