Wide-Oblique-Incident-Angle Stable Polarization-Insensitive Ultra-Wideband Metamaterial Perfect Absorber for Visible Optical Wavelength Applications

Metamaterial absorbers are very attractive due to their significant absorption behavior at optical wavelengths, which can be implemented for energy harvesting, plasmonic sensors, imaging, optical modulators, photovoltaic detectors, etc. This paper presents a numerical study of an ultra-wide-band double square ring (DSR) metamaterial absorber (MMA) for the complete visible optical wavelength region, which is designed with a three-layer (tungsten-silicon dioxide-tungsten) substrate material. Due to the symmetricity, a polarization-insensitive absorption is obtained for both transverse electric (TE) and transverse magnetic (TM) modes by simulation. An absorption above 92.2% and an average absorption of 97% are achieved in the visible optical wavelength region. A peak absorption of 99.99% is achieved at 521.83 nm. A wide range of oblique incident angle stabilities is found for stable absorption properties. A similar absorption is found for different banding angles, which may occur due to external forces during the installation of the absorber. The absorption is calculated by the interference theory (IT) model, and the polarization conversion ratio (PCR) is also validated to verify the perfect MMA. The electric field and magnetic field of the structure analysis are performed to understand the absorption property of the MMA. The presented MMA may be used in various applications such as solar cells, light detection, the biomedical field, sensors, and imaging.


Introduction
Veselago first demonstrated the unusual electromagnetic properties such as negative permittivity or permeability of complex materials [1], which is known as metamaterials (MMs) [2]. The characteristics of MMs depend on its physical structure, not its chemical composition, which makes it suitable for use in various applications such as absorbers [3], filters [4], antennae [5], imaging [6], invisible clocks [7], sensors [8], polarization converters [9], and lenses [10]. Landy first introduced the metamaterial perfect absorber (MPA) [11], Figure 1 shows the three-layer metamaterial absorber structure, where the middle layer uses silicon dioxide (SiO 2 ) as the dielectric substrate and tungsten (W) as a patch and ground on both sides of the SiO 2 . The material properties of SiO 2 and W were considered from [28]. The purpose of using W in the patch and ground of this design was due to the high intrinsic loss and impedance matching feature with free space in optical wavelength applications. Therefore, this leads to an increased absorption behavior by offering low reflection and transmission properties [29]. The SiO 2 was selected as a dielectric substrate due to its lossless property at optical wavelengths and large negative part of permittivity in the whole optical region instated of the large imaginary part. These features of SiO 2 provide a good impedance matching and lead to a high and more comprehensive absorption bandwidth [30]. The proposed DSR metamaterial absorber was designed and simulated using Frequency Integration Technique (FIT)-based Computer Simulation Technology (CST) [31]. The overall design parameters are listed in Table 1. Unit cell boundary conditions on the x-and y-axis and open add space on the z-axis were applied for the simulation. The absorption property A(ω) was calculated by using Equation (1) [8].

