Método Híbrido ANN-Z para Modelagem de Superfícies Inteligentes Reconfiguráveis Baseadas em Nanotubos de Carbono para Direcionamento de Feixe em Terahertz

Resumo

Este protocolo apresenta um método híbrido de rede neural artificial e transformada Z para modelagem eletromagnética precisa de superfícies inteligentes reconfiguráveis baseadas em nanotubos de carbono de parede simples, operando na faixa de terahertz (0,5-30 THz) para aplicações de comunicação sem fio 6G, alcançando aceleração computacional de 180× com sintonização da fase de reflexão superior a 310° e alcance de direcionamento de feixe de ±45°.

Resumo

Superfícies inteligentes reconfiguráveis por meio de luz baseadas em nanotubos de carbono de parede simples oferecem soluções promissoras para direcionamento de feixes em terahertz e manipulação de ondas fotônicas em futuros sistemas sem fio 6G. No entanto, a modelagem precisa dessas estruturas permanece desafiadora devido a efeitos de transporte quântico, indutância cinética e comportamento excitônico multiressonante em amplas faixas de frequência. Este protocolo descreve uma estrutura computacional híbrida que integra a modelagem de condutividade quântica utilizando a formalização de Kubo, suavização de dados baseada em regressão polinomial e análise de função de transferência no domínio Z para caracterização precisa de células unitárias de superfícies inteligentes reconfiguráveis por luz baseadas em nanotubos de carbono de parede simples. O método inicia com o projeto de ressoadores em nano-fitas de nanotubos de carbono de parede simples entrelaçados sobre um substrato de quartzo com quiralidade (10,5) (diâmetro 0,60 nm, banda proibida 1,762 eV), seguido por simulação eletromagnética de onda completa no CST Microwave Studio na faixa de 0,5 a 30 THz. Um modelo de regressão polinomial de ordem 8 processa os parâmetros S extraídos para remover flutuações numéricas e prever respostas eletromagnéticas suavizadas. Uma função de transferência discreta H(z) com ordem do numerador 6 e ordem do denominador 7 é então ajustada usando otimização por mínimos quadrados com decomposição QR, permitindo análise de estabilidade por polos e zeros e verificação de passividade. O protocolo incorpora ainda o ajuste da condutividade quântica por meio da modulação do potencial químico para otimização do direcionamento de feixe. Resultados representativos demonstram sintonização da fase de reflexão superior a 310°, aumento da absorção até 92,3%, alcance de direcionamento de feixe de ±45° com níveis de lóbulo secundário abaixo de -12 dB e aceleração computacional de 180× em comparação com métodos convencionais de otimização por onda completa. A regressão polinomial alcançou erro quadrático médio de teste de 0,0688 com coeficiente R² de 0,994, enquanto o ajuste de H(z) alcançou erro quadrático médio de 0,89 dB. A análise de estabilidade confirmou que todos os polos estão dentro do círculo unitário. Este protocolo fornece um caminho eficiente e reprodutível para o projeto de metassuperfícies fotônicas programáveis e sistemas inteligentes de comunicação em terahertz para 6G e além.

Introdução

A rápida evolução rumo aos sistemas sem fio de sexta geração (6G) acelerou a exploração das faixas de frequência terahertz (THz) e óptica para alcançar taxas de dados ultraelevadas superiores a 1 Tbps, sensores inteligentes, formação de feixe holográfica e engenharia adaptativa de frente de onda^1,2,3. A faixa terahertz (0,1–30 THz) oferece largura de banda abundante, mas sofre com severa atenuação no espaço livre (aproximadamente 20 dB·km^-1 em 1 THz), absorção atmosférica molecular devido ao vapor d'água em 557 GHz, 752 GHz, 988 GHz e 1,13 THz, além de extrema sensibilidade ao bloqueio por partículas atmosféricas e chuva⁴.

As superfícies inteligentes reconfiguráveis surgiram como uma tecnologia habilitadora capaz de manipular dinamicamente frentes de onda eletromagnéticas por meio de reflexão programável, absorção, controle de fase e focagem de ondas^5,6,7. Quando estendidas a frequências ópticas, essas superfícies oferecem um controle sem precedentes sobre a propagação da luz, permitindo aplicações como LiDAR, comunicações ópticas em espaço livre, displays holográficos e direcionamento adaptativo de feixes⁸.

