## Floquet Engineering of Polaritonic Amplification in Dispersive Photonic Time Crystals

Mustafa Goksu Ozlu, ∗ Vahagn Mkhitaryan, ∗ Colton B. Fruhling, Alexandra Boltasseva, and Vladimir M. Shalaev Elmore Family School of Electrical and Computer Engineering Birck Nanotechnology Center, Purdue University, West Lafayette, IN 47907, USA

(Dated: August 2, 2024)

In this study, we investigate the dynamics of dispersive photonic time crystals (PTCs) and their potential applications for controlling light-matter interaction. By employing the Lorentz-Drude model, we analyze theoretically and via numerical simulation the effects of periodic modulation of dispersion parameters, revealing the emergence of hybrid bandgaps from interaction of polaritonic branches with unique characteristics. Our study demonstrates that dispersive PTCs offer novel excitation channels and amplification possibilities, that require lower modulation frequencies compared to non-dispersive systems thus alleviating experimental challenges for the realization of PTCs in the optical regime. These findings pave the way for advancements in polaritonic lasing and resonant Raman scattering.

## I. INTRODUCTION

tion while the probe beam propagates inside the medium, relaxing the typical phase-matching conditions [26].

The study of time-varying photonic media has a history spanning several decades [1-4]. Nevertheless, recent years have witnessed a resurgence of interest in this field revealing numerous new physical phenomena [5] such as the creation of squeezed and entangled photon states [69], spatiotemporal beam steering [10, 11], and the breaking of Lorentz reciprocity [12, 13]. One of the most notable outcomes of engineering temporal dynamics is the concept of photonic time crystals (PTCs)[4, 14], where the permittivity of the host material undergoes a periodic modulation with time, analogous to spatial photonic crystals which are periodically patterned in real space [15, 16]. Thanks to the non-Hermitian properties of PTCs these systems exhibit unique light-matter interactions, leading to intriguing phenomena such as generation and amplification of radiation fields from free electrons and dipoles [17, 18].

Abrupt and strong changes in material's optical properties over time result in reflections and refractions at temporal interfaces [19, 20], akin to spatial boundaries. However, unlike spatial boundaries, temporal ones require dynamical modulation of optical properties, creating an open system and allowing time-reflected and refracted waves to have different energies than the incident beam. PTCs exploit this to create and amplify electromagnetic radiation in their momentum bandgap region, where time-reflected and refracted waves constructively interfere [17, 21]. This amplification resembles parametric amplification [1, 3, 22, 23] but differs in two distinct features. Firstly, the presence of strong modulation of material parameters, such as strongly modulating permittivity around the epsilon near zero point[24, 25], breaks the applicability limits of the conventional perturbative nonlinear treatment. Secondly, it is assumed that the entire medium undergoes simultaneous modula-

∗ These authors contributed equally to this work

Amplification within the bandgap of PTCs fundamentally arises from the transfer of energy from the modulation to the probe. This process necessitates interactions among various frequency components within the medium, leading to the creation of new components at the probe frequency. In a non-dispersive medium, where only a single frequency mode is supported for each momentum, effective amplification requires that the modulation frequency be twice that of the probe beam. Thus, accessing the momentum bandgap in non-dispersive media demands modulation at twice the probe beam frequency, which poses significant challenges for the realization of PTCs at optical or near-infrared (NIR) frequencies. This is why time reflection and exponential growth in the momentum bandgaps have so far only been observed experimentally in the microwave domain [27, 28]. However, recent advances in optical materials that can be strongly modulated on an femtosecond time scale [24, 29-31] create a potential for experimental realization of PTCs also in the optical domain [32, 33]. These materials are typically strongly dispersive and lossy. The material dispersion enforces electric field to evolve continuously in time, and potentially makes it more challenging to observe unique PTC behavior[34]. Additionally, the dispersive nature of materials couples energy from the light to mechanical motion in the material forming hybrid states of light and matter known as polaritons, which are new quasiparticles of the coupled system with well-defined energy and momentum. Building on the nondispersive foundations of recent PTC studies [2, 17, 18, 21, 35], we explore the inclusion of material dispersion in our analysis which opens up the possibility of investigating the interaction between matter and radiation field in the presence of optical modulation. We further discuss and examine more realistic experimental scenarios exploiting the full potential of such materials.

Recent studies for this exciting prospect have begun by investigating plane wave propagation in a dispersive time-varying media [36], a two-level model with periodic

modulation of resonance frequency [37, 38] and widening of momentum bandgap near material or structural resonance [39]. A broader investigation of polaritonic interactions and its effect on PTC amplification is still needed. In this paper we focus on analyzing the overall PTC dispersion, also referred to as optical band-diagram, in the presence of Lorentz-Drude type material dispersion. Specifically, we explore the modulation of the coupling between different polaritonic branches. We observe the emergence of hybrid bandgaps which stems from interaction between upper and lower polaritonic branches. These hybrid bandgaps have experimentally favorable conditions.