Absorber Design
where S 11 and S 21 are the reflection coefficient and transmission coefficient, respectively. R and T are the reflectance and transmittance, respectively. Due to the thick back layer, the transmission of electromagnetic waves blocks and the transmission coefficient S 21 become zero.
Materials 2022, 15, x FOR PEER REVIEW 3 of 14 absorption bandwidth [30]. The proposed DSR metamaterial absorber was designed and simulated using Frequency Integration Technique (FIT)-based Computer Simulation Technology (CST) [31]. The overall design parameters are listed in Table 1. Unit cell boundary conditions on the x-and y-axis and open add space on the z-axis were applied for the simulation. The absorption property A(ω) was calculated by using Equation (1) [8].
where S11 and S21 are the reflection coefficient and transmission coefficient, respectively. R and T are the reflectance and transmittance, respectively. Due to the thick back layer, the transmission of electromagnetic waves blocks and the transmission coefficient S21 become zero.   Figure 2a presents the absorption properties for the different layers, where no absorption was observed only for the SiO2 dielectric layer. Due to all incident waves transmitted through the dielectric material, an absorption around 20% was realized by adding a patch on top of the dielectric, but the maximum incident wave still transmitted through the material. After that, the W base ground plane was added to prevent the transmission of the incident wave and a 97.1% average absorption level was attained in the entire visible optical wavelength region. The geometrical design of the patch also played a crucial role in achieving higher and wider absorption levels, which can be understood from the design evolution of the proposed MMA. The absorption property of the proposed MMA of various design evolutions is presented in Figure 2b. The peak absorption reached 99% at 553.63 nm, and an average absorption of 92.46% was achieved for design 1, where only a square resonator was used. In design two, a square slot was made in the patch of design 1, which increased the average absorption level to 94.46%, and the peak absorption shifted toward 495 nm. By adding another square patch inside the square ring in design 2, the   Figure 2a presents the absorption properties for the different layers, where no absorption was observed only for the SiO 2 dielectric layer. Due to all incident waves transmitted through the dielectric material, an absorption around 20% was realized by adding a patch on top of the dielectric, but the maximum incident wave still transmitted through the material. After that, the W base ground plane was added to prevent the transmission of the incident wave and a 97.1% average absorption level was attained in the entire visible optical wavelength region. The geometrical design of the patch also played a crucial role in achieving higher and wider absorption levels, which can be understood from the design evolution of the proposed MMA. The absorption property of the proposed MMA of various design evolutions is presented in Figure 2b. The peak absorption reached 99% at 553.63 nm, and an average absorption of 92.46% was achieved for design 1, where only a square resonator was used. In design two, a square slot was made in the patch of design 1, which increased the average absorption level to 94.46%, and the peak absorption shifted toward 495 nm. By adding another square patch inside the square ring in design 2, the average absorption increased up to 96.70% in design 3. Finally, an average absorption of 97.1% was attained by making a square slot at the center of the square plane in design 3. (a) (b) Figure 2. Absorption of (a) different layers and (b) design evolution.
The parametric analysis was also very important for understanding the absorption properties of the proposed MMA. The parameter h1 represents the thickness of W at the ground plane. The absorption property was investigated from 20 to 140 nm with a 20 nm interval for h1 and is presented in Figure 3a. The absorption level increased with h1 to 80 nm and displayed an almost unique absorption curve for higher thicknesses. The MMA design chose the h1 = 120 nm parametric value, which completely blocked the transmission and provided a higher absorption level. The dielectric thickness (SiO2) also significantly influenced higher absorption levels. The parameters h2 represents the dielectric thickness, and the absorption at different dielectric thicknesses is presented in Figure 3b. It is shown that the absorption value was higher in the lower optical wavelength region than the upper optical wavelength region for h2 = 40 nm and 50 nm. On the other hand, for higher thicknesses, the absorption value was higher in the upper optical wavelength region than the lower optical wavelength region. The maximum average absorption level was achieved at h2 = 60 nm. The thickness of the MMA patch was also investigated for h3. A constant absorption was found for different h3 values with slight distortion in the upper optical wavelength region. For maximum absorption level, h3 = 15 nm was settled. Figure  3d shows the impact of the outer square ring resonator from 660 to 675 nm with a 5 nm interval, where a stable absorption bandwidth was realized, and very few deformations were found in the upper wavelength region. W1 = 670 nm was chosen for the higher absorption level. The absorption at different values of the inner square ring length (W2) from 310 to 390 nm is presented in Figure 3e, which shows a slight distortion in the absorption curve. Figure 3f demonstrates the effect of inner ring slot length (W3) from 130 to 160 nm, showing a stable absorption for different values. The investigation showed that the optimized top layer parameters (h3, W1, W2, and W3) of the proposed design will provide a sound absorption level.  The parametric analysis was also very important for understanding the absorption properties of the proposed MMA. The parameter h1 represents the thickness of W at the ground plane. The absorption property was investigated from 20 to 140 nm with a 20 nm interval for h1 and is presented in Figure 3a. The absorption level increased with h1 to 80 nm and displayed an almost unique absorption curve for higher thicknesses. The MMA design chose the h1 = 120 nm parametric value, which completely blocked the transmission and provided a higher absorption level. The dielectric thickness (SiO 2 ) also significantly influenced higher absorption levels. The parameters h2 represents the dielectric thickness, and the absorption at different dielectric thicknesses is presented in Figure 3b. It is shown that the absorption value was higher in the lower optical wavelength region than the upper optical wavelength region for h2 = 40 nm and 50 nm. On the other hand, for higher thicknesses, the absorption value was higher in the upper optical wavelength region than the lower optical wavelength region. The maximum average absorption level was achieved at h2 = 60 nm. The thickness of the MMA patch was also investigated for h3. A constant absorption was found for different h3 values with slight distortion in the upper optical wavelength region. For maximum absorption level, h3 = 15 nm was settled. Figure 3d shows the impact of the outer square ring resonator from 660 to 675 nm with a 5 nm interval, where a stable absorption bandwidth was realized, and very few deformations were found in the upper wavelength region. W1 = 670 nm was chosen for the higher absorption level. The absorption at different values of the inner square ring length (W2) from 310 to 390 nm is presented in Figure 3e, which shows a slight distortion in the absorption curve. Figure 3f demonstrates the effect of inner ring slot length (W3) from 130 to 160 nm, showing a stable absorption for different values. The investigation showed that the optimized top layer parameters (h3, W1, W2, and W3) of the proposed design will provide a sound absorption level.
were found in the upper wavelength region. W1 = 670 nm was chosen for the higher absorption level. The absorption at different values of the inner square ring length (W2) from 310 to 390 nm is presented in Figure 3e, which shows a slight distortion in the absorption curve. Figure 3f demonstrates the effect of inner ring slot length (W3) from 130 to 160 nm, showing a stable absorption for different values. The investigation showed that the optimized top layer parameters (h3, W1, W2, and W3) of the proposed design will provide a sound absorption level.

