Floquet Damping Matrix: Concepts & Applications

• Floquet damping matrix is a finite-dimensional operator that captures decay, growth, and oscillatory modes via the spectral analysis of one-period propagators.
• It bridges applications from Floquet-Lindblad open quantum systems to classical non-Markovian and electronic friction models, revealing stroboscopic critical damping and exceptional transitions.
• Its construction via matrix logarithms unifies diverse methodologies, offering precise control over dissipative dynamics and enabling analysis of transient behaviors in driven systems.

Searching arXiv for the specified papers and closely related work on Floquet damping matrices and Floquet-Lindblad/open-system generators. I’ll retrieve arXiv records for the cited works and related papers before composing the article. arXiv search query: Floquet damping matrix open quantum systems Floquet Lindblad propagator exceptional points non-Markovian open system dynamics A Floquet damping matrix is a finite-dimensional object that encodes decay, growth, or dissipative mode structure in a periodically driven system over one driving period. In the cited literature, the term does not have a single universal meaning. It can denote a dissipative Floquet generator obtained from the matrix logarithm of a one-period open-system propagator, a Floquet electronic friction tensor governing Langevin or Fokker–Planck dynamics, or the logarithm of a monodromy matrix in Floquet–Lyapunov theory for linear systems with periodic coefficients (Mickiewicz et al., 11 Nov 2025, Wang et al., 2023, Landa et al., 2012, Gunderson et al., 2020). The common mathematical core is the replacement of time-local periodic dynamics by a stroboscopic operator over one period, followed by spectral analysis of that operator or its logarithm.

1. Terminological scope and formal setting

The cited works use closely related but non-identical constructions. The shared structure is a $T$-periodic evolution law and a one-period propagator whose spectrum determines long-time behavior.

Setting	Object	Definition or role
Periodic linear system	$R$	$R=\tfrac{1}{T}\ln \Phi(T)$
Floquet-Lindblad evolution	$\mathcal L_F$	Defined implicitly by $G(T)=e^{T\mathcal L_F}$
Non-Markovian driven open system	$G$	$Q_F=e^{GT}$
Floquet electronic friction	$\gamma_{\alpha\beta}(\mathbf R)$	Friction or damping tensor

In the linear periodic setting, the fundamental solution $\Phi(t)$ admits a Floquet–Lyapunov factorization

$$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$

with $R$0 the monodromy matrix. The eigenvalues of $R$1 separate growth or damping from oscillation by their real and imaginary parts (Landa et al., 2012).

In open quantum systems, the analogous object is constructed from a stroboscopic propagator on Liouville space or on an enlarged system–auxiliary space. The generator extracted by a matrix logarithm is the dissipative Floquet analogue of a Floquet Hamiltonian, but its eigenvalues generally have negative real parts rather than being purely imaginary (Mickiewicz et al., 11 Nov 2025).

A central terminological caution follows immediately. In the Floquet-Lindblad qubit study, the analysis is performed entirely in the time domain through the one-period propagator and an implicitly defined $R$2; that work does not introduce an explicit infinite-dimensional harmonic-space Floquet matrix or a truncated Fourier-space “damping matrix” $R$3 (Gunderson et al., 2020).

2. Stroboscopic generators and matrix logarithms

For a periodically driven open system in Lindblad form, vectorization of the density matrix yields

$R$4

where $R$5 is an $R$6 Liouvillian matrix,

$R$7

If $R$8, the one-period propagator is

$R$9

and the Floquet problem is formulated through

$R=\tfrac{1}{T}\ln \Phi(T)$0

The steady state satisfies $R=\tfrac{1}{T}\ln \Phi(T)$1, while the remaining $R=\tfrac{1}{T}\ln \Phi(T)$2 control decay toward it (Gunderson et al., 2020).

The same stroboscopic logic appears in exact non-Markovian open-system dynamics, but there the propagator is built on a combined system–bond space rather than directly from a Lindbladian. One defines a finite-dimensional $R=\tfrac{1}{T}\ln \Phi(T)$3 over one period and then sets

$R=\tfrac{1}{T}\ln \Phi(T)$4

where $R=\tfrac{1}{T}\ln \Phi(T)$5 are the eigenvalues of $R=\tfrac{1}{T}\ln \Phi(T)$6. The real parts of $R=\tfrac{1}{T}\ln \Phi(T)$7 are damping rates and the imaginary parts are Floquet quasi-frequencies modulo $R=\tfrac{1}{T}\ln \Phi(T)$8 (Mickiewicz et al., 11 Nov 2025).