Estruturas convencionais de superfícies inteligentes reconfiguráveis baseadas em ressoadores metálicos enfrentam limitações significativas nas frequências terahertz e ópticas devido ao aumento das perdas ôhmicas provenientes da condutividade de Drude, amortecimento plasmônico, restrições de fabricação com tamanhos de características abaixo de 100 nm e sintonização limitada de propriedades dielétricas fixas^9,10,11. Metassuperfícies baseadas em nanomateriais, particularmente aquelas que empregam nanotubos de carbono de parede simples (SWCNTs), oferecem alternativas promissoras devido às suas excepcionais propriedades eletromagnéticas quânticas: dimensões em escala nanométrica com diâmetros de 0,7-3 nm, condutividade quântica sintonizável controlada pelo potencial químico, mobilidade extraordinária de portadores superior a 100.000 cm²·(V·s^-1), alta condutividade térmica de aproximadamente 3000 W·(m·K^-1) e fortes ressonâncias excitônicas nas regiões visível e infravermelha próxima^12,13,14,15. Os SWCNTs exibem fortes interações eletromagnéticas nas faixas terahertz e fotônica por meio do acoplamento excitão-fóton, permitindo ressoadores ópticos compactos com respostas eletromagnéticas programáveis¹⁶. Essas propriedades tornam os SWCNTs candidatos atrativos para superfícies inteligentes reconfiguráveis ópticas, onde tamanho compacto, resposta sintonizável e manipulação eficiente de ondas são requisitos importantes de projeto.

Apesar dessas vantagens, a modelagem precisa de superfícies inteligentes reconfiguráveis baseadas em SWCNT permanece desafiadora devido aos efeitos de transporte quântico que exigem abordagens baseadas na função de Green fora do equilíbrio, à indutância cinética dominante em frequências acima de 100 GHz, ao comportamento excitônico multirressonante com energias de ligação de até 0,4 eV, aos fenômenos de dispersão não linear decorrentes do acoplamento elétron-fônon, às propriedades ópticas dependentes da quiralidade, que exigem a especificação dos índices quirais (n, m), e à condutividade dependente da temperatura, que exige modelagem térmica autoconsistente^17,18,19,20. Abordagens tradicionais baseadas em circuitos equivalentes e métodos convencionais de ajuste eletromagnético, como aproximação racional e ajuste vetorial, frequentemente falham em representar com precisão essas interações complexas em largas faixas de frequência que abrangem 0,5–30 THz, uma largura de banda de 60×^21,22,23. Essas limitações podem reduzir a precisão da modelagem e aumentar a complexidade computacional durante o projeto de superfícies inteligentes reconfiguráveis em terahertz, destacando a necessidade de abordagens eficientes de modelagem que preservem a interpretabilidade física enquanto capturam com precisão o comportamento eletromagnético em larga banda.

Este protocolo apresenta uma estrutura híbrida de modelagem que combina de forma única a análise de condutividade quântica utilizando a formalização de Kubo com o cálculo de bandgap dependente da quiralidade, regressão polinomial para redução de ruído e suavização da resposta, análise da função de transferência por transformada Z para estabilidade de polos e zeros e interpretabilidade física, e simulação eletromagnética de onda completa de células unitárias baseadas em SWCNT. O fluxo de trabalho fornece um procedimento passo a passo para construir, analisar e otimizar superfícies inteligentes reconfiguráveis ​​ópticas baseadas em SWCNT na faixa de frequência terahertz. O método destina-se a pesquisadores e engenheiros que trabalham com superfícies eletromagnéticas avançadas, nanofotônica e sistemas de comunicação sem fio de próxima geração. O protocolo permite aos usuários gerar modelos estáveis e fisicamente interpretáveis ​​de superfícies inteligentes reconfiguráveis ​​ópticas baseadas em SWCNT para aplicações de direcionamento de feixe e controle de ondas eletromagnéticas. O protocolo pressupõe familiaridade com conceitos de simulação eletromagnética, mas fornece etapas detalhadas para replicação por pesquisadores novos na área.

Protocolo

1. Projeto da Célula Unitária Óptica RIS de SWCNT

