Correspondence Approach to Time-Dependent Decaying Quantized Fields

1. Introduction

The decay of electromagnetic (EM) fields is a central theme in electrodynamics, with manifestations ranging from charge relaxation in conductors and magnetic diffusion to radiation damping and cavity losses [1,2,3]. Classically, the functional form of the decay; exponential, diffusive, power-law, oscillatory damped, stretched exponential, or multi-rate; is fixed by the material response, boundary conditions, and geometry [4]. In many applications the decay envelope is not merely a nuisance parameter; it is an experimentally accessible fingerprint of loss channels, transport mechanisms, and mode structure. A unified way of translating such classical decay laws into a quantum description is therefore useful both conceptually and as a modeling tool.

When the field is quantized, dissipation is most often modeled by coupling the mode of interest to a thermal reservoir of harmonic oscillators, including the $T=0$ limit [5,6,7]. Under the Markovian (memoryless) approximation this construction yields a Lindblad master equation with a constant damping rate $\gamma$, a structure that guarantees complete positivity and underpins a large fraction of textbook quantum-optical models. This widely used limit is, however, tightly linked to an exponentially decaying classical amplitude, $\Gamma \left(t\right)\propto exp(-\gamma t/2)$; it effectively builds the exponential envelope into the generator. As a result, describing decay laws that are time-dependent or non-exponential typically requires either a microscopic non-Markovian treatment or a phenomenological generalization whose physical interpretation and consistency must be checked.

Many realistic systems deviate strongly from a single exponential. Disordered and complex dielectrics often display stretched-exponential Kohlrausch–Williams–Watts relaxation [8,9]; open resonators and nanophotonic structures can exhibit multi-mode quasi-normal-mode decay with several competing rates [10,11]; and transient radiation fields may develop power-law tails at late times [1]. These examples motivate a practical and conceptual question that we address here: given a prescribed classical envelope $\Gamma \left(t\right)$ for the field amplitude, what is the quantum master equation that reproduces it, and under what conditions does the resulting evolution remain of (time-local) Lindblad type?

We answer this question using a correspondence-principle strategy tailored to bosonic modes. The starting point is that a coherent state $|\alpha \rangle$ provides the closest quantum analogue of a classical complex amplitude; we therefore represent a decaying classical mode by the coherent-state trajectory $|\alpha \Gamma (t)\rangle$. By imposing that the density operator $\rho \left(t\right)=|\alpha \Gamma (t)\rangle \langle \alpha \Gamma (t\left)\right|$ is an exact solution of the quantum dynamics, we obtain a time-local master equation whose rate is fixed directly by the chosen envelope $\Gamma \left(t\right)$. This construction recovers the standard amplitude-damping Lindblad master equation when $\Gamma \left(t\right)$ is exponential, and it yields a transparent generalization for time-dependent and non-exponential damping scenarios.

The manuscript is organized as follows. Section 2 summarizes representative classical decay regimes and the corresponding envelopes $\Gamma \left(t\right)$. Section 3 derives the associated time-local master equation, relates it to the GKSL form and CP-divisibility conditions, and notes the limitations that arise when $\Gamma \left(t\right)$ crosses zero. Section 4 discusses how the correspondence rule can be used as a practical dictionary, and outlines extensions to finite temperature and multi-mode settings. Section 5 presents the conclusions.

2. Classical Decay Laws of Electromagnetic Fields

The temporal evolution of electromagnetic fields is governed by Maxwell’s equations [1,3]. The functional form of any observed decay, however, is not universal; it depends on the constitutive response of the medium, the relevant length and time scales, boundary conditions, and geometry [1,2]. For the correspondence approach pursued in the next section, the key classical object is an effective amplitude envelope $\Gamma \left(t\right)$ that captures how the complex field amplitude of a given mode decreases in time. The purpose of this section is therefore not to provide an exhaustive review, but to identify representative decay laws and express them in a form that can be summarized by a mode envelope $\Gamma \left(t\right)$.

