Strong Coupling, Hyperbolic Metamaterials and Optical Tamm States in Layered Dielectric-Plasmonic Media

Thin films of noble metals with thickness smaller than the wavelength of light constitute one of the most investigated structures in plasmonics. The fact that surface plasmon modes can be excited in these films by different ways and the simplicity of fabrication offer ideal conditions for applications in nanophotonics. The generation of optical modes in coupled Fabry-Pérot planar cavities and their migration to hyperbolic metamaterials is investigated. Coupled Fabry-Pérot cavities behave as simple coupled resonators. When the intra-cavity media have different refractive indices in two or more coupled cavities resonance anti-crossings arise. The application of this kind of strong coupling in sensing is foreseen. Beyond the cavity modes excited by propagating waves, also long range plasmonic guided modes can be excited using emitters or evanescent waves. A periodic structure made by multiple plasmonic films and dielectrica supports bulk plasmons, of large propagation constant and increasing field amplitude. The optical response of these structures approaches that of the hyperbolic metamaterial predicted by the effective medium theory. Light can propagate with full transmission in a structure made of a photonic crystal based on quarter wavelength layers and a second photonic crystal with an overlapping forbidden band, but presenting a non-trivial topological phase achieved by band inversion. This is due to excitation of optical Tamm states at the boundary between both crystals. The extension to multiple optical Tamm states using dielectric and plasmonic materials and the symmetries of the edge states is investigated.


INTRODUCTION
Since the discovery of surface plasmons in thin metal films excited by electron beams, by Ritchie (1957), the investigation of plasmonic resonances on layered structures based on noble metals and their applications has a steady increase of attention by the scientific community. This is largely due the geometrical simplicity of the one-dimensional structures and the properties of the nearfields generated by the excitation of surface plasmon modes. Over the last decades planar layered stacks of noble metal films combined with dielectrics have been used in applications as refractive index sensors (Homola, 2008), structures for fluorescence decay engineering (Amos and Barnes, 1997;Barnes, 1998;Liebermann and Knoll, 2000), hyperbolic metamaterials (Cortes et al., 2012;Poddubny et al., 2013;Guo et al., 2020), strong-coupling between emitters and plasmon modes (Bellessa et al., 2004;Törmä and Barnes, 2014;Chikkaraddy et al., 2016), and the one-dimensional and twodimensional photonic topological insulators (Su et al., 1979;Kaliteevski et al., 2007;Wang et al., 2008;Obana et al., 2019;Henriques et al., 2020).
Despite the intensive research already done, the motivation to study some of the optical properties of layered media and the wide range of applications of the plasmon resonances did not fade. Layered dielectric media play a very important role in optics, as in lasers, optical filters, and in non-linear optics. The insertion of plasmonic films into dielectric layered media not only adds mathematical complexity and increases the degrees of freedom for optical manipulation, but allows to confine and propagate optical fields to a level never achieved by photonic media. This degree of optical confinement is at the base of the modern applications of classical and quantum physics.
In this article several less investigated properties of surface plasmon modes in layered media are addressed. The main features of the near-fields and propagation properties are compared among different structures. Indeed, small changes in the layer thickness leads to completely distinct optical regimes, but some characteristics of the optical modes are preserved. This includes the application of planar cavity Fabry-Pérot resonators to generate strong coupling. The analysis of coupled surface plasmon modes in finite sized metal-dielectric stacks to obtain hyperbolic metamaterials and optical Tamm states follows.
To my knowledge the application of Fabry-Pérot cavities to demonstrate classical strong coupling between light modes and its potential on sensing applications was not yet explored. The strong coupling between small plasmonic cavities (Ameling and Giessen, 2012) and beyond the light line (evanescent waves) has been investigated (Menghrajani and Barnes, 2020).
The investigation of hyperbolic metamaterials based on periodic metal-dielectric bilayers is often based on the effective medium theory and on the Kronig-Penney model. In this article the plasmonic modes excited in finite layered media are analyzed without approximation models. The plasmonic modes exhibit rich symmetry properties, not possible to investigate using the theoretical models based on infinite layered media. On the other hand, the investigation of the optical Tamm states closely related to photonic topological insulators in one dimension has been limited to a single interface. In this article the effects arising for double edge states are discussed.
Samples were fabricated for the experimental characterization of Fabry-Pérot coupled optical cavities and the experimental demonstration of strong coupling between light modes. Samples for experimental verification of the coupled optical Tamm states were also fabricated and characterized.
The article is structured in the following way. A short description of the simulation methods used is given. The theory of plasmon modes in single and double metal-dielectric interfaces in thin films is summarized. In more complex structures made of multiple layers, where the solutions of Maxwell's equations hardly can be found analytically, the transfer matrix method and the scattering matrix method are used to obtain the mode dispersion and the field profiles. The fabrication techniques and the optical characterization is presented in the following section.
Three examples of layered structures and the application of the corresponding surface plasmon modes are discussed. In first place, the properties of Fabry-Pérot modes in planar cavities based on silver mirrors and their application in the generation of strong coupling between cavities is discussed. These resonances are due to cavity modes and not to guided plasmonic modes, but both can coexist. These structures offer a significant advantage in sensing over the classical plasmonic resonators because they provide optical resonances with a quality factor or the order of Q ∼ 100 and simultaneously, the silver surface can be isolated from an analyte.
When several identical Fabry-Pérot resonators are stacked together and the layer dimensions become much smaller than the wavelength another effect arises. Fabry-Pérot resonances cease to exist, but a new family of guided plasmon modes arises: bulk plasmons with large k-vectors. These modes are representative of the hyperbolic dispersion arising for an infinite number of bilayers. To understand this uniaxial dispersion the effective medium theory based on the homogenization of the dielectric function in two orthogonal directions has been used. However, this theory cannot be used to investigate single plasmon modes and their propagation constants. Finite length plasmonicdielectric structures have a number of modes identical to the number of interfaces and their properties can be analized using the transfer matrix method.
The third class of structures is based on the optical Tamm states, arising in coupled one-dimensional photonic crystals with a non-trivial Zak phase and between a photonic crystal and plasmonic film. The discussion addresses cases beyond the single interface.

THEORY AND SIMULATION METHODS OF PLASMONIC MODES
In this section the properties of plasmonic modes in one and two interfaces are summarized. The application of the transfer matrix method and the scattering matrix method in multiple layers in order to determine the dispersion relation is presented.

