---
title: Eigenpulses in Time-Varying Media
url: https://www.emergentmind.com/topics/eigenpulses
type: topic
---

# Eigenpulses in Time-Varying Media

Eigenpulses are incident waveforms whose spectra are preserved under scattering by a time-varying medium, up to multiplication by a scalar reflection or transmission factor. In the literature on dispersive and ultrafast time-varying media, the term denotes eigenfunctions of scattering operators or of spectrally resolved reflection matrices, so that the output retains the input spectral shape even though the medium exchanges energy with the wave [2208.11778] [2508.12753]. In exact treatments of moving impedance profiles, the same term is used for specially shaped broadband pulses that enter a finite modulated section without reflection, providing a pulse-level analogue of a transmission resonance [2404.09075]. Across these settings, eigenpulses generalize the role of monochromatic modes in static media to temporally modulated systems.

## 1. Definition and core eigenvalue problems

For a single-mode line terminated by a time-varying impedance, the input–output relation can be written in the time domain as
$$
\bar v(t)=\int_{-\infty}^{t} {\rm d}t'\;\bar R\bigl(t,t-t'\bigr)\,\bar u(t').
$$
After Fourier transformation, the reflected spectrum obeys
$$
v(\omega)=\int_{-\infty}^{\infty}\frac{d\omega'}{2\pi}\,R(\omega,\omega')\,u(\omega').
$$
For periodic time modulation with period $T$, this becomes the discrete matrix equation
$$
v_m=\sum_n R_{mn}u_n,\qquad v=Ru,
$$
where the sampled frequencies are separated by $2\pi/T$. The eigenpulses are then defined by the spectral integral equation
$$
\int R(\omega,\omega')\,a(\omega')\,d\omega'=\lambda\,a(\omega),
$$
or, in the discrete basis,
$$
\sum_n R_{mn}\,a_n=\lambda\,a_m.
$$
When $R(\omega,\omega')$ is square-integrable, the kernel is Hilbert–Schmidt, so one obtains a countable set of eigenvalues $\lambda_k$ and orthonormal eigenfunctions $a_k(\omega)$; after discretization over $N$ frequency points, this reduces to a standard $N\times N$ complex matrix eigenproblem [2508.12753].

The associated time-domain waveform is
$$
a(t)=\int_{-\infty}^{\infty} a(\omega)\,e^{-i\omega t}\,d\omega.
$$
By construction, injecting $a(\omega)$ yields
$$
v(\omega)=\lambda\,a(\omega),\qquad v(t)=\lambda\,a(t),
$$
so the reflection is spectrally distortion-free: the pulse shape, including its relative spectral phases, is unchanged apart from an overall scale factor [2508.12753].

A parallel formulation appears in operator-based scattering theory for dispersive time-varying media. There, a reflection eigenpulse satisfies
$$
\hat r\,\psi(\omega)=r\,\psi(\omega),
$$
which is equivalent to
$$
\hat Z\,\psi(\omega)=z\,\psi(\omega),\qquad r=\frac{1-z}{1+z},
$$
with normalization
$$
\int d\omega\,|\psi(\omega)|^2=1.
$$
Any such $\psi$ is an eigenpulse because its reflected spectrum is $r\,\psi(\omega)$; the same eigenproblem also determines transmission eigenpulses, with transmitted amplitude $t=2/(1+z)$ [2208.11778].

## 2. Operator, matrix, and exact-solution frameworks

The three main formulations appearing in the current literature are complementary rather than competing. One is operator-theoretic, one is exact for a specific class of space–time-varying impedances, and one is directly experimental through measured reflection matrices.

| Framework | Setting | Eigenpulse characterization |
|---|---|---|
| Operator scattering | Dispersive time-varying media | Eigenvectors of $\hat r$ or $\hat Z$ |
| Exact moving-impedance solution | Constant $n$, impedance $Z(x-vt)$ | Periodic condition $G(w+2d)=G(w)$ |
| Measured reflection matrix | Ultrafast periodically modulated termination | Eigenvectors of $R$ |

In the operator formulation, the medium response is described by a two-time conductivity kernel,
$$
\mathbf{j}(t)=\int_{-\infty}^{\infty}dt'\;\sigma\bigl(t,t-t'\bigr)\,\mathbf E(t'),
$$
which becomes an operator-valued conductivity $\hat\sigma(-i\partial_\omega,\omega)$ in the frequency domain. Maxwell’s equations reduce to an operator Helmholtz equation, and continuity at an interface yields operator-valued Fresnel relations such as
$$
\hat r_s=(1-\hat Z_s)(1+\hat Z_s)^{-1},\qquad \hat t_s=2(1+\hat Z_s)^{-1}.
$$
The transverse-magnetic case has the same form with $\hat Z_p$ in place of $\hat Z_s$. For a finite slab, one similarly obtains explicit operator expressions for $\hat r_{\rm slab}$ and $\hat t_{\rm slab}$, after which diagonalization or singular-value decomposition identifies reflection eigenpulses and maximal-transmission pulses [2208.11778].

In the exact moving-impedance treatment, the refractive index is constant while the impedance depends on the co-moving coordinate $u=x-vt$. Introducing
$$
u=x-vt,\qquad w=x+vt,
$$
reduces Maxwell’s equations to a form admitting closed-form solutions in terms of arbitrary functions $g(w)$ and $h(u)$. For a finite slab $0\le x\le d$, matching to vacuum at both interfaces completely fixes the reflected and transmitted waves in terms of the incident pulse [2404.09075].

In the experimental matrix formulation, the medium is characterized by a measured spectrally resolved reflection matrix $R_{mn}$, and eigenpulses are obtained numerically from the matrix eigenproblem $Ra_k=\lambda_k a_k$ using standard linear-algebra routines such as MATLAB `eig` or Python `numpy.linalg.eig` [2508.12753].

## 3. Broadband reflectionless pulses and the moving-impedance picture

For a finite section of a medium with constant refractive index and impedance profile $Z(x-vt)$, Li and Horsley identified exact broadband reflectionless eigenpulses. The key condition is that the auxiliary function $G(w)$ be periodic with period $2d$,
$$
G(w+2d)=G(w).
$$
Under this condition, the net reflection from the two interfaces cancels exactly, yielding eigenpulses that enter without reflection. The paper characterizes this as an analogue of a transmission resonance, but at the level of pulse shaping rather than discrete monochromatic frequencies [2404.09075].

This pulse-level resonance generalizes the familiar Fabry–Pérot condition. In a uniform slab, zero reflection occurs at discrete frequencies satisfying $2\omega_mnd/c=2\pi m$. In the moving-impedance problem, the analogous condition is not a single frequency constraint but the broadband pulse-shaping condition $G(w+2d)=G(w)$, which makes the field on the front face at time $t$ identical to that on the back face at the delayed time $t+2d/v$. The two partial reflections therefore cancel for all time [2404.09075].

Any periodic profile $G'(w)$ with period $2d$ admits the Fourier series
$$
G'(w)=\sum_{m\in\mathbb Z} g_m\,e^{-\,i\,\frac{m\pi}{d}\,w},
$$
so the incident eigenpulse acquires a comb-like spectrum with lines at
$$
\omega_m=\frac{m\pi v}{d},
$$
weighted by the Fourier transform of the impedance-envelope prefactor. By contrast, a non-eigen incident pulse excites both co- and counter-propagating waves, and the rapidly varying impedance mixes frequencies in the reflected field. For weak, rapid modulations of the form
$$
Z(u)=\eta_0\bigl[1+\delta\eta\cos(\kappa u)\bigr],\qquad |\delta\eta|\ll 1,
$$
the leading reflected amplitude carries a factor $\delta\eta'\sim\kappa\delta\eta$, corresponding to an up-shift of the spectrum by $\omega\to\omega\pm\kappa v$ [2404.09075].

A common misconception is to equate eigenpulses with ordinary continuous waves. The moving-impedance results make clear that broadband pulse trains can be the exact reflectionless solutions, and that their defining property is a spatiotemporal matching condition rather than monochromaticity [2404.09075].

## 4. Experimental characterization in ultrafast time-varying media

An experimental realization of eigenpulses in an ultrafast time-varying medium was reported using a ring resonator on a PCB with two varactor diodes in opposite gaps, inductively coupled to a coaxial transmission line via a small loop antenna. Reverse-biasing the diodes sets a static resonance, while fast modulation of that bias produces a response function $R(t,\tau)$ with strong coupling between many frequencies [2508.12753].

Two classes of control-voltage modulation were used: a sinusoidal modulation at $1\,{\rm MHz}$, and a “random” modulation formed as a sum of 100 tones from $1\to100\,{\rm MHz}$ with random phases and amplitudes. All modulations were periodic with $T=1\,\mu{\rm s}$, giving a frequency grid $\Delta\omega=2\pi\cdot1\,{\rm MHz}$ [2508.12753].

The reflection matrix was measured sequentially. Continuous-wave tones were injected at $f_n=1,2,\ldots,1000\,{\rm MHz}$ in $1\,{\rm MHz}$ steps using an AWG. For each injected tone, the time trace of the reflected signal was recorded over one period $T$ using a fast oscilloscope. The resulting matrix $R_{t,n}$ was corrected by subtracting a $\pi$-phase-shifted measurement to remove spurious bias-modulation leakage, and a Fourier transform along the time axis produced the spectrally resolved matrix $R_{m,n}\equiv R(\omega_m,\omega_n)$ [2508.12753].

Under the $1\,{\rm MHz}$ sinusoidal modulation, the measured $R$ displayed a near-diagonal “ladder” of sidebands corresponding to first-, second-, and higher-order scattering. The eigenvalues formed a quasi-continuous band, and each eigenvector $a_k(\omega)$ spanned several sidebands. When three representative eigenpulses were injected, the measured reflected spectra perfectly overlapped $\lambda_k a_k(\omega)$, and the time-domain waveforms $a_k(t)$ and $\lambda_k a_k(t)$ were identical up to scale. Under the “random” modulation, $R$ became a fully dense matrix, yet two sample eigenpulses still reflected without spectral distortion [2508.12753].

These measurements establish that the eigenpulse concept does not depend on a sparse sideband structure. A plausible implication is that the relevant object is the full frequency-coupling operator itself, whether it is near-diagonal or fully dense.

## 5. Power extremization, absorption, and spectral focusing

The measured reflection matrix also supports optimization problems beyond distortion-free reflection. To find pulses that extremize the total reflected power, one forms $R^\dagger R$ and diagonalizes it. Its right singular vectors are the relevant incident waveforms, with the top singular vector maximizing total reflected power and the bottom singular vector minimizing it. In the reported experiment, the minimally reflected pulse yielded broadband absorption greater than $98\%$ [2508.12753].

Selective spectral concentration is obtained by partitioning the reflected spectrum into masked bands $A$ and $B$, defining submatrices $R_A$ and $R_B$, and constructing the contrast operator
$$
C=(R_B^\dagger R_B)^{-1}R_A^\dagger R_A.
$$
The top eigenvector of $C$ injects a pulse whose reflected energy is concentrated into the user-defined bands $A$; swapping $A$ and $B$ concentrates the reflection into $B$ instead. The reported contrast exceeded $10^4$ [2508.12753].

Within the terminology of the paper, these optimized inputs are frequency eigenchannels of the ultrafast time-varying medium. The eigenpulse problem identifies waveforms that are invariant in spectral shape, while the singular-vector and contrast-operator problems identify waveforms that extremize global reflected power or redistribute that power into selected spectral windows. This suggests that eigenpulse analysis is part of a broader linear-algebraic control framework for dynamic scattering [2508.12753].

## 6. Physical interpretation, bound states, and prospective applications

The physical intuition offered across the literature is that eigenpulses are “temporal wavefronts.” In a static medium, the natural distortion-free modes are monochromatic continuous waves. In a time-varying medium, by contrast, the incident pulse spectrum must be pre-shaped so that temporal modulation and dispersion cancel. This is explicitly compared to spatial wavefront shaping and eigenchannels in multimode fibres [2508.12753].

The corresponding eigenvalues encode qualitatively distinct regimes. In the experimental reflection-matrix picture, $|\lambda_k|\le 1$ when there is loss, interpreted as absorption into unguided channels. The same framework also identifies waveforms for which $|\lambda|>1$, associated with parametric amplification, and waveforms with $|\lambda|\ll1$, associated with broadband absorption [2508.12753]. In the operator theory, reflection and transmission amplitudes are set by the impedance eigenvalue through $r=(1-z)/(1+z)$ and $t=2/(1+z)$, so eigenpulse behavior is directly linked to the spectral properties of $\hat Z$ [2208.11778].

Operator poles introduce a distinct but related phenomenon. In dispersive time-varying media, the poles of the reflection or transmission operators represent non-time-harmonic bound states. Since
$$
\hat r=(1-\hat Z)(1+\hat Z)^{-1},
$$
poles occur when
$$
\det(1+\hat Z)=0,
$$
equivalently when there exists a nontrivial solution of
$$
(1+\hat Z)\phi(\omega)=0,\qquad \hat Z\phi=-\phi.
$$
These $\phi(\omega)$ are described as surface-plasmon-like eigenpulses confined to the interface and decaying evanescently away from it [2208.11778].

Potential applications listed in the experimental work include distortion-free communications through dynamic channels by encoding in the eigenbasis, dynamic spectral filtering and focusing, temporal dispersion compensation in optics and microwaves, secure key-exchange keyed by knowledge of the time-varying reflection matrix, and inverse design of time-varying media to realize a target set of eigenpulses [2508.12753]. In the moving-impedance setting, the strong asymmetry between co-propagating and counter-propagating waves further suggests one-way, spectrum-reshaping mirrors, while the exact zero-reflection solutions provide a broadband analogue of resonance-based impedance matching [2404.09075].

Taken together, these results place eigenpulses at the intersection of operator scattering theory, exact space–time electromagnetics, and experimentally measured frequency-eigenchannel control. They are not merely pulses that survive propagation unchanged in a trivial sense; rather, they are the specific spectral states selected by a time-varying scattering operator, a moving-impedance boundary-value problem, or a measured ultrafast reflection matrix, and they therefore provide a natural modal language for dynamic wave–matter interactions.

Source: https://www.emergentmind.com/topics/eigenpulses