Metodo ibrido ANN-Z per la modellazione di superfici intelligenti riconfigurabili basate su nanotubi di carbonio per l'orientamento del fascio in frequenze terahertz

Sommario

Questo protocollo presenta un metodo ibrido basato su reti neurali artificiali e trasformata Z per la modellizzazione elettromagnetica accurata di superfici intelligenti riconfigurabili basate su nanotubi di carbonio a parete singola, operanti nella banda terahertz (0,5-30 THz) per applicazioni di comunicazione wireless 6G, ottenendo un'accelerazione computazionale di 180× con una sintonizzabilità della fase di riflessione superiore a 310° e un intervallo di puntamento del fascio di ±45°.

Abstract

Le superfici intelligenti riconfigurabili otticamente basate su nanotubi di carbonio a parete singola offrono soluzioni promettenti per la deflessione del fascio terahertz e la manipolazione delle onde fotoniche nei futuri sistemi wireless 6G. Tuttavia, la modellizzazione accurata di queste strutture rimane impegnativa a causa degli effetti di trasporto quantistico, dell'induttanza cinetica e del comportamento eccitonico multi-risonante su ampie gamme di frequenza. Questo protocollo descrive un framework computazionale ibrido che integra la modellizzazione della conduttività quantistica mediante la formalizzazione di Kubo, una regressione polinomiale per la regolarizzazione dei dati e un'analisi della funzione di trasferimento nel dominio Z per una caratterizzazione accurata delle celle unitarie delle superfici intelligenti riconfigurabili otticamente a base di nanotubi di carbonio a parete singola. Il metodo inizia con la progettazione di risonatori a nanostriscia incrociati di nanotubi di carbonio a parete singola su un substrato di quarzo con chiralità (10,5) (diametro 0,60 nm, bandgap 1,762 eV), seguita da una simulazione elettromagnetica full-wave in CST Microwave Studio nella banda 0,5-30 THz. Un modello di regressione polinomiale di ordine 8 elabora i parametri S estratti per rimuovere le fluttuazioni numeriche e prevedere risposte elettromagnetiche regolarizzate. Viene quindi adattata una funzione di trasferimento discreta H(z) con ordine del numeratore pari a 6 e ordine del denominatore pari a 7, utilizzando un'ottimizzazione ai minimi quadrati con decomposizione QR, consentendo l'analisi di stabilità tramite poli e zeri e la verifica della passività. Il protocollo include inoltre la regolazione della conduttività quantistica mediante modulazione del potenziale chimico per l'ottimizzazione della deflessione del fascio. Risultati rappresentativi dimostrano una sintonizzabilità della fase di riflessione superiore a 310°, un aumento dell'assorbimento fino al 92,3%, un intervallo di deflessione del fascio di ±45° con livelli dei lobi laterali inferiori a -12 dB e un'accelerazione computazionale di 180× rispetto ai metodi convenzionali di ottimizzazione full-wave. La regressione polinomiale ha ottenuto un errore quadratico medio sul test di 0,0688 con coefficiente R² di 0,994, mentre l'adattamento di H(z) ha raggiunto un errore quadratico medio di 0,89 dB. L'analisi di stabilità ha confermato che tutti i poli sono all'interno del cerchio unitario. Questo protocollo fornisce un percorso efficiente e riproducibile per la progettazione di metasuperfici fotoniche programmabili e di sistemi di comunicazione terahertz intelligenti per la 6G e oltre.

Introduzione

L'evoluzione rapida verso i sistemi wireless di sesta generazione (6G) ha accelerato l'esplorazione delle bande di frequenza terahertz (THz) e ottiche per raggiungere velocità di trasmissione ultra-elevate superiori a 1 Tbps, rilevamento intelligente, formazione del fascio olografica e ingegnerizzazione adattiva del fronte d'onda^1,2,3. La banda terahertz (0,1-30 THz) offre una larghezza di banda abbondante, ma soffre di un forte attenuazione nello spazio libero (circa 20 dB·km^-1 a 1 THz), di assorbimento atmosferico molecolare dovuto al vapore acqueo a 557 GHz, 752 GHz, 988 GHz e 1,13 THz, e di un'elevata sensibilità all'oscuramento causato da particelle atmosferiche e pioggia⁴.

Le superfici intelligenti riconfigurabili si sono affermate come una tecnologia abilitante in grado di manipolare dinamicamente i fronti d'onda elettromagnetici attraverso riflessione programmabile, assorbimento, controllo della fase e focalizzazione delle onde^5,6,7. Quando estese a frequenze ottiche, queste superfici offrono un controllo senza precedenti sulla propagazione della luce, permettendo applicazioni come LiDAR, comunicazioni ottiche nello spazio libero, display olografici e puntamento adattivo del fascio⁸.