Absorption Characteristics
The absorption characteristics depend on the impedance matching of the structure. The relative impedance (Z) of the proposed three-layer sandwich model was calculated using Equation (2) [32], and Figure 4a presents the relative impedance of both TE and TM modes. The near-unity value of the real part and near null value of the imaginary part shows that the proper effective impedance of the structure matches with the free space impedance, which offers a high absorption level. Thus, the free space impedance is 0 The reflectance (R) of the TE and TM mode can be calculated by Equations (3) and (4), respectively [17]:

Absorption Characteristics
The absorption characteristics depend on the impedance matching of the structure. The relative impedance (Z) of the proposed three-layer sandwich model was calculated using Equation (2) [32], and Figure 4a presents the relative impedance of both TE and TM modes. The near-unity value of the real part and near null value of the imaginary part shows that the proper effective impedance of the structure matches with the free space impedance, which offers a high absorption level.
where µ = µ r µ 0 and ε = ε r ε are the permeability and permittivity of the MMA, respectively. µ r and ε r are the relative permeability and permittivity of the medium, respectively. µ 0 and ε 0 are the permeability and permittivity of free space, respectively. Thus, the free space impedance is Z 0 = µ 0 /ε 0 . The reflectance (R) of the TE and TM mode can be calculated by Equations (3) and (4), respectively [17]: where n is the refractive index, and θ is the incident angle of the wave for typical incidence. Thus, Equations (3) and (4) become: Materials 2022, 15, x FOR PEER REVIEW 6 of 14 Equation (5) represents the reflectance of the MMA, which is significantly controlled by impedance matching as well as the metamaterial. As the transmission is zero due to the back layer, absorption (A) can be calculated by The reflectance, transmittance, and absorption characteristics of TE and TM modes are illustrated in Figure 4b. The symmetrical design of the proposed MMA provides a unique absorption curve for both TE and TM modes. A peak absorption of 99.99% is realized at a wavelength of 523 nm, and an absorption above 99% appears from 474.36 to 578.75 nm. In addition, an absorption more than 92.2% is obtained for the entire visible optical wavelength region with an average absorption of 97.1%. This near-unity absorption feature of the proposed MMA can be used for solar energy harvesting at visible optical wavelengths, stealth technology, sensors, and micro-imaging technology [21].  Absorption characteristics can also be calculated by applying the Interference Theory (IT) model, which also explains the underlying physics of the metamaterial absorber [33,34]. The patch (layer 1) with a particular metallic pattern acts as a partial reflection surface. This layer modifies complex transmission and reflection coefficients, and the ground plane (layer 3) works as a perfect reflector. Figure 5a shows the patch between air and a dielectric substrate. The incident EM wave partially transmits and reflects the at-  Equation (5) represents the reflectance of the MMA, which is significantly controlled by impedance matching as well as the metamaterial. As the transmission is zero due to the back layer, absorption (A) can be calculated by A = 1 − R. The reflectance, transmittance, and absorption characteristics of TE and TM modes are illustrated in Figure 4b. The symmetrical design of the proposed MMA provides a unique absorption curve for both TE and TM modes. A peak absorption of 99.99% is realized at a wavelength of 523 nm, and an absorption above 99% appears from 474.36 to 578.75 nm. In addition, an absorption more than 92.2% is obtained for the entire visible optical wavelength region with an average absorption of 97.1%. This near-unity absorption feature of the proposed MMA can be used for solar energy harvesting at visible optical wavelengths, stealth technology, sensors, and micro-imaging technology [21].
Absorption characteristics can also be calculated by applying the Interference Theory (IT) model, which also explains the underlying physics of the metamaterial absorber [33,34]. The patch (layer 1) with a particular metallic pattern acts as a partial reflection surface. This layer modifies complex transmission and reflection coefficients, and the ground plane (layer 3) works as a perfect reflector. Figure 5a shows the patch between air and a dielectric substrate. The incident EM wave partially transmits and reflects the atmosphere with transmission coefficient t 12 (ω) = t 12 (ω)e iθ 12 (ω) and reflection coefficient r 12 (ω) = r 12 (ω)e iϕ 12 (ω) . The transmitted wave constantly propagates until reflected from the ground plane. The complex propagation constant of the substrate layer can be written by β = β 1 + iβ 2 = √ ε d k 0 d, where K 0 is the free space wavenumber, β 1 is the propagation phase, absorption of the dielectric is represented by β 2 , d is the substrate thickness, and ε d is the substrate dielectric constant. The entire wave is reflected due to the blocking of transmission by the ground layer after propagation phase delay β. The wave propagates toward the patch (layer 1) and is partially transmitted and reflected with transmission coefficient t 21 (ω) = t 21 (ω)e iθ 21 (ω) and reflection coefficient r 21 (ω) = r 21 (ω)e iϕ 21 (ω) , respectively. Multiple reflection and transmission occur with the substrate, and total reflection is the superposition of all orders, which can be calculated by Equation (6). Figure 5b shows sound agreement between the IT model and CST-calculated absorption.