Our main findings reveal that when the parameters of dispersive materials are modulated, additional bandgaps with lower modulation speeds emerge in the system. These bandgaps are accessible from the modes of the unmodulated medium and enable excitation channels for coherently interfering time-reflected and refracted waves, establishing exponentially growing fields. This relaxes the fast modulation requirement of nondispersive materials and suggests more accessible experimental pathways for observing such effects with much slower modulation frequencies.

As amplification in PTCs stems from non-perturbative nonlinear interactions, the amplification of waves within the polaritonic PTCs is intricately connected to the nonlinear interaction among polaritons [40]. Thus, by examining polaritonic PTCs we may gain valuable insights into the physics of these strong, non-perturbative interactions among polaritons that pave the way forward for enhanced resonant Raman scattering and polaritonic lasing [41, 42].

## II. POLARITONIC BANDS WITH MODULATION

When the optical properties of a medium are periodically modulated with a frequency Ω, its properties remain time-translation invariant under a temporal shift of T = 1 Ω, i.e. / Ψ( t + T ) = Ψ( ), where Ψ is a physical t quantity such as electromagnetic field amplitude or polarization in the medium and t is the time. This temporal periodicity causes the ω -k dispersion of the supported modes to also become periodic in frequency with period Ω. Therefore, the dispersion of the system is completely characterized by the frequency range of 0 to Ω, known as the 1st Brillouin zone. Owing to this periodicity, the modes of the medium can be described using the Floquet theory, by expanding the relevant quantities in the form Ψ ( ) = t e iω k t ∑ ∞ n = -∞ Ψ n e in Ω t ,with ω k +2 πn = ω k , where n is an integer.

In the case of a stationary, nondispersive optical medium, the dispersion relation is given by the expression ω = kc/n , where c is the speed of light in a vacuum, and n is the refractive index. This relation defines a single line for positive ω and k shown in Figure 1 by black dashed lines. The effect of periodic modulation of the dielectric function on the dispersion relation for such nondispersive materials can be visualized as dispersion lines of the unmodulated medium being translated along the frequency axis by integer multiples of the modulation frequency Ω . In this case the dispersion for negative frequencies also becomes relevant as they get shifted to the positive frequency domain and mix with other dispersion lines. The translation of the dispersion lines causes the branches in different Brillouin zones to cross each other leading to complex frequency eigenvalues, ω k , with degenerate real parts. This contrasts with the unmodulated, lossless medium, which has only real eigenfrequencies. Following the analogy of the frequency band gaps in spatial photonic crystals, these regions with imaginary eigenfrequencies have been called momentum bandgaps with the characteristic signature of exponentially growing photon fields [17, 18]. The emergence of the momentum bandgaps and complex eigenvalues are illustrated for a nondispersive medium in Figure 1. The modes, as they are linear, cross each other at frequency Ω / 2 resulting in a bandgap at that frequency.

FIG. 1. (a) The band diagram of a non-dispersive photonic time crystal with electric permittivity ε = 1 + 0 3 sin (Ω ) . t the dashed green line shows the dispersion relation of the unmodulated medium with ε = 1. The blue lines indicate the real parts of eigenvalues while the red indicates the imaginary parts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-d2a932e8-f078-4583-8082-2b0adc11ddda.jpg)

In the case of a dispersive medium, the energy no longer remains purely in the fields. The electromagnetic waves hybridize with material excitations and form polaritons. The response of the materials that possess bound charges can be described classically with Lorentz type frequency-dependent dielectric function.

$$\varepsilon ( \omega ) = \epsilon _ { \infty } + \frac { \omega _ { p } ^ { 2 } } { \omega _ { 0 } ^ { 2 } - \omega ^ { 2 } - i \gamma \omega } \quad \quad ( 1 )$$

where, ω 0 , γ , and ε ∞ are respectively the central frequency, Drude damping and high-frequency permittivity

FIG. 2. (a) The band diagram of a non-dispersive photonic time crystal with electric permittivity ε = 1+0 3sin(Ω ) the dashed . t green line shows the dispersion relation of the unmodulated medium (b) The dispersion diagram of a Lorentz type material with ω p = 1 6Ω . ω 0 = 0 6Ω . ε ∞ = 1 (c) PTC band diagram of the dispersive material presented in part b. The plasma frequency is modulated with ω p ( ) = 1 6 + 0 8 sin (Ω ). t . . t The blue lines show the real and the red lines show the imaginary part of the frequency eigenvalues. The dispersion branches of the unmodulated medium is presented with dashed lines (d) A close up of the region where two types of momentum bandgaps are highlighted (e) A close up of the region where the avoided crossing is highlighted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-13e0dc2e-b5d1-4907-a073-18c0cc16e789.jpg)

