Multi-Time Quantum-Field Correlations

• Multi-Time Quantum-Field Correlations are defined by evaluating quantum field operators at distinct times with combined normal and time ordering, capturing nonclassical temporal processes.
• The methodology uses multitime characteristic functions and Fourier-based regularization to reconstruct singular quasiprobabilities in practical experiments.
• Time ordering, nonequal-time commutators, and open-system dynamics are key for understanding memory effects and interference patterns in diverse quantum applications.

Searching arXiv for recent and foundational papers on multi-time quantum-field correlations. Multi-time quantum-field correlations are correlations of quantum fields evaluated at distinct times, typically through normally and time-ordered products of Heisenberg-picture field operators or through detector probabilities built from unequal-time correlation functions. In quantum optics and related areas, the problem is not merely to characterize the field state at a single instant, but to determine whether the full temporal process generated by operators such as $\{\hat a(t_i)\}$ can be represented by a classical stochastic process, how time ordering and nonequal-time commutators modify that process, and how such correlations can be reconstructed or regularized when the underlying phase-space objects are singular. A central line of work formulates these questions through a time-dependent $P$ functional and its multitime characteristic function, yielding necessary-and-sufficient nonclassicality criteria for arbitrary numbers of times (Krumm et al., 2016). Closely related developments regularize the singular multitime $P$ functional into a smooth quasiprobability accessible by correlated homodyne measurements (Krumm et al., 2017). Beyond quantum optics, the same subject interfaces with detector-based formulations of quantum field theory in curved spacetime, process-tensor approaches to non-Markovian open systems, and multi-time wave-function formulations of QFT (Anastopoulos et al., 2019, Jørgensen et al., 2020, Petrat et al., 2013).

1. Multitime quantum fields as temporal stochastic processes

The modern phase-space formulation begins by asking whether the family of field amplitudes sampled at several times can be understood as a classical stochastic process. In the single-time case, nonclassicality is tied to the Glauber–Sudarshan distribution

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$

whose failure to be a genuine probability density signals nonclassicality. The multitime problem generalizes this from a state property at one instant to a property of the whole temporal process $\{\hat a(t_i)\}_{i=1}^k$, with physically relevant observables given by normally and time-ordered products rather than equal-time moments alone (Krumm et al., 2016).

For $k$ times $t_1,\dots,t_k$, the multitime $P$ functional is defined as

$$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$

Here the ordering symbol denotes combined normal and time ordering, with creation operators sorted with increasing time arguments from left to right and annihilation operators with decreasing time arguments, as in photocounting theory (Krumm et al., 2016). This object plays the role of a formal joint quasiprobability for coherent amplitudes $\alpha_i\equiv \alpha(t_i)$ at several times (Krumm et al., 2017).

Its operational meaning is fixed by

$P$0

Thus the multitime $P$1 functional encodes the full hierarchy of multitime normally ordered field correlations. Nonclassicality of the temporal process is then defined by the existence of an operator function $P$2 such that

$P$3

which excludes any classical stochastic-process interpretation (Krumm et al., 2016).

This formulation is specific to quantum fields in several ways. The basic observables are bosonic field operators, the ordering prescription is the normal-and-time ordering relevant to photodetection theory, and nonequal-time commutators and time-ordered unitary evolution enter essentially. The same papers stress that this is not merely a temporal analogue of equal-time multimode optics; the multitime setting can generate singularities stronger than those known from equal-time multimode phase-space descriptions (Krumm et al., 2017).

2. Characteristic functions, Bochner hierarchies, and nonclassicality tests

Because the multitime $P$4 functional is typically highly singular, the technically central object is its Fourier transform, the multitime characteristic function

$P$5

with

$P$6

For ordered times $P$7, it can be written as

$P$8

The decisive practical point is that $P$9 remains well behaved even when the $P$0 functional itself is singular, and its derivatives at the origin generate normally and time-ordered moments (Krumm et al., 2016, Krumm et al., 2017).