1. Seleção da Quiralidade do SWCNT
   1. Selecione a quiralidade do SWCNT (10,5) com base na análise de condutividade quântica.
   2. Calcule o diâmetro do nanotubo usando a fórmula mencionada abaixo
      
      em que a_cc=0,142 nm é o comprimento da ligação carbono-carbono. A quiralidade (10,5) resulta em um diâmetro de 0,60 nm e uma banda proibida de 1,762 eV, ideal para operação na faixa de terahertz.
   3. Calcule o ângulo quiral usando
2. Definição da Geometria da Célula Unitária
   1. Projete ressoadores em nano-fitas de SWCNT cruzadas sobre um substrato de quartzo com os seguintes parâmetros: permissividade do substrato de 3,8, espessura do substrato de 500 nm, dimensões da célula unitária de 1,5 × 1,5 µm², largura da fita de SWCNT de 50 nm, comprimento da fita de SWCNT de 700 nm e espaçamento entre fitas de 120 nm. Essas dimensões garantem operação em escala sub-ondulatória na faixa de 0,5-30 THz.
3. Configuração da Simulação no CST Microwave Studio
   1. Inicie o CST Microwave Studio e crie um novo projeto utilizando o resolvedor no domínio da frequência para micro-ondas e RF.
   2. Construa a geometria da célula unitária usando os parâmetros definidos. Aplique condições de contorno de célula unitária nas direções x e y com fronteiras periódicas. Defina portas de Floquet ao longo da direção z para excitação por onda plana.
   3. Defina a faixa de varredura de frequência de 0,5 THz a 30 THz com um passo de 0,05 THz, gerando 590 pontos de frequência. Configure o resolvedor no domínio do tempo com refinamento adaptativo da malha e defina o alvo de convergência como -40 dB.
   4. Execute a simulação eletromagnética de onda completa. Exporte o coeficiente complexo de reflexão S₁₁, o coeficiente de transmissão S₂₁, a fase da reflexão e os espectros de absorção como arquivos CSV.

2. Modelagem da Condutividade Quântica Utilizando o Formalismo de Kubo

1. Inicialização das Constantes Físicas
   1. Inicialize as constantes físicas: carga elementar e=1.602×10⁻¹⁹C, constante de Planck reduzida ħ=1.0546×10⁻³⁴J·s, constante de Boltzmann k_B=1.3806×10⁻²³J·K^-1, temperatura T=300K, velocidade de Fermi v_F=8×10⁵m·s^-1 e impedância do espaço livre η₀=377Ω.
2. Cálculo da Condutividade Intra-banda
   1. Calcule a condutividade intra-banda para cada ponto de frequência usando:
      
      em que τ=0,5ps é o tempo de relaxação dos portadores e E_F é a energia de Fermi (variada de 0,1eV a 0,4eV).
3. Cálculo da Condutividade Inter-banda
   1. Calcule a condutividade inter-banda usando:
      
   2. Calcule a condutividade superficial total como σ_total(ω) = σ_intra(ω) + σ_inter(ω).
4. Cálculo da Resposta Óptica
   1. Calcule o atraso de fase do substrato para quartzo (ε_r=3,8, espessura=500nm) usando φ_sub=2πf(n_sub)d_sub/c, onde n_sub=.
   2. Calcule o coeficiente de reflexão a partir da condutividade usando
      
   3. Calcule a refletância como R(ω)=|S₁₁(ω)|². Calcule a absorção como A(ω)=1-R(ω)-T(ω).

3. Geração e Pré-processamento de Dados CST

1. Geração de Dados CST Emulados
   1. Gere resultados simulados de onda completa CST emulados adicionando ressonâncias realistas ao S₁₁ quântico.
   2. Defina cinco ressonâncias: exciton E₁₁ em 2,8 THz (magnitude -15,2 dB, fase 45°), exciton E₂₂ em 5,6 THz (-22,8 dB, -120°), ressonância plasmônica em 12,4 THz (-8,5 dB, 60°), modo de cavidade em 18,9 THz (-12,1 dB, -30°) e ressonância assistida por fônon em 24,7 THz (-6,8 dB, 15°).
   3. Adicione ruído numérico com nível de ruído 0,03 (SNR=30,5 dB) para simular flutuações numéricas do CST.