Polarization Conversion Ratio (PCR)
The PCR value of the proposed MMA calculated for verifying the proposed design is not a polarization converter but an absorber. Figure 6a shows the co-and cross-polarization as Txx and Tyx of x-polarization of the incident wave, respectively. The unique coand cross-polarization of the x-and y-polarization of the incident wave were also achieved due to the structural symmetricity. The PCR value was calculated by Equations (7) and (8) for y and x polarization, respectively [35]. The PCR value was calculated as near-zero because of the very small cross-polarization Tyx = Txy < −90 dB. Figure 6b shows the PCR value of x-and y-polarized waves, which verify that the proposed design is an absorber, not a polarization converter.
X to y polarized wave ( ) The above calculation was performed for the normal incident angle. For oblique incident angle θ, the propagation phase delay was modified because the propagation length increases, d = d/cos θ , and the phase delay becomes β = √ ε d k 0 d ; the refractive angle θ' can be calculated by Snell's law √ ε d sin θ = sin θ.

Polarization Conversion Ratio (PCR)
The PCR value of the proposed MMA calculated for verifying the proposed design is not a polarization converter but an absorber. Figure 6a shows the co-and cross-polarization as Txx and Tyx of x-polarization of the incident wave, respectively. The unique co-and cross-polarization of the x-and y-polarization of the incident wave were also achieved due to the structural symmetricity. The PCR value was calculated by Equations (7) and (8) for y and x polarization, respectively [35]. The PCR value was calculated as near-zero because of the very small cross-polarization Tyx = Txy < −90 dB. Figure 6b shows the PCR value of x-and y-polarized waves, which verify that the proposed design is an absorber, not a polarization converter.
Y to x polarized wave, PCR y = T 2 xy / T 2 yy + T 2 xy (7) X to y polarized wave PCR  Figure 7a,b present the metamaterial properties of the proposed MMA for TE and TM modes, respectively. In Figure 7a, the permittivity is negative from 427 to 541 nm and from 605 to 750 nm. The permeability is negative from 380 to 394 nm and from 506 to 750 nm. Both permeability and permittivity are negative from 505 to 540 nm and from 602 to 750 nm. A negative value of refractive index appears from 472 to 750 nm. As a result of this negative value of metamaterial property, the proposed MMA exhibits high-level broadband absorption properties in the TM mode. The TE mode alternatively achieves the negative value of metamaterial property, where a negative value of permeability is reached from 514 to 750 nm. A unique absorption curve is completed due to the structural symmetricity shown in Figure 4.