of the material and ω 2 p = 4 πe N/m 2 ∗ is the plasma frequency of the material defined by free carrier concentration N and effective mass m ∗ . This form of dielectric function follows from the harmonically oscillating polarization density of the material excited by electromagnetic waves of the same frequency.

$$\frac { d ^ { 2 } P } { d t ^ { 2 } } + \gamma \frac { d P } { d t } + \omega _ { 0 } ^ { 2 } P = \frac { N e ^ { 2 } } { m ^ { * } } E \, \quad \quad ( 2 ) \, \quad \, \text{an} \,,$$

Polaritons in such media are formed when electromagnetic waves couple to the mechanical oscillations of these bound charges through induced polarization created due to their harmonic motion. The dispersion relation for such material follows from Maxwell's equations and is given by k 2 = ω ε ω /c 2 ( ) 2 . The solution of this quadratic equation leads to two dispersion branches above and below characteristic frequency ω / 0 √ ε ∞ , which are conventionally called upper and lower polaritonic branches. These two branches arise from the coupling of phononic and photonic dispersion lines. Due to the coupling, the branches bend outward and form avoided crossing with an energy difference ∆ E known as the Rabi splitting. An example of such dispersion curves is shown in Figure 2.a calculated for normalized oscillator parameters ω p = 1 6Ω, . ω 0 = 0 6Ω, and . ε ∞ = 1 .

We now incorporate temporal changes in material properties using dispersion parameters. Although our model includes dispersive effects the response of disper-

sion parameters to external stimuli is still assumed to be instantaneous. Experimentally the most accessible modulation is through optical pumping of semiconductors and changing conduction band electron density or effective mass. These two parameters directly affect the plasma frequency, therefore in our model, we assume to be modulating the plasma frequency.

$$\omega _ { p } ^ { 2 } ( t ) = \omega _ { p 0 } ^ { 2 } [ 1 + \Delta ^ { 2 } \cos { ( \Omega t ) } ] \quad \quad ( 3 ) \quad \underset { \sin } { \text{sim} }$$

Due to this periodicity, the fields inside the medium can be written as Floquet expansion and the following eigenvalue equation is obtained from Maxwell's equations.

$$\omega _ { k } \Psi _ { p } = \sum _ { n = - \infty } ^ { \infty } ( M _ { p - n } - n \Omega I _ { 4 \times 4 } \delta _ { p n } ) \Psi _ { n } \quad \ \ ( 4 ) \quad \text{ditto} \\ \text{bran}$$

where Ψ k p = [ E xp , H yp , J xp , P xp ] T , with superscript T meaning transpose operation and J xp = ∂ P t xp , and matrices ¯ M m are given by

$$M _ { m } = \left [ \begin{array} { c c c c } 0 & - c k \delta _ { m 0 } / \epsilon _ { \infty } & 4 \pi i \delta _ { m 0 } & 0 \\ - c k \delta _ { m 0 } & 0 & 0 & 0 \\ - i f _ { m } & 0 & i \gamma _ { m } & i \omega _ { o m } ^ { 2 } \\ 0 & 0 & - i \delta _ { m 0 } & 0 \end{array} \right ] \quad ( 5 ) & \text{if $i$ to} \\ \text{branc} \\ \text{italic}$$

The summation in Eq.4 runs over all possible Floquet modes, however in numerical implementation, we consider only finite number modes and truncate this infinite system of equations into generalized eigenvalue problem:

$$\omega _ { k } \Psi ^ { ( k ) } = \bar { A } \Psi ^ { ( k ) } & & ( 6 ) & & \text{loc}$$

with Ψ = [ Ψ ( k ) 1 , Ψ ( k ) 2 , ..., Ψ ( k ) ] and ¯ A is a matrix formed by elements from Eq.5. Here for an isotropic medium, without losing the generality we assumed the waves propagate along the z-axis. The equations for waves with orthogonal polarization Ψ = [ E , H y x , J y , P y ] T have very similar form, hence we will focus only on one of them for our discussions. The details of the calculation can be found in the Supplementary Information(SI).

Solving the eigenvalue Eq.6 results in the band diagram for a dispersive photonic time crystal along with the corresponding eigenfunctions. Introducing modulation to a dispersive medium, like in the nondispersive case, causes the dispersion branches to be translated with the modulation frequency. The overlapping polaritonic branches interact and give rise to bandgaps or avoided crossing based on the type of dispersion line as shown in Figure 2.b. The dashed lines show the translated dispersion lines of the stationary medium and the blue lines show the dispersion line of the modulated medium. The dispersion parameters and eigenfrequency of the Floquet modes are expressed relative to the modulation frequency Ω.