Le strutture convenzionali di superfici intelligenti riconfigurabili basate su risonatori metallici incontrano notevoli limitazioni alle frequenze terahertz e ottiche a causa dell'aumento delle perdite ohmiche dovute alla conducibilità di Drude, all'attenuazione plasmonica, ai vincoli di fabbricazione con dimensioni strutturali inferiori a 100 nm e alla limitata sintonizzabilità delle proprietà dielettriche fisse^9,10,11. Le metasuperfici basate su nanomateriali, in particolare quelle che impiegano nanotubi di carbonio a parete singola (SWCNT), offrono alternative promettenti grazie alle loro eccezionali proprietà quantistiche ed elettromagnetiche: dimensioni nanometriche con diametri compresi tra 0,7 e 3 nm, conducibilità quantistica sintonizzabile controllata dal potenziale chimico, mobilità dei portatori straordinaria superiore a 100.000 cm²·(V·s^-1), elevata conducibilità termica di circa 3000 W·(m·K^-1) e forti risonanze eccitoniche nei regimi visibile e vicino infrarosso^12,13,14,15. Gli SWCNT mostrano forti interazioni elettromagnetiche nei regimi terahertz e fotonico grazie al coupling eccitone-fotone, consentendo risonatori ottici compatti con risposte elettromagnetiche programmabili¹⁶. Queste proprietà rendono gli SWCNT candidati interessanti per superfici intelligenti riconfigurabili ottiche, in cui dimensioni ridotte, risposta sintonizzabile e manipolazione efficiente delle onde sono requisiti progettuali importanti.

Nonostante questi vantaggi, una modellizzazione accurata delle superfici intelligenti riconfigurabili basate su nanotubi di carbonio a parete singola (SWCNT) rimane una sfida a causa degli effetti di trasporto quantistico che richiedono approcci basati sulla funzione di Green fuori equilibrio, dell'induttanza cinetica dominante a frequenze superiori a 100 GHz, del comportamento eccitonico multi-risonante con energie di legame fino a 0,4 eV, dei fenomeni di dispersione non lineare derivanti dall'accoppiamento elettrone-fonone, delle proprietà ottiche dipendenti dalla chiralità, che richiedono la specificazione degli indici chirali (n, m), e della conducibilità dipendente dalla temperatura, che richiede una modellizzazione termica autoconsistente^17,18,19,20. Gli approcci tradizionali basati su circuiti equivalenti e i metodi convenzionali di adattamento elettromagnetico, come l'approssimazione razionale e il vector fitting, spesso non riescono a rappresentare con precisione queste interazioni complesse su ampie bande di frequenza che coprono da 0,5 a 30 THz, un intervallo di 60×^21,22,23. Queste limitazioni possono ridurre l'accuratezza della modellizzazione e aumentare la complessità computazionale durante la progettazione di superfici intelligenti riconfigurabili nella banda terahertz, evidenziando la necessità di approcci di modellizzazione efficienti che preservino l'interpretabilità fisica pur catturando con precisione il comportamento elettromagnetico su larga banda.

Questo protocollo presenta un framework ibrido di modellizzazione che combina in modo unico l'analisi della conducibilità quantistica mediante la formalizzazione di Kubo con il calcolo del bandgap dipendente dalla chiralità, la regressione polinomiale per la riduzione del rumore e la regolarizzazione della risposta, l'analisi della funzione di trasferimento tramite trasformata Z per la stabilità poli-zeri e l'interpretazione fisica, e la simulazione elettromagnetica in onda intera di celle unitarie basate su SWCNT. Il flusso di lavoro fornisce una procedura passo-passo per costruire, analizzare e ottimizzare superfici intelligenti ottiche riconfigurabili basate su SWCNT nell'intervallo di frequenze terahertz. Il metodo è pensato per ricercatori e ingegneri che operano su superfici elettromagnetiche avanzate, nanofotonica e sistemi di comunicazione wireless di nuova generazione. Il protocollo consente agli utenti di generare modelli stabili e fisicamente interpretabili di superfici intelligenti ottiche riconfigurabili basate su SWCNT per applicazioni di controllo del fascio e delle onde elettromagnetiche. Si presuppone una conoscenza dei concetti di simulazione elettromagnetica, ma vengono forniti passaggi dettagliati per consentire la replicazione del metodo anche a ricercatori nuovi del settore.

Protocollo

1. Progettazione della cella unitaria ottica RIS per SWCNT