Polarization Insensitivity and Oblique Incident Angle Stability
The polarization insensitivity of the proposed MMA was investigated to verify its absorption efficiency. Figure 8 shows the absorption property of different polarization incident angles (ф) for TE and TM modes. The TE mode z-axis is the direction of wave propagation, and the magnetic field vector Hz is along the z-axis. The electric field vector (Ex) and magnetic field vector (Hy) are along the x-and y-axis, respectively. On the other hand, in the TM mode, the electric field vector (Ez) is along the direction of wave propagation, and the magnetic (Hx) and electric (Ey) field vectors are along the x-and y-axis, respectively. Due to the axial and rotational symmetricity, the proposed MMA achieves unique absorption properties for polarization incident angles (ф) up to 180°. In all the above discussions, all results have been obtained from a normal incident angle (θ = 0°), but the EM wave occurs at an oblique incident angle on the MMA. Therefore, discussion of the  Figure 7a,b present the metamaterial properties of the proposed MMA for TE and TM modes, respectively. In Figure 7a, the permittivity is negative from 427 to 541 nm and from 605 to 750 nm. The permeability is negative from 380 to 394 nm and from 506 to 750 nm. Both permeability and permittivity are negative from 505 to 540 nm and from 602 to 750 nm. A negative value of refractive index appears from 472 to 750 nm. As a result of this negative value of metamaterial property, the proposed MMA exhibits high-level broadband absorption properties in the TM mode. The TE mode alternatively achieves the negative value of metamaterial property, where a negative value of permeability is reached from 514 to 750 nm. A unique absorption curve is completed due to the structural symmetricity shown in Figure 4.  Figure 7a,b present the metamaterial properties of the proposed MMA for TE and TM modes, respectively. In Figure 7a, the permittivity is negative from 427 to 541 nm and from 605 to 750 nm. The permeability is negative from 380 to 394 nm and from 506 to 750 nm. Both permeability and permittivity are negative from 505 to 540 nm and from 602 to 750 nm. A negative value of refractive index appears from 472 to 750 nm. As a result of this negative value of metamaterial property, the proposed MMA exhibits high-level broadband absorption properties in the TM mode. The TE mode alternatively achieves the negative value of metamaterial property, where a negative value of permeability is reached from 514 to 750 nm. A unique absorption curve is completed due to the structural symmetricity shown in Figure 4.

Polarization Insensitivity and Oblique Incident Angle Stability
The polarization insensitivity of the proposed MMA was investigated to verify its absorption efficiency. Figure 8 shows the absorption property of different polarization incident angles (ф) for TE and TM modes. The TE mode z-axis is the direction of wave propagation, and the magnetic field vector Hz is along the z-axis. The electric field vector (Ex) and magnetic field vector (Hy) are along the x-and y-axis, respectively. On the other hand, in the TM mode, the electric field vector (Ez) is along the direction of wave propagation, and the magnetic (Hx) and electric (Ey) field vectors are along the x-and y-axis, respectively. Due to the axial and rotational symmetricity, the proposed MMA achieves unique absorption properties for polarization incident angles (ф) up to 180°. In all the above discussions, all results have been obtained from a normal incident angle (θ = 0°), but the EM wave occurs at an oblique incident angle on the MMA. Therefore, discussion of the

