مادة ميتامادية جرافين غير متجانسة في تيراهيرتز متعددة النطاقات تعتمد على تكوين شريطي مزدوج متعامد
Anisotropic terahertz graphene metamaterial multi-band absorber based on an orthogonal dual-ribbon configuration

المجلة: Optical and Quantum Electronics، المجلد: 58، العدد: 2
DOI: https://doi.org/10.1007/s11082-026-08673-y
تاريخ النشر: 2026-01-09
المؤلف: Somayyeh Asgari وآخرون
الموضوع الرئيسي: تطبيقات المواد الميتامادية والأسطح الميتامادية

نظرة عامة

تبحث هذه الدراسة في مادة امتصاصية مصنوعة من الجرافين غير المتجانسة القابلة لإعادة التكوين ديناميكيًا، والتي تعمل ضمن نطاق تردد 0.5-4.7 تيراهيرتز. يتميز تصميم وحدة الخلية بوجود زوجين متعامدين من شرائط الجرافين، مما يسهل الاستجابات الكهرومغناطيسية المعتمدة على الاتجاه والاستقطاب. من خلال المحاكاة العددية التي أجريت باستخدام CST Studio Suite ونموذج الدائرة المكافئ (ECM) الذي تم تحليله في MATLAB، تحدد الدراسة الرنينات الكهربائية المستعرضة (TE) والرنينات المغناطيسية المستعرضة (TM).

عند طاقة فيرمي تبلغ 0.6 إلكترون فولت، يظهر الممتص قابلية ضبط كبيرة وثنائية انكسار خطية عالية (LD) تصل إلى 99.1%. تدعم البنية أربعة رنينات عالية الامتصاص من نوع TE ورنينتين من نوع TM، محققة معدلات امتصاص تتجاوز 94.5% لـ TE و82% لـ TM، مع قيم امتصاص قصوى تصل إلى 99.68% (TE) و99.92% (TM). تشير هذه الخصائص الطيفية الفريدة لـ TE/TM إلى تطبيقات واعدة في أجهزة THz الحساسة للاستقطاب، والمنصات الضوئية، وتقنيات التعرف الآمن المشفر بالتردد (FESID).

مقدمة

تناقش المقدمة خصائص وتطبيقات المواد غير المتجانسة، وخاصة تلك المعتمدة على الجرافين، والتي تظهر استجابات تعتمد على الاتجاه والاستقطاب لموجات الكهرومغناطيسية. يتم تعريف المواد غير المتجانسة من خلال هندستها المعتمدة على الاتجاه، وخصائص المواد الجوهرية، أو المجالات المغناطيسية الخارجية، مما يؤدي إلى وظائف فريدة مثل ثنائية الانكسار الخطية (LD). أظهرت التطورات الأخيرة في المواد غير المتجانسة المعتمدة على الجرافين قدرتها على دعم استجابات تعتمد على الاستقطاب قابلة للضبط، مما يمكّن من تصميم ممتصات تيراهيرتز (THz) متعددة النطاقات مع رنينات كهربائية مستعرضة (TE) ورنينات مغناطيسية مستعرضة (TM) متميزة. هذه الميزات مفيدة بشكل خاص للتطبيقات الضوئية الحساسة للاستقطاب، بما في ذلك الاستشعار عالي الدقة والترشيح الكهرومغناطيسي.

على الرغم من مزايا التصاميم غير الحساسة للاستقطاب، إلا أنها تحد من الوصول إلى الوظائف المعتمدة على الاستقطاب. يسلط البحث الضوء على الحاجة إلى ممتصات جرافين مضغوطة وقابلة للضبط كهربائيًا تجمع بين عدم التجانس القوي، والتشغيل متعدد النطاقات، وكفاءة امتصاص عالية. يتميز التصميم المقترح بهندسة مزدوجة شرائط متعامدة تدعم ستة أوضاع رنين متميزة – أربعة TE واثنان TM – يمكن ضبط مواقعها الطيفية ديناميكيًا عن طريق تعديل طاقة فيرمي للجرافين. لا تعزز هذه البنية المبتكرة كفاءة الامتصاص وثنائية الانكسار الخطية فحسب، بل تمنع أيضًا تداخل الطيف بين أوضاع TE وTM، مما يمكّن من الضبط المستقل. تشير النتائج إلى أن هذا الممتص يمكن أن يعمل كمنصة قوية لأنظمة THz الضوئية المتقدمة وتطبيقات مثل التعرف الآمن المشفر بالتردد (FESID).

