Stroboscopic Dynamics

Stroboscopic dynamics observes a periodically driven system at one fixed phase of every drive cycle:

$$t_n=t_0+nT, \qquad n\in\mathbb Z.$$

At those times, continuous evolution collapses to repeated application of the one-period unitary

$$U_F(t_0) = U(t_0+T,t_0).$$

This reduction is exact at the sampling times. It is also incomplete: it discards motion within each period. Understanding both statements is essential when interpreting Floquet effective Hamiltonians, numerical maps, kicked systems, pulse cycles, and digital quantum simulations.

Floquet Operators owns the one-period spectrum, logarithm branches, and numerical construction. Quantum Maps and Discrete-Time Evolution owns abstract repeated unitary maps. This page owns the sampling viewpoint and the information it preserves or loses.

Stroboscopic Map

For a $T$-periodic Hamiltonian,

$$H(t+T)=H(t),$$

propagator covariance gives

$$U(t_0+nT,t_0) = U_F(t_0)^n.$$

The sampled pure states are

$$\lvert\psi_n\rangle = U_F^n \lvert\psi_0\rangle.$$

Density operators obey the discrete unitary channel

$$\rho_n = U_F^n \rho_0 \left( U_F^* \right)^n.$$

It is useful to name the one-step superoperator:

$$\mathcal U_F(\rho) = U_F\rho U_F^*,$$

so that

$$\rho_n = \mathcal U_F^n(\rho_0).$$

For a closed system, $\mathcal U_F$ preserves trace, positivity, entropy, and purity. It does not generically drive states toward one attractor.

The sampling phase $t_0$ is part of the definition. Sampling at a different drive phase gives a unitarily conjugate Floquet operator, but the state and measurement representatives at that phase also change.

Stroboscopic Observables

In the discrete Heisenberg description,

$$A_n = \left( U_F^* \right)^n A U_F^n.$$

Expectation values agree:

$$\operatorname{Tr}(\rho_nA) = \operatorname{Tr}(\rho_0A_n).$$

The same construction defines sampled correlation functions. For example,

$$C_{AB}(n) = \operatorname{Tr} \left[ \rho_0A_nB \right].$$

If an observable is fixed by the one-period map,

$$U_F^*AU_F=A,$$

then its stroboscopic expectation is conserved. This is weaker than continuous conservation. The observable may vary strongly during each period and return only at the sampling phase.

Likewise, a density operator satisfying

$$U_F\rho U_F^*=\rho$$

is stroboscopically stationary. It need not commute with $H(t)$ at every instant.

Spectral Form of Sampled Motion

In finite dimension, choose Floquet eigenvectors

$$U_F \lvert\phi_\alpha\rangle = e^{-i\theta_\alpha} \lvert\phi_\alpha\rangle.$$

Expand

$$\lvert\psi_0\rangle = \sum_\alpha c_\alpha \lvert\phi_\alpha\rangle.$$

Then

$$\lvert\psi_n\rangle = \sum_\alpha c_\alpha e^{-in\theta_\alpha} \lvert\phi_\alpha\rangle.$$

For an observable $A$,

$$\langle A\rangle_n = \sum_{\alpha,\beta} c_\alpha^*c_\beta e^{in(\theta_\alpha-\theta_\beta)} A_{\alpha\beta}.$$

Sampled oscillations are controlled by eigenphase differences modulo $2\pi$. Degenerate eigenphases give stationary coherences; nondegenerate differences give discrete-time interference and recurrences.

An eigenstate of $U_F$ returns to the same ray after each period. A superposition generally does not, because its components accumulate different eigenphases.

Micromotion Versus Stroboscopic Evolution

Let

$$0\leq\tau\lt T.$$

At a fixed intra-period offset, exact composition gives

$$\begin{aligned} U(t_0+nT+\tau,t_0) &= U(t_0+nT+\tau,t_0+nT) U_F^n\\ &= U(t_0+\tau,t_0) U_F^n. \end{aligned}$$