The presence of polaritonic bands and their modulation introduce novel phenomena into the physics of time-varying medium. One of the most prominent effects for this case is the interaction between lower and upper polariton branches. The crossing of lower and upper branches belonging to positive and negative frequency quadrants causes the real part of the frequency eigenvalues to split and give rise to a frequency bandgap highlighted in Fig 2.c which is otherwise not observed in nondispersive case.

The interaction of polaritonic branches causes the emergence of bandgaps as well. First, a bandgap with similar characteristics to non-dispersive PTC bandgaps arises from the interaction between two lower branches, highlighted in Fig 2.d as same branch bandgap. These bandgaps occur in the middle of the unit cell at frequency Ω 2, with a flat dispersion line. / A distinguishing factor for these bandgaps from the non-dispersive case is the additional excitation channels made available by the upper branch with matching momentum.

Furthermore, the crossing of upper and lower polaritonic branches leads to the formation of a novel bandgap which we term as hybrid bandgaps with distinct features. First, hybrid bandgaps does not occur on a flat line resulting in nonzero group velocity. Therefore, the energy is amplified directionally. Secondly, the modulation frequency required to create this bandgap is not necessarily twice the frequency of the unmodulated branch, but it is the energy difference between the upper and lower branches. Therefore, the requirement of having the modulation speed twice the probe beam's frequency does not apply.

The hybrid bandgap's size and amplification heavily rely on the coupling strength between the upper and lower branches, and is maximized when the group velocities of the branches match to allow a more profound overlap. This alignment occurs when the hybrid bandgap coincides with the momentum region where Rabi splitting is present between the unmodulated branches. We accentuate this effect using dispersion parameters of a well-known phononic optical material hexagonal boron nitride (h-BN)(see Table I)[43]. The band diagram for this case is given in Fig 3.a.

For a phononic medium, the minimum separation of the lower and upper branches is equal to ω / p √ ε ∞ which is proportional to the difference between the transverse phonon (TO) and longitudinal phonon (LO) frequencies such that ∆ ω = √ ω 2 LO -ω 2 TO [36]. The modulation frequency has been chosen to be half of the Rabi splitting (∆ ω ) between the upper and lower branches. This allows maximal overlap coupling between lower and upper branches where their group velocities match, as shown in Fig 3.b, creating profound bandgaps. We investigate

TABLE I. The dispersion parameters of hBN in tangential direction. Modulation frequency Ω = 12 THz

| ω p           | ω 0         | γ        | ε ∞   |
|---------------|-------------|----------|-------|
| 58 THz        | 41 THz      | 0.2 THz  | 4.9   |
| 1920 cm - 1 √ | 1360 cm - 1 | 7 cm - 1 |       |
| 2 ε ∞ Ω       | 3.40 Ω      | 0.017 Ω  | 4.9   |

FIG. 3. (a) PTC band diagram of the hBN like dispersive material presented in Table 1. The plasma frequency is modulated with ω p ( ) = 0 61+0 3 sin (Ω ) where Ω = 13THz. The t . . t blue lines show the real and the red lines show the imaginary part of the frequency eigenvalues (b) Group velocity of the upper and lower branches for the unmodulated hBN dispersion

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-8027f95a-ea83-465a-a000-361c2acb73b2.jpg)

the characteristics of wave propagation in these bandgaps and the possibility of obtaining amplification with slower modulations in the subsequent section.

## III. EXCITATION OF MODES

For a given medium, the eigen solutions and their dispersions discussed in the previous section characterize the natural modes supported in that medium. When this medium is excited by external sources, the induced field can be represented analytically as time-dependent superpositions of these eigenmodes. Such analytical formulations give great flexibility to study the dynamics of the modes in the system. In particular, it allows us to observe how individual modes are excited, and how they evolve and interact among themselves. To demonstrate this,

FIG. 4. (a) The excitation spectrum of the dispersive PTC discussed in Figure 2.c excited with a dipole source at frequency 0 75Ω in log scale. The blue dashed lines show the real . part of the band-diagram (b) Time-momentum graph showing the evolution of the fields in the excited medium in log scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-9169cd9f-560a-47b0-b29b-18876284e8de.jpg)

we discuss the dynamics of a dispersive, time-varying medium under dipole and plane wave excitations.

We first illustrate the system dynamics using a dipole excitation. This study builds on previous discussions regarding a non-dispersive medium [17]. Specifically, we examine the excitation of the system whose band diagram is depicted in Figure 2.b. The system is, excited by a time-harmonic dipole source with an oscillation frequency ω d = 0 7Ω, chosen to fall outside of any bandgaps. . Fig 4.a shows the excitation spectrum as density plots superimposed by the band diagram (green dashed lines).