طرق

تستعرض هذه القسم الطرق المستخدمة في تصميم وتحليل ممتص مادة تيراهيرتز غير متجانسة معتمدة على الجرافين. تتكون البنية من ثلاث طبقات: مصفوفة رنانات الجرافين، وركيزة PTFE، وطبقة خلفية من الذهب، مع طبقة هلام أيوني لتوجيه الجرافين. يحتوي الهلام الأيوني على سماحية قدرها $\epsilon_{ig} = 2.0164$ وسمك يبلغ 150 نانومتر، بينما طبقة الذهب، التي تعمل كقاعدة موصلة، لها سمك يبلغ 0.5 ميكرومتر. تم تعيين طاقة فيرمي للجرافين بين $E_F = 0.6 – 0.8$ إلكترون فولت، متوافقة مع القيم التجريبية، ويفترض أن أوقات استرخاء الحاملات تتراوح بين 1-3 بيكو ثانية. يهدف التصميم إلى زيادة الامتصاص لكل من الاستقطابات TE وTM وتعزيز ثنائية الانكسار الخطية عبر نطاق ترددات التيراهيرتز.

تستخدم التحليل العددي طريقة العناصر المحدودة (FEM) ضمن CST Microwave Studio 2018، مع حدود دورية في الاتجاهين x وy وظروف حدود امتصاصية في الاتجاه z. تم تعيين فترات وحدة الخلية للمادة غير المتجانسة عند $P_x = P_y = 21$ ميكرومتر، وهي أصغر من الطول الموجي الأدنى البالغ 63.83 ميكرومتر، لتقليل أوضاع فلوكيت من الدرجة الأعلى. تم تطوير نموذج الدائرة المكافئ (ECM) لطبقة رنان الجرافين لكل من أوضاع TE وTM، مع حسابات محددة للمقاومة والسعة بناءً على اتجاه المجال الكهربائي. تؤكد الدراسة على نهج منهجي لتحسين المعلمات دون تعديلات حدسية، مما يضمن جدوى التصميم لتحقيق التنفيذ التجريبي.

نتائج

تظهر النتائج أن تكوين الشرائط المزدوجة المتعامدة للممتص غير المتجانس المقترح يمكّن من تفعيل كامل للأوضاع وضبط مستقل للرنينات في كل من أوضاع الكهرباء المستعرضة (TE) والمغناطيسية المستعرضة (TM). تكشف المحاكاة أن الممتص يعمل بفعالية عبر نطاق 0.5-4.7 تيراهيرتز، محققًا مستويات امتصاص قصوى تبلغ 99.92% لأوضاع TE و99.68% لأوضاع TM، مع ثنائية انكسار خطية (LD) تبلغ 99.1%. تبرز بصمات الطيف المتميزة، التي تتميز بأربعة رنينات TE ورنينتين TM، إمكانيات الجهاز لتطبيقات التعرف الآمن المشفر بالتردد (FESID).

تستكشف الدراسة أيضًا قابلية الضبط الديناميكي للممتص من خلال تغيير طاقة فيرمي للجرافين، مما يؤدي إلى انزياح أزرق لترددات الرنين. يحدث أقوى امتصاص عند طاقة فيرمي تبلغ 0.6 إلكترون فولت، حيث تصل LD إلى ذروتها عند 99.1%. تشير الدراسات المعلماتية الشاملة إلى أن الممتص يحافظ على خصائص رنين مستقرة وأداء امتصاص عالي على الرغم من التغيرات في ت tolerances التصنيع، مما يبرز قوته للتطبيقات العملية. بشكل عام، تضع النتائج الممتص غير المتجانس المقترح كمرشح واعد لأنظمة THz المتقدمة، خاصة في الأجهزة الانتقائية للاستقطاب والمصادقة متعددة الترددات الآمنة.