2.1. Exponential Decay: Conductive Relaxation

A paradigmatic example is Ohmic relaxation in a homogeneous conductor. For finite conductivity $\sigma$ and permittivity $\epsilon$, free charge relaxes exponentially and the associated fields decay on the time scale $\tau :=\epsilon /\sigma$ [1,2,3]:

$$\mathbf{E}\left(t\right)=\mathbf{E}\left(0\right)e^{-t/\tau },\tau :={\displaystyle \frac{\epsilon }{\sigma }}.$$

(1)

In terms of our envelope language, a single-mode amplitude may be idealized as $\Gamma \left(t\right)=exp(-t/\tau )$, or, in quantum-optics convention for energy damping, $\Gamma \left(t\right)=exp(-\gamma t/2)$. Exponential decay also arises as a good approximation for a long-lived mode of a high-Q resonator, where Q is the quality factor, when radiation leakage can be characterized by a single rate [10].

2.2. Diffusive Decay: Magnetic Field Penetration

When magnetic fields penetrate a conductor, the relevant dynamics is often diffusive. In a simple geometry one obtains ${\partial }_t\mathbf{B}=D{\nabla }^2\mathbf{B}$ with $D:={\left(\mu \sigma \right)}^{-1}$; for a step-like boundary excitation, the field at depth x behaves as [1,2]

$$B(x,t)\sim {\displaystyle \frac{1}{\sqrt{t}}}{e}^{-x^2/\left(4Dt\right)}.$$

(2)

At fixed x this implies an algebraic envelope $\Gamma \left(t\right)\propto t^{-1/2}$ for the peak amplitude, rather than a constant-rate exponential. Such diffusive decay underlies the skin effect and magnetic field penetration in conductors and superconductors [12].

2.3. Oscillatory Damped Decay: Lossy Dielectrics and Conductors

Lossy media generically support damped waves with an oscillatory carrier and a decaying envelope. For a monochromatic plane wave the field can be written as

$$\mathbf{E}(\mathbf{r},t)={\mathbf{E}}_0{e}^{-\kappa z}cos(kz-\omega t),$$

(3)

where $\kappa$ is the attenuation coefficient. In the time domain, a localized wave packet therefore exhibits an exponentially decaying envelope multiplying an oscillatory carrier. In the present envelope language one may identify $\Gamma \left(t\right)=exp(-\gamma t)$ for the amplitude damping, while the oscillation $cos(\omega t+\varphi )$ represents a phase rotation of the underlying complex amplitude. A minimal time-domain model is

$$E\left(t\right)\sim e^{-\gamma t}cos(\omega t+\varphi ).$$

(4)

This form is relevant for pulses in lossy dielectrics and conductors, and for under-damped cavity modes [1,13].

2.4. Power-Law Decay: Radiation from Transient Sources

Transient sources can generate late-time radiation tails that decay algebraically. For a localized charge-current distribution that is switched on or off over a finite time, the radiation-zone field at late times often follows a power law; a representative case is [1,14]

$$E\left(t\right)\sim {\displaystyle \frac{1}{t}}(t\to \infty ).$$

(5)

More generally, the exponent is set by the lowest non-vanishing multipole moment of the transient source; in envelope language this corresponds to $\Gamma \left(t\right)\propto t^{-p}$ with $p>0$ [1]. Power-law tails also appear in related wave and quantum-decay contexts, including photonic-band-gap emission and late-time open-system dynamics [15,16,17].

2.5. Stretched Exponential: Disordered and Complex Media

In many disordered materials; glasses, polymers, porous media, and biological tissues; relaxation is broadly distributed in time and is not well described by a single Debye relaxation time. A common empirical description is the Kohlrausch–Williams–Watts (KWW) function for the susceptibility [8,9,18,19]

$$\chi \left(t\right)\propto e^{-{(t/\tau )}^{\beta }},0<\beta <1,$$

(6)