The excitation frequency is represented by the black line. Initially, the light couples most strongly to the branches of the dispersion diagram closest to the excitation frequency. However, the entire band is also excited because the radiated light from a point source has all the momentum components. Additionally, higher-order bands are excited through modulation-induced interaction. Similar to the non-dispersive case, waves excited at the momentum bandgap grow exponentially and appear with the strongest amplitude.

Exponential growth can also be observed from the time-momentum graph in Figure 42b. Initially, the mo-

mentum values corresponding to the supported modes of the system at the source frequency are excited. As time progresses, the modes with the correct momentum range to couple into the bandgaps grow exponentially, eventually surpassing fields at any other wavenumber.

Excitation via dipole sources displays the range of modes. However, experimental excitation is usually confined to propagating or evanescent waves inside the medium, which are bound within the supported momentum range. To investigate excitation with a wave we employ FDTD simulations with a probe beam propagating within a modulated medium. The details of the calculation are provided in the SI.

Before modulation, the probe beam can have momentum and frequency residing only on one of the dispersion branches of the non-modulated medium. Once the modulation starts, the beam's momentum remains constant while its frequency couples to the Floquet modes. When the probe beam momentum is within the bandgap, the beam couples into the momentum gap mode and starts growing. For a nondispersive medium, the bandgap always occurs at k = Ω (2 / n avg c ) [17, 18]. In such a medium, for a wave to have this momentum, its frequency should always be Ω / 2 Hence, exciting the bandgap requires very fast and strong modulation of the optical properties, making it challenging to possibly observe photonic time crystals in the optical regime. However, the situation is dramatically different for dispersive materials and these strict conditions can be relaxed as there are multiple features, including new types of band gaps and avoided crossings in the band diagram that lead to different excitation regimes.

First, we investigate the same-branch bandgaps. We again use the dispersion parameters of the medium in Figure 2.c. In this case the bandgap occurs at 2.5 k = Ω /c . Additionally, contrary to the non-dispersive case there are two frequencies where propagation with this momentum is allowed for the non-modulated medium, on the upper and lower branches. We excite this medium from the upper branch with ω = 3Ω as highlighted in Figure 5.a. Once the modulation starts some of the energy remains at the original branch and continues to propagate without a change in the amplitude. This is a result of the original branch being almost identically present in the dispersion of the modulated medium and the bandgap does not occur on this branch. As shown in figure 5.b some portion of the energy of the incident beam couples into the bandgap created by the lower branches occurring at Ω / 2 and at every integer multiple of Ω / 2. The light in this branch shows the same characteristics as the bandgap mode in a nondispersive medium. As the modulation starts, the beam coupled to the bandgap is equally separated into forward and backward propagating modes and grows exponentially. The bandgap occurs on a flat band, therefore, the group velocity of the beam is zero and the wave remains localized in space until modulation ends. It is observed that the new modes are excited in both quadrants equally, which means the forward and backward propagating beams have equal amplitudes. The video of the temporal evolution of the beam for this excitation is given in Supplementary Video 1. Higher order frequencies are also visible at 3Ω / 2 as a result of harmonic generation with the modulation. Although this case demonstrates transferring energy between polaritons of the medium, the bandgap mode still occurs at half the frequency of the modulation.

This is not the scenario for the hybrid bandgap. We investigate the band diagram shown in Figure 3.a. The hybrid bandgap occurs around k = 7Ω /c where the upper and lower branches cross each other. To couple into the bandgap mode, we excite from the upper branch at ω = 4 23Ω and . k = 6 7Ω . /c as shown in Figure 6.a. After modulation is introduced, the bandgap occurs on the branch where the beam was already propagating before the modulation, therefore all the energy directly couples to the bandgap mode resulting in amplification. Notably, the bandgap in this case is not a flat band as in the same-branch case but is instead slanted. Therefore, the resulting amplified beam propagates with non-zero group velocity. This is a distinct difference between the dispersionless cases investigated previously. Figure 6.b shows the ω -k relationship of the beam inside the medium during the modulation where one can observe the growth in the field amplitude and the excitation of the lower branches Ω away in frequency from the probe beam. Supplementary Video 2 shows the temporal evolution of this excitation.

One notable observation of this structure is the requirement of modulation speed is more than 8 times lower than that of a non-dispersive medium. The interaction between the polaritonic branches allows for PTC bandgaplike amplification to happen at significantly slower modulation. In the example shown a modulation at 13 THz gives rise to directional amplification of a beam at 55 THz. This feature is a clear path forward for observing PTC bandgaps in the infrared (IR) regime and possibly in the optics regime with appropriate selection of material and modulation properties.