The matrix logarithm is therefore not a formal ornament but the step that converts stroboscopic evolution into a generator-like object. The branch choice of $R=\tfrac{1}{T}\ln \Phi(T)$9 matters because quasi-frequencies are defined only modulo $\mathcal L_F$0 (Mickiewicz et al., 11 Nov 2025).

3. Floquet-Lindblad decay, Bloch reduction, and exceptional contours

A concrete Floquet-Lindblad realization is provided by a periodically driven qubit with a single dissipator, treated either as mode-selective spontaneous emission or as pure phase damping. In Bloch-vector form,

$\mathcal L_F$1

For mode-selective spontaneous emission with $\mathcal L_F$2, $\mathcal L_F$3, $\mathcal L_F$4, and

$\mathcal L_F$5

with $\mathcal L_F$6 in pure $\mathcal L_F$7 dissipation. For pure phase damping with $\mathcal L_F$8, $\mathcal L_F$9, $G(T)=e^{T\mathcal L_F}$0, and

$G(T)=e^{T\mathcal L_F}$1

with $G(T)=e^{T\mathcal L_F}$2 (Gunderson et al., 2020).

Exceptional points are detected from the eigenvectors of the one-period propagator rather than from a Fourier-expanded Floquet matrix. The criterion is based on the maximum inner product among the three trace-zero eigenmatrices,

$G(T)=e^{T\mathcal L_F}$3

with an exceptional point identified when $G(T)=e^{T\mathcal L_F}$4. By scanning drive or dissipation strength together with $G(T)=e^{T\mathcal L_F}$5, one obtains contours of exceptional points in parameter space. For sinusoidal $G(T)=e^{T\mathcal L_F}$6, the contours emerge periodically at subharmonic resonance $G(T)=e^{T\mathcal L_F}$7; for modulated $G(T)=e^{T\mathcal L_F}$8, they appear at $G(T)=e^{T\mathcal L_F}$9 (Gunderson et al., 2020).

The physical interpretation is stroboscopic critical damping. In the static Lindblad problem, an exceptional point marks the transition between underdamped and overdamped relaxation toward the steady state. In the periodically driven case, the same distinction is read off from the eigenvalues $G$0 of the one-period map: real $G$1 correspond to overdamped decay, while a complex-conjugate pair corresponds to underdamped oscillatory decay in the stroboscopic dynamics. The reported Floquet exceptional points persist down to vanishingly small dissipation when $G$2 is tuned to a multiphoton resonance, a feature absent in the static Lindblad spectrum (Gunderson et al., 2020).

4. Exact non-Markovian construction on system–bond space

In the strong-coupling, non-Markovian setting, the starting point is a small driven system with $G$3-periodic Hamiltonian $G$4, coupled bilinearly to a bosonic environment. Time is discretized as $G$5. The environment is traced out through the Feynman–Vernon influence functional, which is represented as a periodic matrix product operator. For a Gaussian bath, the undriven influence functional is compressed into matrices $G$6 of bond dimension $G$7, together with boundary vectors $G$8 (Mickiewicz et al., 11 Nov 2025).

Including the local system drive by a Trotter breakup yields micro-step dissipative propagators $G$9 acting on system-$Q_F=e^{GT}$0-bond space: $Q_F=e^{GT}$1 Periodic driving gives $Q_F=e^{GT}$2, so the one-period stroboscopic propagator is

$Q_F=e^{GT}$3

By construction, $Q_F=e^{GT}$4 is a finite-dimensional matrix of size $Q_F=e^{GT}$5, where $Q_F=e^{GT}$6, and it gives the exact stroboscopic open-system evolution including non-Markovian memory (Mickiewicz et al., 11 Nov 2025).

The Floquet damping matrix is then defined as

$Q_F=e^{GT}$7

or equivalently through the spectral decomposition

$Q_F=e^{GT}$8

with $Q_F=e^{GT}$9. Typically $\gamma_{\alpha\beta}(\mathbf R)$0 except for one steady-state mode with $\gamma_{\alpha\beta}(\mathbf R)$1, so the corresponding exponents $\gamma_{\alpha\beta}(\mathbf R)$2 have negative real parts. Their imaginary parts, taken modulo $\gamma_{\alpha\beta}(\mathbf R)$3, are the Floquet quasi-frequencies of damped modes (Mickiewicz et al., 11 Nov 2025).

This construction reduces continuously to the closed-system Floquet problem. In the limit $\gamma_{\alpha\beta}(\mathbf R)$4, the tensors $\gamma_{\alpha\beta}(\mathbf R)$5 become $\gamma_{\alpha\beta}(\mathbf R)$6, the $\gamma_{\alpha\beta}(\mathbf R)$7 become pure-unitary superoperators, and $\gamma_{\alpha\beta}(\mathbf R)$8 with

$\gamma_{\alpha\beta}(\mathbf R)$9