1. Selezione della chiralità del SWCNT
   1. Selezionare la chiralità del SWCNT (10,5) in base all'analisi della conducibilità quantistica.
   2. Calcolare il diametro del nanotubo utilizzando la formula riportata di seguito
      
      dove a_cc=0.142 nm è la lunghezza del legame carbonio-carbonio. La chiralità (10,5) produce un diametro di 0,60 nm e un bandgap di 1,762 eV, ottimali per il funzionamento in banda terahertz.
   3. Calcolare l'angolo chirale utilizzando
2. Definizione della geometria della cella unitaria
   1. Progettare risonatori a nanostriscia di SWCNT incrociati su un substrato di quarzo con i seguenti parametri: permittività del substrato pari a 3,8, spessore del substrato di 500 nm, dimensioni della cella unitaria di 1,5 × 1,5 µm², larghezza della striscia di SWCNT di 50 nm, lunghezza della striscia di SWCNT di 700 nm e distanza tra le strisce di 120 nm. Queste dimensioni garantiscono un funzionamento sublunghezza d'onda nell'intervallo 0,5-30 THz.
3. Configurazione della simulazione in CST Microwave Studio
   1. Avviare CST Microwave Studio e creare un nuovo progetto utilizzando il risolutore nel dominio della frequenza per microonde e radiofrequenze.
   2. Costruire la geometria della cella unitaria utilizzando i parametri definiti. Applicare condizioni al contorno di cella unitaria nelle direzioni x e y con contorni periodici. Definire porte di Floquet lungo la direzione z per l'eccitazione dell'onda piana.
   3. Impostare l'intervallo di scansione della frequenza da 0,5 THz a 30 THz con un passo di 0,05 THz, generando 590 punti frequenziali. Configurare il risolutore nel dominio del tempo con raffinamento adattivo della mesh e impostare l'obiettivo di convergenza a -40 dB.
   4. Eseguire la simulazione elettromagnetica in campo completo. Esportare il coefficiente di riflessione complesso S₁₁, il coefficiente di trasmissione S₂₁, la fase di riflessione e gli spettri di assorbimento come file CSV.

2. Modellizzazione della conducibilità quantistica mediante formalismo di Kubo

1. Inizializzazione delle costanti fisiche
   1. Inizializzare le costanti fisiche: carica elementare e=1,602×10⁻¹⁹C, costante di Planck ridotta ħ=1,0546×10⁻³⁴J·s, costante di Boltzmann k_B=1,3806×10⁻²³J·K^-1, temperatura T=300K, velocità di Fermi v_F=8×10⁵m·s^-1 e impedenza del vuoto η₀=377Ω.
2. Calcolo della conducibilità intra-banda
   1. Calcolare la conducibilità intra-banda per ogni punto di frequenza utilizzando:
      
      dove τ=0,5ps è il tempo di rilassamento dei portatori ed E_F è l'energia di Fermi (variata da 0,1eV a 0,4eV).
3. Calcolo della conducibilità inter-banda
   1. Calcolare la conducibilità inter-banda utilizzando:
      
   2. Calcolare la conducibilità superficiale totale come σ_totale(ω) = σ_intra(ω) + σ_inter(ω).
4. Calcolo della risposta ottica
   1. Calcolare il ritardo di fase del substrato per il quarzo (ε_r=3,8, spessore=500nm) utilizzando φ_sub=2πf(n_sub)d_sub/c, dove n_sub=.
   2. Calcolare il coefficiente di riflessione a partire dalla conducibilità utilizzando
      
   3. Calcolare la riflettanza come R(ω)=|S₁₁(ω)|². Calcolare l'assorbanza come A(ω)=1-R(ω)-T(ω).

3. Generazione e preelaborazione dei dati CST

1. Generazione di dati CST emulati
   1. Generare risultati simulati di onde complete CST emulati aggiungendo risonanze realistiche al S₁₁ quantistico.
   2. Definire cinque risonanze: eccitone E₁₁ a 2,8 THz (ampiezza -15,2 dB, fase 45°), eccitone E₂₂ a 5,6 THz (-22,8 dB, -120°), risonanza plasmonica a 12,4 THz (-8,5 dB, 60°), modo della cavità a 18,9 THz (-12,1 dB, -30°) e risonanza assistita da fonone a 24,7 THz (-6,8 dB, 15°).
   3. Aggiungere rumore numerico con livello di rumore pari a 0,03 (SNR=30,5 dB) per simulare le fluttuazioni numeriche CST.