Define the intra-period propagator

$$M(\tau;t_0) = U(t_0+\tau,t_0).$$

Then

$$\lvert\psi(t_0+nT+\tau)\rangle = M(\tau;t_0) U_F^n \lvert\psi_0\rangle.$$

The Floquet operator controls cycle-to-cycle evolution; $M(\tau;t_0)$ reconstructs the motion inside a cycle. Knowing $U_F$ alone does not determine $M(\tau;t_0)$.

This is the operational meaning of micromotion. Measurements locked to $\tau=0$ see only $U_F^n$. Measurements at another fixed phase see a conjugated or dressed observable. Measurements with timing jitter can mix micromotion into what was intended to be a stroboscopic signal.

Same Map, Different Micromotion

Consider two protocols on a period $T$.

The first has

$$H_a(t)=0,$$

so

$$U_{F,a}=I.$$

The second has

$$H_b(t) = \begin{cases} A, & 0\leq t\lt T/2, \\ -A, & T/2\leq t\lt T. \end{cases}$$

Because the two segments commute,

$$U_{F,b} = e^{iAT/(2\hbar)} e^{-iAT/(2\hbar)} =I.$$

Both protocols are identical at integer periods. Inside the period they are different: the first state is constant, while the second evolves away from its initial value and then returns.

This example shows why a stroboscopic effective Hamiltonian is not a full reconstruction of the laboratory drive.

Frequency Aliasing

Suppose a continuous signal contains

$$e^{-i\nu t}.$$

At the stroboscopic times,

$$e^{-i\nu t_n} = e^{-i\nu t_0} e^{-in\nu T}.$$

Replacing $\nu$ by

$$\nu'=\nu+m\Omega, \qquad m\in\mathbb Z,$$

does not change the sampled sequence because

$$e^{-in(\nu+m\Omega)T} = e^{-in\nu T} e^{-i2\pi mn} = e^{-in\nu T}.$$

Sampling once per period therefore identifies frequencies modulo $\Omega$. This is temporal aliasing.

Quasienergy modularity is the spectral version of the same fact:

$$\varepsilon \sim \varepsilon+m\hbar\Omega.$$

A stroboscopic record cannot determine which representative generated an eigenphase without additional continuity, bandwidth, or microscopic information.

Sampling at multiple intra-period phases can recover some micromotion harmonics. It does not remove all inverse-problem ambiguities automatically.

Effective Hamiltonian

Choose a logarithm branch and write

$$U_F = e^{-iH_FT/\hbar}.$$

Then

$$U_F^n = e^{-iH_FnT/\hbar},$$

so $H_F$ exactly generates the selected stroboscopic sequence.

The logarithm is multi-valued. If

$$U_F = \sum_\alpha e^{-i\theta_\alpha} \lvert\phi_\alpha\rangle \langle\phi_\alpha\rvert,$$

then

$$H_F = \sum_\alpha \frac{\hbar(\theta_\alpha+2\pi m_\alpha)}{T} \lvert\phi_\alpha\rangle \langle\phi_\alpha\rvert$$

generates the same $U_F$ for any integers $m_\alpha$.

Additional criteria may select a useful branch:

• continuity as the drive is turned on;
• locality or quasilocality;
• a chosen quasienergy zone;
• matching to a high-frequency expansion;
• agreement with known weak-drive limits.

Even after a branch is chosen, $H_F$ determines only the stroboscopic part. A micromotion operator is still required for general times.

Sampling Phase

Let

$$U_F(t_0) = U(t_0+T,t_0).$$

At another phase $t_1$,

$$U_F(t_1) = U(t_1,t_0) U_F(t_0) U(t_1,t_0)^*.$$

Thus the eigenphases are unchanged, but the eigenvectors and effective-Hamiltonian representative are dressed by micromotion.

A laboratory stroboscopic protocol should specify:

1. the drive phase used as $t_0$;
2. the measurement offset $\tau$;
3. whether state preparation is repeated at that phase;
4. timing jitter and finite measurement duration;
5. whether a reported observable is in the laboratory or a rotating frame.