The multitime nonclassicality criteria are built by expanding a general test operator in a finite Fourier series of multitime displacement operators,

$P$1

which yields the quadratic form

$P$2

where $P$3 is a Bochner-type matrix whose entries are characteristic-function values at phase-space differences (Krumm et al., 2016). By Bochner’s theorem and Sylvester’s criterion, the temporal process is nonclassical exactly when some finite principal minor becomes negative: $P$4 This gives an infinite hierarchy of necessary and sufficient nonclassicality probes for temporal quantum correlations (Krumm et al., 2016).

At the lowest nontrivial order $P$5, the determinant condition reduces to

$P$6

hence the simple sufficient witness

$P$7

This is the direct multitime analogue of the single-time condition $P$8 (Krumm et al., 2016). The formal significance is that multitime nonclassicality becomes an experimentally accessible property of a regular characteristic function rather than a singular quasiprobability.

A distinct but related detector-based route appears in the deferred-measurement formulation of time correlators. There, intermediate von Neumann projections are replaced by entanglement with quantum memories and a final readout. In the weak-measurement limit, the irreducible detector probabilities yield detector-weighted combinations of

$P$9

while in the strong-measurement limit the projected correlator is recovered exactly (Oehri et al., 2015). This reinforces a general point: multitime quantum correlations are not uniquely specified without an ordering and measurement prescription.

3. Time ordering, nonequal-time commutators, and dynamical structure

A defining feature of multi-time quantum-field correlations is that time ordering is not optional. Since Heisenberg operators $$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$0 at different times need not commute, the field evolution itself depends on the time-ordered unitary

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$1

whenever

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$2

This dependence propagates directly into multitime characteristic functions and into the temporal evolution of nonclassicality (Krumm et al., 2016).

To isolate time-ordering corrections, the Magnus expansion is used: $$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$3 with

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$4

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$5 is the evolution without time ordering, while all higher terms are time-ordering corrections (Krumm et al., 2016). In explicitly time-dependent interactions, those corrections can change the rise, suppression, or oscillation of multitime nonclassicality.

The optical benchmark is degenerate parametric down-conversion with frequency mismatch,

$$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$6

For $$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$7, time ordering is irrelevant; for $$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$8, it becomes nontrivial (Krumm et al., 2016, Krumm et al., 2017). The Heisenberg equations can be reduced to a $$P(\alpha)=\langle \mathord{:} \hat \delta (\hat a - \alpha) \mathord{:} \rangle,$$9 matrix problem, and the Magnus series converges if

$\{\hat a(t_i)\}_{i=1}^k$0

For this model, $\{\hat a(t_i)\}_{i=1}^k$1, so with $\{\hat a(t_i)\}_{i=1}^k$2, convergence is guaranteed for $\{\hat a(t_i)\}_{i=1}^k$3 (Krumm et al., 2016).

A multitime-specific indicator of the underlying dynamical algebra is the nonequal-time commutator. In the parametric model it takes the form

$\{\hat a(t_i)\}_{i=1}^k$4

and is generally nonzero for time-dependent Hamiltonians (Krumm et al., 2016). In the regularization analysis, the corresponding central commutator produces an extra time-ordering factor in the two-time characteristic function,

$\{\hat a(t_i)\}_{i=1}^k$5

which can yield factors proportional to $\{\hat a(t_i)\}_{i=1}^k$6 whose strength grows with $\{\hat a(t_i)\}_{i=1}^k$7 (Krumm et al., 2017). This shows explicitly that time ordering can change the singular structure of multitime phase-space objects in ways without an equal-time analogue.

A different unequal-time phenomenon appears in quantum electrodynamics of source fields. For a stationary dipole, the Heisenberg-picture source field itself is retarded, but certain unequal-time vacuum-subtracted field correlations contain advanced-wave-like contributions through mixed vacuum-source commutators. These advanced terms occur in two-time correlators, remain consistent with Einstein causality, and can affect radiative-force statistics even though they do not significantly contribute to ordinary photodetection amplitudes in the vacuum state (Stokes, 2017). The broader implication is that causal field propagation does not by itself determine the structure of unequal-time observables.

