Key claim

A finite temporal bond dimension χ bounds the Schmidt rank of an influence matrix across a time cut. In a converged calculation it can also impose an operational infrared timescale, but χ alone neither defines that timescale nor turns numerical truncation into a physical bath.

Near an impurity quantum critical point, this distinction becomes visible. The physical relaxation time diverges, while finite temperature and finite temporal compression eventually stop the numerical growth of the fitted relaxation time.

The two effects may produce the same leading cutoff in one observable. They need not produce the same multi-time process.

Theme

Temporal tensor networks, critical slowing down, and cutoff-limited Kibble–Zurek dynamics.

Guiding question

When a Keldysh influence functional is compressed into a finite-χ temporal matrix-product state, in what sense does χ act as an infrared regulator, and where does the analogy with temperature fail?

Setup

Consider the two-impurity Anderson model used by Lotem and collaborators. Two interacting impurity orbitals couple to independent fermionic baths and to each other through an exchange K.

At a critical value Kc, the model reaches a non-Fermi-liquid fixed point in the two-channel-Kondo universality class. Write the tuning field as

g=K−Kc.

The leading symmetry-allowed perturbation has scaling dimension

Δ=12.

Its RG eigenvalue is therefore

y=1−Δ=12,

and close to the fixed point

dgdℓ=yg+O(g2).

Running the flow until g(ℓ∗)∼1 gives

eℓ∗∼|g|−1/y.

The corresponding crossover energy and time obey

T∗∼|g|1/y=|g|2,τphys∼(T∗)−1=A|g|−2.

Here A is a nonuniversal amplitude with the units required by the definition of g and time. The exponent 2 is the universal datum.

Universality class

The two-impurity critical point is in the 2CK universality class under the stated symmetries. Calling it simply “the 2CK model” would erase the microscopic distinction.

Analysis

1. What temporal compression represents

Trace out a bath after discretizing the forward and backward Keldysh contours. Its complete action on an impurity history is encoded in an influence object

I[s+,s−].

After folding the two contours, time steps become sites of a one-dimensional tensor network. A temporal cut separates an earlier history from a later history, so the exact influence object admits a Schmidt decomposition

I=∑a=1RsaIa[<]⊗Ia[>].

Approximating this state by a temporal MPS of bond dimension χ retains at most χ Schmidt channels across each cut:

Rretained≤χ.

This is the exact kinematic statement. It says how much temporal information can cross the cut in the chosen representation.

It does not say that correlations vanish after a fixed number of time steps. A finite-dimensional transfer operator can support long exponential tails, degeneracies, or oscillatory modes. Finite bond dimension, finite Markov order, and finite memory duration are different properties.

Rank is not range

Bond dimension limits the number of transmitted temporal modes. It does not put a hard upper bound on how long the slowest retained mode can persist.

For a time-translation-invariant bath, a semigroup influence matrix uses a single repeated tensor. Suppose its normalized transfer operator has eigenvalues

|λ1|=1>|λ2|≥|λ3|≥⋯.

If λ1 is isolated and the observable overlaps with the λ2 sector, the transfer spectrum defines a correlation time

τtr=−δtlog⁡|λ2/λ1|.

This formula is meaningful only after specifying the transfer operator, time step, normalization, spectral sector, and observable overlap. A near-degeneracy or a vanishing overlap changes the relevant time.

An operational cutoff τχ is therefore best defined through evidence: convergence of chosen observables with χ, convergence of the relevant transfer spectrum, or a scaling collapse controlled by one extracted time. It should not be inferred from χ alone.

2. What the 2026 calculation actually finds

Lotem et al. combine separate semigroup influence matrices for four fermionic baths with a spatial MPS for the coupled auxiliary spaces and impurity orbitals. They apply the construction to sudden quenches and ramps through the two-impurity critical point.

Away from Kc, the late-time interimpurity correlation is fitted by

⟨S1⋅S2⟩(t)−⟨S1⋅S2⟩eq∼e−t/τ.

Both temperature T and finite SGIM bond dimension limit the fitted τ. The authors introduce the empirical leading-rate ansatz

τ−1(K;T,χ)≃τ−1(K)+c(T,χ).

Because the zero-temperature physical rate vanishes at criticality, they estimate

c(T,χ)=τ−1(Kc;T,χ)

and recover

τ−1(K)≃τ−1(K;T,χ)−τ−1(Kc;T,χ).