which interpolates between an exponential ($\beta =1$) and a slow, broad relaxation ($\beta <1$). For step-like driving, the field or polarization can be modeled with an envelope of the form $\Gamma \left(t\right)=exp[-{(t/\tau )}^{\beta }/2]$ [8,20]. The KWW law can be viewed as a superposition of Debye-type exponentials with a distribution of relaxation times, reflecting heterogeneous energy landscapes and anomalous transport [21,22].

2.6. Multi-Mode Exponential Sum: Quasi-Normal Modes

Open systems; leaky cavities, antennas, and plasmonic resonators; are often inherently multi-modal. In that setting, the field can be expanded in quasi-normal modes (QNMs) with complex eigenfrequencies ${\tilde{\omega }}_n:={\omega }_n-i{\gamma }_n/2$ [10,11,23]:

$$\mathbf{E}\left(t\right)=\sum _n{c}_n{e}^{-i{\tilde{\omega }}_{n}t}=\sum _n{c}_n{e}^{-{\gamma }_{n}t/2}{e}^{-i{\omega }_{n}t}.$$

(7)

Each term has an exponential envelope, but the overall decay is a superposition of rates; equivalently, one may view the effective envelope $\Gamma \left(t\right)$ as a multi-rate sum that can become non-exponential when several QNMs contribute appreciably [23,24].

Taken together, these examples show that the envelope $\Gamma \left(t\right)$ can be exponential, algebraic, stretched exponential, oscillatory damped, or a multi-rate sum. In the next section we use $\Gamma \left(t\right)$ as the classical input and derive a quantum master equation whose dissipative rate is determined directly by the chosen envelope.

3. Correspondence Principle and the Quantum Master Equation

In this section we turn the classical envelope $\Gamma \left(t\right)$ into a quantum dynamical law for the density operator [5,6,25]. The goal is to construct a time-local master equation whose solution reproduces the correspondence ansatz $\rho \left(t\right)=|\alpha \Gamma (t)\rangle \langle \alpha \Gamma (t\left)\right|$ for an initially coherent field [26,27], and to express the resulting dissipative rate directly in terms of $\Gamma \left(t\right)$. We also highlight the connection to Lindblad generators and the conditions under which the dynamics retains a Lindblad-like structure [28,29].

At the classical level, a single field mode is conveniently characterized by a complex amplitude $\alpha \left(t\right)$ whose modulus and argument encode, respectively, the field strength and phase. We parameterize damping by introducing an envelope $\Gamma \left(t\right)$ normalized at the initial time, $\Gamma \left(0\right):=1$, so that the amplitude can be written as

$$\alpha \left(t\right)=\alpha \left(0\right)\Gamma \left(t\right)\equiv \alpha \Gamma \left(t\right),$$

(8)

where $\alpha :=\alpha \left(0\right)$ may be complex. In smooth dissipative scenarios $\Gamma \left(t\right)$ is typically positive and decaying, which leads to a regular rate $\gamma \left(t\right)$ and a GKSL generator; in oscillatory damped cases $\Gamma \left(t\right)$ may cross zero and the resulting time-local rate must be interpreted piecewise. In this parametrization, phase rotations reside in the complex phase of $\alpha \left(t\right)$, while $\Gamma \left(t\right)$ controls the damping entering the time-local generator derived below. Representative envelopes include:

• Exponential relaxation: $\Gamma \left(t\right)=exp(-\gamma t/2)$;
• Algebraic decay: $\Gamma \left(t\right)={(1+t/{\tau }_0)}^{-1}$;
• Stretched-exponential relaxation: $\Gamma \left(t\right)=exp[-{(t/\tau )}^{\beta }/2]$.

To pass from this classical description to a quantum one, we invoke the correspondence principle in its quantum-optical form: a coherent state $|\alpha \rangle$ provides the closest quantum analogue of a classical mode with complex amplitude $\alpha$ [26,27]. We therefore postulate that the quantized counterpart of a classically decaying mode follows the coherent-state trajectory $|\alpha \Gamma (t)\rangle$, and we require the associated density operator