Surfaces Plasmons in Single Metal Layer
A p-polarized plane electromagnetic shinnes on a metal film, of thickness d ≫ λ and parallel to the XY-plane, with complex dielectric function ǫ m (ω) = ǫ m,r + ǫ m,i . The upper medium has permittivity ǫ 1 . The magnetic field near the interface is given by 1 where β = k r + ik i and electric field components are 1 The upper and lower media can be analyzed even more generally than dielectric and metal as in Yang et al. (1991).
The function f (z) is an exponential decaying function with form with k 0 = ω/c. The dispersion relation is obtained from the continuity of tangential components of the electric field Then the real part of β reads This solution is real and therefore, there is a surface wave propagating at the interface metal-dielectric. However, the kvector of the incoming light wave k 0 n 1 cos θ never equals k r , independently of the angle of incidence θ . That is, any plane wave cannot excite the surface plasmons at a single interface. Simultaneously, a s-polarized wave, where the electric field oscillates in the XY-plane, never excites surface plasmon in flat layers.
When the thickness of the metal film d decreases to values compared to its skin depth and the film is sandwiched between two dielectric media of constants ǫ 1 and ǫ 2 , there are two metaldielectric interfaces whose surface modes interact. In that case the dispersion relation is tanh(α 2 ) = − ǫ 2 α 2 (ǫ 1 α 3 + ǫ 3 α 1 ) ǫ 1 ǫ 3 α 2 2 + α 1 α 3 ǫ 2 2 (7) where α j = β 2 − k 2 0 ǫ j , with j = 1, 2, 3. When ǫ 1 = ǫ 3 the equation above has two solutions. One correspondends to an antisymmetric electric field distribution on both sides of the metal film (Smith et al., 2008), and the dispersion tanh The other has a symmetric electric field distribution tanh Only for p-polarized light these modes can be excited. For spolarized light there are no solutions. The imaginary part of surface wave propagation constant k i = ℑ(β) of the antisymmetric mode is munch smaller than that of the symmetric mode. For this reason the surface plasmon associated with an antisymmetric oscillation of the electric field in both sides of the metal is called long range surface plasmon (LRSPP). The symmetric mode, appearing at longer wavelength for fixed β, has larger damping due to larger k i value and is therefore called short range surface plasmon (SRSPP) (Sarid, 1981;Yang et al., 1991;Berini, 2009;Konopsky and Alieva, 2009). The the excitation of a single surface plasmon mode in a thin metal film is easily achieved using different dielectric media (asymmetrical structure dielectric-metal-dielectric, or DMD also called insulator-metal-insulator, or IMI) in both sides of the metal (Kretschmann and Raether, 1968;Otto, 1968;Raether, 1988). The excitation of both modes is more challenging, because the k-vector of the light is always smaller than k SPP . For that it is necessary to use evanescent wave excitation. The propagation length of a surface plasmon mode is defined by Raether (1988), A large imaginary part of β is reflected in a short propagation length. Long propagation length is desired in applications and plasmonic waveguides and long range optical interactions (Berini et al., 2007).

Matrix Formalisms for Layered Media
The calculation of the reflectance, transmittance, absorbance, and fields in the layered structure are based on the transfer matrix method (Yeh, 2005). A layered medium is made of N optical layers, each with refractive index n j and thickness d j (i = 1, 2, . . . , N), surrounded by semi-infinite media of refractive index n 0 and n s , respectively. The layered medium is illuminated with a plane electromagnetic wave of p-polarization (TM), or spolarization (TE). Using the amplitudes A j and B j in each layer, associated with the forward and the backward propagation waves we can use the transfer-matrix to obtain the amplitudes of the fields in each layer. The field (either the electric field component E y for the s-polarization, or the H y magnetic field component for the p-polarization) in the layer j, defined by the interval between z j and z j+1 , is described by a sum of plane waves traveling in the opposite direction of the z-axis. It reads with β = n 2 j k 2 0 − k 2 z,j . The relationship between amplitudes for two consecutive layers j and j + 1 reads for j = 1, 2, . . . , N. By iterating for the full layered structure where Frontiers in Nanotechnology | www.frontiersin.org and n j cos θ j −n j cos θ j for s-polarization cos θ j cos θ j n j −n j for p-polarization . (15) From the matrix M we can obtain the reflectance and transmittance using and T = n s cos θ s n 0 cos θ 0 respectively. An alternative way to obtain the fields in the layers and the measured quantities as the reflectance and transmittance is by using the scattering-matrix formalism, which has better numerical stability for evanescent waves (Whittaker and Culshaw, 1999;Defrance et al., 2016;Orfanidis, 2016). The scattering-matrix was introduced by Wheeler in the context of quantum physics (Wheeler, 1937), but has been also applied in solid state physics (Ko and Inkson, 1988), electrical engineering (Orfanidis, 2016), and optics (Whittaker and Culshaw, 1999;Yuffa and Scales, 2012;Defrance et al., 2016). The S-matrix relates the forward and backward field amplitudes in the layer l and in the layer l ′ in an entirely different way. The outgoing field amplitudes are determined from the incoming field amplitudes by The elements of the M-matrix and the S-matrix in an arbitrary layer p are related by Li (1996) or equivalently In the calculations done in this article two matrix-based codes were used. (a) An adapted version of the Python TMM code based on the T-matrix calculation (Byrnes, 2019) was used for the reflectance, transmittance and fields in the layered structure. (b) An adapted version on the MATLAB Moosh code was used to obtain the guided modes and field profiles of plasmonic-dielectric layered structures (Defrance et al., 2016). The fields inside the structure and the reflectance and transmittance spectra of the planar photonic crystals used in the generation of optical Tamm states were verified using COMSOL Multiphysics.