Their low-T, sufficiently large-χ data collapse supports this subtraction for that fitted observable and parameter range. The paper presents additivity as an assumption tested numerically, not as a theorem about temporal MPS truncation.

A Matthiessen-type ansatz

Independent rates often add when distinct slow processes produce simple poles. Operator mixing, K-dependent truncation error, or nonexponential tails can invalidate that picture.

At K=Kc, the same work reports relaxation consistent with t−3/2. It also states that no established theoretical prediction for this exponent is known. The observation is a numerical result of the study, not a general 2CK law.

This qualification matters because an algebraic tail has no unique exponential relaxation time. The fitted off-critical τ, the critical power-law window, and the SGIM transfer time are three distinct quantities that may agree only in a controlled scaling regime.

3. Temperature and finite χ stop the calculation differently

RegulatorWhat changesPhysical statusWhat a matched relaxation time proves
T>0The bath state and its correlation functions; KMS periodicity introduces a thermal timescalePhysical control parameterThe selected observable has the same leading cutoff time
Physical decoherenceThe system–environment dynamics or reduced generatorPhysical processThe selected decay rate can match
Finite χThe variational representation of the multi-time influence objectNumerical approximationThe calculation resolves the selected observable only up to a comparable time
Finite observation timeThe available fitting window and frequency resolutionMeasurement or analysis limitationOnly that the inferred scale cannot exceed the window reliably

Temperature acts on every correlation constrained by the thermal state and KMS relations. A finite-χ approximation discards particular temporal Schmidt channels selected by the compression procedure.

Matching one scalar time does not imply equality of response functions, counting statistics, out-of-time-order correlations, or interventions at several times. Those are properties of the full process tensor.

The safe analogy is structural:

finite spatial size↔finite accessible correlation length,
finite temporal accuracy↔finite accessible scaling time.

Spatial finite-entanglement scaling supplies a useful precedent: a finite-χ MPS often generates an effective correlation length near a critical state. Transferring that idea to a temporal influence matrix is a hypothesis to test through the temporal transfer spectrum and observable collapse, not an automatic identity.

4. Cutoff-limited Kibble–Zurek scaling

Ramp the tuning field linearly through the critical point:

g(t)=vt.

Write the physical critical relaxation law in general form as

τ(g)=A|g|−a,a=11−Δ.

The Kibble–Zurek time follows by equating the remaining time to the critical point with the instantaneous relaxation time:

tQ=A(vtQ)−a.

Therefore

tQ1+a=Av−a,

and

tQ=A1−Δ2−Δv−12−Δ.

For Δ=1/2,

tQ=A1/3v−2/3.

Now suppose temperature, compression, or observation time limits the resolved relaxation to an operational ceiling τIR. The zero-cutoff KZ construction applies only while

tQ≲τIR.

Setting tQ(vIR)=τIR gives

vIR∼A1−ΔτIR−(2−Δ).

Thus, for the 2CK value,

vIR∼A1/2τIR−3/2.

The exponent 3/2 follows from critical scaling. The prefactor depends on the normalization of g, the relaxation amplitude, and the operational definition of the cutoff.

Lotem et al. define the impurity dissipated work for a K ramp by

⟨Wd⟩=∫K0Kf[⟨S1⋅S2⟩(t)−⟨S1⋅S2⟩eq]dK.

For this protocol, their scaling ansatz and simulations give three velocity regimes:

  1. Fast ramps are dominated by microscopic transients and are nonuniversal.

  2. Intermediate ramps approach the impurity KZ law

    ⟨Wd⟩∼v1/(2−Δ)=v2/3.
  3. At the slowest ramps, finite temperature and finite-χ resolution cut off the critical growth. The reported response crosses to

    ⟨Wd⟩∝v.

The last statement belongs to this definition of integrated dissipated work and this ramp protocol. Linear response can produce different powers for a local rate, a total work, or a protocol whose duration is held fixed.

Slower can be worse

Reducing v increases the ideal freeze-out time. Once that time exceeds the numerical or thermal window, a slower ramp probes the regulator more strongly than the critical fixed point.

False claim to diagnose

If finite-T and finite-χ relaxation curves collapse after subtracting one constant rate, then finite bond dimension is physically equivalent to raising the bath temperature.

The collapse establishes a narrower result: the dominant cutoff in one fitted relaxation rate can be parameterized by the same scalar correction over the tested range.