4. Regularized quasiprobabilities and experimental reconstruction

The central obstacle to direct phase-space visualization is that the multitime $\{\hat a(t_i)\}_{i=1}^k$8 functional can be more singular than the single-time multimode $\{\hat a(t_i)\}_{i=1}^k$9 function. In equal-time multimode optics one has the standard bound

$k$0

but in the multitime case non-equal-time commutators may be operator-valued and noncentral, BCH simplifications can fail, and time ordering can induce stronger-than-Gaussian growth (Krumm et al., 2017). This motivates a universal regularization.

The solution is a compact-support filter applied directly to the multitime characteristic function: $k$1 leading to a smooth regularized quasiprobability

$k$2

The filter is chosen as

$k$3

with triangular function

$k$4

Because the filter has compact support, it suppresses any growth of $k$5, including stronger-than-Gaussian growth, and its Fourier transform is a bona fide probability density, so filtering does not create negativities (Krumm et al., 2017). Negativities of $k$6 remain necessary and sufficient indicators of genuine temporal nonclassicality.

For the parametric process with mismatch, the filtered two-time quasiprobability is smooth and exhibits negativities, providing a direct phase-space visualization of temporal quantum correlations (Krumm et al., 2017). The same paper also gives a concrete reconstruction scheme using two correlated balanced homodyne detection setups, one probing the field at time $k$7 and the other at $k$8. In the strong-local-oscillator limit, the two-dimensional Fourier transform of the correlated difference statistics gives direct access to the multitime characteristic function: $k$9 Using the product filter, the regularized $t_1,\dots,t_k$0 is then sampled by pattern functions and estimated empirically from measured quadrature data (Krumm et al., 2017).

This reconstruction program is significant because it moves multitime nonclassicality from a purely formal criterion in Fourier space to a directly samplable quasiprobability in phase space. It also clarifies a common misconception: while equal-time quantities such as $t_1,\dots,t_k$1 isolate selected aspects of temporal structure, the regularized multitime $t_1,\dots,t_k$2 functional is designed to encode the full nonclassical temporal process rather than one chosen moment hierarchy (Krumm et al., 2017).

5. Open systems, continuous monitoring, and operational temporal processes

In open quantum systems, multitime correlations are not determined in general by single-time reduced dynamics. A two-time correlator already depends on how the environment stores memory across interventions, which is why the ordinary quantum regression theorem can fail in non-Markovian regimes (Jørgensen et al., 2020). The process-tensor formalism addresses this by encoding all multi-time statistics in a Choi-state object $t_1,\dots,t_k$3, with probabilities of intervention records given by a spatio-temporal Born rule

$t_1,\dots,t_k$4

Within this framework, generalized transfer tensors $t_1,\dots,t_k$5 provide a discrete memory-kernel decomposition of the full multi-time dynamics, and three-time transfer tensors $t_1,\dots,t_k$6 capture irreducible memory across an inserted operation (Jørgensen et al., 2020).

A recent development shows that for time-independent Hamiltonians and finite memory times $t_1,\dots,t_k$7, there is an exact factorization rule for multitime correlations in non-Markovian open quantum systems. Higher-order correlations reduce to products of lower-order correlations whenever adjacent operator insertions are separated by more than $t_1,\dots,t_k$8, and all information required to reconstruct $t_1,\dots,t_k$9-time correlations is contained in a temporal volume of $P$0 (Bracht et al., 21 May 2026). In the example of quantum dots coupled to phonons, this yields exact agreement with brute-force calculations while avoiding the breakdown of the standard QRT (Bracht et al., 21 May 2026).

Continuous-measurement theory provides a complementary operational perspective. For a qubit monitored continuously by linear detectors, the output currents

$P$1

define multitime correlators

$P$2

