INTRODUCTION
One-dimensional photonic crystals (1D-PCs) have been a hot topic in optical physics over the past two decades [1-5].  The 1D-PCs are a periodic structure created from distinct refractive indices of dielectrics. Refractive index periodicity leads to the formation of forbidden bands, called photonic band gaps (PBGs), controlling the propagation of light in specific wavelength ranges [6-8]. Also, the alternation of the forbidden band can be broken by entering a defect in PC [9], namely, the generation gap. Defects make the light pass through narrowband wavelengths in the forbidden band, which is known as the defect mode. In comparison with two-dimensional (2D) photonic crystals [10-12], one-dimensional (1D) photonic crystals are the most straightforward in their geometry and handling, allowing for controlling the properties of these structures and studying the influence of various physical parameters such as the layers order, number of periods, thickness of layers, defect layer thickness and the angle of incident [13-15]. Recently, there has been an increasing interest in 1D-PCs due to their capability to control and localize light [16-17], providing important consequences for designing waveguides, optical filters and omnidirectional mirrors [18-21]. An optical filter that operates as a Bragg mirror can pass a range of wavelengths of light, while also blocking all other wavelengths. The laser eye protection is one of the most important optical filters used [22-24]. There are multiple ways for designing such filters including the use of polymer coatings [25], absorptive glass [26] and 1D-PC [22-24]. Many theoretical and numerical methods have been developed to understand and design 1D-PC structures. Among them, the plane wave expansion method (PWEM) [27], finite-difference time-domain (FDTD) [28-29], and transfer-matrix method (TMM) [30-35] have become popular. 
In this paper, 1D-PCs consisting of [Si/SiO2]6 stacks on the Si3N4 substrate are designed, where SiO2 and Si layers act as the low and high refractive index media, respectively. This structure can induce constructive and destructive interference due to the difference in refractive index among the materials. The SiO2, Si and Si3N4 thin films are chosen in a way to provide minimum absorption in the visible wavelength region [36-37, 39]. Since the number of layers is selected in a long wavelength region (thus preventing the crystal from passing the light), the resulting filter has a low transmission, and is suitable to be used as an eye protection filter. A method of transferring matrix is then employed to determine the impact of defects (e.g., TiO2), functional wavelength, defect layer thickness, and angle of incident on the optimization of 1D-PCs. 

THEORETICAL FORMALISM
Calculating the reflectance and transmittance of the periodic multilayered structure is accomplished using the transfer matrix method. In this regard, the characteristic matrix corresponding to one stack is initially calculated and then distributed to a PC. The incident beam on the interface is split into two waves: a transmitted wave proceeding into the second layer and a reflected wave propagating back into the first layer (Fig. 1).

Where Ei,j-1 is the incident wave, Er,j-1is the reflected wave from the dielectric of j-1th layer, Et,j-1is the transmitted wave, and Er,j is the reflected wave from the dielectric of jth layer. The parameters a and c are amplitude of the incident and transmitted components, respectively. The parameters b and d indicate amplitude of the reflected components. Moreover, Kj is the wavenumber of the jth layer Kj=2πnj/?. According to the boundary conditions (Eq. 3), when the free charge density or current density in the dielectric is 0, the tangential components of the resultant electric field are continuous across the interface as given by Eq. (4):

Similarly, all tangential magnetic field components are also continuous, so that one obtains the following equations from |B|= n/c |E|:

                               
These formulations relate the field components of the first boundary to those at the next boundary, and can be written in the matrix form as follows:

Accordingly, the transfer matrix of the jth layer, Aj, is expressed as follows for TE and TM modes by Eqs. (9) and (10), respectively:

The effect of propagation through a layer of index nj, and thickness Lj is captured by the propagation matrix:

Each layer of a multilayer has its own transfer matrix, and the total transfer matrix is the product of individual transfer matrices, which can be deduced as follows:

where air is denoted as 0, and the substrate is S. When the electromagnetic field propagates, it will produce a different phase φ between the reflected or transmitted waves, leading to interference effects, while also forming the PBG or TP. The reflectance and transmittance coefficients of the PC are provided with the elements of the transfer matrix as follows:

Finally, the reflection and transmission parameters of a 1D-PC can be calculated using the following formula [30-35]:

Fig. 2 shows the schematic of 1D-photonic crystal structure with Si/SiO2 as stacks and TiO2 as the defect. Here, Si and SiO2 denote the high index and low index materials, respectively. Also, the number of the repetition for the basic structure is 6. The thicknesses of Si and SiO2 thin films are L1 and L2, respectively.

RESULTS AND DISCUSSION
A 1D-PC consisting of [Si/SiO2]m/Si3N4 stacks is simulated where Si3N4 is the substrate, and Si and SiO2 are thin films with high and low refractive indexes, respectively. At the start of the simulation, a functional wavelength (FW) is selected to be 511 nm (in conformity with Copper Vapor laser) and layer thicknesses are then calculated. The wavelength is changed from 400 to 1000 nm in order to compute the transmission spectrum. Complex refractive indexes of thin films as a function of wavelength are extracted based on the literature [36-39] and are illustrated in Fig. 3. The thicknesses of all layers are obtained from L= λF/(4nj), where λF is the FW, being initially set to 511 nm. Since the refractive indexes of Si and TiO2 are complex, their real parts are only taken into account for selecting the thicknesses. Therefore, thicknesses of Si3N4, Si, SiO2 and TiO2 thin films are selected to be 64.21, 37.37, 84.54 and 45.99 nm, respectively. 
Fig. 4(a) displays the transmittance spectrum of the [Si/SiO2]m/Si3N4 structure. In this notation, the bracket demonstrates the layer pairs (the stack) and m denotes the number of stacks. The increase in the period number m results in the formation and growth of the width of the photonic forbidden band or PBG, blocking the transmission of photons or light energy in certain frequencies (transmittance less than 2). Fig. 4(b) represents the transmission spectra of [Si/SiO2]6/Si3N4 and [SiO2/Si]6 /Si3N4 1D-PC structures, indicating the importance and effect of the arrangement of Si and SiO2 layers in the stacks. Since the only difference between the two structures is in the order of Si and SiO2 thin films, obtaining similar transmission spectra is predictable and reasonable. Therefore, [Si/SiO2]6/Si3N4 structure is chosen arbitrarily. 
In the following, TiO2 layer is inserted as a defect layer, and sandwiched between three stacks of [Si/SiO2], resulting in [Si/SiO2]3/TiO2/[Si/SiO2]3/Si3N4 structure. It is known that inserting a defect layer in a 1D-PC leads to the creation of a defect mode in its PBG [40-41]. Fig. 5 illustrates transmission spectra of the [Si/SiO2]3/TiO2/[Si/SiO2]3/Si3N4 structure, which is calculated for FWs of 461 nm (Rhodamine 6G laser), 511 nm (Copper Vapor laser), 594 nm (Red He-Ne laser), 612 nm (Yellow He-Ne laser), 711 nm (NIR Laser Diodes), 811 nm (NIR Laser Diodes), and 911 nm (NIR Laser Diodes). 
Fig. 5(a) indicates the occurrence of one transmission peak (TP) for the different FWs. Also, it is found that the TP shifts to longer wavelengths (i.e., a red shift) with increasing FW value. In order to investigate the effect of higher FWs on TP shift, transmission spectrum of [Si/SiO2]3/TiO2/[Si/SiO2]3/Si3N4 PC is calculated for FWs of 711, 811 and 911 nm, and the results are shown in Fig. 5(b). Note that, all layer thicknesses are calculated from the following formula: L= λF/(4n), positioning the TP location exactly on the selected FW. However, in addition to a TP at 711 nm, there is a transmission region for FW = 711 nm whose average transmission is up to 80%, according to Fig. 5(b). For FW = 911 nm, there is also a TP at 911 nm and a transmission region at the wavelength ranges between 450-684 nm, as is shown with the dotted arrow. As a result, one can find that the location of the TP, the number of TPs, the PBG and the intensity of the transmitted light in Fig. 5(b) can be completely engineered and controlled. Table 1 shows that, as the FW increases, the maximum transmission and the PBGs width increase, whereas the width of the right PBG shows a noticeable change with the increase in FW. 
According to the results of Figs. 5(a) and 5(b), one can change the PBG position with changing the thicknesses of all layers. However, we aim to change the position of PBG in a simpler manner, due to the defect layer thickness. To this end, the defect layer thickness is only changed from 20 to 100 nm, and the results obtained are shown in Fig. 6.
As can be seen, when the FWs of all layer thicknesses are fixed at 511 nm, and the thickness of the defect layer varies from 20 to 100 nm, the location of TP changes and can be divided into the two following regions: thicknesses between 20-70 nm where the position of TP is located between wavelengths of 438-590 nm, and thicknesses between 71-100 nm where the two TPs are located in the wavelength ranges of 618-668 nm and 421-451 nm. In this way, we have to define two different equations to understand the relationship between the thickness of defect layer and the TP location. For thicknesses between 20-100 nm, the relationship between the selected thickness of defect layer and the location of TP can be calculated from the following equation #1: Y = 380.97 + 2.93 X. On the other hand, for thicknesses between 71-100 nm, the relationship between the selected thickness and the location of the second TP can be calculated from the following equation #2: Y = 300.67 + 1.5 X. Using equations #1 and #2, it is possible to set the TP location on the desired thickness. For example, by using equation #1 and selecting the thickness of the defect layer to be 70 nm, it is possible to locate the TP at 590 nm. Therefore, this structure allows for locating the TP on any desired thickness. It should be noted that this ability is just possible for the structures whose stacks are the same on both sides of the defect layer.
Next, we analyze the effects of the incident angle on performance of the PCs. With increasing the incident angle from 0° to 60°, the TP location shifts to regions of longer wavelengths, and the TP height decreases. According to the results shown in Fig. 7, the TP height of the TE mode is 0.76 at 0°, 0.56 at 30°, and 0.1 at 60°. As well, the TP height of the TM mode is 0.57 at 30°, and 0.23 at 60°. This demonstrates that the incident angle affects the TP height, so that the TP height is maximum when the incident angle θ is 0°. It is also interesting to note that there is another TP at TM mode from 41-60 degrees, decreasing the TP height with the increase in the angle of incidence. In this case, the dependence of TP location and TP height on the incident angle is shown in Fig. 8.
From Figs. 5(a) and 5(b), putting the two structures together may lead to an optical filter with two or three TPs, leading to a so-called multi-optical filter. This arises from a transmission gap and a forbidden band for transmission in each structure. To verify this idea, two consecutive structures ([Si/SiO2]3/TiO2/[Si/SiO2]3 structure with FW of 511 nm and [Si/SiO2]3/TiO2/[Si/SiO2]3/Si3N4 structure with FW of 911 nm) along with the transmission of a merged structure as a function of wavelength are taken into consideration, and the calculated results are presented in Fig. 9. Evidently, the transmission is found to be possible only for the two aforementioned wavelengths, resulting from the combination of the two different transmissions of the consecutive structures. Note that the transmission is about zero for the rest of the spectrum. Hence, since the transmission of the two structures is below unity (< 1) at each FW, it is reasonable to have an aggregate transmission of less than unity for their combined structure at each specific wavelength. In this case, the aggregate transmission for all TPs is less than 1. 
On the other hand, three TPs are observed in wavelengths of 461, 694 and 949 nm for the combination of FWs = 461 and 950 nm. Also, three TPs are observed in wavelengths of 511, 776 and 1063 nm for the combination of FWs = 511 and 1064 nm. From an experimental standpoint, one can design a consecutive structure of [Si/SiO2] stacks to obtain the results of our simulations, thus forming an optical filter with three TPs.

CONCLUSION
The transfer matrix method was used to simulate the transmission spectrum in the 1D-PC consisting of [Si/SiO2] stacks. The results showed that the forbidden band can be broken, creating a TP by inserting a TiO2 thin film as a defect. The simulation calculations indicated that the TP shifted toward longer wavelengths with increasing FW, positioning the TP location exactly on the selected FW. By fixing the FW of all layer stacks and changing thickness for the defect layer, the location of TP was shown to be divided into two regions, thus allowing for locating the TP at each desired thickness. Moreover, the relationship between selected thickness and the TP location was discussed. The incident angle affected the TP location, and the TP height. Finally, a filter with two or three TPs was obtained by putting two stacked structures together, being applicable for designing and fabricating optical filters.

CONFLICT OF INTEREST
The authors declare that there is no conflict of interests regarding the publication of this manuscript.