Polarization Insensitivity and Oblique Incident Angle Stability
The polarization insensitivity of the proposed MMA was investigated to verify its absorption efficiency. Figure 8 shows the absorption property of different polarization incident angles (ф) for TE and TM modes. The TE mode z-axis is the direction of wave propagation, and the magnetic field vector Hz is along the z-axis. The electric field vector (Ex) and magnetic field vector (Hy) are along the x-and y-axis, respectively. On the other hand, in the TM mode, the electric field vector (Ez) is along the direction of wave propagation, and the magnetic (Hx) and electric (Ey) field vectors are along the x-and y-axis, respectively. Due to the axial and rotational symmetricity, the proposed MMA achieves unique absorption properties for polarization incident angles (ф) up to 180 • . In all the above discussions, all results have been obtained from a normal incident angle (θ = 0 • ), but the EM wave occurs at an oblique incident angle on the MMA. Therefore, discussion of the absorption behavior of oblique incidents is also vital. The oblique incident angle (θ) is made between the direction of wave propagation and the z-axis of the MMA structure. Figure 9a,b show the absorption curve of both TE and TM modes for the oblique incident angles (θ) up to 60 • and 70 • , respectively, where an operating frequency above 70% absorption is found. Absorption properties of various oblique incident angles are listed in Table 2.
Materials 2022, 15, x FOR PEER REVIEW 9 of 14 absorption behavior of oblique incidents is also vital. The oblique incident angle (θ) is made between the direction of wave propagation and the z-axis of the MMA structure. Figure 9a,b show the absorption curve of both TE and TM modes for the oblique incident angles (θ) up to 60° and 70°, respectively, where an operating frequency above 70% absorption is found. Absorption properties of various oblique incident angles are listed in Table 2.

Bendable Property
The MMA sheet may be curved due to the external effect during installation [36]. The banding effect on absorption behavior was investigated for both TE and TM modes. Figure 10 shows the absorption curve of banding angles θ = 0° to 18° with 3° intervals. In the TE mode, the absorption increases in the lower wavelength region and slightly decreases absorption behavior of oblique incidents is also vital. The oblique incident angle (θ) is made between the direction of wave propagation and the z-axis of the MMA structure. Figure 9a,b show the absorption curve of both TE and TM modes for the oblique incident angles (θ) up to 60° and 70°, respectively, where an operating frequency above 70% absorption is found. Absorption properties of various oblique incident angles are listed in Table 2.

Bendable Property
The MMA sheet may be curved due to the external effect during installation [36]. The banding effect on absorption behavior was investigated for both TE and TM modes. Figure 10 shows the absorption curve of banding angles θ = 0° to 18° with 3° intervals. In the TE mode, the absorption increases in the lower wavelength region and slightly decreases

Bendable Property
The MMA sheet may be curved due to the external effect during installation [36]. The banding effect on absorption behavior was investigated for both TE and TM modes. Figure 10 shows the absorption curve of banding angles θ = 0 • to 18 • with 3 • intervals. In the TE mode, the absorption increases in the lower wavelength region and slightly decreases in the upper wavelength region with the increment in θ. The same absorption behavior is also shown in the TM mode. A unique absorption property is achieved for both TE and TM modes of various banding angles, which offers a bendable property of the proposed absorber. in the upper wavelength region with the increment in θ. The same absorption behavior is also shown in the TM mode. A unique absorption property is achieved for both TE and TM modes of various banding angles, which offers a bendable property of the proposed absorber.
(a) (b)  Figures 11 and 12 show the e-field (V/m) and h-field (A/m) distribution of the designed MMA for the normal incident angle θ = 0° of the TE and TM modes, respectively; the absorption mechanism can also be analyzed from this electromagnetic behavior. The electromagnetic property is enhanced and agitated at a different area of the absorber.  Figures 11 and 12 show the e-field (V/m) and h-field (A/m) distribution of the designed MMA for the normal incident angle θ = 0 • of the TE and TM modes, respectively; the absorption mechanism can also be analyzed from this electromagnetic behavior. The electromagnetic property is enhanced and agitated at a different area of the absorber.