نقاش

في هذا القسم، يستنتج المؤلفون المعادلات الرئيسية المتعلقة بالأطوال الكهربائية والمقاومات لهلام أيوني وركيزة عازلة داخل ممتص مادة غير متجانسة، تحديدًا في أوضاع الكهرباء المستعرضة (TE) والمغناطيسية المستعرضة (TM). يتم حساب الأطوال الكهربائية، التي يشار إليها بـ $\theta_{ig}$ و$\theta_d$، باستخدام العلاقة بين زاوية السقوط $\theta_{in}$ وخصائص العازل، مع توفير المعادلات (11) إلى (15) الصيغ الخاصة بالمقاومة لكل من الأوضاع.

يتم صياغة مصفوفة T العامة لممتص المادة غير المتجانسة من خلال دمج مصفوفات T الفردية لهلام الأيون، وطبقة الذهب، والركيزة العازلة، كما هو معبر عنه في المعادلة (16). يتم اشتقاق المقاومة المدخلة للمادة غير المتجانسة في المعادلة (17)، مما يبسط إلى شكل يبرز المساهمات من عناصر المصفوفة. علاوة على ذلك، يتم حساب معاملات الانعكاس لأوضاع TE وTM باستخدام المعادلات (18) و(19)، مما يؤدي إلى معاملات الامتصاص، كما هو موضح في المعادلة (20). توضح هذه النتائج أداء المادة غير المتجانسة من حيث الانعكاس والامتصاص، وهو أمر حاسم للتطبيقات في معالجة الموجات الكهرومغناطيسية.

Journal: Optical and Quantum Electronics, Volume: 58, Issue: 2
DOI: https://doi.org/10.1007/s11082-026-08673-y
Publication Date: 2026-01-09
Author(s): Somayyeh Asgari et al.
Primary Topic: Metamaterials and Metasurfaces Applications

Overview

This research investigates a dynamically reconfigurable anisotropic graphene-based metamaterial absorber that operates within the 0.5-4.7 THz frequency range. The unit-cell design features two orthogonal pairs of graphene ribbons, which facilitate direction- and polarization-dependent electromagnetic responses. Through numerical simulations conducted with CST Studio Suite and an equivalent circuit model (ECM) analyzed in MATLAB, the study identifies transverse electric (TE) and transverse magnetic (TM) resonances.

At a Fermi energy of 0.6 eV, the absorber exhibits significant tunability and a high linear dichroism (LD) of up to 99.1%. The structure supports four TE and two TM high-absorption resonances, achieving absorption rates exceeding 94.5% for TE and 82% for TM, with peak absorption values reaching 99.68% (TE) and 99.92% (TM). These unique TE/TM spectral characteristics indicate promising applications in polarization-sensitive THz devices, photonic platforms, and Frequency-Encoded Secure Identification (FESID) technologies.

Introduction

The introduction discusses the characteristics and potential applications of anisotropic metamaterials, particularly those based on graphene, which exhibit direction- and polarization-dependent responses to electromagnetic waves. Anisotropic metamaterials are defined by their orientation-dependent geometry, intrinsic material properties, or external magnetic fields, leading to unique functionalities such as linear dichroism (LD). Recent advancements in graphene-based anisotropic metamaterials have demonstrated their capability to support tunable polarization-dependent responses, enabling the design of multi-band terahertz (THz) absorbers with distinct transverse electric (TE) and transverse magnetic (TM) resonances. These features are particularly beneficial for polarization-sensitive photonic applications, including high-resolution sensing and electromagnetic filtering.