Accordingly,

$\Phi(t)$0

all eigenvalues lie on the unit circle, and the damping rates vanish (Mickiewicz et al., 11 Nov 2025).

5. Spectral information, heating, and transient entanglement

The utility of a Floquet damping matrix lies in what its spectrum organizes. In the driven spin–boson model,

$\Phi(t)$1

the steady-state eigenvector $\Phi(t)$2 of $\Phi(t)$3 with $\Phi(t)$4 gives the Floquet stationary density matrix $\Phi(t)$5. The period-averaged heat current density

$\Phi(t)$6

is obtained from the End–Matter formula for $\Phi(t)$7, with two-point correlations $\Phi(t)$8 reconstructed by interleaving the micro-motion operators $\Phi(t)$9 with the system-state vector in bond space. The spectral decomposition of $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$0 also yields the long-time decay rates of correlations, which enter the frequency dependence of the heat current (Mickiewicz et al., 11 Nov 2025).

For two driven qubits with

$$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$1

the spectrum of $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$2 contains one steady mode $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$3 and a ladder of damped modes with nonzero $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$4 and negative $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$5. In the undriven limit, the auxiliary propagator $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$6 alone has eigenmodes governing quench dynamics,

$$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$7

One transient mode $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$8 has an associated $$\Phi(t)=P(t)e^{Rt}, \qquad P(t+T)=P(t), \qquad R=\tfrac1T\ln \Phi(T),$$9 with high concurrence. Choosing $R$00 resonantly stabilizes this transient entangled mode in the Floquet steady state, and numerically the asymptotic concurrence from the leading eigenvector of $R$01 reaches a maximum $R$02 near $R$03 (Mickiewicz et al., 11 Nov 2025).

In the Floquet-Lindblad qubit problem, the spectrum of the one-period propagator plays a different but related role. There the relevant structure is not entanglement stabilization or heat transport but the boundary between underdamped and overdamped approach to the Lindblad steady state. The exceptional contours of the stroboscopic map therefore act as critical-damping manifolds in parameter space (Gunderson et al., 2020).

6. Floquet electronic friction, multidimensional damping tensors, and conceptual distinctions

A distinct usage of “Floquet damping matrix” appears in nonadiabatic dynamics near metal surfaces, where periodic driving of an Anderson–Holstein model leads first to a Floquet classical master equation and, in the fast-drive and fast-electron limits, to a Floquet Fokker–Planck equation. For the one-dimensional nuclear problem, the resulting effective equation is

$R$04

with mean force

$R$05

electronic friction coefficient

$R$06

and diffusion strength

$R$07

The equivalent Langevin dynamics is

$R$08

with $R$09 (Wang et al., 2023).

In multiple nuclear dimensions, the scalar friction becomes a tensor $R$10. The explicit Floquet expression is

$R$11

where $R$12, $R$13 are Bessel weights with $R$14, and the trace over electronic indices is implicit in the matrix product (Wang et al., 2023).

The physical distinction between this tensor and the stroboscopic logarithm $R$15 is substantial. The former is a friction kernel entering an effective nuclear Langevin description; the latter is a generator of one-period dissipative evolution for an open quantum system. Both encode damping, but they do so in different state spaces and at different levels of coarse graining.

Several common conceptual confusions are resolved by the cited literature. First, a Floquet damping matrix is not necessarily a harmonic-space block-Toeplitz construction. The Floquet-Lindblad qubit analysis extracts all spectral information directly from the finite one-period propagator and does not construct an infinite-dimensional harmonic-space matrix or any truncated $R$16 (Gunderson et al., 2020). Second, a Floquet damping matrix is not necessarily Markovian: the system–bond-space construction produces an exact non-Markovian $R$17 and generator $R$18 (Mickiewicz et al., 11 Nov 2025). Third, periodic driving can invalidate equilibrium fluctuation–dissipation structure. In the Floquet electronic friction formulation, $R$19 in general, so the second fluctuation–dissipation theorem is violated, and the driven nuclear degree of freedom heats up to an effective temperature $R$20 in the long-time limit (Wang et al., 2023).

From a computational standpoint, the required numerics depend on which object is meant. In Floquet electronic friction, one truncates sidebands $R$21, builds the Floquet Hamiltonian in extended space, computes $R$22 and $R$23, evaluates the double sum and energy integral, and assembles $R$24 at each geometry. In the exact non-Markovian quantum construction, one compresses the influence functional with an MPO method such as uniTEMPO, builds the periodic micro-step tensors $R$25, forms $R$26, and diagonalizes it or takes its logarithm to obtain the damping spectrum (Wang et al., 2023, Mickiewicz et al., 11 Nov 2025).