## IV. CONCLUSION

In this work, we have investigated the effects of periodic temporal modulation on dispersive mediums. The dispersion lines of the material start translating and mixing via modulation. The resulting band diagrams exhibit bandgaps that emanate from the interaction of polaritonic branches. Due to dispersive materials being able to support several modes at a given momentum, there are more possibilities for interactions. This allows for the occurrence of new bandgaps and new channels for excitation.

One significant feature observed in the polaritonic PTCs is the hybrid bandgap that occurs through the interaction of a lower and an upper dispersion branch. The region where the two branches have the closest en-

FIG. 5. (a) Frequency-momentum relation of the excitation wave inside the dispersive medium without modulation. The blue lines show the dispersion relation of the unmodulated medium. (b) The frequency-momentum relation of the field inside the dispersive PTC. The blue lines show the band-diagram of the PTC

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-96db904f-0612-4734-9f9f-a623e92927f0.jpg)

FIG. 6. (a) Frequency-momentum relation of the excitation wave inside the medium with hBN like dispersion without modulation. The blue lines show the dispersion relation of the unmodulated medium. (b) The frequency-momentum relation of the field inside the dispersive PTC. The blue lines show the band-diagram of the PTC.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00552v1-figure-f9a489d2-c33a-4f30-af1c-293de3412de4.jpg)

ergy and exhibit avoided crossing is where they are most strongly coupled. When driven by their energy difference in this regime, it is possible to create a strong interaction between polaritons and the modulation, that allows for energy to be transferred from modulation to the probe beam. This is achieved with much lower modulation frequencies than what is necessary for non-dispersive PTCs. Allowing for experimental observations in higher energy regimes of the electromagnetic spectrum.

When modulated with a frequency proportional to the energy separation of the LO and TO phonons, phononic materials can be used for polaritonic amplification. Although creating strong optically induced changes in these materials are challenging, it has been shown that it is possible to drive nonlinear response in phononic materials through Born effective charge and dielectric screening of the electric field by electrons. This nonlinearity has been used to observe parametric amplification of optical

phonons.[23]

Apart from employing the phononic material dispersion one other possible realization of polaritonic amplification is through constructing the desired dispersion lines through structured materials. Bulk and surface plasmons of thin films hybridize near the epsilon near zero region of the thin films to give rise to Berreman and ENZ modes, which can be coupled through relatively low frequency modulation to give rise to amplified polaritons. The interaction of polaritonic modes in dispersive photonic time crystals (PTCs) has the potential for amplification across a broad electromagnetic spectrum. This enhancement

- [1] W. H. Louisell, A. Yariv, and A. E. Siegman, Quantum Fluctuations and Noise in Parametric Processes. I., Physical Review 124 , 1646 (1961).
- [2] E. S. Cassedy and A. A. Oliner, Dispersion Relations in Time-Space Periodic Media: Part I-Stable Interactions, Proceedings of the IEEE 51 , 1342 (1963).
- [3] B. R. Mollow and R. J. Glauber, Quantum theory of parametric amplification. i, Physical Review 160 , 1076 (1967).
- [4] J. R. Zurita-S´nchez, a P. Halevi, and J. C. CervantesGonz´ alez, Reflection and transmission of a wave incident on a slab with a time-periodic dielectric function (t), Physical Review A - Atomic, Molecular, and Optical Physics 79 , 1 (2009).
- [5] E. Galiffi, R. Tirole, S. Yin, H. Li, S. Vezzoli, P. A. Huidobro, M. G. Silveirinha, R. Sapienza, A. Al`, u and J. B. Pendry, Photonics of time-varying media, https://doi.org/10.1117/1.AP.4.1.014002 4 , 014002 (2022).
- [6] J. T. Mendon¸ ca, A. Guerreiro, and A. M. Martins, Quantum theory of time refraction, Phys. Rev. A 62 , 033805 (2000).
- [7] V. V. Dodonov, A. B. Klimov, and D. E. Nikonov, Quantum phenomena in nonstationary media, Physical Review A 47 , 4422 (1993).
- [8] V. V. Dodonov, M. A. Marchiolli, Y. A. Korennoy, V. I. Man'ko, and Y. A. Moukhin, Dynamical squeezing of photon-added coherent states, Physical Review A 58 , 4087 (1998).
- [9] J. T. Mendon¸ ca, G. Brodin, and M. Marklund, Vacuum effects in a vibrating cavity: Time refraction, dynamical Casimir effect, and effective Unruh acceleration, Physics Letters, Section A: General, Atomic and Solid State Physics 372 , 5621 (2008).
- [10] A. Shaltout, A. Kildishev, and V. Shalaev, Time-varying metasurfaces and lorentz nonreciprocity, Optical Materials Express 10.1364/OME.5.002459 (2015).
- [11] X. Wang, V. S. Asadchy, S. Fan, and S. A. Tretyakov, Space-Time Metasurfaces for Perfect Power Combining of Waves , Tech. Rep.
- [12] D. L. Sounas and A. Al`, Non-reciprocal photonics based u on time modulation, Nature Photonics 2017 11:12 11 , 774 (2017).
- [13] N. Chamanara, S. Taravati, Z.-L. Deck-L´ eger, and C. Caloz, Optical isolation based on space-time engineered asymmetric photonic band gaps, PHYSICAL RE-

