Eq. (40) — the partial trace over the site intertwines conjugation with conjugation
ProvedNonAbelianDMNRG.Recursion.eq_40p2o-batch-pfp2ap2o-gran-per-chapterp2o-plan-paperp2o-v1partial-tracesymmetry
Let and be finite index sets and consider operators on , where the second factor plays the role of the last NRG site. Write for the partial trace over the second factor. For every matrix on , every unitary matrix on (that is, ) and every operator on ,
In the paper the symmetry transformation on the first sites factors as , with acting on site and on sites ; this identity is the step that moves the symmetry through the trace over site .
Formalization Note The statement holds for an arbitrary ; only needs to be unitary. Operators are indexed with the block factor first and the site factor second; the partial trace is the published WildeQIT.partialTraceRight.
Preamble
import Mathlib import Definitions.Def_WildeQIT_partialTrace open Matrix open scoped Kronecker
Formal statement
namespace NonAbelianDMNRG.Recursion
/-- **Eq. (40), p. 7.** Tracing out the right factor (the site) intertwines conjugation by
`V ⊗ L` with conjugation by `V`, for any `V` and any unitary `L`:
`Tr_site{(V ⊗ L) R (V ⊗ L)⁺} = V Tr_site{R} V⁺`. -/
theorem eq_40 {a b : Type} [Fintype a] [Fintype b] [DecidableEq b]
(V : Matrix a a ℂ) (L : Matrix b b ℂ) (hL : Lᴴ * L = 1)
(R : Matrix (a × b) (a × b) ℂ) :
WildeQIT.partialTraceRight ((V ⊗ₖ L) * R * (V ⊗ₖ L)ᴴ) =
V * WildeQIT.partialTraceRight R * Vᴴ := by sorry
end NonAbelianDMNRG.Recursion
Source
Tóth, Moca, Legeza, Zaránd, Density matrix numerical renormalization group for non-Abelian symmetries, arXiv:0802.4332v1, p. 7, eq. (40)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.