2. Preprocessing dei dati
   1. Normalizzare l'asse della frequenza su [0,1] utilizzando
      
   2. Estrarre le parti reale e immaginaria di S₁₁ e l'ampiezza di S₂₁.
   3. Suddividere il dataset in training (80%), validazione (10%) e test (10%) mediante permutazione casuale. Normalizzare gli ingressi a media zero e varianza unitaria.

4. Regressione polinomiale per la levigatura dei dati

1. Applicare un polinomio di 8^° ordine alla parte reale di S₁₁ utilizzando polyfit: P_real = polyfit(X_train, Y_train(:,1), 8). Applicare un polinomio di 8^° ordine alla parte immaginaria: P_imag = polyfit(X_train, Y_train(:,2), 8). Applicare un polinomio di 8^° ordine al modulo di S₂₁: P_mag = polyfit(X_train, Y_train(:,3), 8). Generare previsioni lisce utilizzando polyval sull'intero intervallo di frequenze. Calcolare S₁₁ liscia come S₁₁_ann = S11_real_ann + i × S11_imag_ann.

5. Adattamento della Funzione di Trasferimento nel Dominio Z

1. Impostare la frequenza di campionamento F_s=60THz (2× la frequenza massima secondo il criterio di Nyquist). Mappare le frequenze nel dominio Z utilizzando z=e^{i2πf/Fs}.
2. Definire la funzione di trasferimento discreta nel dominio Z, espressa come:
   
   dove l'ordine del numeratore è n=6 e l'ordine del denominatore è m=7.
3. Per l'ottimizzazione ai minimi quadrati, impostare le equazioni lineari A × x = B, dove A contiene i termini del numeratore e del denominatore e B contiene i valori di S₁₁ predetti dalla rete neurale (ANN). Risolvere utilizzando la decomposizione QR per garantire stabilità numerica:
   
4. Analisi di poli, zeri e stabilità
   1. Estrarre i poli risolvendo il polinomio del denominatore. Assicurare la stabilità proiettando i poli con |polo| ≥ 1 all'interno del cerchio unitario utilizzando poli(p) = poli(p) / (|poli(p)| + 0,1).
   2. Ricostruire il denominatore a partire dai poli stabilizzati. Estrarre gli zeri risolvendo il polinomio del numeratore. Valutare H(z) su una griglia di frequenze e calcolare l'errore quadratico medio.

6. Ottimizzazione della Deviazione del Fascio

1. Il fattore di antenna normalizzato per un array lineare a scansione elettronica è dato da:
   ​
   dove I_n = 1 per eccitazione uniforme, k = 2π/λ, d = λ/2 (distanza tra elementi), β_n è lo sfasamento progressivo. Per angoli di puntamento di ±45°, ±30° e 0°, calcolare lo sfasamento tra elementi adiacenti come .EQUATION
   1. Calcolare i diagrammi del campo lontano e valutare i livelli dei lobi laterali e l'ampiezza del fascio a mezza potenza.

Risultati

Selezione della chiralità dei nanotubi di carbonio singolo parete (SWCNT)
Il protocollo descritto è stato implementato per una cella unitaria RIS ottica basata su SWCNT (10,5) operante nella banda 0,5-30 THz. I risultati rappresentativi dimostrano l'efficacia dell'approccio modellistico ibrido polinomiale-Z per la caratterizzazione elettromagnetica accurata e l'ottimizzazione della deviazione del fascio.

Analisi della Conduttività Quantistica
Il formal...

Discussione

Le fasi critiche del protocollo richiedono un'attenzione particolare per garantire un'implementazione corretta. In primo luogo, la selezione accurata degli indici chirali dei nanotubi di carbonio a parete singola (SWCNT) è essenziale poiché il bandgap e la risposta ottica dipendono fortemente dalla chiralità. La chiralità (10,5) indicata in questo protocollo fornisce un bandgap ottimale di 1,762 eV per il funzionamento in banda terahertz, ma gli utenti che mirano a bande di frequenza diverse dovrebbero calcolare la ...

Materiali

Elenco dei materiali utilizzati in questo articolo

Nome	Azienda	Numero di catalogo	Commenti
CST Microwave Studio	Dassault Systèmes	N/A	Versione 2024, risolutore nel dominio della frequenza
MATLAB	MathWorks	N/A	Versione R2014a o successiva
Substrato in quarzo	University Wafer	4526	Spessore di 500 nm, ε_r = 3.8
Chiralità SWCNT (10,5)	NanoIntegris	SWCNT-105	Diametro di 0,60 nm, >90% semiconduttivo
Computer personale	N/A	N/A	32 GB di RAM, minimo 8 core CPU