Physical equivalence would require agreement of the complete multi-time process, including thermal consistency conditions and several independent observables. The reported collapse does not supply that evidence.

What follows — and what does not

StatementStatus
A temporal MPS with bond dimension χ has Schmidt rank at most χ across a temporal cut.Exact representation statement
Every finite-χ temporal MPS has a hard memory range.False
A gapped, normalized uniform transfer operator defines exponential correlation times.True under spectral and overlap assumptions
τχ is a universal function of χ alone.Unsupported
The 2IAM critical point used here has Δ=1/2 and $\tau\simK-K_c
Finite T and finite χ support an additive correction to the fitted rate in the reported data.Numerical result plus an empirical ansatz
Additive rates hold for every observable and temporal truncation.False
The observed t−3/2 critical decay is an established universal theorem.False; the paper labels it an unexplained numerical observation
The KZ time scales as v−1/(2−Δ).Follows from the stated relaxation law and linear ramp
The amplitude-free formula vIR=τIR−(2−Δ) is exact.False unless units set A=1
Matching τIR makes finite temperature and finite χ the same physical perturbation.False

Exercise

Assume

τ(g)=A|g|−1/(1−Δ)

until an operational ceiling τIR is reached, and ramp g(t)=vt.

  1. Derive tQ(v) without setting A=1.
  2. Derive the crossover velocity vIR from tQ(vIR)=τIR.
  3. Evaluate both exponents at Δ=1/2.
  4. Suppose two simulations have the same τIR, one from finite T and one from finite χ. Name one one-time observable and one multi-time diagnostic you would compare before calling the regulators equivalent.
  5. A uniform SGIM transfer operator has λ1=1 and λ2=0.99 at time step δt. Compute its leading transfer time. State one reason this eigenvalue might not control the measured impurity observable.
Hint 1

With a=1/(1−Δ), solve tQ=A(vtQ)−a by collecting all powers of tQ on the left.

Hint 2

An eigenmode contributes only if the boundary conditions and the observable have nonzero overlap with its left and right eigenvectors.

Oral check 1. What does finite χ bound exactly?

Oral check 2. Why is a critical power law incompatible with assigning a unique exponential relaxation pole?

Solution

Let a=1/(1−Δ). Then

tQ1+a=Av−a.

Because

a1+a=12−Δ,11+a=1−Δ2−Δ,

we obtain

tQ=A1−Δ2−Δv−12−Δ.

At crossover,

τIR=A1−Δ2−ΔvIR−12−Δ,

so

vIR=A1−ΔτIR−(2−Δ).

For Δ=1/2,

tQ=A1/3v−2/3,vIR=A1/2τIR−3/2.

A useful one-time comparison is the full relaxation curve of ⟨S1⋅S2⟩(t) at several K. A stronger multi-time comparison could use a two-time response function, full counting statistics, or interventions reconstructed as a process tensor. Agreement of the fitted τ alone is insufficient.

For the transfer spectrum,

τtr=−δtlog⁡0.99≈99.5δt.

This mode will not control an observable if the observable has zero overlap with its symmetry sector. Another sector, a nonnormal transient, or a power-law window may then set the apparent timescale.

Check your understanding

You increase χ and find that the fitted relaxation time doubles, while the second transfer-matrix eigenvalue and a two-time response function have not converged.

Which statement is justified: “the physical memory doubled,” “the accessible relaxation window increased,” or “the critical exponent changed”? Explain what further convergence test would distinguish them.

You may also reply with “deeper,” “too easy,” “too hard,” or your derivation.

Further Reading

Connections and next step

  • RG layer: the relevant field with Δ=1/2 gives T∗∼g2 and τ∼g−2.
  • Keldysh layer: tracing out the baths produces a multi-time influence object without a Born–Markov approximation.
  • Tensor-network layer: χ bounds temporal Schmidt rank; the transfer spectrum, not χ alone, supplies candidate correlation times.
  • Numerical layer: rate subtraction is a tested collapse ansatz for one observable, not a universal law of truncation.
  • Dynamical layer: an infrared ceiling removes the slow end of the KZ window and can expose cutoff-limited linear response.
  • Next step: extract the relevant τχ from symmetry-resolved SGIM transfer eigenvalues and compare it with observable-by-observable convergence.
  • Revisit: contrast this numerical regulator with the physical KMS and fluctuation–dissipation constraints in the July 31 entry.