Without this information, two nominally identical stroboscopic datasets may correspond to different sampled observables.

Examples

Slow dynamics under a fast carrier

Near-resonant two-level control often separates a rapid carrier rotation from a slow Rabi envelope. Sampling once per carrier period suppresses the visible carrier phase and reveals the slow cycle-to-cycle motion. A rotating frame explains the same simplification continuously.

The two descriptions are complementary. A rotating frame retains a continuous slow trajectory; stroboscopic sampling retains only selected points of the laboratory trajectory.

Repeated pulse cycle

A pulse sequence with one-cycle unitary

$$U_{\rm cyc} = U_m\cdots U_2U_1$$

has stroboscopic evolution $U_{\rm cyc}^n$. Reordering pulses generally changes the cycle map when the generators do not commute.

Cycle-boundary measurements can diagnose accumulated errors, while within-cycle measurements reveal transient populations and toggling-frame motion.

Delta-kicked rotor

For a rotor receiving an impulsive potential once per period, the cycle map factorizes exactly into a kick and a free-rotation unitary. Kicked Rotor Preview uses this map to compare classical momentum diffusion with quantum interference, dynamical localization, and resonant ballistic growth.

Subharmonic return

If

$$U_F^q\lvert\psi_0\rangle = e^{i\chi}\lvert\psi_0\rangle$$

for some integer $q\gt1$, the ray returns after $qT$ rather than $T$. Such subharmonic behavior can occur in finite systems by spectral commensurability. Robust many-body subharmonic response requires additional physics and should not be inferred from a short finite-time oscillation alone.

Quantum Simulation Connection

Digital quantum simulation naturally produces a repeated cycle:

$$U_{\rm step} = U_m\cdots U_1.$$

After $n$ cycles,

$$\lvert\psi_n\rangle = U_{\rm step}^n \lvert\psi_0\rangle.$$

This is stroboscopic dynamics whether the cycle represents:

• one Trotter step approximating a target Hamiltonian;
• an engineered Floquet period;
• a compiled gate layer;
• a discrete-time quantum walk;
• a kicked-map model.

The interpretation matters. For Trotterization, the within-step gate sequence is an implementation detail and the cycle map approximates a target continuous propagator. For Floquet engineering, the micromotion may be physical and intentionally designed.

Sampling only at cycle boundaries can conceal coherent within-cycle errors. Conversely, boundary sampling can isolate the effective long-time map from large but reversible micromotion.

The error bounds for product formulas belong to Trotter Product Formula. The direct one-period series belongs to Floquet–Magnus Expansion, while the broader inverse-frequency and prethermal analysis belongs to High-Frequency Expansions.

Scope and Limitations

The unitary stroboscopic framework does not directly describe:

• periodic open-system channels or Liouvillians;
• measurement-conditioned maps;
• noisy cycle-to-cycle timing or control fluctuations;
• quasiperiodic drives with no common period;
• long-time many-body heating without further analysis;
• continuous observables reconstructed from one sampled phase alone.

For a periodic open system, one can still define a one-period quantum channel. Its spectrum can contain decay rates and attractors, unlike a finite-dimensional unitary Floquet map. That belongs to open-system dynamics.

Common Mistakes

• Treating stroboscopic equality as equality throughout a period.
• Omitting the sampling phase $t_0$.
• Assuming $U_F$ uniquely determines micromotion.
• Assuming a logarithm of $U_F$ is unique.
• Interpreting aliased frequencies or quasienergies as absolute.
• Calling a stroboscopically fixed observable continuously conserved.
• Ignoring timing jitter or finite measurement duration.
• Confusing a physical Floquet cycle with a numerical Trotter step.
• Inferring robust subharmonic order from a short recurrence.
• Applying unitary-map intuition to a dissipative one-period channel.

References