For unital ensemble evolution and no phase backaction, these factorize into sequential two-time blocks: $P$3 and similarly with an extra one-time mean for odd $P$4 (Atalaya et al., 2017). The paper stresses that this is not a Wick theorem or Gaussian factorization, but a special structural property of the monitored output record in the unital informational-backaction regime (Atalaya et al., 2017).

At a more general operational level, sequential invasive measurements can be reformulated without intermediate collapse by deferring all readout to final-time measurements on co-evolving ancillas. This yields a unified description of weak and strong measurements and reproduces projected correlators or symmetrized/antisymmetrized correlators in the corresponding limits (Oehri et al., 2015). A plausible implication is that multitime quantum-field correlations always require an explicit statement of the intervention model whenever they are given an operational meaning.

6. Curved spacetime, black holes, and field-theoretic multi-time amplitudes

In quantum field theory in curved spacetime, the same distinction between single-time thermality and multitime structure becomes especially sharp. Detector-based analyses of Hawking radiation show that one-time observables depend on the Wightman two-point function, but two-time coincidence probabilities are governed by connected closed-time-path four-point functions of the form

$P$5

and, more generally, $P$6-time detection probabilities depend on $P$7-point closed-time-path correlators (Anastopoulos et al., 2019). For a scalar field in Schwarzschild spacetime, the resulting multi-time correlations depend explicitly on angular variables, transmission and reflection amplitudes, and scattering phase delays, even though the one-time detector response is the standard thermal Hawking rate with greybody factor (Anastopoulos et al., 2019).

This logic is sharpened in the collapse-and-evaporation analysis of black holes. In a $P$8-dimensional collapsing-shell model, the in-vacuum Wightman function is exactly computable, and the one-event detector probability is thermal: $P$9 But the two-event coherence function

$$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$0

splits into an uncorrelated term, an Unruh-vacuum contribution, a non-Unruh correction, and a memory term $$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$1 whose phase depends explicitly on the collapse parameter $$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$2 (Xenos et al., 11 Jun 2026). The paper’s claim is that late-time multi-time correlations are not fully reproduced by the Unruh vacuum and can retain pre-collapse information even when all single-time observables are thermal (Xenos et al., 11 Jun 2026). This does not settle the information paradox, but it shows that formulations based solely on single-time reduced states are incomplete.

A different field-theoretic line replaces correlators by spacetime amplitudes. In multi-time wave-function formulations of QFT, the basic object is a multi-time Fock function $$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$3 defined on spacelike configurations, with one time variable per particle. In the emission–absorption model, $$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$4 satisfies coupled multi-time PDEs across particle-number sectors and is related on spacelike configurations to vacuum-to-state matrix elements of Heisenberg fields, schematically

$$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$5

Under natural assumptions, this multi-time formulation, the Tomonaga–Schwinger picture, and the Heisenberg-picture field representation are equivalent (Petrat et al., 2013, Lienert et al., 2017). These objects are not standard Wightman or Feynman correlators, but they provide a particle-position representation of off-diagonal field amplitudes across spacetime points.

The field-theoretic landscape also includes operational proposals that access otherwise unavailable spacetime correlations. Quantum-controlled superpositions of detector activation times make interference terms such as

$$P\big[\{ \alpha(t_i) ; t_i\}_{i=1}^k\big]=\Big\langle \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \prod_{i=1}^k\hat \delta (\hat a(t_i) - \alpha(t_i)) \begin{smallmatrix} \circ \ \circ \end{smallmatrix} \Big\rangle .$$6

appear in the reduced detector state; these terms are built from vacuum Wightman functions evaluated between spacetime regions belonging to different branches of a temporal superposition (Henderson et al., 2020). This indicates that some multi-time field correlations become operationally accessible only when the temporal structure of the probe itself is coherent.

Taken together, these developments define multi-time quantum-field correlations as a subject at the intersection of phase-space optics, detector theory, open-system memory, and spacetime QFT. The recurring themes are the nontrivial role of time ordering, the insufficiency of single-time reduced descriptions, the need for explicit operational prescriptions, and the emergence of temporal structures—nonclassicality, memory, interference, or information retention—that are invisible at equal time.