Electric and Magnetic Field
in the upper wavelength region with the increment in θ. The same absorption behavior is also shown in the TM mode. A unique absorption property is achieved for both TE and TM modes of various banding angles, which offers a bendable property of the proposed absorber.  Figures 11 and 12 show the e-field (V/m) and h-field (A/m) distribution of the designed MMA for the normal incident angle θ = 0° of the TE and TM modes, respectively; the absorption mechanism can also be analyzed from this electromagnetic behavior. The electromagnetic property is enhanced and agitated at a different area of the absorber.  The e-field and h-field properties of 380 nm, 521.83 nm, and 750 nm are presented. The metamaterial property also significantly influences the electromagnetic behavior. The relationship between metamaterial properties and electromagnetic behavior can be understood from Equations (9) and (10).

Electric and Magnetic Field
where D is the electric flux density, B is the magnetic flux density, eff ε = effective permittivity, 0 ε = permittivity of free space, eff μ = effective permeability, 0 μ = permeability of free space, and E and H are the electric and magnetic field intensities, respectively. In the TE mode, at a wavelength of 380 nm, the maximum e-field is assimilated in the bounding area of the patch, where a h-field maximum intensity appears in the metal region of the patch. A high e-field intensity is also visible in the nonmetallic area of the dielectric medium, and the high magnetic field appears in the left and right arms of the outer ring. The e-field intensity is reduced at 750 nm compared with 350 nm, which shows the side view of the absorber and causes a lower 92.2% absorption at 750 nm compared with 94% absorption at 350 nm. The peak of 99.99% absorption is achieved at 521.83 nm, the e-field intensity is less than those at 380 and 750 nm, but a high h-field is created in the dielectric substrate, shown in the side view at 521.83 nm wavelength, and leads to The e-field and h-field properties of 380 nm, 521.83 nm, and 750 nm are presented. The metamaterial property also significantly influences the electromagnetic behavior. The relationship between metamaterial properties and electromagnetic behavior can be understood from Equations (9) and (10).
where D is the electric flux density, B is the magnetic flux density, ε e f f = effective permittivity, ε 0 = permittivity of free space, µ e f f = effective permeability, µ 0 = permeability of free space, and E and H are the electric and magnetic field intensities, respectively. In the TE mode, at a wavelength of 380 nm, the maximum e-field is assimilated in the bounding area of the patch, where a h-field maximum intensity appears in the metal region of the patch.
A high e-field intensity is also visible in the nonmetallic area of the dielectric medium, and the high magnetic field appears in the left and right arms of the outer ring. The e-field intensity is reduced at 750 nm compared with 350 nm, which shows the side view of the absorber and causes a lower 92.2% absorption at 750 nm compared with 94% absorption at 350 nm. The peak of 99.99% absorption is achieved at 521.83 nm, the e-field intensity is less than those at 380 and 750 nm, but a high h-field is created in the dielectric substrate, shown in the side view at 521.83 nm wavelength, and leads to an increased absorption peak. Figure 12 shows the electromagnetic behavior of the TM mode, where the field distribution is similar to that of the TE mode, but it rotates vertically to horizontal in the TM mode. A comparison with the existing absorbers is presented in Table 3. The proposed SiO 2 and W-based metamaterial absorber achieves 70 • stability for an absorption above 70% with a wider operational bandwidth. An absorption above 92.2% and a peak absorption of 99.99% are also achieved. These excellent absorption results prioritize the proposed MMA over the listed absorber in Table 3.

Conclusions
This paper proposed a numerical analysis of a three-layer (W-SiO 2 -W) ultra-wideband MMA for visible optical wavelength application. The projected design evaluation and geometrical parameters were investigated for achieving near-unity absorption, which was also validated by Interference Theory (IT). Analysis of the PCR value proved that the proposed model is an absorber, not a polarization converter. The absorption calculation of TE and TM modes was simultaneously presented, showing uniform absorption and increasing the acceptability of the proposed MMA. A peak absorption of 99.99% was achieved at 521.83 nm, and an absorption above 92.2% and an average absorption of 97.1% were achieved in the entire visible operational wavelength region. Both TE and TM modes showed polarization insensitivity, wide oblique incident angle stability, and sound absorption of 18 • bending effects. These properties make the proposed MMA a potential candidate for various applications such as solar energy harvesting, light trapping, and optical sensors. Finally, a detailed comparison with existing MMAs showed the achieved unique features of the proposed MMA over the existing MMAs.

Data Availability Statement:
The data presented in this study are presented in this article.