$$\rho \left(t\right)=|\alpha \Gamma (t)\rangle \langle \alpha \Gamma (t\left)\right|$$

(9)

to be an exact solution of a time-local master equation. The problem is then: determine the generator that produces this evolution for an arbitrary admissible envelope $\Gamma \left(t\right)$. With this correspondence ansatz in place, we now derive the master equation explicitly.

3.1. Derivation

Our goal is to determine a time-local generator ${\mathcal{L}}_t$ such that $\dot{\rho }\left(t\right)={\mathcal{L}}_t\left[\rho \left(t\right)\right]$ reproduces the correspondence trajectory $\rho \left(t\right)=|\alpha \Gamma (t)\rangle \langle \alpha \Gamma (t\left)\right|$. We will show that ${\mathcal{L}}_t$ can be cast in GKSL form with a single jump operator a, i.e., ${\mathcal{L}}_t\left[\rho \right]=\gamma \left(t\right)\mathcal{D}\left[a\right]\rho$, where $\mathcal{D}\left[a\right]\rho :=a\rho a^{†}-\frac{1}{2}\{a^{†}a,\rho \}$. We proceed constructively by computing the first-order increment $\rho (t+h)-\rho \left(t\right)$ and then taking the limit $h\to 0$.

We start from the coherent-state expansion in the Fock basis,

$$|\alpha \Gamma \left(t\right)\rangle =e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)/2}{\displaystyle \sum _{n=0}^{\infty }}{\displaystyle \frac{{[\alpha \Gamma \left(t\right)]}^n}{\sqrt{n!}}}|n\rangle ,$$

(10)

which yields the density operator

$$\rho \left(t\right)=e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}{\displaystyle \sum _{n,m=0}^{\infty }}|n\rangle \langle m|{\displaystyle \frac{{\alpha }^n{\alpha }^{*m}{\Gamma }^{n+m}\left(t\right)}{\sqrt{n!m!}}}.$$

(11)

To compute $\rho (t+h)$ to first order in h, we expand each factor that depends on t. Using $\Gamma (t+h)=\Gamma \left(t\right)+h\dot{\Gamma }\left(t\right)+\mathcal{O}\left(h^2\right)$ and retaining only $\mathcal{O}\left(h\right)$ terms gives the set of expansions

$$\begin{array}{cc}\hfill \Gamma (t+h)& \approx \Gamma \left(t\right)+h\dot{\Gamma }\left(t\right),\hfill \end{array}$$

(12a)

$$\begin{array}{cc}\hfill {\Gamma }^{n+m}(t+h)& \approx {\Gamma }^{n+m}\left(t\right)+(n+m)h{\Gamma }^{n+m}\left(t\right){\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}},\hfill \\ \hfill e^{-{\left|\alpha \right|}^2{\Gamma }^2(t+h)}& \approx e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}\left[1-{2h\left|\alpha \right|}^2\Gamma \left(t\right)\dot{\Gamma }\left(t\right)\right]\hfill \end{array}$$

(12b)