SAMPLES FABRICATION AND OPTICAL CHARACTERIZATION
All samples were fabricated by successive electron-beam evaporations of metal and dielectrics on glass cover slide substrates. The vacuum pressure before the evaporation was 10 −6 mbar. Samples fabricated for layered hyperbolic metamaterials were based on coatings of Al 2 O 3 and silver. Aluminum oxide has a major advantage over silicon dioxide when some of the layers are of noble metals. It offers much better inter-layer adhesion, which prevents cracking and detachment from the substrate.
In the fabrication of coupled Fabry-Pérot cavities and hyperbolic metamaterials layers of silver, Al 2 O 3 and MgF 2 were deposited on glass substrates. For the samples based on photonic crystals for the generation of optical Tamm states, other oxides were purchased from Kurt J. Lesker: TiO 2 , ZrO 2 , and Ta 2 O 5 . The purity of TiO 2 and ZrO 2 is of 99.9%. The purity of Ta 2 O 5 is 99.95%. The pieces for evaporation of Ag, MgF 2 , Al 2 O 3 , and SiO 2 purchased have higher purity, reaching 99.99%. The optical constants of these materials for the visible and nearinfrared were obtained from Johnson and Christy (1972) Wood and Nassau (1982) (ZrO 2 ), DeVore (1951) (using the Sellmeier equation of the ordinary ray for TiO 2 ), and Dodge (1984) (using the Sellmeier equation of the ordinary ray for MgF 2 ). MgF 2 in the crystalline solid form is one of the optical materials with the lowest refractive index (n ∼ 1.37 in the visible), although in thin films the refractive index reaches 1.4 (de Marcos et al., 2017). It presents transparency in a wide range of wavelengths between 150 nm and 7 µm and is often used in distributed Bragg reflectors and anti-reflection coatings. All the oxides have high transparency in the visible and NIR and are therefore suitable for photonic applications were low loss is desired. The experimental band gaps of SiO 2 , Al 2 O 3 , ZrO 2 , Ta 2 O 5 , and TiO 2 are 9, 8.8, 5.8, 4.4, and 3.5 eV, respectively (Robertson, 2004). For wavelengths in the UV and above 10 µm optical absorption emerges, as for Ta 2 O 5 (Franke et al., 2000;Bright et al., 2013). The effective refractive index of each coating depends on the evaporation conditions and small deviations from those values in the literature are expected.
The samples fabricated were inspected by optical microscope. For samples placed in the vertical direction to the evaporation source, the film thickness is in good agreement with the measured value by the quartz balance. For some films, a comparison was made by measuring the topography of a film stripe by AFM (atomic force microscopy). In substrates placed off the vertical direction the film thickness decreases with increasing angle between evaporation source and substrate measured from the vertical. Further experimental characterization of the dielectric films, could be made using ellipsometry. The profile of multilayers can be imaged by scanning electron microscopy, but was not done so far.
In order to achieve an optical wide bandgap in a planar photonic crystal based on two dielectric media with a lower refractive index n l and a higher refractive index n h a large difference in the refractive index of the both materials is required (Yeh, 2005). Thus, a natural choice is to employ TiO 2 and MgF 2 . However, the film growth conditions required for a good adhesion and constant stoichiometry could not be fulfilled with the evaporation system available. The substrates could not be heated and the temperature controlled, though this is a recommendation of the materials suppliers. Better films quality and adhesion properties was achieved using the pair ZrO 2 and SiO 2 . However, multiple PCs deposited in two or more steps lead to widespread cracking. The cracking is absent in a single PC up to 14 bilayers, but becomes common if the number of layers duplicate, or the thickness of the ZrO 2 increases. The best results were achieved using Ta 2 O 5 and SiO 2 . The adhesion of these films is very good, and almost no cracking was visible, despite to modest conditions of the evaporation. Namely, the recommended partial pressure of O 2 could not be used and no ion assisted deposition was employed in the electron-beam evaporation process. Despite that, the reflectance experiments show that these films can be prepared to obtain optical Tamm states and verify the theoretical predictions.
Binary layers of Ta 2 O 5 and SiO 2 have been used in the highly reflective optical coatings of the test masses of the LIGO experiment (Abbott et al., 2016;Pinard et al., 2016;Steinlechner et al., 2018). A advantageous property of Ta 2 O 5 is its low thermal noise and thermal stability, comparing with other oxides or large refractive index. The state-of-the-art of these anti-reflection coatings was discussed in Granata et al. (2020). All the materials referred above were selected because of their optical properties in the visible and NIR spectral range. For longer wavelengths other materials of larger refractive index and low absorption can be used.
The reflectance of the samples was measured using a WITec SNOM/Confocal microscope equipped with an Zeiss Achrostigmat (5x, NA = 0.12) for illumination and light collection. A halogen white light source was used for excitation and an Avantes AvaSpec-ULS-TEC spectrometer was used for the spectral characterization. The angle-resolved reflectance of the coupled Fabry-Pérot cavities presenting strong coupling were characterized using a collimated halogen light source (Ocean Optics HL-2000-FHSA) and an optical spectrometer (Ocean Optics S2000). The samples were mounted in a rotating stage. Measurements of the reflectance were done for angles between 10 • and 80 • to the vertical in steps of 2 • .

STRONG COUPLING IN PLASMONIC LAYERED MEDIA
Layers of a metal separated by a dielectric medium of refractive index n and thickness d ≥ λ/2 form a Fabry-Pérot (FP) resonator. Despite the penetration of field into the metal the cavity length can be approximated to the half wavelength. In a stack formed by 3 metal layers and two intercalated dielectric layers of the same thickness two resonances in the reflectance and transmittance spctrum arise, with a wavelength separation λ dependent on the thickness of the central metal layer. Thus, a stack of two Fabry-Pérot cavities coupled by a middle metal layer corresponds to the optical analog of the classical coupled oscillators. When the thickness of last metal layer largely exceeds the light penetration depth (i.e., 200 nm for gold and silver), almost total absorption is reached at the resonances of the single or double Fabry-Pérot resonators. The calculated and experimental Q-factor of the cavity resonances can reach values above 100. Larger values can be reached in the IR, where the metal absorption of gold is low.
The half-wavelength resonance is the first mode. Higher modes are excited for shorter wavelengths. The resonances at vertical incidence show very low sensitivity to changes in the refractive index of the first and last media. However, for an angle of incidence > 30 • the resonances start to have large sensitivity to a variation of the refractive index in the first medium.
When two FP cavities based on two different intra-cavity dielectric media are coupled, an anti-crossing in the dispersion of their resonances occurs. This corresponds to a exchange of energy between the cavities leading to a splitting in the resonance, as in the case of identical intra-cavity medium. However, this anticrossing can be tuned by varying the angle of incidence. This anticrossing of resonant modes corresponds to the classical analog of the strong coupling regime between a quantum emitter and an optical cavity (Bellessa et al., 2004;Novotny, 2010;Ameling and Giessen, 2012;Törmä and Barnes, 2014;Chikkaraddy et al., 2016).
The classical analog of the strong coupling can occur not only in cavity modes excited by propagating waves, rather between surface plasmons modes, when β > k 0 (Ameling and Giessen, 2012;Menghrajani and Barnes, 2020). The only requirement is the design of dispersion curves that would cross, if considered individually. The results presented is this article are however restricted to the case of excitation by propagating waves. The advantage of the strong coupling in coupled cavities, using two different intra-cavity dielectrics, e.g., Al 2 O 3 with n h ∼ 1.76 and MgF 2 with n l ∼ 1.38 is the large sensitivity to the refractive index of the first medium, for angles of incidence close to the angle where anti-crossing in the dispersion is reached.
Silver films offer the best plasmonic properties, namely long propagation length of surface plasmons, low optical absorption, but are not adequate in biological sensing applications because of the chemical activity of silver. Sensing variations of the refractive index of a dielectric medium in the Kretschmanm-Raether setup, requires either a solution, an inert solid material, or any material of low chemical activity in contact with the silver, or gold film (Homola, 2008). Coupled Fabry-Pérot cavities presenting anti-crossing can be sandwiched by protective layers of Al 2 O 3 and be used in refractive index sensing, even for chemically active substances. Figure 1A presents the dispersion of two coupled Fabry-Pérot cavities using silver films and Al 2 O 3 as a intra-cavity medium. In Figure 1B are presented experimental spectra obtained at normal incidence for single and coupled cavities, indicating moderate Q-factors and narrow resonance bandwidth in the visible.
In Figure 2 are presented calculated and experimental results of the dispersion for two an three coupled cavities. In the calculated results (Figures 2A,B) the first and the last medium are the same. In the experimental results the illumination if from air and the substrate is glass. In the three coupled cavities there are two anti-crossings. The effect of the variation of the refrative index of the first medium is illustrated in Figure 3. For angles above 30 • the lower curve deflects rapidly with an increase of the refractive index n 1 .