can impact the scope and applications of polaritonic lasing, resonant Raman scattering, and polariton nonlinearities.

## ACKNOWLEDGMENTS

Authors acknowledge support by the U.S. Department of Energy under Award DE-SC0017717 (Control of Light), and the Air Force Office of Scientific Research under Award FA9550-21-1-0299

## VIEW B 96 , 155409 (2017).

- [14] J. R. Reyes-Ayona and P. Halevi, Observation of genuine wave vector (k or β ) gap in a dynamic transmission line and temporal photonic crystals, Applied Physics Letters 107 , 074101 (2015), https://pubs.aip.org/aip/apl/articlepdf/doi/10.1063/1.4928659/14106877/074101 1 online.pdf.
- [15] S. John, Strong localization of photons in certain disordered dielectric superlattices, Physical Review Letters 58 , 2486 (1987).
- [16] E. Yablonovitch, Inhibited Spontaneous Emission in Solid-State Physics and Electronics, Physical Review Letters 58 (1987).
- [17] M. Lyubarov, Y. Lumer, A. Dikopoltsev, E. Lustig, Y. Sharabi, and M. Segev, Amplified emission and lasing in photonic time crystals, SCIENCE 377 , 425 (2022).
- [18] A. Dikopoltsev, Y. Sharabi, M. Lyubarov, Y. Lumer, S. Tsesses, E. Lustig, I. Kaminer, and M. Segev, Light emission by free electrons in photonic timecrystals, Proceedings of the National Academy of Sciences of the United States of America 119 , 10.1073/PNAS.2119705119 (2022).
- [19] J. T. Mendon¸ ca, A. M. Martins, and A. Guerreiro, Temporal beam splitter and temporal interference, Phys. Rev. A 68 , 043801 (2003).
- [20] A. M. Shaltout, J. Fang, A. V. Kildishev, and V. M. Shalaev, Photonic Time-Crystals and Momentum BandGaps, Conference on Lasers and Electro-Optics (2016), paper FM1D.4 , FM1D.4 (2016).
- [21] J. B. Pendry, E. Galiffi, and P. A. Huidobro, Gain in time-dependent media-a new mechanism, Journal of the Optical Society 38 , 10.1364/JOSAB.427682 (2021).
- [22] S. Sugiura, E. A. Demler, M. Lukin, and D. Podolsky, Resonantly enhanced polariton wave mixing and parametric instability in a Floquet medium, The Journal of Chemical Physics 156 , 174110 (2022), https://pubs.aip.org/aip/jcp/articlepdf/doi/10.1063/5.0091718/16542313/174110 1 online.pdf.
- [23] A. Cartella, T. F. Nova, M. Fechner, R. Merlin, and A. Cavalleri, Parametric amplification of optical phonons, Proceedings of the National Academy of Sciences of the United States of America 115 , 12148 (2018).
- [24] N. Kinsey, C. DeVault, J. Kim, M. Ferrera, V. M. Shalaev, and A. Boltasseva, Epsilon-near-zero Al-doped ZnO for ultrafast switching at telecom wavelengths, Optica 2 , 616 (2015).