2. Pré-processamento dos Dados
   1. Normalize o eixo de frequência para [0,1] usando
      
   2. Extraia as partes real e imaginária de S₁₁ e a magnitude de S₂₁.
   3. Divida o conjunto de dados em treinamento (80%), validação (10%) e teste (10%) usando permutação aleatória. Normalize as entradas para média zero e variância unitária.

4. Regressão Polinomial para Suavização de Dados

1. Adapte um 8^th ajuste um polinômio de ordem 8 à parte real de S₁₁ usando polyfit: P_real = polyfit(X_train, Y_train(:,1), 8).^th polinômio de ordem à parte imaginária: P_imag = polyfit(X_train, Y_train(:,2), 8). Ajuste um polinômio de 8º grau^th polinômio de ordem para a magnitude de S₂₁: P_mag = polyfit(X_train, Y_train(:,3), 8). Gerar previsões suavizadas utilizando polyval em toda a faixa de frequência. Calcular S₁₁ suavizado como S₁₁_ann = S11_real_ann + i × S11_imag_ann.

5. Ajuste da Função de Transferência do Domínio Z

1. Defina a frequência de amostragem F_s=60THz (2×frequência máxima segundo o critério de Nyquist). Mapeie as frequências para o domínio Z usando z=e^{i2πf/Fs}.
2. Defina a função de transferência discreta no domínio Z, expressa como:
   
   em que a ordem do numerador é n=6 e a ordem do denominador é m=7.
3. Para a Otimização por Mínimos Quadrados, monte as equações lineares A × x = B, onde A contém os termos do numerador e do denominador e B contém S₁₁ previsto pela RNA. Resolva utilizando decomposição QR para garantir estabilidade numérica:
   
4. Análise de Polos e Zeros e Estabilidade
   1. Extraia os polos resolvendo o polinômio do denominador. Garanta a estabilidade projetando polos com |polo| ≥ 1 para dentro do círculo unitário usando polos(p) = polos(p) / (|polos(p)| + 0,1).
   2. Reconstrua o denominador a partir dos polos estabilizados. Extraia os zeros resolvendo o polinômio do numerador. Avalie H(z) em uma grade de frequências e calcule o erro quadrático médio.

6. Otimização da Direção do Feixe

1. O fator de arranjo normalizado para um arranjo linear com fase controlada é dado por:
   ​
   em que I_n = 1 para excitação uniforme, k = 2π/λ, d = λ/2 de espaçamento, β_n é o deslocamento de fase progressivo. Para ângulos de direcionamento alvo de ±45°, ±30° e 0°, calcule o deslocamento de fase entre elementos adjacentes como .EQUATION
   1. Calcule os padrões de campo distante e avalie os níveis dos lobos laterais e a largura do feixe a meia potência.

Resultados

Seleção da Quiralidade do SWCNT
O protocolo descrito foi implementado para uma célula unitária óptica de SWCNT (10,5) operando na faixa de 0,5-30 THz. Os resultados representativos demonstram a eficácia da abordagem híbrida de modelagem polinomial-Z para a caracterização eletromagnética precisa e a otimização do direcionamento de feixe.

Análise da Condutividade Quântica
O formalismo de Kubo revelou que o SWCNT (10,5) exibe uma condutividade superficial...

Discussão

Etapas críticas do protocolo exigem atenção cuidadosa para garantir a implementação bem-sucedida. Primeiro, a seleção precisa dos índices quiral de SWCNT é essencial, pois a banda proibida e a resposta óptica são altamente dependentes da quiralidade. A quiralidade (10,5) especificada neste protocolo fornece uma banda proibida ótima de 1,762 eV para operação em terahertz, mas usuários que visam bandas de frequência diferentes devem calcular a quiralidade correspondente usando a fórmula da banda proibida E_g = 2ħv_F/d...

Materiais

Lista de materiais utilizados neste artigo

Nome	Empresa	Número de catálogo	Comentários
CST Microwave Studio	Dassault Systèmes	N/A	Versão 2024, resolvedor no domínio da frequência
MATLAB	MathWorks	N/A	Versão R2014a ou posterior
Substrato de quartzo	University Wafer	4526	Espessura de 500 nm, ε_r = 3.8
Quiralidade SWCNT (10,5)	NanoIntegris	SWCNT-105	Diâmetro de 0,60 nm, >90% semicondutor
Computador pessoal	N/A	N/A	Mínimo de 32 GB de RAM, 8 núcleos de CPU