Despite the advantages of polarization-insensitive designs, they limit access to polarization-dependent functionalities. The paper highlights the need for compact, electrically tunable graphene-based absorbers that combine strong anisotropy, multi-band operation, and high absorption efficiency. The proposed design features an orthogonal dual-ribbon geometry that supports six distinct resonance modes—four TE and two TM—whose spectral positions can be dynamically tuned by modulating the graphene Fermi energy. This innovative structure not only enhances absorption efficiencies and linear dichroism but also prevents spectral overlap between TE and TM modes, thereby enabling independent tunability. The findings suggest that this absorber could serve as a robust platform for advanced THz photonic systems and applications such as Frequency-Encoded Secure Identification (FESID).

Methods

The section outlines the methods used to design and analyze an anisotropic terahertz graphene-based metamaterial absorber. The structure comprises three layers: a graphene resonator array, a PTFE substrate, and a gold backing layer, with an ion-gel layer for biasing the graphene. The ion gel has a permittivity of $\epsilon_{ig} = 2.0164$ and a thickness of 150 nm, while the gold layer, serving as a conductive base, has a thickness of 0.5 μm. The graphene’s Fermi energy is set between $E_F = 0.6 – 0.8$ eV, consistent with experimental values, and the carrier relaxation times are assumed to be between 1-3 ps. The design aims to maximize absorption for both TE and TM polarizations and enhance linear dichroism across the terahertz frequency range.

The numerical analysis employs the finite element method (FEM) within CST Microwave Studio 2018, with periodic boundaries in the x and y directions and an absorbing boundary condition in the z direction. The metamaterial’s unit-cell periods are set at $P_x = P_y = 21$ μm, smaller than the minimum wavelength of 63.83 μm, to suppress higher-order Floquet modes. The equivalent circuit model (ECM) for the graphene resonator layer is developed for both TE and TM modes, with specific impedance and capacitance calculations based on the orientation of the electric field. The study emphasizes a systematic approach to parameter optimization without heuristic adjustments, ensuring the design’s feasibility for experimental realization.

Results

The results demonstrate that the orthogonal dual-ribbon configuration of the proposed anisotropic metamaterial absorber enables full mode activation and independent tunability of resonances in both transverse electric (TE) and transverse magnetic (TM) modes. Simulations reveal that the absorber operates effectively across the 0.5-4.7 THz range, achieving peak absorption levels of 99.92% for TE and 99.68% for TM modes, with a linear dichroism (LD) of 99.1%. The distinct spectral fingerprints, characterized by four TE resonances and two TM resonances, highlight the device’s potential for Frequency-Encoded Secure Identification (FESID) applications.

The study further explores the dynamic tunability of the absorber by varying the graphene Fermi energy, which results in a blueshift of resonance frequencies. The strongest absorption occurs at a Fermi energy of 0.6 eV, where the LD peaks at 99.1%. Extensive parametric studies indicate that the absorber maintains stable resonance characteristics and high absorption performance despite variations in fabrication tolerances, underscoring its robustness for practical applications. Overall, the findings position the proposed metamaterial absorber as a promising candidate for advanced THz systems, particularly in polarization-selective devices and secure multi-frequency authentication.

Discussion

In this section, the authors derive key equations related to the electrical lengths and impedances of an ion gel and dielectric substrate within a metamaterial absorber, specifically in transverse electric (TE) and transverse magnetic (TM) modes. The electrical lengths, denoted as $\theta_{ig}$ and $\theta_d$, are calculated using the relationship between the incident angle $\theta_{in}$ and the dielectric properties, with equations (11) to (15) providing the respective impedance formulas for both modes.

The overall T-matrix of the metamaterial absorber is formulated by combining the individual T-matrices of the ion gel, gold layer, and dielectric substrate, as expressed in equation (16). The input impedance for the metamaterial is derived in equation (17), simplifying to a form that highlights the contributions from the matrix elements. Furthermore, the reflection coefficients for TE and TM modes are calculated using equations (18) and (19), leading to the absorption coefficients, as shown in equation (20). These findings illustrate the metamaterial’s performance in terms of reflection and absorption, which are critical for applications in electromagnetic wave manipulation.