• J. H. Shirley, “Solution of the Schrödinger equation with a Hamiltonian periodic in time,” Physical Review 138, B979–B987, 1965.
• H. Sambe, “Steady states and quasienergies of a quantum-mechanical system in an oscillating field,” Physical Review A 7, 2203–2213, 1973.
• F. Haake, Quantum Signatures of Chaos, 3rd ed., Springer, 2010.
• M. Bukov, L. D’Alessio, and A. Polkovnikov, “Universal high-frequency behavior of periodically driven systems,” Advances in Physics 64, 139–226, 2015.
• A. Eckardt, “Colloquium: Atomic quantum gases in periodically driven optical lattices,” Reviews of Modern Physics 89, 011004, 2017.
• M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.

Exercises

1. Derive the exact reconstruction formula at a fixed intra-period offset.

Solution

Use composition:

$$\begin{aligned} U(t_0+nT+\tau,t_0) &= U(t_0+nT+\tau,t_0+nT) U(t_0+nT,t_0). \end{aligned}$$

Periodicity gives

$$U(t_0+nT+\tau,t_0+nT) = U(t_0+\tau,t_0),$$

and stroboscopic evolution gives

$$U(t_0+nT,t_0) = U_F^n.$$

Therefore

$$U(t_0+nT+\tau,t_0) = M(\tau;t_0)U_F^n.$$

2. Derive the sampled expectation value in the Floquet eigenbasis.

Solution

With

$$\lvert\psi_0\rangle = \sum_\alpha c_\alpha \lvert\phi_\alpha\rangle$$

and

$$U_F\lvert\phi_\alpha\rangle = e^{-i\theta_\alpha} \lvert\phi_\alpha\rangle,$$

the sampled state is

$$\lvert\psi_n\rangle = \sum_\alpha c_\alpha e^{-in\theta_\alpha} \lvert\phi_\alpha\rangle.$$

Thus

$$\begin{aligned} \langle A\rangle_n &= \langle\psi_n\vert A\vert\psi_n\rangle\\ &= \sum_{\alpha,\beta} c_\alpha^*c_\beta e^{in(\theta_\alpha-\theta_\beta)} A_{\alpha\beta}. \end{aligned}$$

3. Prove that sampling once per period aliases frequencies modulo $\Omega$.

Solution

At $t_n=t_0+nT$,

$$e^{-i(\nu+m\Omega)t_n} = e^{-i(\nu+m\Omega)t_0} e^{-in(\nu+m\Omega)T}.$$

The cycle-to-cycle factor is

$$e^{-in(\nu+m\Omega)T} = e^{-in\nu T} e^{-i2\pi mn} = e^{-in\nu T}.$$

The overall prefactor at $t_0$ depends on the sampling phase, but the discrete evolution cannot distinguish $\nu$ from $\nu+m\Omega$.

4. Construct two protocols with the same Floquet operator and different micromotion.

Solution

Take $H_a(t)=0$, so $U_{F,a}=I$ and the state never moves.

For the second protocol, use

$$H_b(t) = \begin{cases} A, & 0\leq t\lt T/2, \\ -A, & T/2\leq t\lt T. \end{cases}$$

Then

$$U_{F,b} = e^{iAT/(2\hbar)} e^{-iAT/(2\hbar)} =I.$$

At $T/2$, however,

$$U_b(T/2,0) = e^{-iAT/(2\hbar)},$$

which is generally not the identity. The two protocols agree stroboscopically and differ inside the period.

5. Explain the difference between a Floquet cycle and a Trotter step.

Solution

A Floquet cycle is the physical period of a time-periodic Hamiltonian. Its intra-period micromotion can be observable and may be part of the intended dynamics.

A Trotter step is usually a numerical or digital approximation to evolution under a target Hamiltonian. The within-step gate ordering creates approximation error relative to that target. Repeating the step produces a stroboscopic approximation:

$$U_{\rm target}(n\Delta t) \approx U_{\rm step}^n.$$

The mathematics of repeated unitaries is shared, but the physical interpretation and error question differ.