Dans le domaine de l'optique, les métasurfaces exploitent des structures sub-longueur d'onde pour manipuler la lumière de manière avancée. Les métasurfaces chérales se distinguent par leur capacité à répondre de façon asymétrique à la lumière polarisée circulairement de différentes hélicités. Ce document décrit la méthodologie de simulation de la réponse chérale d'une métasurface à nanopores en utilisant le logiciel Comsol Multiphysics.
- Principes des métasurfaces à nanopores et de la réponse chérale
La chiralité désigne une propriété géométrique où un objet et son image miroir ne sont pas superposables, analogue à la main gauche et droite. Dans les métasurfaces à nanopores, la conception soignée des paramètres tels que la forme, l'agencement et les dimensions des pores induit une réponse optique chérale intense. Lorsqu'une lumière polarisée circulairement interagit avec une telle structure, les composantes droite et gauche subissent des interactions distinctes, entraînant des phénomènes comme la dichroïsme circulaire.
- Configuration de la simulation dans Comsol
2.1 Construction du modèle
Lancer Comsol Multiphysics et sélectionner l'interface « Ondes électromagnétiques, fréquence » dans le module RF, car la lumière est une onde électromagnétique et l'analyse fréquentielle permet d'évaluer la réponse chérale. Construire une géométrie périodique représentant la métasurface, par exemple un substrat carré avec des nanopores disposés en motif hélicoïdal pour amplifier la chiralité.
// Illustration conceptuelle de la création du modèle
// Initialisation d'un modèle d'onde électromagnétique en 2D
model_initial = setupModel('dimension', 2, 'physics', 'emwave', 'study', 'frequency');
2.2 Définition des propriétés des matériaux
Pour la métasurface, des métaux comme l'or ou des diélectriques sont couramment utilisés. En prenant l'exemple de l'or, sa permittivité diélectrique est modélisée par l'équation de Drude-Lorentz dans Comsol :
// Expression de la permittivité selon le modèle de Drude-Lorentz
permittivity = eps_inf - (plasma_freq^2) / (freq * (freq + 1i * damping_coeff));
Ici, eps_inf représente la permittivité relative à haute fréquence, plasma_freq la fréquence plasmon, freq la fréquence angulaire de la lumière incidente, et damping_coeff le coefficient d'amortissement. Cette formulation assure une modélisation précise du comportement optique de l'or sur une plage de fréquences.
2.3 Conditions aux limites et excitation
Pour simuler une métasurface périodique infinie, appliquer des conditions périodiques sur les bords opposés dans Comsol, assurant ainsi la continuité latérale de la structure. L'excitation est définie comme une lumière polarisée circulairement en réglant la différence de phase entre les composantes x et y du champ électrique à ±π/2, correspondant aux polarisations droite ou gauche.
// Définition du champ électrique pour une polarisation circulaire gauche
Ex = amp * cos(ang_freq * time);
Ey = amp * sin(ang_freq * time);
- Résultats et analyse de la simulation
L'exécution de la simulation révèle la distribution du champ électrique autour des nanopores, avec une localisation et une déformation intenses, indicative d'une réponse chérale. Le dichroïsme circulaire (CD) est calculé en intégrant la différence d'absorption entre les polarisations droite et gauche sur une bande de fréquences spécifiée :
// Calcul du dichroïsme circulaire
dc_spectrum = numerical_integration(abs(absorb_RCP - absorb_LCP), freq_low, freq_high);
Les variables absorb_RCP et absorb_LCP dénotent l'absorption des polarisations droite et gauche, tandis que freq_low et freq_high définissent l'intervalle de fréquences d'intérêt. Une valeur non nulle du CD confirme la présence d'une réponse chérale. L'analyse du spectre CD permet d'identifier les résonances et les caractéristiques clés de la métasurface.
Cette approche simulée offre un outil puissant pour investiguer les propriétés optiques chérales à l'échelle nanométrique, facilitant la conception de dispositifs photoniques innovants.