--- title: "Design of 1D photonic crystals sustaining optical surface modes" authors: ["Valery Konopsky 1"] affiliation: "% Institute of Spectroscopy, Fizicheskaya, 5, Troitsk, Moscow, 108840, Russia; konopsky@isan.troitsk.ru" journal: "Coatings" year: 2022 volume: "" issue: "" article_number: "" pages: "" doi: "10.3390/coatings12101489" type: journal-article site_group: "" url_abstract: "" url_pdf: "https://valery.konopsky.com/kvnlocal/coatings-12-01489-v3.pdf" language: en source_tex: "Z:\\ValeryData\\Valery_New\\my_articles\\Coatings2022\\Coatings\\manuscript\\Konopsky2Coatings.tex" source_pdf: "Z:\\ValeryData\\Valery_New\\my_articles\\Coatings2022\\Coatings\\published\\coatings-12-01489-v3.pdf" --- ## Abstract An impedance approach has been implemented to design truncated 1D photonic crystals, sustaining optical surface modes, with any predetermined wavelength and wavevector. The implementation is realized as a free Windows program that calculates both the thicknesses of the double layers and the thickness of the final truncated layer at given refractive indices of the layers. The dispersion of the refractive indices can be given in the form of the Sellmeier/Drude formulas or in the form of a wavelength-n-k table. For mixed layers, the Maxwell Garnett theory can be used. This approach is suitable for studying and visualizing the field distribution inside photonic crystals, dispersion, and other aspects of the designed structures that sustain optical surface modes. Therefore, this program should promote scientific development and implementation of practical applications in this area. # Introduction Optical surface waves are excitations of electromagnetic modes that exist near the interface between two media. Photonic Crystal Surface Modes (PC SMs), which are also called ‘photonic band-gap surface modes’, ‘modes of (asymmetric) planar Bragg waveguide’, ‘surface waves in periodic layered medium’, ‘photonic crystal surface waves’, ‘optical Bloch surface waves’, and ‘surface waves in multilayer coating’ are modes that are bound to the external surface of an one-dimensional (1D) photonic crystal. ![Excitation of optical surface waves in a Kretschmann-like scheme.](media/Coatings2022/Fig1.eps) *Excitation of optical surface waves in a Kretschmann-like scheme.* Photonic crystals are materials that possess a periodic modulation of their refraction indices on the scale of the wavelength of light [41]. Such materials can exhibit photonic band gaps that are very much like the electronic band gaps for electron waves traveling in the periodic potential of the crystal. In both cases, frequency intervals exist in which wave propagation is forbidden. This analogy may be extended [30] to include surface levels, which can exist in band gaps of electronic crystals. In PCs, they correspond to optical surface modes with dispersion curves located inside the photonic band gap. The one-dimensional photonic crystal (1D PC) is a simple periodic multi-layer stack. Optical surface modes in 1D PCs were studied in the 1970s, both theoretically [3, 44] and experimentally [43]. Twenty years later, the excitation of optical surface modes in a Kretschmann-like configuration was demonstrated [37, 38]. A scheme of the Kretschmann-like excitation of PC SMs is presented in Fig. 1. In recent years, the PC SMs have been used in ever-widening applications in the fields of optical sensors [26, 14, 28, 16, 2, 32, 15], optical biosensors [25, 13, 29, 35, 20, 18, 39, 17, 33] and in other fields [10, 27, 12, 19, 6, 21, 31]. To excite the PC SM at any predetermined wavelength and at any predetermined wavevector (i.e., at any predetermined angle in the Kretschmann scheme), it is necessary to calculate the thickness of the last dielectric layer, which depends on the thicknesses of the double layers and their refractive indices (RIs), as well as RI of external environment. The thickness of the double layer should be pre-calculated in advance in order to maximize a photonic band gap at these predetermined wavelength and wavevector, since the confinement of the PC SM near the interface is the result of the photonic band gap on one side and the total internal reflection (TIR) on another side of the interface. Various approaches and strategies have been proposed to find the optimal parameters of the 1D PC structure and the thickness of the terminated layer [36, 9, 11, 7]. Here, to solve this problem, we present a free Windows program [1] based on an impedance approach, which was used in [22, 23]. This software permits both modeling and visualization of parameters of an 1D PC structure before submitting it for fabrication. # Materials and Methods ## Theoretical background: impedance approach All calculations in the program are performed using the impedance approach, which provides a unified theoretical description of multilayer systems for s- and p-polarizations by the same equations. With this approach, the ‘normal impedance’ $Z$ is the ratio of the tangential components of the electric field to the magnetic field, which has the following forms (in a $j$ th layer) for *s*- and *p*-polarizations [5, 8, 22]: ``` math \begin{eqnarray} Z_{s(j)}= &\displaystyle\frac{1}{n_j \cos(\theta_j)} &=\displaystyle\frac{1}{n_j\sqrt{1-(\rho/n_j)^2}} \[8pt] Z_{p(j)}= &\displaystyle\frac{\cos(\theta_j)}{n_j} &=\displaystyle\frac{\sqrt{1-(\rho/n_j)^2}}{n_j} \; . \end{eqnarray} ``` Here and hereafter the numerical aperture $\rho=n\sin(\theta)$ will be used as an angle variable in a planar multilayer system. It is a unified angle variable for all layers, since according to Snell’s law $\rho=n_j\sin(\theta_j)$, for any $j$. ## Reflection and transmission for a multilayer with N layers One of the advantages of using the impedance approach is the ability to represent the impedances of *several* layers as a *single* ‘apparent’ input impedance through a recursive relation. For example, if the multilayer is made up of $N$ plane-parallel, homogeneous, isotropic layers (with refractive indices $n_j$ and geometrical thicknesses $d_j$, where $j=1,2, \ldots, N$) between semi-infinite incident $_{(0)}$ and external $_{(e)}$ media (see Fig. 2), the apparent input impedance $Z^\mathrm{into}_{(j)}$ of a semi-infinite external medium $_{(e)}$ and layers from $N$ to $j$ may be calculated by the following *recursion relation* [22, 23, 5, 8]: ``` math \begin{equation} Z_{(j)}^\mathrm{into}=Z_{(j)}{\frac {{ Z^\mathrm{into}_{(j+1)}}-iZ_{(j)}\tan(\alpha_{j})}{Z_{(j)}-i{ Z^\mathrm{into}_{(j+1)}}\tan(\alpha_{j})}} \, , \end{equation} ``` where $\alpha_j=k_{z(j)}\,d_j=(2\pi/\lambda)n_j\cos(\theta_j)\,d_j$; $j=N, N-1,\ldots, 2,1$ and $Z_{(N+1)}^\mathrm{into}=Z_{(N+1)}=Z_{(e)}$, $n_{N+1}=n_{e}$ while $d_{N+1}=d_{e}=0$ by definition. This procedure should be continued recursively from the last layer (where, on the first iteration, $Z_{(e)}$ would be replaced by $Z_{(N)}^ \mathrm{into}$) – to the first layer until $Z_{ (1)}^ \mathrm{into}$ is obtained. The equation for reflection coefficients of *s*- or *p*-polarized waves from any complex multilayer (see Fig. 2) has a very simple form in the impedance terms: ``` math \begin{equation} R=\displaystyle\frac{ Z_{(1)}^\mathrm{into}-{Z}_{(0)}}{ Z_{(1)}^\mathrm{into}+{Z}_{(0)}} \, , \end{equation} ``` where $Z_{(1)}^\mathrm{into}$ is an apparent input impedance for a multilayer, i.e., it is the impedance that is seen by an incoming wave as it approaches to the interface. In the absence of the multilayer, it reduces to the standard Fresnel’s formula for *s*- or *p*-polarization, respectively. ![Reflection and transmission for a multilayer with N layers in terms of impedance.](media/Coatings2022/Fig2.eps) *Reflection and transmission for a multilayer with N layers in terms of impedance.* Fresnel’s formulas for multilayer transmission coefficients are as follows: ``` math \begin{eqnarray} &T_s=\prod\limits^{j=N}_{j=0} T_{s{j+1\choose j}}\, , \quad \mbox{with}\nonumber \[6pt] &T_{s{j+1\choose j}} = -\displaystyle{\frac {\left (Z_{s\,(j+1)}^\mathrm{into}+Z_{s\,(j+1)}\right )} {\left (Z_{s\,(j+1)}^\mathrm{into}+Z_{s\,(j)}\right )}} {e^{i\alpha_{j+1}}} \end{eqnarray} ``` and: ``` math \begin{eqnarray} &T_p=\prod\limits^{j=N}_{j=0} T_{p\,{j+1\choose j}} \, , \quad \mbox{with}\nonumber \[6pt] &T_{p\,{j+1\choose j}}=-\displaystyle{\frac {n_{j} Z_{p\,(j)}}{n_{j+1} {Z_{p\,{(j+1)}}}}} \displaystyle{\frac {\left( Z_{p\,(j+1)}^\mathrm{into} + Z_{p\,(j+1)}\right)}{\left (Z_{p\,(j+1)}^\mathrm{into}+Z_{p\,(j)}\right )}}{e^{i\alpha_{j+1}}} \quad \end{eqnarray} ``` where ${T}_{j+1\choose j}$ are transmission coefficients at an interface between the $j$ layer and the $j+1$ layer. ## Input impedance of a semi-infinite 1D PC A strategy for finding multilayer thicknesses and the thickness of the last, truncated layer that sustain the propagation of optical surface waves, is still a matter of discussion [9]. Our strategy in the first step is to choose a double layer thickness of the PC multilayer structure to maximize the bandgap extinction at a pre-selected wavelength and wavevector, as this ensures maximum confinement of the PC SM from the PC side of the interface. Then we solve the dispersion equation PC SM to find the thickness of the truncated (last) dielectric layer. Both operations employ the input impedances of a semi-infinite 1D PC and then apply the obtained thicknesses to a practical finite 1D PC. This approach shows an excellent match of the resonance parameters with a pre-selected wavelength and wavevector. The impedance of a semi-infinite 1D PC was derived in [22]: ``` math \begin{equation} Z^\mathrm{into}_{(PC)}=-\frac{i}{2} {\frac{\left( ( Z_{(2)}^{2}- Z_{(1)}^{2} )\tan (\alpha_{2})\tan(\alpha_{1})\pm\sqrt {{s}}\right)} {{ Z_{(2)}} \tan(\alpha_{1})+Z_{(1)}\tan(\alpha_{2})}} , \end{equation} ``` where ``` math \begin{eqnarray} s=-4 { Z_{(1)}} { Z_{(2)}} \left ({ Z_{(2)}} \tan({ \alpha_1})+{ Z_{(1)}} \tan({ \alpha_2})\right ) \left ({ Z_{(1)}} \tan({ \alpha_1})+{ Z_{(2)}}\tan({ \alpha_2})\right )+ \nonumber \\ \qquad \left [\left (Z_{(2)}^2-Z_{(1)}^{2} \right )\tan({ \alpha_1})\tan({ \alpha_2})\right ]^2 \nonumber \end{eqnarray} ``` and $\alpha_j=k_{z(j)}\,d_j=\pm(2\pi/\lambda)n_j d_j \sqrt{1-(\rho/n_j)^2}$. ## Band gap maximum extinction per length As a rule, in practical applications, we have values of two RIs ($n_{1}$ and $n_{2}$) of alternative media in the 1D PC, and the purpose is to find the thickness of each alternative layer, which provides the maximum extinction (or the maximum extinction per length) at given RIs, wavelength and angle. The first method for maximizing extinction at pre-selected wavelength and wavevector is to use ‘quarter-wave-length’ layer’s thicknesses for the double layer, which are: ``` math \begin{equation} d_j=\frac{\lambda}{4n_j \cos(\theta_j )} = \frac{\lambda}{4 \sqrt{ n^2_j - \rho^2 }} , \end{equation} ``` where $j=$ 1 or 2. The second method is more optimal, especially at large incident (grazing) angles (i.e. at $\cos(\theta_j)\rightarrow 0$), and it permits one to find the desired values of the thicknesses $d_1=d_{1\mathrm{max}}$ and $d_2=d_{2\mathrm{max}}$ which maximize expression (20) in [22] and provide the maximum extinction per length $d_1 + d_2$. In the first step of the program, a user may choose either of these two methods to determine the double layer thicknesses ($d_{1}$ and $d_{2}$). ## Dispersion relation for PC SM and its solution for the truncated layer thickness There are several programs and algorithms for calculating electromagnetic wave propagation through planar stratified media [45, 34]. In contrast to them, our program additionally calculates the thickness of the truncated layer $d_3$ making propagation of PC SM possible. The last (or penultimate - see below) layer with impedance $Z_{(3)}$ must be truncated to ensure the propagation of surface waves along the planar interface of stratified media at a given wavelength and wavevector. A general condition for the existence of a surface mode between two media with impedances $Z_{\mathrm{left}}$ and $Z_{\mathrm{right}}$ is: ``` math \begin{equation} Z_{\mathrm{left}}+Z_{\mathrm{right}}=0\, . \end{equation} ``` In our case, this condition takes the form: ``` math \begin{equation} Z_{(PC\& Z_{(3)})}^\mathrm{into}+Z_{(e)}=0\, . \end{equation} ``` The input impedance of a 1D PC structure plus the truncated layer, according to (2.3), will be: ``` math \begin{equation} Z_{(PC\& Z_{(3)})}^\mathrm{into}=Z_{(3)}{\frac {{ Z^\mathrm{into}_{(PC)}}-iZ_{(3)}\tan(\alpha_{3})}{Z_{(3)}-i{ Z^\mathrm{into}_{(PC)}}\tan(\alpha_{3})}} \, . \end{equation} ``` By solving equations (2.10) and (2.11), we obtain the dispersion relation for the optical surface waves in the 1D PC: ``` math \begin{equation} \alpha_3\equiv\left[ k_{z(3)} d_3 \right] = \pi {M}+\arctan\left({\frac{-i\left({ Z^\mathrm{into}_{(PC)}}+{ Z_{(e)}} \right){Z_{(3)}}}{Z_{(3)}^{2}+{ Z^\mathrm{into}_{(PC)}} { Z_{(e)}}}}\right) \, , \end{equation} ``` where $M$ is a whole number. If the last dielectric layer to be truncated is followed by a thin metal layer, sustaining long-range surface plasmons [24], then one can use equations (2.12), where the external medium impedance $Z_{(e)}$ is replaced by the input impedance of the external medium plus an impedance of the metal films with a pre-selected thickness $d_\mathrm{m}$, using (2.3): $Z^\mathrm{into}_{(e)\& (\mathrm{m})}$. In this case, the thickness of the last dielectric layer (or the penultimate layer, if we are counting the metal layer) can ensure the propagation of long-range surface plasmons along the external surface if the condition for minimizing the electromagnetic field inside the metal film is satisfied. In our program, this conditions (see (31) in [22]): ``` math \begin{equation} \rho_{1/2}= n_e + \frac{n_e^{3}}{2}\,\left[{\pi}{\frac {{d_\mathrm{m}}}{{\lambda}}}\right]^{2}\, , \end{equation} ``` can be selected by the appropriate checkbox. # Results ## Practical implementation in the program An implementation of the impedance approach described in this paper is available as a free Windows program at: . Both the program interface and numerical calculations were implemented in the C# programming language within .NET 6 framework. The program is distributed as a self-contained single file, which contains all components of the application, including the .NET libraries and target runtime libraries. The program is isolated from other .NET applications and does not use a locally installed shared runtime. The executable file `1DPC4all.exe` can be run on any 64-bit Windows above Windows 7. The user of the program is not required to download and install any versions of .NET. ## Refractive indices data The refractive indices data of the layers and their dispersion can be represented in the form of the Sellmeier formulas for dielectric layers and the Drude formulas for metal layers. Alternatively, the RI dispersion can be presented as a set of experimental data in the form of a wavelength-n-k table. Users can add their RI data to the existing ones in all of these forms. ![Step 1 program’s window for selection of double layer materials, wavelength and angle.](media/Coatings2022/Fig3.eps) *Step 1 program’s window for selection of double layer materials, wavelength and angle.* ### Sellmeier formula The seven coefficients for the Sellmeier formula may be written in an ASCII file with the extension `.slmr` in subfolder ‘\1DPC4all\Resources\calcRI\\’. These coefficients ($c_i$) will be substituted into the Sellmeier formula in the following form: ``` math \begin{equation} n^2= c_0 + \frac{c_1\lambda^2}{\lambda^2-c_2} + \frac{c_3\lambda^2}{\lambda^2-c_4} + \frac{c_5\lambda^2}{\lambda^2-c_6}\, , \end{equation} ``` where $\lambda$ is in $\mu m$. ### Drude formula The coefficients for the Drude formula may be also written in subfolder ‘\1DPC4all\Resources\calcRI\\’ as an ASCII file with the extension `.drd`. These nine coefficients ($c_i$) will be substituted into the Drude formula in the following form: ``` math \begin{equation} n^2 = c_0- \frac{c_1^2}{\omega^2+i \omega c_2 }+ \frac{c_3^2}{c_4^2-\omega^2-i \omega c_5 }+ \frac{c_6^2}{c_7^2-\omega^2-i \omega c_8 } \, , \end{equation} ``` where $\omega=1/\lambda$ in $1/\mu m$. In both cases, the user can insert $c_i =0$ for the trailing $i$ if the shorter form of the Sellmeier/Drude equation is used. ### Experimental n-k dataset Users can also represent the real and imaginary parts of RI as a tab/space-separated table, where each line has the form ‘${\lambda}$[Å] n k’. This table must be saved as an ASCII file with the extension `.nk` and added to the zip file ‘\1DPC4all\Resources\allRI.zip’, where 188 `.nk` files (representing various materials) are already stored. The values of n and k between the wavelengths presented in the table will be found by linear interpolation. ### Maxwell Garnett approximation for mixed layers The effective medium theory, in the form of the Maxwell Garnett approximation [4], can also be used for RI of mixed layers. The RI of a mixed layer is specified in an ASCII file with the extension `.gnt` in the subfolder ‘\1DPC4all\Resources\calcRI\\’. Each line in the file must represent the matrix medium and one (or two) inclusion(s) followed by their volume percentage. For example, a file with lines: ‘SiO2.nk 97 Au.drd 2 air.slmr 1 ’ will provide an effective RI for a layer consisting of 97% $SiO_2$ as a matrix material (given as `.nk` dataset ‘${\lambda}$ n k’), 2% Au inclusion (represented by the Drude formula) and 1% of air bubbles (represented by Sellmeier formula). Such file can also be automatically created from two or three selected layers in Step 2 of the program, where the volume percentage will be given by the thicknesses of the layers (hints in the program will provide all the details). Note that the Maxwell Garnett theory gives a good approximation when one of the materials (the matrix) prevails over the inclusions in terms of volume percentage. ![Step 2 program’s window for choosing the number of layers and final adjustment of the structure.](media/Coatings2022/Fig4.eps) *Step 2 program’s window for choosing the number of layers and final adjustment of the structure.* ## Step 1: selection of double layer materials, wavelength and angle The design of 1D PC structures in the program is divided into three simple steps. After starting the program, the window of Step 1 appears (see Fig. 3). In this window, the user must select the double layer materials ($n_1$ and $n_2$) and set the (central) wavelength and angle of incidence of the optical waves. If the angle of incidence is greater than the TIR angle ($\rho>n_0$), the program will also calculate the final layer thickness $d_3$ (with RI $n_3$) that sustains the surface optical modes. But if the user sets $\rho < n_0$ (no TIR), the program will not calculate the final layer (sustaining surface optical modes), and in the following steps it will only be possible to calculate reflection and transmission through 1D PC (without PC SM). After selecting the quarter-wavelength thickness (or the optimal thickness) checkbox, the button to Step 2 will be enabled, as indicated by the green dotted arrows in Fig. 3. The program will simultaneously calculate the thicknesses of the double layers $d_1$, $d_2$ and the possible thicknesses of the third layer, $d_3(M)$, which will sustain surface waves at a given wavelength and at a given angle. ![Reflection and transmission in Step 3 program’s window.](media/Coatings2022/Fig5.eps) *Reflection and transmission in Step 3 program’s window.* ## Step 2: choosing the number of layers and final adjustment of the structure In the Step 2 window shown in Fig. 4, the user must select the number of layers for the 1D PC structure. In the Kretschmann scheme, the optimal number of layers (for given losses in the structure) can be checked by the resonant dip of the reflection coefficient, which should drop to zero when the number is optimal. Also, at this step, any final adjustments to thicknesses, RIs, and materials of structure may be made. ![The spatial field distribution inside 1D PC.](media/Coatings2022/Fig6.eps) *The spatial field distribution inside 1D PC.* ## Step 3: presentation and analysis At Step 3 of the program, one can visualize the reflectance and transmission coefficients of the structure as a function of the wavelength $\lambda$ or the angle $\rho$ (see Fig. 5). It is also possible to choose any point of the angular dependence $T(\rho)$ to visualize the spatial distribution of the field (see Fig. 6) at the selected angle $\rho$. One can see in Fig. 6 that the electromagnetic fields reach their minimum in the center of the metal film (layer \#14), which minimizes total losses in the metal film when the condition (2.13) is satisfied. This condition for $\rho_{1/2}$ can be selected in the Step 1 by an appropriate checkbox. The dispersion of the structure can also be calculated and visualized in the form of transmission field enhancement near the interface $\log_{10}T(\lambda,\rho)$, as shown in Fig. 7. ## Additional features The program also facilitates the calculation of some special 1D PC and their features. ### 1D PC structures with two metal nanolayers It has recently been shown [23, 21], that it is possible to design multilayer systems, in which *two* metal nanolayers – one on each side of a thin dielectric film – sustain the long-range propagation of SPs (LRSPs). It this case, the correct thickness of the metal-bounded dielectric film is important, which ensures that the optical electric field has minima inside *both* metal nanolayers, and LRSP propagation will be possible in structures containing *two* metal nanolayers. In Step 2, the user can calculate and insert a dielectric film of the proper thickness, sandwiched between two metal layers, into a structure sustaining LRSP (if the checkbox for $\rho_{1/2}$ was preliminary selected in the Step 1). ![Optical dispersion in Step 3 program’s window.](media/Coatings2022/Fig7.eps) *Optical dispersion in Step 3 program’s window.* ### Luminescence from 1D PC structures One of the possible applications of such structures is the electroluminescence of an active (sandwiched) film in PC [21]. The program can calculate the dispersion of luminescence from layers that RI values in Step 2 marked with ‘+-’, ‘- -’ or ‘++’. These calculations are performed using the reciprocity theorem [42, 40]. In this case, it means that the integral optical electric field generated in the far zone by the dipoles located in the luminescence layer is the same as the electric field from the far zone dipoles (i.e. from plane optical waves) generated in this (marked) layer. # Discussion We have presented a software for calculations of one-dimensional PC structures. To the best of our knowledge, this is the first free program that can not only calculate the reflection and transmission of optical waves through a multilayer coating, but also calculate parameters for the excitation of surface optical waves propagating along the interface. The use of the impedance approach made it possible to develop software that calculates the thickness of the last (or penultimate) layer, as well as the thickness of the active layer (sandwiched between two metal nanolayers) that sustain surface waves. We hope that this program should facilitate scientific development and implementation of practical applications in this growing field. -0cm 999 Yablonovitch, E. Photonic band-gap structures. **1993**, *10*, 283–295. Kossel, D. Analogies between thin-film optics and electron band theory of solids. **1966**, *56*, 1434–1434. Arnaud, J.A.; Saleh, A.A.M. Guidance of surface waves by multilayer coatings. **1974**, *13*, 2343–2345. Yeh, P.; Yariv, A.; Hong, C.S. lectromagnetic propagation in periodic stratified media. I. General theory. **1977**, *67*, 423–438. Yeh, P.; Yariv, A.; Cho, A.Y. Optical surface waves in periodic layered media. **1978**, *32*, 104–105. Robertson, W.M.; May, M.S. Surface electromagnetic waves on one-dimensional photonic band gap arrays. **1999**, *74*, 1800–1802. Shinn, A.; Robertson, W. Surface plasmon-like sensor based on surface electromagnetic waves in a photonic band-gap material. **2005**, *105*, 360–364. Konopsky, V.N.; Alieva, E.V. Long-range plasmons in lossy metal films on photonic crystal surfaces. **2009**, *34*, 479–481. Hamidi, S.; Ramezani, R.; Bananej, A. Hydrogen gas sensor based on long-range surface plasmons in lossy palladium film placed on photonic crystal stack. **2016**, *53*, 201–208. Konopsky, V.N.; Basmanov, D.V.; Alieva, E.V.; Sekatskii, S.K.; Dietler, G. Size-dependent hydrogen uptake behavior of Pd nanoparticles revealed by photonic crystal surface waves. **2012**, *100*, 083 108. Ignatyeva, D.O.; Knyazev, G.A.; Kapralov, P.O.; Dietler, G.; Sekatskii, S.K.; Belotelov, V.I. Magneto-optical plasmonic heterostructure with ultranarrow resonance for sensing applications. **2016**, *6*, 28077. Alieva, E.; Konopsky, V.; Basmanov, D.; Sekatskii, S.; Dietler, G. Blue surface plasmon propagation along thin gold film–gas interface and its use for sensitive nitrogen dioxide detection. **2013**, *309*, 148–152. Li, J.; Tang, T.; Zhang, Y.; Luo, L.; Sun, P. Magneto-plasmonic sensor with one dimensional photonic crystal for methane detection. **2018**, *155*, 74–80. Ignatyeva, D.; Kapralov, P.; Golovko, P.; Shilina, P.; Khramova, A.; Sekatskii, S.; Nur-E-Alam, M.; Alameh, K.; Vasiliev, M.; Kalish, A.; et al. Sensing of surface and bulk refractive index using magnetophotonic crystal with hybrid magneto-optical response. **2021**, *21*, 1984. Konopsky, V.N.; Alieva, E.V. Photonic crystal surface waves for optical biosensors. **2007**, *79*, 4729–4735. Guo, Y.; Ye, J.Y.; Divin, C.; Huang, B.; Thomas, T.P.; Baker, Jr., J.R.; Norris, T.B. Real-time biomolecular binding detection using a sensitive photonic crystal biosensor. **2010**, *82*, 5211–5218. Konopsky, V.N.; Karakouz, T.; Alieva, E.V.; Vicario, C.; Sekatskii, S.K.; Dietler, G. Photonic crystal biosensor based on optical surface waves. **2013**, *13*, 2566–2578. Rivolo, P.; Michelotti, F.; Frascella, F.; Digregorio, G.; Mandracci, P.; Dominici, L.; Giorgis, F.; Descrovi, E. Real time secondary antibody detection by means of silicon-based multilayers sustaining Bloch surface waves. **2012**, *161*, 1046–1052. Konopsky, V.; Mitko, T.; Aldarov, K.; Alieva, E.; Basmanov, D.; Moskalets, A.; Matveeva, A.; Morozova, O.; Klinov, D. Photonic crystal surface mode imaging for multiplexed and high-throughput label-free biosensing. **2020**, *168*, 112575. Khodami, M.; Hirbodvash, Z.; Krupin, O.; Wong, W.R.; Lisicka-Skrzek, E.; Northfield, H.; Hahn, C.; Berini, P. Fabrication of Bloch Long Range Surface Plasmon Waveguides Integrating Counter Electrodes and Microfluidic Channels for Multimodal Biosensing. **2021**, *30*, 686–695. Sizova, S.; Shakurov, R.; Mitko, T.; Shirshikov, F.; Solovyeva, D.; Konopsky, V.; Alieva, E.; Klinov, D.; Bespyatykh, J.; Basmanov, D. The Elaboration of Effective Coatings for Photonic Crystal Chips in Optical Biosensors. **2021**, *14*, 152. Kalas, B.; Ferencz, K.; Saftics, A.; Czigany, Z.; Fried, M.; Petrik, P. Bloch surface waves biosensing in the ultraviolet wavelength range – Bragg structure design for investigating protein adsorption by in situ Kretschmann-Raether ellipsometry. **2021**, *536*, 147869. Petrova, I.; Konopsky, V.; Nabiev, I.; Sukhanova, A. Label-Free Flow Multiplex Biosensing via Photonic Crystal Surface Mode Detection. **2019**, *9*, 8745. Delfan, A.; Liscidini, M.; Sipe, J.E. Surface enhanced Raman scattering in the presence of multilayer dielectric structures. **2012**, *29*, 1863–1874. Konopsky, V.N.; Alieva, E.V.; Alyatkin, S.Y.; Melnikov, A.A.; Chekalin, S.V.; Agranovich, V.M. Phase-matched third-harmonic generation via doubly resonant optical surface modes in 1D photonic crystals. **2016**, *5*, e16168. Fong, N.R.; Menotti, M.; Lisicka-Skrzek, E.; Northfield, H.; Olivieri, A.; Tait, N.; Liscidini, M.; Berini, P. Bloch long-range surface plasmon polaritons on metal stripe waveguides on a multilayer substrate. **2017**, *4*, 593–599. Konopsky, V. Long-range surface plasmon amplification with current injection on a one-dimensional photonic crystal surface. **2015**, *40*, 2261–2264. Degli-Eredi, I.; Sipe, J.; Vermeulen, N. -polarized graphene modes sustained by photonic crystal structures. **2015**, *40*, 2076–2079. Konopsky, V.; Prokhorov, V.; Lypenko, D.; Dmitriev, A.; Alieva, E.; Dietler, G.; Sekatskii, S. Electrical excitation of long-range surface plasmons in PC/OLED structure with two metal nanolayers. **2020**, *12*, 1–8. Kovalevich, T.; Belharet, D.; Robert, L.; Ulliac, G.; Kim, M.S.; Herzig, H.P.; Grosjean, T.; Bernal, M.P. Bloch surface waves at the telecommunication wavelength with lithium niobate as the top layer for integrated optics. **2019**, *58*, 1757–1762. Rizzo, R.; Danz, N.; Michelotti, F.; Maillart, E.; Anopchenko, A.; Wächter, C. Optimization of angularly resolved Bloch surface wave biosensors. **2014**, *22*, 23202–23214. Delfan, A.; Degli-Eredi, I.; Sipe, J. Long-range surface plasmons in multilayer structures. **2015**, *32*, 1615–1623. Fong, N.; Menotti, M.; Lisicka-Skrzek, E.; Northfield, H.; Olivieri, A.; Tait, N.; Liscidini, M.; Berini, P. Guided Bloch long-range surface plasmon polaritons. In Proceedings of the 2017 19th International Conference on Transparent Optical Networks (ICTON). IEEE, 2017, pp. 1–4. Degli-Eredi, I.; Sipe, J.; Vermeulen, N. Power-flow-based design strategy for Bloch surface wave biosensors. **2018**, *43*, 1095–1098. Konopsky, V., 2022. , Version 1.0.8134.15784. Konopsky, V.N. Plasmon-polariton waves in nanofilms on one-dimensional photonic crystal surfaces. **2010**, *12*, 093 006. Konopsky, V.N. Long-range surface plasmons on duplex metal nanolayers. **2020**, *39*, 100788. Brekhovskikh, L. ; Academic: New-York, 1980. Delano, E.; Pegis, R. Methods of syntesis for dielectric multilayer filters. In *Progress in Optics*; Wolf, E., Ed.; North-Holland: Amsterdam, 1969; Vol. VII, chapter 2, pp. 67–137 (see pp.77, 130). Yuffa, A.J.; Scales, J.A. Object-oriented electrodynamic S-matrix code with modern applications. **2012**, *231*, 4823–4835. Reiser, P. Calculation of lossy dielectric multilayer filter response. **2005**, *25*, 499–513. Konopsky, V.N.; Alieva, E.V. Long-range propagation of plasmon polaritons in a thin metal film on a one-dimensional photonic crystal surface. **2006**, *97*, 253 904. Bohren, C.; Huffman, D. ; Wiley & Sons: New-York, 1983. Yang, H.; Alexopoulos, N. Gain enhancement methods for printed circuit antennas through multiple superstrates. **1987**, *35*, 860–863. Wu, X.H.; Kishk, A.A.; Glisson, A.W. A transmission line method to compute the far-field radiation of arbitrarily directed Hertzian dipoles in a multilayer dielectric structure: theory and applications. **2006**, *54*, 2731–2741. # Short Biography of Authors **Valery Konopsky** received the Ph.D. and M.S. degrees from the Moscow Institute of Physics and Technology in 1996 and 1993, respectively. During 1990–1993 he worked at Lebedev Physical Institute in Moscow and then, during 1993–1996, he worked at the Institute of Spectroscopy in Troitsk. Since 1996 he is a staff member of the Institute of Spectroscopy and currently he is a Senior Scientific Researcher, headed a group that is active in optical investigation of surfaces and surface nanostructures. His present research interests include surface optical waves, optical chemical sensors and biosensors. ## References 1. %Author 1, T. The title of the cited article. em Journal Abbreviation bf 2008, em 10, 142--149. %% Reference 2 %. 2. %Author 2, L. The title of the cited contribution. In em The Book Title; Editor 1, F., Editor 2, A., Eds.; Publishing House: City, Country, 2007; pp. 32--58. %% Reference 3 %. 3. %Author 1, A.; Author 2, B. Book Title, 3rd ed.; Publisher: Publisher Location, Country, 2008; pp. 154--196. %% Reference 4 %. 4. %Author 1, A.B.; Author 2, C. Title of Unpublished Work. Abbreviated Journal Name year, phrase indicating stage of publication (submitted; accepted; in press). %% Reference 5 %. 5. %Author 1, A.B. (University, City, State, Country); Author 2, C. (Institute, City, State, Country). Personal communication, 2012. %% Reference 6 %. 6. %Author 1, A.B.; Author 2, C.D.; Author 3, E.F. Title of presentation. In Proceedings of the Name of the Conference, Location of Conference, Country, Date of Conference (Day Month Year); Abstract Number (optional), Pagination (optional). %% Reference 7 %. 7. %Author 1, A.B. Title of Thesis. Level of Thesis, Degree-Granting University, Location of University, Date of Completion. %% Reference 8 %. 8. %Title of Site. Available online: URL (accessed on Day Month Year). %thebibliography. --- Generated: 2026-08-30 from Z:\ValeryData\Valery_New\my_articles\Coatings2022\Coatings\manuscript\Konopsky2Coatings.tex, validated against Z:\ValeryData\Valery_New\my_articles\Coatings2022\Coatings\published\coatings-12-01489-v3.pdf.