- [25] O. Reshef, I. De Leon, M. Z. Alam, and R. W. Boyd, Nonlinear optical effects in epsilon-near-zero media, Nature Reviews Materials 2019 4:8 4 , 535 (2019).
- [26] J. B. Khurgin, Photonic time crystals and parametric amplification: similarity and distinction (2024), arXiv:2305.15243 [physics.optics].
- [27] X. Wang, M. S. Mirmoosa, V. S. Asadchy, C. Rockstuhl, S. Fan, and S. A. Tretyakov, Metasurface-based realization of photonic time crystals, Science Advances 9 (2023).
- [28] H. Moussa, G. Xu, S. Yin, E. Galiffi, Y. Ra'di, and A. Al`, Observation of temporal reflection u and broadband frequency translation at photonic time interfaces, Nature Physics 2023 19:6 19 , 863 (2023).
- [29] M. Z. Alam, I. De Leon, and R. W. Boyd, Large optical nonlinearity of indium tin oxide in its epsilon-near-zero region, Science 352 , 795 (2016).
- [30] E. Lustig, O. Segal, S. Saha, C. Fruhling, V. M. Shalaev, A. Boltasseva, and M. Segev, Photonic time-crystals -fundamental concepts, Opt. Express 31 , 9165 (2023).
- [31] S. Saha, O. Segal, C. Fruhling, E. Lustig, M. Segev, A. Boltasseva, and V. M. Shalaev, Photonic time crystals: a materials perspective, Opt. Express 31 , 8267 (2023).
- [32] E. Lustig, O. Segal, S. Saha, E. Bordo, S. N. Chowdhury, Y. Sharabi, A. Fleischer, A. Boltasseva, O. Cohen, V. M. Shalaev, and M. Segev, Time-refraction optics with single cycle modulation, Nanophotonics 12 , 2221 (2023).
- [33] E. Lustig, S. Saha, E. Bordo, C. DeVault, S. N. Chowdhury, Y. Sharabi, A. Boltasseva, O. Cohen, V. M. Shalaev, and M. Segev, Towards photonic time-crystals: observation of a femtosecond time-boundary in the refractive index, in Conference on Lasers and Electro-Optics (Optica Publishing Group, 2021) p. FF2H.1.
- [34] Z. Hayran, J. B. Khurgin, and F. Monticone, ℏ ω ; versus ℏ k: dispersion and energy constraints on time-varying photonic materials and time crystals, Opt. Mater. Express 12 , 3904 (2022).
- [35] L. Felsen and G. Whitman, Wave Propagation in TimeVarying Media, IEEE Transactions on Antennas and Propagation Ap-18 , 242 (1970).
- [36] D. M. Sol´ ıs, R. Kastner, R. Kastner, and N. Engheta, Time-varying materials in the presence of dispersion: plane-wave propagation in a Lorentzian medium with temporal discontinuity, Photonics Research, Vol. 9, Issue 9, pp. 1842-1853 9 , 1842 (2021).
- [37] J. Sloan, N. Rivera, and M. Soljacic, Dispersion in photonic time crystals, Optics InfoBase Conference Papers Part F182-, 4 (2020).
- [38] E. Arqu´ e L´pez, V. Di Giulio, and F. J. Garc´ ıa de Abajo, o Atomic floquet physics revealed by free electrons, Phys. Rev. Res. 4 , 013241 (2022).
- [39] X. Wang, P. Garg, M. S. Mirmoosa, A. G. Lamprianidis, C. Rockstuhl, and V. S. Asadchy, Unleashing infinite momentum bandgap using resonant material systems, Arxiv , 1 (2023).
- [40] D. L. Mills and E. Burstein, Polaritons: The electromagnetic modes of media, Reports on Progress in Physics 37 , 817 (1974).
- [41] S. K´na-Cohen and S. R. Forrest, Room-temperature poe lariton lasing in an organic single-crystal microcavity, Nature Photonics 2010 4:6 4 , 371 (2010).
- [42] V. S. Gorelik and M. M. Sushchinski˘ ı, RAMAN SCATTERING OF LIGHT IN CRYSTALS, Soviet Physics Uspekhi 12 , 399 (1969).
- [43] J. D. Caldwell, A. V. Kretinin, Y. Chen, V. Giannini, M. M. Fogler, Y. Francescato, C. T. Ellis, J. G. Tischler, C. R. Woods, A. J. Giles, M. Hong, K. Watanabe, T. Taniguchi, S. A. Maier, and K. S. Novoselov, Sub-diffractional volume-confined polaritons in the natural hyperbolic material hexagonal boron nitride, Nature Communications 2014 5:1 5 , 1 (2014).
Mustafa Goksu Ozlu Vahagn Mkhitaryan Colton B. Fruhling Alexandra Boltasseva Vladimir M. Shalaev
Floquet Engineering of Polaritonic Amplification in Dispersive Photonic Time Crystals
In this study, we investigate the dynamics of dispersive photonic time
crystals (PTCs) and their potential applications for controlling light-matter
interaction. By employing the Lorentz-Drude model, we analyze theoretically and
via numerical simulation the effects of periodic modulation of dispersion
parameters, revealing the emergence of hybrid bandgaps from interaction of
polaritonic branches with unique characteristics. Our study demonstrates that
dispersive PTCs offer novel excitation channels and amplification
possibilities, that require lower modulation frequencies compared to
non-dispersive systems thus alleviating experimental challenges for the
realization of PTCs in the optical regime. These findings pave the way for
advancements in polaritonic lasing and resonant Raman scattering.