1. Introduction
The periodicity of crystals and gratings exhibits universal scattering patterns, which are manifested by the well-known grating formula and grating’s spectral resolution [
1].
In general, the grating’s scattering is much more complicated and it depends on the grating’s material, profile, and grooving method [
1,
2,
3]. In particular, plasmon excitations in the metal grating are responsible for anomalous scattering [
4,
5,
6]. Moreover, it was found that even dielectric gratings exhibit anomalous scattering behavior [
7,
8]. These results were in clear contrast to the traditional belief that narrow dielectric gratings can have only a small effect on the incident light beam [
9]. It was shown that even when the dielectric grating has weak modulation, the reflection coefficient can be very high [
10,
11]. The modulation can be either in the refraction index [
12] or in the grating’s grooves [
7]. These results show that in principle the reflectance at specific incident angles can be arbitrarily close to 100%. It should be stressed that these angles (with 100% reflectance) have nothing to do with the grating formula.
This anomalous conduct is related to a specific case of Fano resonances [
8,
13,
14,
15]. Recently, it was shown that in theory perfect reflection (zero transmission) can occur for arbitrarily narrow dielectric grating layer [
16]. The layer thickness can be arbitrarily narrow and the grating modulation amplitude can be arbitrarily low and still the effect of 100% reflectance occurs.
This conduct resembles a similar effect in Quantum Mechanics, where the transmission via an oscillating barrier can be suppressed completely at certain incoming particles’ energies. Even though the quantum phenomenon was discovered after its optical equivalent, it was thoroughly investigated for both scattering scenarios [
17,
18,
19,
20,
21,
22,
23] as well as for transition between adjacent quantum wells [
24,
25,
26].
This effect occurs when the particle’s energy approximately equals the oscillations’ quanta. However, this is a relatively simple quantum system. When the system becomes more complex the dynamics are much less intuitive. In particular, it has been shown that if the oscillating potential is a well, which is located in an opaque barrier, then this system exhibits complex behavior. The current is extremely sensitive to the system’s parameter, it can either increase substantially (activation and elevation) or almost vanish (complete suppression) [
27,
28,
29,
30].
When the barrier is absent this complex behavior vanishes. The reason is that the system’s behavior is determined by the ratio between the oscillating potential’s amplitude and the wavenumber in its vicinity. If this ratio is imaginary no complex behavior appears.
The barrier’s presence makes sure that the wavenumber is approximately real, and since the potential amplitude is real then the ratio is approximately real as well.
However, if the barrier is absent then the ratio can still be real provided the potential amplitude is imaginary. In that case, some of this complex behavior reappears.
While the physical application of an imaginary barrier is not entirely clear, the optical counterpart, i.e., imaginary wavenumber and imaginary refractive index are very common.
These conclusions can be applied to optical equivalent systems where the wavefunction is replaced with the electromagnetic field, and the oscillating potential is replaced with a spatial dielectric (grating as in Ref. [
16] but for complex dielectric coefficient).
Then, similarly, the complex structure’s dynamics can reappear by either putting the grating in a low index of refraction environment or by replacing the dielectric grating with an amplifying/absorbing one. In what follows we will focus on the latter: the scattering pattern of an amplifying/absorbing grating is going to be investigated.
Modulated amplification can be done by changing gradually the pumping over the grating causing a gradual change in the population inversion.
However, it is much simpler to create a uniform amplifying layer and attach to it a modulated absorber. If both layers are considerably narrower than the wavelength, then the two-layer system can be regarded as a single layer with modulated amplification/absorption.
2. The system under study
The system under study is a plane wave’s scattering from a grating. This scenario is presented in
Figure 1. A TE-polarized plane-wave beam hits an active grating at an incident angle
. The grating is located on the
x=0 plane. The gratings variations occur in the z-direction, in which case the electromagnetic Maxwell’s wave equation can be written
where the electric field is polarized in the y-direction (TE polarization)
, and
satisfies
where the grating’s frequency is the reciprocal of the groove spacing
(the pitch), i.e.
In (2) we followed Ref. [
16] and chose a delta function to represent the narrow grating’s layer. As was explained in [
16] (and references therein) a delta function is an excellent approximation for layers, which are shorter than the electromagnetic field’s wavelength. To be more specific, a dielectric with the following properties
where is the narrow layer’s thickness, can be approximated with great accuracy by . Since in what follows we focus on the limit of extremely narrow layers, this is an excellent approximation. Moreover, the validity of this approximation is improved when the incident angle increases, since the x-component of the beam’s wavenumber decreases. We will see that it is easier to measure the effect in the large incident angle domain. When the incident angle approaches π/2 the approximation is accurate.
3. Exact Numerical Solution
Let
be the incident angle, therefore the wavenumber is
(see
Figure 1), and the electromagnetic field can be written in all space as
where
and
, i.e.,
and
. The boundary criterions on the grating requires
and using the dimensionless parameters
the problem reduces to the following difference equation
which can be rewritten in a more compact matrix form
where
The numerical solution of this set of equations is presented in
Figure 2.
Eq.(7) illustrates the fact, which was mentioned in the introduction, that the dynamics is governed by the ratio between oscillations amplitude and . When this ratio is imaginary (note that is the imaginary part of the oscillations’ amplitude), the spectrum cannot have a complex structure.
As can be seen from this figure, there are incident angles in which the transmission coefficient of the non-negative, i.e.
, modes vanish. A small change in the incoming angle (~1.4% in
Figure 2) reduces the positive modes’ transmission probability from maximum value to zero. On the other hand, the incoming angle has a small effect on the negative modes (see
Figure 2). Unlike the positive modes regime, where the differences between mode’s amplitudes is small, on the negative regime, the spectrum consists of ridges and valleys. This spectrum structure and behavior is similar to the quantum system of Refs. [
29].
In
Figure 3 the grating’s scattering pattern is presented for these two incident angles. When the incident angle is
then the e-m field is scattered only in the
angles. Nothing is scattered at angles, which exceed the incident angle
. Below the incident angle, the e-m field is scattered at only specific angles.
On the other hand when the incident angle is the e-m is scattered at all angles even beyond the incident angle . Below the incident angle the scattering pattern is similar – in this regime, the e-m field is scattered at only specific angles.
Since all the
modes vanish simultaneously, one can focus on the central mode
. In
Figure 4 the amplitude of the central mode is presented as a function of the incident angle. As can be seen, beyond a certain angle, the transmission coefficient vanishes for several angles.
When the modulation depth
varies, so does the number of incident angles, in which the non-negative modes vanish. As can be seen in
Figure 5, the larger the modulation depth, the larger the number of incident angles, in which the non-negative modes vanish.
4. Approximate Analytical Derivation
The numerical studies shows that this effect occurs when the grating constant is considerably larger than the incident beam’s wavelength. In this regime, Eq.(7) can be simplified to allow for analytical derivation. In this regime,
which follows that Eq.(7) can be simplified to
The solution of this difference equation can be solved with Bessel functions. Using Eq.(9.1.27) of Ref. [
31]
the solution of (11) can be written
where
,
,
are the Bessel functions of the first kind [
31] and the coefficients
A and
C can be calculated by matching the solution and substituting the solution in (12), i.e.,
The solution therefore reads
Now we have an approximate solution of the entire spatial spectrum, from which two main results can be derived: (A) the incident angles in which the scattered e-m is bounded to the lower angles, i.e., the scattered angles are lower than the incident angles. (B) The scattering angles below the incident one.
5. Calculating the Incident Angles for Which Scattering Is Bounded
From (12) and (14) if follows that suppression of high angles scattering occurs for
, i.e. for
and in which case all the
modes vanish as well.
Since this approximation was derived for the limit of
then both relations
and
hold. Therefore, using Eq.(9.3.3) of Ref. [
31],
where in this case
For a solution, in which the argument of the cosine function of (16) is small whilst still
, the parameter
must be accordingly small. Therefore, we seek a solution around
in which case, (15) and (16) yields
Solving (19) for
, we finally obtain the condition for high angles suppression
or equivalently, one can find the incident angles
for which high angle scattering does not occur
In
Figure 6 the approximation (20) is plotted on top of the exact numerical simulation (
Figure 5). As can be seen, this approximation shows high agreement with numerical results.
6. Scattering Angles
The scattered beam takes the following form
where
A is the constant (13).
To calculate the scattering angles, one can use the fact that
. In this case the summation can be approximated by an integral
Using the method of stationary phase approximation, this integral can be evaluated at
Note that (24) is the classical grating formula.
Therefore, the scattering amplitude at the scattering angle
is approximately
Again, since
one can use Eq.(9.3.3) of Ref. [
31], in which case
the angles which yields the maximum value of (25) must satisfy
Substituting (26) in (27) and solving the angle, while keeping in mind that
one finally get the grating formula
As can be seen in
Figure 7 the scattered angles
for
, and
shows high agreement with the approximation (28). It should be stressed that the grating formula (28) appears on top of the classical one (24).
7. Physical Realization
The realization of this grating requires simultaneous changes in both the real and imaginary components of the refractive index, however, they do not change in the same manner.
To simplify the model, the harmonic change occurs in the imaginary component of the dielectric constant
. Therefore, the real and imaginary parts of the refraction index are
These components are plotted in
Figure 8 for
. As can be seen, for this value, the real part of the refractive index changes between 1 and 1.7, while the imaginary part varies between -1.5 to +1.5. Then, if the layer thickness is
then
, which is sufficient to measure several angles which shows this effect (according to
Figure 5). It should be noted that if we seek the effect at high, i.e., grazing, angles, the layers width can be increased substantially, while keeping the narrow layer approximation still. In this case multiple angles should be detected. Moreover, Most of this work concentrated, for simplicity, on the
scenario, however, by varying
the effect will appear even for very weak
, in which case, it will be easier to measure the effect. Nevertheless, a full study of this effect for different values of
and different modulations is beyond the scope of this paper.
8. Summary
The scattering from an active grating layer was analyzed. In this model, the imaginary component of the dielectric coefficient varies harmonically. The grating’s layer width is shorter than the optical beam wavelength, while the grating’s wavelength is longer than the beam’s wavelength. A new grating formula is formulated, which is the profile on top of the classical grating formula. An approximate analytical expression for this formula was derived. Besides the scattering orders, a special scattering pattern appears for certain incident angles.
It was found that in some cases, the scattering angles cannot exceed the incident angle , i.e., for specific the scattered angles obey .
References
- Wilcox, C.H. Scattering Theory for Diffraction Gratings; Springer-Verlag: New York, NY, USA, 1984. [Google Scholar]
-
Gratings: Theory and Numeric Applications; Popov, E. (Ed.) Presses universitaires de Provence; PUP; 2012. [Google Scholar]
- Palmer, C.; Loewen, E. Diffraction Grating Handbook; Newport Corporation, 2005. [Google Scholar]
- Piper, J.R.; Fan, S. Total absorption in a Graphene Monolayer in the Optical Regime by Critical Coupling with a Photonic Crystal Guided Resonance. ACS Photonics 2014, 1, 347–353. [Google Scholar] [CrossRef]
- Hutley, M.C.; Maystre, D. Total absorption of light by a diffraction grating. Opt. Commun. 1976, 19, 431–436. [Google Scholar] [CrossRef]
- Jory, M.J.; Vukusic, P.S.; Sambles, J.R. Development of a prototype gas sensor using surface plasmon resonance on gratings. Sens. Actuators B 1994, 17, 203–209. [Google Scholar] [CrossRef]
- Mashev, L.; Popov, E. Zero order anomaly of dielectric coated gratings. Optics. Commun. 1985, 55, 377. [Google Scholar] [CrossRef]
- Popov, E.; Mashev, L.; Maystre, D. Theoretical study of the anomalies of coated dielectric gratings. Optica Acta 1986, 33, 607–619. [Google Scholar] [CrossRef]
- Burckhardt, C.B. Efficiency of a Dielectric Grating. J. Opt. Soc. Am. 1967, 57, 601. [Google Scholar] [CrossRef]
- Tibuleac, S.; Magnusson, R. Reflection and Transmission guided-mode resonance filters. J. Opt. Soc. Am. A 1997, 14, 1617. [Google Scholar] [CrossRef]
- Tibuleac, S.; Magnusson, R. Narrow-linewidth bandpass filters with diffractive thin-film layers. Opt. Lett. 2001, 26, 584. [Google Scholar] [CrossRef] [PubMed]
- Wang, S.S.; Magnusson, R.; Bagby, J.S.; Moharam, M.G. Guided-mode resonance in planar dielectric-layer diffraction gratings. J. Opt. Soc. Am. A 1990, 7, 1470. [Google Scholar] [CrossRef]
- Fano, U. The theory of anomalous diffraction gratings and of quasi-stationary waves on metallic surfaces (Sommerfield’s waves). J. Opt. Soc. Am. 1941, 31, 213–222. [Google Scholar] [CrossRef]
- Hessel, A.; Oliner, A.A. A new theory of Wood’s anomalies on optical gratings. Appl. Opt. 1965, 4, 1275–1297. [Google Scholar] [CrossRef]
- Fan, S.; Suh, W.; Joannopoulos, J.D. Temporal coupled-mode theory for Fano resonant mode in optical resonators. J. Opt. Soc. Am. A 2003, 20, 569–572. [Google Scholar] [CrossRef] [PubMed]
- Granot, E. Derivation of analytical expressions for anomalous reflection in the limit of zero thickness and weakly modulated dielectric grating. J. Opt. Soc. Am. A 2022, 39, 2205–2213. [Google Scholar] [CrossRef] [PubMed]
- Cota, E.; Jose, J.V.; Rojas, F. Electronic transmission through time-periodic oscillating barriers. Nanostructured Mater. 1993, 3, 349. [Google Scholar] [CrossRef]
- Bagwell, P.F.; Lake, R.K. Resonances in transmission through an oscillating barrier. Phys. Rev. B 1992, 46, 15329. [Google Scholar] [CrossRef] [PubMed]
- Li, W.; Reichl, L.E. Floquet scattering through a time-periodic potential. Phys. Rev. B 1999, 60, 15732. [Google Scholar] [CrossRef]
- Saraga, D.S.; Sassoli-de-Bianchi, M. On the One-Dimensional Scattering by Time-Periodic Potentials: General Theory and Application to Two Specific Models. Helv. Phys. Acta 1997, 70, 751. [Google Scholar]
- Martinez, D.F.; Reichl, L.E. Transmission Properties of the oscillating delta-function potential. Phys. Rev. B 2001, 64. [Google Scholar]
- Haavig, D.L.; Reifenberger, R. Dynamic transmission and reflection phenomena for a time-dependent rectangular potential. Phys. Rev. B 1982, 26, 6408. [Google Scholar] [CrossRef]
- Cai, W.; Hu, P.; Zheng, T.F.; Yudanin, B. Resonance of the one-dimensional electron transmission above a quantum well with dissipation. Phys. Rev. B 1990, 41, 3513. [Google Scholar] [CrossRef]
- Kayanuma, Y. Role of phase coherence in the transition dynamics of a periodically driven two-level system. Phys. Rev. A 1994, 50, 843. [Google Scholar] [CrossRef] [PubMed]
- Kayanuma, Y. Stokes phase and geometrical phase in a driven two-level system. Phys. Rev. A 1997, 55, 2495. [Google Scholar] [CrossRef]
- del Campo, A.; Garcia-Calderon, G.; Muga, J.G. Quantum transients. Phys.Rep. 2009, 476, 1. [Google Scholar] [CrossRef]
- Azbel, M.Y. Eigenstate assisted activation. Phys. Rev. Lett 1992, 68, 98. [Google Scholar] [CrossRef] [PubMed]
- Zangwill, G.; Granot, E. Eigenstate suppressed activation. Phys. B: Condens. Matter 2015, 461, 140. [Google Scholar] [CrossRef]
- Granot, E. Selected elevation in quantum tunneling. Europhys. Lett. 2003, 61, 817. [Google Scholar] [CrossRef]
- Granot, E. The Tunneling Current through Oscillating Resonance and the Sisyphus Effect. Adv. Cond. Matter Phys. 2017, 2017, 2435857. [Google Scholar]
- Abramowitz, M.; Stegun, A. Handbook of Mathematical Functions; Dover Publications: New York, NY, USA, 1965. [Google Scholar]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).