Guided Modes in Plasmonic-Dielectric Layered Media
The cavity modes discussed in the previous section are excited by propagating waves. Evanescent wave modes also exist for coupled Fabry-Pérot cavities. They propagate along any direction in the plane of the layers with constant β. For single and double metals layers forming IMI and IMIMI structures analytic and closed form solutions of the dispersion relation of surface plasmons exist (Economou, 1969;Raether, 1988). In structures with more layers of dispersive materials, the roots of the complex equations are difficult to obtain and other semi-analytical methods are more convenient.
The calculation of the dispersion of surface plasmon is done using two alternative methods. The first is based on the transfer matrix method, using a point dipole source located near the first dielectric-metal. The power density spectrum of the light radiated by the dipole is either transmitted through the stack, absorbed by the materials, or coupled into guided modes. The theory of the power loss of a dipole near a layered medium was developed by several researchers in the 70's and 80's of the last century (Chance et al., 1978;Sipe, 1981;Ford and Weber, 1984) and applied in problems of radiative decay engineering of fluorescent dyes near interfaces (Barnes, 1998(Barnes, , 1999Novotny and Hecht, 2012).
The dispersion of plasmonic modes is achieved by evaluation of the integrand functions of the power loss of a harmonic oscillating dipole very close to a metal surface. According to the semi-classical theory, the decay rate of a dipole in a homogeneous medium of dielectric permittivity ǫ 1 is given by γ = P/(hω), with where µ is the dipole moment. Near an interface separating the medium where the dipole is located from other semi-infinite medium, or layered media, the power loss is dependent on the reflection coefficients of these media and the distance to the dipole. There is a radiating contribution and a non-radiating contribution to the decay rate, depending on the k-vector of the radiating field. The radiating contribution to the power spectrum FIGURE 3 | Reflectance of coupled planar cavities presenting strong coupling with dependence on the refractive index of the first medium. At vertical incidence (θ = 0) the reflectance is independent of the refractive index of the first medium, but for angles larger then 30 • fast deflection of the lower curve arises.
is given by where k 1 is the k-vector of the medium where the dipole is located. The non-radiating contribution includes guided modes, namely surface plasmon modes associated with non-propagating k-vectors larger than k 1 . These modes can be calculated based on the Fresnel reflection coefficients from the boundary for the s-and p-polarization. The quantum efficiency of the radiation is given by the ratio The decay rate can be decomposed into a parallel (in-plane) and a perpendicular (out-of-plane) term. They affect differently the decay rate, because they rely on different contributions of the Fresnel reflection coefficients for p-and s-polarization. The density of power loss of an isotropic dipole, as it is expected from a dye molecule embedded in a dielectric transparent medium, is given by For a dipole near in a medium of dielectric permittivity ǫ 1 at a distance d to a second medium of permittivity ǫ 2 (or an arbitrary layered medium) the decay rate for the parallel and perpendicular components is given by Ford and Weber (1984) and Barnes (1998) and k 0 = ω/c, u = k 1 /k 0 and l i = −i ǫ j /ǫ 1 − u 2 1/2 with j = 1, 2. However, in order to obtain the dispersion relation these integrals have not to be evaluated. If the second medium is substituted by a layered medium of dielectric and plasmonic layers, the corresponding reflection coefficients are obtained from the transfer matrix calculation.
The second method to obtain the dispersion relation relies on the scattering matrix calculation. It is a semi-analytical calculation method. The main advantage of the scattering matrix comparing to the transfer matrix is its numerical stability, in particular when absorbing media are involved. The dispersion is found by solving the roots of the following equation, by the steepest descent method (Defrance et al., 2016) f (β, λ) = 0 (27) and f is determined using the scattering matrix for the layered structure. The MATLAB library MOOSH was employed for this purpose (Defrance et al., 2016). The dispersion of the IMIMI stack calculated using the scattering matrix formalism is presented in Figure 4. The density of power loss of an isotropic dipole localized 5 nm above the first silver layer is presented in Figure 5. Both dispersions relations are equivalent, but the normalized k-vector is different. Hence, the differences in the values of the modes in the absissa. In Figure 4 k x was normalized by the vacuum k 0 , whereas in Figure 5 k x the normalization constant is k 1 = k 0 n 1 . The multiple plasmon modes arise as result of a hybridization between the antisymmetric and the symmetric surface plasmon modes of a single film and the cavity mode between the metal layers. They form two branches: the lower energy branch includes both fundamental modes and a number of modes of larger β identical to the number of coupled cavities. For a single Fabry-Pérot cavity is one. At higher energy levels the second branch has two modes, the antisymmetric mode and a volume mode. Both branches converge to a flat dispersion for large propagation constant around 3 eV (see Figure 5). The modes splitting and saturation for long propagation constant is a feature shared by layered hyperbolic metamaterials (Cortes et al., 2012;, discussed in the next section. The upper branch corresponds to a Type I hyperbolic metamaterial and the lower branch of modes corresponds to Type II, respectively (Cortes et al., 2012;. The plasmon modes of large k-vector propagate inside the layered medium and are therefore called volume plasmons.
Samples were fabricated using bilayers of Al 2 O 3 and silver with thickness d ≪ λ. The experimental investigation of these samples, namely the properties of propagating plasmons and the fluorescence lifetime of emitters near the first silver layer can provide more information for their potential application in quantum optics experiments.

HYPERBOLIC METAMATERIALS
A hyperbolic metamaterial is a material which has a dispersion relation described by a hyperbolic equation . The isofrequency surfaces of isotropic materials have spherical shape. When a dielectric material is birefringent in a single axial direction it is called uniaxial and its dispersion equation has isofrequency surfaces of spheroidal shape, whereas biaxial materials have ellipsoidal isofrequency surfaces. The effective medium theory (EMT) (Bruggeman, 1935;Rytov, 1956) is a method of averaging the dielectric function according to the symmetry axes of the lattice. An early application of the EMT goes back to investigation of the colors of colloidal particles in glasses (Garnett, 1904). The homogenization of the EMT requires a unit cell size much smaller than the wavelength of light in the frequency range considered. The resulting dielectric function is calculated based on geometrical parameters and on the dielectric functions of the constituent materials and their lengths along the symmetry axes of the periodic structure. For a layered medium made of binary layers of materials with dielectric functions ǫ 1 (ω) and ǫ 2 (ω) and thicknesses d 1 and d 2 , respectively, the EMT dielectric function results into two components: one for any axis parallel to the layers and one for the axis perpendicular to the layers. Assuming that the layers are parallel to XY-plane, then ǫ xx = ǫ yy = ǫ and ǫ zz = ǫ ⊥ and using k 0 = ω/c the in-plane and the out-of-plane dielectric functions read The corresponding dispersion equation takes the form (Guo et al., 2020) ( The first term corresponds to the ordinary TE wave and the corresponding isofrequency surfaces (ω = const.) are spherical. The second term corresponds to the extraordinary TM waves. When ǫ > 0 and ǫ ⊥ < 0 both modes exist and the second term leads to a dispersion whose isofrequency have double sheet hyperboloid shape (Type I). For ǫ < 0 and ǫ ⊥ > 0 only the TM waves can propagate and isofrequency surfaces have single sheet hyperboloid shape (Type II). Type II hyperboloids have negative Gaussian curvature everywhere, whereas Type I hyperboloids have positive Gaussian curvature. A hyperbolic material of Type II dispersion requires one constituent layer allowing large anisotropy, for example a material of negative dispersion (plasmonic). By other hand, an array of dielectric rods embedded in a distinct dielectric medium achieve a EMT dielectric function made of two parts: a spheroidal dielectric function for small k-vectors and a Type I hyperboloidal function for large k-vectors . However, according to the EMT a transition between Type II and type I can be achieved, depending on the filling factor of the plasmonic material in the unit cell and on the wavelength range (Cortes et al., 2012).
Array of rods of dielectric function ǫ m (ω) embedded in a isotropic medium of dielectric function ǫ d form other example of hyperbolic metamaterial, where the EMT homogenization can also be applied. Rods of diameter r and lattice constant L arranged in a squared lattice and oriented along the z-axis have the following EMT dielectric functions : The filling factor ρ = a/A, is defined as the ratio between the area of the cross-section of the rods a and the area of the unit cell of the lattice A.
The effective dielectric function components of a hyperbolic material made of an infinite number of silver (10 nm thickness) and Al 2 O 3 (20 nm) layers are presented in Figure 6. The FIGURE 4 | Dispersion relation of the surface plasmon modes supported in a pair of silver films separated by a layer of aluminum oxide. The semi-infinite media embedding the three layers is also aluminum oxide. The dielectric constants for silver were taken from Johnson and Christy (1972). The dielectric function used for the aluminum oxide was based on the Sellmeier dispersion equation found by Malitson (1962). The tabulated optical constant values were extracted from the online database https://refractiveindex.info (Polyanskiy, 2020). The dark regions correspond to the plasmon guided modes. Plasmon modes have always a k-vector larger then k 1 = k 0 n 1 , the propagation constant of the first medium. The bright lines in the left region denote the cavity modes for p-polarized light.
parallel and perpendicular components were calculated using the effective medium theory equations.

Bulk Plasmon Modes in Finite Hyperbolic Stacks
In order to avoid the limitations of the effective medium theory (Kidwai et al., 2012) and the Kronig-Penney model for infinite number of layers (Li and Khurgin, 2016) the calculations of the field profiles in layered media of finite length were based on the scattering matrix method. Magnetic field profiles of surface modes and the bulk plasmon modes for two finite stacks of dielectric and metal layers are presented in Figure 7 and Figure 8. Both structures have bilayers of the same thickness. The difference is the number silver films: 5 in Figure 7 and 6 in Figure 8. The wavelengths at which the profiles were calculated are in the range of the Type II hyperbolic modes branch.
The following properties can be inferred from the profiles of the H y field component: (a) the first mode with a symmetric profile of the magnetic field along the full structure has the smallest k-vector (normalized to k 0 ). It corresponds to the antisymmetric electric field mode. (b) The second mode with a slightly larger propagation constant has an antisymmetric magnetic field distribution, considering the symmetry plane of the full structure. (c) The volume plasmon modes have successive increasing values of the propagation constants and simultaneously increasing amplitudes inside the layers. (d) The imaginary part of the propagation constant ℑ[β] = k i increases with the order of the mode and the respective real part k r . This means that the propagation length of the bulk plasmons decreases with an increase of the propagation constant. Therefore, an arbitrary large number of layers and modes may be counterproductive in applications as imaging (Li et al., 2017) and long-range interactions (Biehs et al., 2016;Newman et al., 2018). Due to the large field strength inside some layers and the wellseparated propagation constants of the volume modes it would be of high relevance the selective excitation of single plasmon modes.

PHOTONIC TOPOLOGICAL MATERIALS IN ONE-DIMENSION
The emergence of the topological order in condensed matter physics has revolutionize materials science, namely in its solid state and optical domains. Starting with the fractional quantum Hall effect and the quantum spin Hall effect, where the electron current is replaced by a spin current, an increasing number of effects has been discovered to posses topological properties as the topological insulators (Hasan and Kane, 2010). Moreover, the concepts of topological phase transition and topological insulator have been also extended into photonics (Ozawa et al., 2019).
It is also noteworthy that an early discovery of electron states in truncated crystals, the edge states, are now investigated from the perspective of topological insulators, giving them a much large degree of generality. The impact of topological materials and topological insulators has achieved a large scale attraction, from effects arising in classical waves to states of quantum matter (Hasan and Kane, 2010;Haldane, 2017). Topology properties are more general than geometry properties as symmetry. A topological space preserves its properties by a continuous deformation, whereas an arbitrary space deformation often means a break of symmetry. Thus, properties of physical systems with topological invariance are robust against small deformations and defects.
One of the most relevant examples of topological properties is the Berry phase (Berry, 1984). It generalizes the concept of geometric phase introduced by Pancharatnam (1956) and it is related to the parallel transport in topological spaces. The close relationship of the Berry phase, arising in a cyclic path of a quantum system described by a Hamiltonian under an adiabatic evolution, with topological properties as the Chern number and other topological invariants arising in condensed matter physics was due to Simon (1983). It is remarkable that the publication date of the article of B. Simon preceded that of the article of M. Berry and coined the name of this phase. Reviews of phenomena in condensed matter physics manifesting Berry phase are available (Resta, 2000;Xiao et al., 2010).
Zak investigated the application to the Berry phase is solid state materials characterized by a periodic lattice in real space and by the corresponding Brillouin zone (BZ) in the momentum space (Zak, 1989). Due to the periodicity of BZ the evolution of FIGURE 6 | Dielectric functions of a plasmonic-dielectric layered medium calculated using the effective medium theory. Each bilayer is composed of silver (10 nm) and aluminum oxide (20 nm). The permittivity of silver was obtained from the optical constants from Johnson and Christy (1972). The spectral regions labeled with Type I and Type II are the regions where the real part of perpendicular and parallel parts of the dielectric function become negative, respectively.
a Bloch wave along a band of the BZ is equivalent to a closed loop. Zak also found that bands of a periodic lattice, in particular in one-dimension, manifest certain symmetry properties when crossing. These properties are required for the generation of edge states (Zak, 1984(Zak, , 1985. The most well-known examples of edge states of electrons in solids are the Tamm states (Tamm, 1932a,b) and Schockley states (Shockley, 1939). Both are variants of the same kind of edges states (Vinogradov et al., 2010). The advent of photonic crystals brought out the investigation of the optical analogues of the edge states known for electrons in solids (Stȩślicka et al., 1990;Tikhodeev, 1991;Kavokin et al., 2005;Kaliteevski et al., 2007;Vinogradov et al., 2010) as well as in two-dimensional materials (Nakada et al., 1996;Delplace et al., 2011). Also related to formation of edge states in one-dimension is the model of Su-Schrieffer-Heeger describing the excitation of solitons in polyacetylene (Su et al., 1979) and the Peierls transition. The electron edges states of one-dimensional atom lattices have an analog in photonics, substituting the periodic potential by a layered structure with periodic refractive index.
Topological invariants of the photonic band structure in photonic crystals and in plasmonic materials coupled with photonic crystals, as the Zak phase play an important role in the edge state properties (Xiao et al., 2014Gao et al., 2015;Henriques et al., 2020). Moreover, a requirement for the generation of edge states in crystals is the band inversion (Zak, 1989;Xiao et al., 2014Xiao et al., , 2015Esmann et al., 2018;Wang et al., 2019). Despite the fact that the edge states of photonic crystals are only analogs of the edge electronic states discovered by Tamm, for historical reasons it is reasonable the call edge states arising in truncated photonic crystals optical Tamm states (OTS).
Presently many publications are dedicated to the study of edge states and the associated topological properties in two-and three-dimensional photonic crystals and in arrays of plasmonic particles (Wang et al., 2008(Wang et al., , 2016van Miert et al., 2016;Proctor et al., 2020). In this article only one-dimensional edge states and their properties are discussed.

Design of Single and Multiple Optical Tamm States
Up to date most of the optical Tamm states investigated are based either on structures made of two planar photonic crystals (Kavokin et al., 2005;Vinogradov et al., 2010;Xiao et al., 2014;Esmann et al., 2018), where one of them presents band inversion, or structures based on a photonic crystal and plasmonic materials (Durach and Rusina, 2012;Henriques et al., 2020), or an excitonic material (Núñez-Sánchez et al., 2016). Moreover, the analogy between photonic layered media and electronic lattices has been also extended to systems supporting mechanical waves Yang et al., 2016;Ma et al., 2019). The design of edge states arising in stacks of photonic crystals relies on the definition of a bulk bandgap, in an infinite photonic crystal (PC1), and by establishing a correspondence between the bulk and the edge state. This requires that the second photonic crystal (PC2) should have an overlapping bandgap with the first after operating a band inversion. This excludes the second photonic crystal to be identical to the first. In photonic crystals with binary layers the band inversion can be achieved by changing either the values of refractive indices of the materials of the layers, or by changing their thicknesses. Using the commercial available oxides with transparency in the visible and NIR spectral regions, it is more practical to change the thickness of each layer. Semiconductors of different materials can also be used either in phononic Tamm states (Esmann et al., 2018), or in electronic Tamm states in superlattices (Steslicka, 1995). The band inversion of one-dimensional photonic crystals is the same kind of band transformation in two-dimensional crystals, in order to obtain topological insulators (Parappurath et al., 2020).
In a photonic crystal of indefinite length the bulk Bloch modes (bands) for illumination at vertical incidence (θ = 0) are found from the dispersion relation of the layered medium, which reads (Yeh, 2005) cos(K B ) = cos(k 1 a) cos(k 2 b)− 1 2 n 2 n 1 + n 1 n 2 sin(k 1 a) sin(k 2 b), with k 1 = (ω/c)n 1 and k 2 = (ω/c)n 2 . The thicknesses of the layers are d 1 = a and d 2 = b and the refractive indices n 1 and n 2 , respectively. K B is the Bloch k-vector and = a + b. A similar equation arises in the dispersion of the Kronig-Penney model of electrons in solids (Kronig and Penney, 1931). The same method was used by I. Tamm to find the edge states of electrons in a onedimensional lattice, but substituting the propagation constants k 1 and k 2 by constants related to the periodic potential function.
Forbidden gaps in the energy occur for ℜ[K B ] = nπ/ , with n = 1, 2, . . . . In the forbidden bands K B = π ± ix, where x is a function of the refractive indices n 1 and n 2 (Yeh, 2005). Bandgaps centered in the wavelength λ = λ g arise when the following condition is verified k 0 n 1 d 1 = k 0 n 2 d 2 = jπ/2, with j = 1, 2, 3, . . . or equivalently, The band edges are also dependent on n 1 and n 2 . They determine the width of each bandgap. By varying a and b opening or closing of bandgaps can be achieved. When two semi-infinite photonic crystals are stacked together and truncated at the same inversion point of the unit cell of the lattice, with two forbidden bands overlap, but the Zak phases of the respective edge bands reverses an edge state arises. This permits total transmission at a very narrow bandwidth inside the bandgap. The tuning of the parameters of the gaps and the respective Zak phases was discussed in detail in Xiao et al. (2014). I should be noted that the definition of the unit cell centered around an inversion point is necessary for a well-defined value of the Zak phase (0 oi π), but it is not a requirement for the generation of an edge state.
Using the treatment of optical impedance leads to the conditions necessary to achieve full transmission inside the gap. This is Z 1 + Z 2 = 0 where Z 1 and Z 2 are the input optical impedances of the first and the second semiinfinite crystals, respectively. Konopsky found analytic solutions for input impedance for s-and p-polarized light combining photonic crystals and dispersive media (Konopsky, 2010). The formulas obtained permit application in quite general situations accounting for angular dependence of the incoming wave. The estimation of optimal conditions of excitation on LRSPP in planar layered media also benefit from this calculation method (Konopsky and Alieva, 2009;Delfan et al., 2015).
Despite the simplicity of the impedance condition, the practical design of optical Tamm states for an arbitrary forbidden band centered around λ g is more cumbersome. In practice we deal with finite crystals where the Kronig-Penney model cannot be applied. A practical way to tailor the thickness of the layers for which the edge states could be found is based on the phase of the reflection coefficient at vertical incidence (Kavokin et al., 2005;Xiao et al., 2014;Gao et al., 2015). For that we calculate the total reflection coefficient of a finite binary photonic crystal using the transfer matrix method. The phase of the reflection coefficient can be made dependent on a parameter δ defined in the following way. We set the central wavelength of the forbidden band to be λ g and determine the quarter wavelength thicknesses to be a 0 and b 0 for the higher and lower refractive index layers, respectively. The reflectance and transmittance and respective reflection and transmission coefficients, as the phase of the reflection coefficient can be obtained for any combination of a and b such that a = a 0 (1 + δ) and b = b 0 (1 − δ), with −1 < δ < 1. Thus, the extreme pairs are a = 2a 0 and b = 0, or a = 0 and b = 2b 0 .
In Figure 9 are shown results of the transmittance and phase of the reflection coefficient of a photonic crystal made of SiO 2 and ZrO 2 with 7 units cells. a 0 and b 0 were estimated for two central wavelengths λ g = 600 nm and λ g = 1, 200 nm. The unit cell is constituted by three layers: two lateral layers of thickness a/2 = λ g /(n h 8) and a central layer of thickness b = λ g /(n l 4) . Therefore, the total optical length is λ g /2. A unit cell built by this way has a mirror symmetry (inversion center) at the middle of the layer of lower refractive index and is required for a well-defined value of the Zak phase in the dispersion bands. Arbitrary values of the thickness of each layer, keeping the total optical length constant, would not change the forbidden band for an infinite crystal but change the symmetry properties of the bands in the dispersion relation (Zak, 1984;Xiao et al., 2014). The forbidden bands are the dark regions in the transmittance (Figure 9). The corresponding phase of the reflection coefficient exhibit gradients either positive from 0 to π, or negative from 0 to −π. Outside the forbidden bands the average phase is large and changes abruptly from π to −π.
By stacking two or more finite crystals (see Figure 10) with alternating sign of the phase of the reflection coefficient, edge states arise when the impedances match. The number of solutions of the impedance condition is not necessarily single valued. Indeed multiple edge states can be achieved. Two examples of edge states generated at the first bandgap of quarter wavelength crystals and their counterparts with band inversion are presented in Figure 11. In a stack of two photonic crystals PC1-PC2 using the same materials, where the dimensions of the PC2 were tailored in order that its second band gap has opposite sign of the reflection coefficient, only ore edge state arises ( Figure 11A). This optical Tamm state results in a narrow transmission line in the middle of bandgap. The wavelength of this state at vertical illumination differs from λ g = 600 nm, due to the small variation   of thickness of the layers and the wavelength dependence of the dielectric functions, not considered in the optical length equation.
When three PCs are stacked together, the edge states at both boundaries can couple and lead to double transmission lines ( Figure 11B). Surprisingly, not only two lines appear but 4. The central transmission narrow bands are expected from the frequency splitting due the coupling between edge states of both boundaries between PC1 and PC2. They are located at almost identical distance to the wavelength of the single edge state. But other two narrow bands arise at larger distance approaching the edges of the bandgap. The full-width at half-maximum (FWHM) of the narrow bands is typically smaller than 0.5 nm for the PCs with 7 unit cells, but decreases rapidly with increasing number of units cells. Therefore, it is worth to look more in detail into the field profiles and respective phase for the edge states in a PC1-PC2-PC1 stack (Figures 12, 13). The field component E y (corresponding to s-polarized light at vertical incidence) and its modulus along the transverse direction of stack indicates the type of edge state that is formed. The field inside the forbidden band decay exponentially with the distance, characteristic of an evanescent wave. Thus, the reflection coefficient can be expressed by r = 1.0 exp(iφ). Its phase has either a positive or negative value. At the edge states the field distribution has a fully distinct pattern inside the structure, peaking at the boundary of the PCs and decaying with z in both directions ( Figure 13A). The number of unit cells of the PC2 determines the strength of the coupling and the corresponding frequency splitting width. Therefore, at large number of unit cells the edge state do not interact and both spectral lines degenerate into the single line characteristic of a single edge state. The smallest thickness of the PC2 is a single unit cell, for which the interaction between edge states and the separation between the corresponding transmission lines is maximum.
The field profiles of the lateral spectral lines in the bandgap that arise in the symmetric structure PC1-PC2-PC1 cannot be explained by a equivalence between interacting edge states and narrow transmitted lines. The field distribution has maxima centered in the stack and not at the PCs boundary. However, they have similar narrow bandwidth and appear independently of the thickness of PC2. More intriguingly, the number of lateral bands inside the bandgap increases slowly with the number of unit cells of PC2. This is the behavior expected for a Fabry-Pérot resonator, with the difference that the lateral modes here occur in both sides of the central OTS modes. The resonance bandwidth of the optical Tamm state modes is comparable to that of Fabry-Pérot resonators based on identical Bragg mirrors of the same length. However, the relative distance of the resonances in the forbidden band to the band edges remains almost constant, when the angle of incidence varies from 0 • to 90 • (see Figure 11). By contrast, in a Fabry-Pérot resonator the cavity resonances tend to cross the forbidden band as function of the angle of incidence,   while the free spectral range remains constant (Saleh and Teich, 2007).
The formation of multiple edge states based on binary and quaternary units cells was investigated in Bianchi and Kahn (2020). The interplaying between Fabry-Pérot modes and optical Tamm states in structure formed by a photonic crystal and two metallic mirrors was studied in Durach and Rusina (2012). A deeper analysis of the OTS and lateral modes arising in multiple alternating PCs with band inversion and between PCs and plasmomic layers is outside the scope of the present article and will be addressed in other publication.
In Figure 14 are presented experimental results of the reflectance on two samples measured at quasi-vertical incidence, using an objective of small numerical aperture for illumination and light collection (NA = 0.12). The samples were fabricated by alternated deposition of layers of Ta 2 O 5 and SiO 2 making a stack of three photonic crystals with 7, 3, and 7 units cells, respectively. The layers of the first and third PC have a nominal thickness of 70 nm (Ta 2 O 5 ) and 100 nm ( SiO 2 ). The nominal thickness of the layers of the central PC are 200 nm (Ta 2 O 5 ) and 100 nm ( SiO 2 ). The spectral shift of the bandgap and resonances are due to fabrication. Small changes in the effective layer thickness are due to different distances from each substrate to the crucible and different evaporation angle. Despite the limited conditions of the deposition pronounced dips corresponding to the OTS and lateral states are visible in the reflectance spectrum. The reflectance dips are probably deeper that measured, because the spectrometer used has a spectral resolution of approximately 5 nm.
When the central protonic crystal (PC2) of the former layered structure is substituted by a homogeneous layer of a material with refractive index n and thickness d ≤ λ/(2n) the ordinary Fabry-Pérot resonator is obtained. It must be stressed that in the structure PC1-PC2 the thickness of the layer containing the boundary of the crystals is d 12 = (a 1 /2 + a 2 /2). Using the thickness values of the structure used in Figure 12A we have d 12 = 135 nm. The material of thickness d 12 is ZrO 2 with a refractive index of n ≈ 2.16 at λ = 600 nm. This is approximately half wavelength in the medium for a vacuum wavelength of 600 nm. Thus, the single OTS can be regarded as the single Fabry-Pérot mode in cavity of refractive index n h and length d 12 . Indeed, using the wavelength of the OTS line of Figure 12A, λ OTS ≈ 584 nm at vertical incidence, d 12 = 135 nm is exactly half wavelength in ZrO 2 of vacuum wavelength λ = 584 nm.
We can as well substitute the central PC2 by a plasmonic layer, as silver or gold. It is well-known that a stack of a quarter wavelength photonic crystal attached to a plasmonic layer produces an optical Tamm state (Vinogradov et al., 2006;Kaliteevski et al., 2007;Henriques et al., 2020). We can expect a double OTS for a thin layer of few tens of nanometer of silver at optical wavelengths. Effectively, two modes arise for a PC made of ZrO 2 , or Ta 2 O 5 and SiO 2 . In Figure 15 are presented the reflectance as function of angle on incidence and the profile of the the electric field in the structure for excitation at vertical incidence. The spectral separation of the modes depends on the thickness of the silver layer, decreasing with increasing thickness. The symmetry properties of the plasmon modes are reflected in the field distribution of the edge states inside the bandgap (see Figure 16). The higher energy mode (λ a = λ 1 = 540.3 nm) has an anti-symmetry electric field distribution, using as reference the center of the silver layer. The mode with lower energy (λ s = λ 2 = 592.3 nm) has smaller amplitude and is symmetric. Hence, the symmetric and anti-symmetric modes of a plasmonic thin film are transformed inot a pair of optical Tamm states width symmetric and antisymmetric electric field distribution in the stack. Unlike in the symmetric configuration of single plasmonic film requiring evanescent wave excitation, the excitation of the both optical Tamm states is here from air, whereas a the substrate is silica.

CONCLUSIONS
In this article three types of resonances with narrow bandwidth involving layered plasmonic media were analyzed: resonances arising in coupled Fabry-Pérot cavities of silver mirrors with heterogeneous intra-cavity media, presenting optical cavity strong coupling; hyperbolic layered media based on subwavelength layers of plasmonic and dielectric films, supporting bulk plasmon modes with large propagation constant and large field amplitudes inside the structure; and finally, multiple optical Tamm states, arising in stacks of photonic crystals and in plasmonic films attached to photonic crystals.
The resonance anti-crossing in coupled Fabry-Pérot cavities is not only one of the most simple examples of classical strong coupling but also represents a platform with potential in refractive index sensing applications. Fabry-Pérot cavities based on distributed Bragg reflectors, instead of silver mirrors, could also be used for the strong coupling, but the total thickness of the stack is much larger and the sensitivity to variations of refractive index of external media much weaker.
The migration of Fabry-Pérot cavities into layered hyperbolic metamaterials is put straightforward. The analysis of the bulk plasmon modes inside layered hyperbolic metamaterials of finite size reveals two fundamental properties: the amplitude of the field increases with the real part of the propagation constant, and the value of the imaginary part k i also increases in the same order, which implies a systematic decrease of the plasmon propagation length. Thus, in applications requiring plasmons of long propagation constant, extremely thin layers in large number is counterproductive. A small number of bilayers supporting fewer plasmon modes could achieve better results in terms of averaged length of propagation. On the other hand, the increasing field amplitude with the plasmon mode order inside the stack can benefit applications based on nonlinear effects (Davoyan et al., 2009). The periodic plasmonicdielectric multilayers can also be applied in fluorescence lifetime engineering. Excited nanoemitters on top of the stack can couple to the evanescent modes of the bulk plasmons. Fast fluorescence decay can be achieved due to the large density of optical modes in the stack. The strong near-field confinement allowed by the surface plasmons is related to the large progation constants with relatively low damping of the evanescent waves and thus, permits the optical interaction between distant objects, as in the Förster resonance energy transfer (FRET) (Roth et al., 2018;Gonçalves et al., 2020).
Multiple optical Tamm states arising in ternary photonic crystals, or symmetric structures made of a plasmonic film and two identical Bragg reflectors form a new family of resonances relying on edge states with topological properties. The number of full or high transmission narrow resonances increases with the number of Tamm states. Two kinds of resonances were found: hybridized Tamm modes resulting from the interaction of single optical Tamm states at short distance and lateral modes with apparent Bragg field distribution. The bandwidth of the resonances calculated using the experimental dielectric functions of the materials is typically less than 0.4 nm, for photonic crystal with a length of 7 unit cells. This corresponds to Q-factors of the order of 1,000s in the visible spectrum. Higher values can be reached by increasing the number of unit cells in the external Bragg mirrors. Such narrow resonances can be of high interest in nanolasers and photonic switching systems, if their spectral position could be changed by variation the dielectric function of one of the materials employed. For some geometries one-way propagation can also be investigated by introducing materials with non-reciprocal properties in the structure as in (Khanikaev et al., 2009). Optical Tamm states present resonant modes with a bandwidth comparable to Fabry-Pérot modes, in cavities based on identical distributed Bragg mirrors, but their spectral dependence on the angle of incidence differ considerably. Multiple OTS modes based on ternary photonic crystal stacks can provide a robust alternative to optical filters based on Fabry-Pérot cavities. By other side, OTS modes arising in structures including plasmonic layers as the PC-Ag-PC, are also useful in sensing applications due to the narrow resonances, but can only be restricted to few plasmonic layers owing to the optical absorption. To my knowledge, the excitation of the two OTS modes with symmetric and antisymmetric field profiles was analyzed for the first time.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
The author confirms being the sole contributor of this work and has approved it for publication.

FUNDING
Some materials used in the preparation of samples were funded by the German Federal Ministry of Education and Research (BMBF) Grant 01DK17040.