$$\begin{array}{cc}\hfill & =e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}\left[1-{2h\left|\alpha \right|}^2{\Gamma }^2\left(t\right){\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\right].\hfill \end{array}$$

(12c)

Substituting Eqs. (12a)–(12c) into Eq. (11) and collecting all $\mathcal{O}\left(h\right)$ contributions yields

$$\begin{array}{c}\rho (t+h)\approx e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}\left[1-{2h\left|\alpha \right|}^2{\Gamma }^2\left(t\right){\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\right]\hfill \\ \hfill \times \sum _{n,m}|n\rangle \langle m|{\displaystyle \frac{{\alpha }^n{\alpha }^{*m}{\Gamma }^{n+m}\left(t\right)}{\sqrt{n!m!}}}\left[1+h(n+m){\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\right].\end{array}$$

(13)

The expression (13) is still written componentwise in the Fock basis. To identify a compact operator form, we relate the factors n and m in the sum to number-operator actions. Using $a|n\rangle =\sqrt{n}|n-1\rangle$ and $a^{†}|n\rangle =\sqrt{n+1}|n+1\rangle$, one finds the operator identities

$$\begin{array}{cc}\hfill a\rho \left(t\right)a^{†}& ={\left|\alpha \right|}^2{\Gamma }^2\left(t\right)\rho \left(t\right),\hfill \end{array}$$

(14a)

$$\begin{array}{cc}\hfill a^{†}a\rho \left(t\right)& =e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}\sum _{n,m}n|n\rangle \langle m|{\displaystyle \frac{{\alpha }^n{\alpha }^{*m}{\Gamma }^{n+m}\left(t\right)}{\sqrt{n!m!}}},\hfill \end{array}$$

(14b)

$$\begin{array}{cc}\hfill \rho \left(t\right)a^{†}a& =e^{-{\left|\alpha \right|}^2{\Gamma }^2\left(t\right)}\sum _{n,m}m|n\rangle \langle m|{\displaystyle \frac{{\alpha }^n{\alpha }^{*m}{\Gamma }^{n+m}\left(t\right)}{\sqrt{n!m!}}}.\hfill \end{array}$$

(14c)

In particular, Eqs. (14b) and (14c) show that the combination $(n+m)$ appearing in Eq. (13) corresponds to $a^{†}a\rho +\rho a^{†}a$, while Eq. (14a) accounts for the term proportional to ${\left|\alpha \right|}^2{\Gamma }^2\left(t\right)$.

Using these identities and expanding Eq. (13) to first order, we obtain the finite-difference form

$$\rho (t+h)\approx \rho \left(t\right)+h{\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\left[a^{†}a\rho \left(t\right)+\rho \left(t\right)a^{†}a-2a\rho \left(t\right)a^{†}\right].$$

(15)

Dividing by h and taking $h\to 0$ gives the master equation

$$\dot{\rho }={\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\left(a^{†}a\rho +\rho a^{†}a-2a\rho a^{†}\right).$$

(16)

This expression can be rewritten in the standard dissipator form

$$\dot{\rho }=\gamma \left(t\right)\left(a\rho a^{†}-\frac{1}{2}\{a^{†}a,\rho \}\right),\gamma \left(t\right):=-2{\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}.$$

(17)

Here $\{A,B\}:=AB+BA$ denotes the anticommutator. The right-hand side of Eq. (17) has the standard Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) structure [28,29]: it consists of a single dissipative channel with jump operator a and rate $\gamma \left(t\right)$, i.e., $\dot{\rho }=\gamma \left(t\right)\mathcal{D}\left[a\right]\rho$ with $\mathcal{D}\left[L\right]\rho :=L\rho L^{†}-\frac{1}{2}\{L^{†}L,\rho \}$. When $\gamma \left(t\right)\ge 0$ for all t, the generator is GKSL at each instant and the evolution is completely positive (CP) and trace preserving; in that case the dynamical map is completely-positive-divisible (CP-divisible), often taken as a Markovianity criterion in the divisibility sense, even though $\gamma \left(t\right)$ may be time dependent [25,30]. The condition $\gamma \left(t\right)\ge 0$ is equivalent to $\dot{\Gamma }\left(t\right)\le 0$, namely a monotonically decaying envelope.

For the special case $\Gamma \left(t\right)=exp(-\gamma t/2)$ with constant $\gamma >0$, Eq. (17) reduces to the familiar amplitude-damping (optical master) equation [5,6,7]. If $\gamma \left(t\right)$ becomes negative over some interval, Eq. (17) is no longer of GKSL form at that time and the dynamics is not CP-divisible; physically, this corresponds to a temporary reversal of the instantaneous damping rate, a common signature of non-Markovian behavior [25,30].

3.2. Properties and Consistency Checks

The master equation (17) was engineered to reproduce the coherent-state trajectory $|\alpha \Gamma (t)\rangle$; nevertheless, it is useful to make explicit a few structural properties that any physically meaningful dissipative evolution should satisfy. First, Eq. (16) is trace preserving: taking the trace and using cyclicity, $Tr\left(a\rho a^{†}\right)=Tr\left(a^{†}a\rho \right)$, yields $Tr\dot{\rho }=0$, so $Tr\rho \left(t\right)$ is conserved. Moreover, when $\gamma \left(t\right)\ge 0$ the generator in Eq. (17) has the Gorini–Kossakowski–Sudarshan–Lindblad structure with a single jump operator a, and therefore preserves Hermiticity and defines a completely positive, trace-preserving (CPTP) dynamical map [25,28,29].

The correspondence content is most transparent at the level of field moments. By construction, $\rho \left(t\right)=|\alpha \Gamma (t)\rangle \langle \alpha \Gamma (t\left)\right|$ is an exact solution of Eq. (16) for any differentiable envelope $\Gamma \left(t\right)$, and the first moment follows the classical amplitude law,

$${\langle a\rangle }_t=Tr\left(a\rho \left(t\right)\right)=\alpha \Gamma \left(t\right).$$

(18)

For the same reason, the mean photon number $\overline{n}\left(t\right)=Tr\left(a^{†}a\rho \right)$ obeys a simple rate equation. For an initial coherent state one has $\overline{n}\left(t\right)={\left|\alpha \right|}^2{\Gamma }^2\left(t\right)$, and differentiation gives

$${\displaystyle \frac{d\overline{n}}{dt}}=2{\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}\overline{n}\left(t\right)=-\gamma \left(t\right)\overline{n}\left(t\right),$$

(19)

which matches the classical energy decay ${\left|\alpha \left(t\right)\right|}^2={\left|\alpha \right|}^2{\Gamma }^2\left(t\right)$.

Finally, the time dependence encoded in $\Gamma \left(t\right)$ controls the divisibility properties of the dynamics. When $\Gamma \left(t\right)$ is not a simple exponential, the decay rate $\gamma \left(t\right)$ becomes time dependent. If $\gamma \left(t\right)\ge 0$ for all t, the evolution remains GKSL at each instant and is CP-divisible; if $\gamma \left(t\right)$ becomes negative over some interval, the dynamics is not CP-divisible, which is a standard indicator of non-Markovian behavior [25,30].

3.3. Examples

The correspondence rule provides a direct recipe: once a classical envelope $\Gamma \left(t\right)$ is specified, the quantum dynamics follows from $\gamma \left(t\right):=-2\dot{\Gamma }\left(t\right)/\Gamma \left(t\right)$. For $\gamma \left(t\right)\ge 0$, Eq. (17) generates a bosonic pure-loss (amplitude-damping) channel with transmissivity $\eta \left(t\right):={\Gamma }^2\left(t\right)$ [5,6]. In the standard input–output picture, this channel is equivalent to mixing the system mode with a vacuum environment mode on a beam splitter of transmissivity $\eta \left(t\right)$; consequently, the field amplitude and energy obey ${\langle a\rangle }_t=\Gamma \left(t\right){\langle a\rangle }_0$ and ${\langle a^{†}a\rangle }_t=\eta \left(t\right){\langle a^{†}a\rangle }_0$, with added vacuum noise required by the commutation relations [5,6]. The coherent-state trajectory emphasized in this work is the special case in which the state remains pure and the only effect is the deterministic rescaling $\alpha \mapsto \alpha \Gamma \left(t\right)$. Below we illustrate how different classical envelopes translate into distinct time-dependent damping rates and physical regimes.

3.3.1. Exponential Decay

For $\Gamma \left(t\right)=exp(-\gamma t/2)$ one has $\dot{\Gamma }/\Gamma =-\gamma /2$ and therefore $\gamma \left(t\right)=\gamma$ (constant). Eq. (17) reduces to the standard optical master equation,

$$\dot{\rho }=\gamma \left(a\rho a^{†}-\frac{1}{2}{a}^{†}a\rho -\frac{1}{2}\rho a^{†}a\right),$$

(20)

which is time-homogeneous and defines a dynamical semigroup [28,29]. Physically, this corresponds to memoryless energy loss at a constant rate, such as cavity damping into a broadband continuum [5,6]. In this case $\eta \left(t\right)=exp(-\gamma t)$ and every normally ordered moment decays exponentially, e.g., ${\langle {\left(a^{†}\right)}^m{a}^n\rangle }_t={\eta }^{(m+n)/2}\left(t\right){\langle {\left(a^{†}\right)}^m{a}^n\rangle }_0$. For an initial coherent state $|\alpha \rangle$, the solution remains coherent, $\rho \left(t\right)=|\alpha exp(-\gamma t/2)\rangle \langle \alpha exp(-\gamma t/2\left)\right|$.

3.3.2. Power-Law Decay

For an algebraic envelope $\Gamma \left(t\right)={(1+t/{\tau }_0)}^{-p}$ with $p>0$, the instantaneous damping rate becomes

$$\gamma \left(t\right)={\displaystyle \frac{2p}{{\tau }_0+t}},$$

(21)

which decreases in time. The corresponding transmissivity $\eta \left(t\right)={(1+t/{\tau }_0)}^{-2p}$ has no single characteristic time scale; instead, the decay slows down progressively, producing long tails. Such behavior is consistent with situations in which the effective loss rate is controlled by diffusion-limited transport, broad distributions of coupling strengths, or other mechanisms that invalidate a constant-rate approximation. In this case the coherent amplitude scales as $\alpha \left(t\right)=\alpha {(1+t/{\tau }_0)}^{-p}$ and the mean photon number decays as $\overline{n}\left(t\right)={\left|\alpha \right|}^2{(1+t/{\tau }_0)}^{-2p}$. Related algebraic laws appear in non-exponential quantum decay and anomalous-transport settings [16,17,21].

3.3.3. Stretched Exponential

For the stretched-exponential envelope $\Gamma \left(t\right)=exp[-{(t/\tau )}^{\beta }/2]$ with $0<\beta <1$, one finds

$$\gamma \left(t\right)={\displaystyle \frac{\beta }{\tau }}{\left({\displaystyle \frac{t}{\tau }}\right)}^{\beta -1}.$$

(22)

The divergence of $\gamma \left(t\right)$ as $t\to 0^{+}$ indicates that the decay is initially faster than exponential, a feature commonly associated with a broad spectrum of microscopic relaxation times rather than a single Debye constant. At long times $\gamma \left(t\right)\to 0$, producing an ever-slower tail and an effective loss that weakens with age. The photon number follows $\overline{n}\left(t\right)={\left|\alpha \right|}^{2}exp[-{(t/\tau )}^{\beta }]$, reproducing the KWW relaxation at the level of the quantized field [8,9].

4. Discussion

The time-dependent decay laws of electromagnetic fields are not universal; rather, they emerge from Maxwell’s equations under specific physical conditions and therefore encode information about material response, geometry, and coupling to external degrees of freedom. The central point of this work is that, once such a classical decay law is distilled into an effective envelope $\Gamma \left(t\right)$ for a mode amplitude, it can be promoted to a quantum description by the correspondence rule

$$\gamma \left(t\right):=-2{\displaystyle \frac{\dot{\Gamma }\left(t\right)}{\Gamma \left(t\right)}}.$$

(23)

This rule fixes the time-local generator in Eq. (17); the classical information is encoded in the logarithmic derivative $\dot{\Gamma }\left(t\right)/\Gamma \left(t\right)$. When $\Gamma \left(t\right)$ is smooth and strictly positive, $\gamma \left(t\right)$ is regular; if $\Gamma \left(t\right)$ crosses zero (as in some oscillatory damped signals), $\gamma \left(t\right)$ becomes singular at the nodes and the corresponding time-local description must be interpreted piecewise in time.

Table 1 collects representative examples and illustrates how it should be read. The second column lists a convenient closed-form $\Gamma \left(t\right)$; the third column gives the corresponding instantaneous rate obtained by differentiating and applying Eq. (23). When $\gamma \left(t\right)\ge 0$ the dynamics is of GKSL type at each time and therefore completely positive and trace preserving, whereas negative intervals of $\gamma \left(t\right)$ signal departures from CP-divisibility. The table is meant as a dictionary: one can start from a classical decay observed in experiment or simulation, infer an effective $\Gamma \left(t\right)$, and then write down the associated quantum master equation.

Two remarks are in order. First, in the diffusive and power-law rows the decay rate is not constant but decreases with time, reflecting the absence of a single relaxation time scale. Second, the multi-mode QNM row emphasizes that even if each contribution decays exponentially, interference between several rates can produce an overall envelope that is effectively non-exponential; in that situation the correspondence rule still applies once a single-mode effective $\Gamma \left(t\right)$ has been identified.

Beyond the coherent-state correspondence used for the derivation, Eq. (17) is linear in $\rho$ and therefore governs the evolution of arbitrary initial states; in particular, it applies to nonclassical inputs such as squeezed states, Fock states, or Schrödinger-cat states [6]. Within its domain of validity, the approach provides a compact phenomenological alternative to fully microscopic reservoir models: it translates a measured or simulated classical envelope into a time-local quantum master equation (equivalently, a time-local generator) while keeping complete positivity under explicit control through the sign of $\gamma \left(t\right)$.

Several extensions follow naturally. At finite temperature, the zero-temperature dissipator can be replaced by the standard thermal Born–Markov (weak-coupling, memoryless) master equation, leading to time-dependent absorption and emission terms with rates determined by $\Gamma \left(t\right)$ and the bath occupation [5,6]. For multi-mode fields, a classical superposition $\mathbf{E}(\mathbf{r},t)={\sum }_{\mathbf{k}\lambda }{\mathcal{E}}_{\mathbf{k}\lambda }{\alpha }_{\mathbf{k}\lambda }{\Gamma }_{\mathbf{k}\lambda }\left(t\right)e^{i\mathbf{k}·\mathbf{r}}+\mathrm{complex}\mathrm{conjugate}(\mathrm{c}.\mathrm{c}.)$ suggests mode-resolved envelopes and, in the absence of inter-mode coupling, independent master equations for each mode. When coupling between modes is relevant, $\Gamma \left(t\right)$ becomes matrix-valued and the correspondence procedure points toward multi-channel GKSL generators with several jump operators.

5. Conclusions

We have presented a correspondence-principle framework that promotes classical decay laws of electromagnetic modes to quantum dissipative dynamics. Starting from a classical envelope $\Gamma \left(t\right)$ for the complex field amplitude, we derived a time-local master equation whose generator is fixed by the logarithmic derivative $\gamma \left(t\right)-2\dot{\Gamma }\left(t\right)/\Gamma \left(t\right)$, yielding Eq. (17).

A key outcome is that the classical information is encoded entirely in $\Gamma \left(t\right)$; different envelopes produce different time-dependent damping rates and therefore distinct physical regimes. When $\gamma \left(t\right)\ge 0$, the generator has GKSL structure at each time and the resulting evolution is a bosonic pure-loss channel with transmissivity $\eta \left(t\right)={\Gamma }^2\left(t\right)$; in particular, coherent states remain coherent with amplitude scaled by $\Gamma \left(t\right)$ and the photon number decays as $\overline{n}\left(t\right)={\left|\alpha \right|}^2{\Gamma }^2\left(t\right)$. The examples discussed; exponential, power-law, and stretched-exponential decay; illustrate how non-exponential classical relaxation leads to non-constant quantum loss rates.

More broadly, the correspondence approach provides a compact route from a measured or simulated classical decay envelope to a quantum master equation, complementing microscopic reservoir constructions. It offers a unified language for dissipative conductive media, disordered dielectrics, and open resonators, and it can be extended to finite temperature and multi-mode settings by introducing additional channels and rates consistent with the underlying envelopes.